Ignore:
Timestamp:
12/21/17 15:11:18 (7 years ago)
Author:
dubos
Message:

devel/unstructured : update some test cases

File:
1 edited

Legend:

Unmodified
Added
Removed
  • codes/icosagcm/devel/Python/test/py/slice_GW_NH.py

    r631 r646  
    11from dynamico.meshes import Cartesian_mesh as Mesh 
     2from dynamico import meshes 
    23from dynamico import unstructured as unst 
    34from dynamico import dyn 
     
    1516    caldyn_thermo = unst.thermo_entropy 
    1617    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) 
    1819 
    1920    Sik = m0ik*gas0.s 
     
    9596xx,ll = mesh.xx[:,0,:]/1000, mesh.ll[:,0,:] 
    9697 
    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.) 
     98metric, 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.) 
    99100 
    100101BC.rho_bot = 1e6*BC.rho_bot # large stiffness 
    101102 
    102 mass_bl,mass_dak,mass_dbk = unst.compute_hybrid_coefs(m0ik) 
    103 unst.init_hybrid(mass_bl,mass_dak,mass_dbk) 
     103mass_bl,mass_dak,mass_dbk = meshes.compute_hybrid_coefs(m0ik) 
     104unst.ker.dynamico_init_hybrid(mass_bl,mass_dak,mass_dbk) 
    104105 
    105106run(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.