Changeset 678 for codes/icosagcm/devel

Timestamp:

02/08/18 15:35:25 (6 years ago)

Author:

dubos

Message:

devel/unstructured : fixed scaling factors in DCMIP31

File:

• codes/icosagcm/devel/Python/test/py/NH_3D_DCMIP31.py (modified) (3 diffs)

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

@@ -13 +13 @@
 
 Cpd, Rd, g = 1004.5, 287., 9.81
+u0, dT, Xfactor =  20., 1., 125.
+radius = 6.37122e6/Xfactor
 N, Teq, peq = 0.01, 300., 1e5   # background 
-lonc, d2 = 2.*math.pi/3., 25e6  # perturbation
 N2, g2 = N*N, g*g
 N2_g2, g2_N2 = N2/g2, g2/N2
 G, kappa = g2/(N2*Cpd), Rd/Cpd
+lonc, d2 = 2.*math.pi/3., (625e3/Xfactor)**2  # perturbation
 
 def psTs(lat) : 
@@ -83 +85 @@
 #------------------------ main program -------------------------
 
-ztop, u0, dT, Xfactor =  1e4, 20., 1., 125.
-radius, grid, llm = 6.37122e6/Xfactor, 10242, 20
-T, Nslice, courant = 360., 10, 3.0
+ztop, grid, llm = 1e4, 10242, 20
+T, Nslice, courant = 45000./Xfactor, 10, 3.0
 
 junk, ptop = T0_p0(0.,0.,g*ztop)
@@ -133 +134 @@
 
     mesh.plot_i(w[:,llm/2])
-    plt.title('vertical velocity at t=%d'%(it*T))
+    plt.title('Vertical velocity at t=%d'%(it*T))
     plt.savefig('fig_NH_3D_DCMIP31/w%02d.png'%it)
+    plt.close()
+    mesh.plot_i(z[:,llm/2])
+    plt.title('Altitude at t=%d'%(it*T))
+    plt.savefig('fig_NH_3D_DCMIP31/z%02d.png'%it)
     plt.close()