def rain(x,y,t):
    """Function to calculate "rain"
    input x,y should be considered to be numpy arrays
    abd t a scalar
    """
    if t<=4.0:
        return (x+y)*1.0
    else:
        return 0.0*x