Пример #1
0
def qinit(x,y):
    """
    Gaussian hump:
    """
    x0 = -93.0
    y0 = 28.0
    d = topotools.gcdist(x,y,x0,y0)  # great circle distance on earth
    ze = -d**2 / 2.e9
    amplitude = 1.0
    z = np.where(ze > -100., amplitude * np.exp(ze), 0.0)
    return z
Пример #2
0
def qinit(x,y):
    """
    Gaussian hump:
    """
    from numpy import where,exp
    x0 = -82.0
    y0 = -47.0
    d = topotools.gcdist(x,y,x0,y0)  # great circle distance on earth
    ze = -d**2/2.e9
    z = where(ze>-100., 20.e0*exp(ze), 0.)
    return z