Changeset 797 for codes

Timestamp:

01/10/19 17:37:58 (5 years ago)

Author:

jisesh

Message:

unstructured/Python : command-line arg for ztop in Baroclinic_3D_ullrich.py

File:

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

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

@@ -18 +18 @@
 getargs.add("--kappa_divgrad",  type=float, default=3.0e15)
 getargs.add("--kappa_curlcurl", type=float, default=3.0e15)
+getargs.add("--ztop", type=float, default=30000.)
 
 getargs.defaults(dt=360., llm=50)
@@ -54 +55 @@
     return pmesh
 
-def baroclinic_3D(pmesh, Lx,Ly,llm,ztop=25000.):
+def baroclinic_3D(pmesh, Lx,Ly,llm,ztop):
     Rd    = 287.0      # Gas constant for dryy air (j kg^-1 K^-1)
     T0    = 288.0      # Reference temperature (K) 
@@ -61 +62 @@
     u0    = 35.        # Reference zonal wind speed (m s^-1)
     a     = 6.371229e6 # Radius of the Earth (m)
-    ptop  = 2000.
     y0    = .5*Ly
     Cpd   = 1004.5
@@ -104 +104 @@
     y_ek, alpha_ek = np.meshgrid(mesh.lat_e+y0, alpha_k, indexing='ij')
 
-    print('ztop(ptop) according to Eq. 7:', T0/lap*(1.-(ptop/p0)**(Rd*lap/g))) 
     print(np.shape(alpha_k),np.shape(alpha_l))
     thermo = dyn.Ideal_perfect(Cpd, Rd, p0, T0)
@@ -130 +129 @@
     ptop = p(eta(1.))
     print( 'ptop (Pa) = ', gas.p[0,-1], ptop)
+    print('ztop(ptop) according to Eq. 7:', T0/lap*(1.-(ptop/p0)**(Rd*lap/g))) 
 
     params=dyn.Struct()
@@ -190 +190 @@
     
     pmesh = create_pmesh(nx,ny)
-    thermo, mesh, params, flow0, gas0 = baroclinic_3D(pmesh, Lx,Ly,llm)
+    thermo, mesh, params, flow0, gas0 = baroclinic_3D(pmesh, Lx,Ly,llm, args.ztop)
 
     mass_bl,mass_dak,mass_dbk = meshes.compute_hybrid_coefs(flow0[0])