Changeset 666 for codes/icosagcm/devel
- Timestamp:
- 01/20/18 18:47:41 (6 years ago)
- Location:
- codes/icosagcm/devel/Python
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/dynamico/__init__.py
r663 r666 5 5 import libs 6 6 libs.libicosa.dynamico_show_openmp() 7 libs.libicosa.dynamico_init_trace() -
codes/icosagcm/devel/Python/test/py/RSW2_MPAS_W02.py
r642 r666 1 print 'Loading DYNAMICO modules ...' 2 from dynamico import unstructured as unst 3 from dynamico.meshes import MPAS_Mesh as Mesh 4 from dynamico import time_step 5 print '...Done' 6 1 7 print 'Loading modules ...' 2 8 import math as math 3 9 import matplotlib.pyplot as plt 4 10 import numpy as np 5 print '...Done'6 7 print 'Loading DYNAMICO modules ...'8 from dynamico import unstructured as unst9 from dynamico.meshes import MPAS_Mesh as Mesh10 from dynamico import time_step11 11 print '...Done' 12 12 -
codes/icosagcm/devel/Python/test/py/slice_GW_NH.py
r656 r666 10 10 import numpy as np 11 11 import time 12 13 unst.setvar('nb_threads', 1)14 12 15 13 def run(mesh,metric,thermo,BC,caldyn_eta,u, m0ik,gas0,Phi0_il, T=300., courant=2.9): … … 97 95 98 96 def deform(x,eta): # x and eta go from 0 to 1 99 h0, dd, xi = 2000., 5000., 4000.97 h0, dd, xi = 500., 5000., 4000. 100 98 cos2=lambda x: np.cos(np.pi*x)**2 101 99 gauss=lambda x: np.exp(-x*x) … … 105 103 106 104 #Lx, nx, llm = 3e5, 100, 1e4, 30 107 Lx, nx, ztop, llm = 3e5, 300, 3e4, 80105 Lx, nx, ztop, llm = 3e5, 300, 1e4, 80 108 106 nqdyn, ny, dx = 1, 1, Lx/nx 109 107 mesh = Mesh(nx,ny,llm,nqdyn,nx*dx,ny*dx,0.)
Note: See TracChangeset
for help on using the changeset viewer.