Changeset 797 for codes/icosagcm/devel
- Timestamp:
- 01/10/19 17:37:58 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/test/py/Baroclinic_3D_ullrich.py
r795 r797 18 18 getargs.add("--kappa_divgrad", type=float, default=3.0e15) 19 19 getargs.add("--kappa_curlcurl", type=float, default=3.0e15) 20 getargs.add("--ztop", type=float, default=30000.) 20 21 21 22 getargs.defaults(dt=360., llm=50) … … 54 55 return pmesh 55 56 56 def baroclinic_3D(pmesh, Lx,Ly,llm,ztop =25000.):57 def baroclinic_3D(pmesh, Lx,Ly,llm,ztop): 57 58 Rd = 287.0 # Gas constant for dryy air (j kg^-1 K^-1) 58 59 T0 = 288.0 # Reference temperature (K) … … 61 62 u0 = 35. # Reference zonal wind speed (m s^-1) 62 63 a = 6.371229e6 # Radius of the Earth (m) 63 ptop = 2000.64 64 y0 = .5*Ly 65 65 Cpd = 1004.5 … … 104 104 y_ek, alpha_ek = np.meshgrid(mesh.lat_e+y0, alpha_k, indexing='ij') 105 105 106 print('ztop(ptop) according to Eq. 7:', T0/lap*(1.-(ptop/p0)**(Rd*lap/g)))107 106 print(np.shape(alpha_k),np.shape(alpha_l)) 108 107 thermo = dyn.Ideal_perfect(Cpd, Rd, p0, T0) … … 130 129 ptop = p(eta(1.)) 131 130 print( 'ptop (Pa) = ', gas.p[0,-1], ptop) 131 print('ztop(ptop) according to Eq. 7:', T0/lap*(1.-(ptop/p0)**(Rd*lap/g))) 132 132 133 133 params=dyn.Struct() … … 190 190 191 191 pmesh = create_pmesh(nx,ny) 192 thermo, mesh, params, flow0, gas0 = baroclinic_3D(pmesh, Lx,Ly,llm )192 thermo, mesh, params, flow0, gas0 = baroclinic_3D(pmesh, Lx,Ly,llm, args.ztop) 193 193 194 194 mass_bl,mass_dak,mass_dbk = meshes.compute_hybrid_coefs(flow0[0])
Note: See TracChangeset
for help on using the changeset viewer.