Changeset 930

Timestamp:

07/02/19 12:02:06 (5 years ago)

Author:

dubos

Message:

devel/Python : periodize initial perturbation for Baroclinic test case

File:

• codes/icosagcm/devel/Python/test/py/Baroclinic_3D_ullrich.py (modified) (1 diff)

codes/icosagcm/devel/Python/test/py/Baroclinic_3D_ullrich.py

@@ -83 +83 @@
     def ulon(x,y,eta):
         u = -u0*(sin(pi*y/Ly)**2)*log(eta)*(eta**(-log(eta)/(b*b))) 
-        u = u + up*exp(-(((x-xc)**2+(y-yc)**2)/lp**2))
-        return  u
+        # periodize the envelope of the initial perturbation
+        A = exp(-(((x-xc)**2+(y-yc)**2)/lp**2))
+        A = A + exp(-(((x-xc-Lx)**2+(y-yc)**2)/lp**2))
+        A = A + exp(-(((x-xc+Lx)**2+(y-yc)**2)/lp**2))
+        return  u*A
 
     def tmean(eta) : return T0*eta**(Rd*lap/g)