Changeset 780

Timestamp:

11/20/18 12:44:31 (5 years ago)

Author:

jisesh

Message:

devel/Python : writing z output and fixed issue with theta values in Baroclinic_3D_ullrich.py

File:

• codes/icosagcm/devel/Python/test/py/Baroclinic_3D_ullrich.py (modified) (4 diffs)

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

@@ -64 +64 @@
     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))
+        u = u + up*exp(-(((x-xc)**2+(y-yc)**2)/lp**2))
         return  u
 
@@ -134 +134 @@
 
 def diagnose(Phi,S,m,W):
-    s=S/m ; s=.5*(s+abs(s))
+    s=S/m
     for l in range(llm):
         v[:,l]=(Phi[:,l+1]-Phi[:,l])/(g*m[:,l])
@@ -201 +201 @@
 
     T  = args.T
-    dt = 360.
+    dt = 100.#180.#360.
     dz = flow0[3].max()/(params.g*llm)
     nt = int(math.ceil(T/dt))
@@ -265 +265 @@
             context.send_field_primal('theta', gas.s)
             context.send_field_primal('uz', w)
+            context.send_field_primal('z', z)
 
             print( 'ptop, model top (m) :', unst.getvar('ptop'), Phi.max()/unst.getvar('g'))