Changeset 979 for codes

Timestamp:

10/08/19 16:11:10 (5 years ago)

Author:

dubos

Message:

devel/Python : increased mountain height in NH orographic test case to 750m

File:

• codes/icosagcm/devel/Python/test/py/slice_orographic.py (modified) (5 diffs)

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

@@ -9 +9 @@
 dx, Nx, llm, tau_relax = 500., 400, 60, 25.
 #dx, Nx, llm, tau_relax = 250., 800, 120, 25.
-T1, N1, N2 = 30., 6, 20 # check energy every 60s, plot every 10min, run 10h
+T1, N1, N2 = 60., 6, 10 # check energy every 60s, plot every 6min, run 1h
 #dx, Nx, llm, tau_relax = 20000., 200, 200, 1e3
-#T1, N1, N2 = 1800., 4, 6 # check energy every 60s, plot every 10min, run 2h
+
 
 Lx=Nx*dx
 metric, thermo, BC, m0ik, gas0, Phi0_il, u0_jk = DCMIP.DCMIP21(
-    Lx, Nx, llm, h0=250., ztop=3e4, dd=4000, xi=3200)
+    Lx, Nx, llm, h0=750., ztop=3e4, dd=4000, xi=3200)
 #    Lx, Nx, llm, h0=250., ztop=3e4, dd=10*dx, xi=8*dx)
 Sik = m0ik*gas0.s
@@ -34 +34 @@
 BC.rigid_bottom = False
 rho_bot = BC.rho_bot
-BC.rho_bot = 1e6*rho_bot
+BC.rho_bot = 1e2*rho_bot
 
 def go(model,courant,scheme,start):
@@ -45 +45 @@
         zil = Phi_il/metric.g0/1000
         zjl = metric.mj(zil)
-        plt.figure(figsize=(12,5))
-#        plt.figure(figsize=(6,5))
+        plt.figure(figsize=(8,4))
 #        plt.pcolormesh(xjl, zjl, gas_ik.T-gas0.T)
-#        plt.contourf(xil, zil, wil, np.arange(-.2,.2,.02))
-        plt.contourf(xil, zil, wil, np.arange(-1.,1.,.1))
+        plt.contourf(xil, zil, wil, np.arange(-3.,3.,.1))
         plt.plot(xil,zil[:,llm/10])
         plt.plot(xil,zil[:,llm/5])
@@ -57 +55 @@
         plt.ylabel('z (km)')
 #        plt.ylim((-1,30.)), plt.xlim((-dx/20.,dx/10.))
-        plt.ylim((0.,10.)), plt.xlim((-10.,10.))
+        plt.ylim((0.,13.)), plt.xlim((-20.,20.))
         plt.colorbar()
         plt.title("t = %f s" % (t))
@@ -80 +78 @@
 time_scheme = lambda dt : time_step.ARK2(bwd_fast_slow, dt, a32=0.7)
 go(model, 3.0, time_scheme, start )
+#time_scheme = lambda dt : time_step.RK4(bwd_fast_slow, dt)
+#go(model, 1.0, time_scheme, start )