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
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