- Timestamp:
- 10/08/19 16:11:10 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/test/py/slice_orographic.py
r617 r979 9 9 dx, Nx, llm, tau_relax = 500., 400, 60, 25. 10 10 #dx, Nx, llm, tau_relax = 250., 800, 120, 25. 11 T1, N1, N2 = 30., 6, 20 # check energy every 60s, plot every 10min, run 10h11 T1, N1, N2 = 60., 6, 10 # check energy every 60s, plot every 6min, run 1h 12 12 #dx, Nx, llm, tau_relax = 20000., 200, 200, 1e3 13 #T1, N1, N2 = 1800., 4, 6 # check energy every 60s, plot every 10min, run 2h 13 14 14 15 15 Lx=Nx*dx 16 16 metric, thermo, BC, m0ik, gas0, Phi0_il, u0_jk = DCMIP.DCMIP21( 17 Lx, Nx, llm, h0= 250., ztop=3e4, dd=4000, xi=3200)17 Lx, Nx, llm, h0=750., ztop=3e4, dd=4000, xi=3200) 18 18 # Lx, Nx, llm, h0=250., ztop=3e4, dd=10*dx, xi=8*dx) 19 19 Sik = m0ik*gas0.s … … 34 34 BC.rigid_bottom = False 35 35 rho_bot = BC.rho_bot 36 BC.rho_bot = 1e 6*rho_bot36 BC.rho_bot = 1e2*rho_bot 37 37 38 38 def go(model,courant,scheme,start): … … 45 45 zil = Phi_il/metric.g0/1000 46 46 zjl = metric.mj(zil) 47 plt.figure(figsize=(12,5)) 48 # plt.figure(figsize=(6,5)) 47 plt.figure(figsize=(8,4)) 49 48 # plt.pcolormesh(xjl, zjl, gas_ik.T-gas0.T) 50 # plt.contourf(xil, zil, wil, np.arange(-.2,.2,.02)) 51 plt.contourf(xil, zil, wil, np.arange(-1.,1.,.1)) 49 plt.contourf(xil, zil, wil, np.arange(-3.,3.,.1)) 52 50 plt.plot(xil,zil[:,llm/10]) 53 51 plt.plot(xil,zil[:,llm/5]) … … 57 55 plt.ylabel('z (km)') 58 56 # plt.ylim((-1,30.)), plt.xlim((-dx/20.,dx/10.)) 59 plt.ylim((0.,1 0.)), plt.xlim((-10.,10.))57 plt.ylim((0.,13.)), plt.xlim((-20.,20.)) 60 58 plt.colorbar() 61 59 plt.title("t = %f s" % (t)) … … 80 78 time_scheme = lambda dt : time_step.ARK2(bwd_fast_slow, dt, a32=0.7) 81 79 go(model, 3.0, time_scheme, start ) 80 #time_scheme = lambda dt : time_step.RK4(bwd_fast_slow, dt) 81 #go(model, 1.0, time_scheme, start )
Note: See TracChangeset
for help on using the changeset viewer.