Changeset 678 for codes/icosagcm/devel/Python/test/py/NH_3D_DCMIP31.py
- Timestamp:
- 02/08/18 15:35:25 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/test/py/NH_3D_DCMIP31.py
r676 r678 13 13 14 14 Cpd, Rd, g = 1004.5, 287., 9.81 15 u0, dT, Xfactor = 20., 1., 125. 16 radius = 6.37122e6/Xfactor 15 17 N, Teq, peq = 0.01, 300., 1e5 # background 16 lonc, d2 = 2.*math.pi/3., 25e6 # perturbation17 18 N2, g2 = N*N, g*g 18 19 N2_g2, g2_N2 = N2/g2, g2/N2 19 20 G, kappa = g2/(N2*Cpd), Rd/Cpd 21 lonc, d2 = 2.*math.pi/3., (625e3/Xfactor)**2 # perturbation 20 22 21 23 def psTs(lat) : … … 83 85 #------------------------ main program ------------------------- 84 86 85 ztop, u0, dT, Xfactor = 1e4, 20., 1., 125. 86 radius, grid, llm = 6.37122e6/Xfactor, 10242, 20 87 T, Nslice, courant = 360., 10, 3.0 87 ztop, grid, llm = 1e4, 10242, 20 88 T, Nslice, courant = 45000./Xfactor, 10, 3.0 88 89 89 90 junk, ptop = T0_p0(0.,0.,g*ztop) … … 133 134 134 135 mesh.plot_i(w[:,llm/2]) 135 plt.title(' vertical velocity at t=%d'%(it*T))136 plt.title('Vertical velocity at t=%d'%(it*T)) 136 137 plt.savefig('fig_NH_3D_DCMIP31/w%02d.png'%it) 138 plt.close() 139 mesh.plot_i(z[:,llm/2]) 140 plt.title('Altitude at t=%d'%(it*T)) 141 plt.savefig('fig_NH_3D_DCMIP31/z%02d.png'%it) 137 142 plt.close() 138 143
Note: See TracChangeset
for help on using the changeset viewer.