Changeset 646 for codes/icosagcm/devel/Python/test/py/slice_GW_NH.py
- Timestamp:
- 12/21/17 15:11:18 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/test/py/slice_GW_NH.py
r631 r646 1 1 from dynamico.meshes import Cartesian_mesh as Mesh 2 from dynamico import meshes 2 3 from dynamico import unstructured as unst 3 4 from dynamico import dyn … … 15 16 caldyn_thermo = unst.thermo_entropy 16 17 g = mesh.dx/metric.dx_g0 17 caldyn = unst.Caldyn_NH(caldyn_thermo,caldyn_eta, mesh, metric,thermo,BC,g)18 caldyn = unst.Caldyn_NH(caldyn_thermo,caldyn_eta, mesh,thermo,BC,g) 18 19 19 20 Sik = m0ik*gas0.s … … 95 96 xx,ll = mesh.xx[:,0,:]/1000, mesh.ll[:,0,:] 96 97 97 #metric, thermo, BC, m0ik, gas0, Phi0_il, u0_jk = DCMIP.DCMIP31(Lx, nx, llm, deform=deform, dTheta=1.)98 metric, thermo, BC, m0ik, gas0, Phi0_il, u0_jk = DCMIP.DCMIP21(Lx, nx, llm, deform=deform, h0=0.)98 metric, thermo, BC, m0ik, gas0, Phi0_il, u0_jk = DCMIP.DCMIP31(Lx, nx, llm, deform=deform, dTheta=1.) 99 #metric, thermo, BC, m0ik, gas0, Phi0_il, u0_jk = DCMIP.DCMIP21(Lx, nx, llm, deform=deform, h0=0.) 99 100 100 101 BC.rho_bot = 1e6*BC.rho_bot # large stiffness 101 102 102 mass_bl,mass_dak,mass_dbk = unst.compute_hybrid_coefs(m0ik)103 unst. init_hybrid(mass_bl,mass_dak,mass_dbk)103 mass_bl,mass_dak,mass_dbk = meshes.compute_hybrid_coefs(m0ik) 104 unst.ker.dynamico_init_hybrid(mass_bl,mass_dak,mass_dbk) 104 105 105 106 run(mesh,metric,thermo,BC,unst.eta_mass,20., m0ik,gas0,Phi0_il, courant=2.) # lower courant number required if h0=2000m
Note: See TracChangeset
for help on using the changeset viewer.