Changeset 631 for codes/icosagcm/devel/Python/test/py/test_xios.py

Timestamp:

12/12/17 16:04:23 (6 years ago)

Author:

dubos

Message:

devel/Python : extract pure Python stuff from cython module unstructured.pyx

File:

• codes/icosagcm/devel/Python/test/py/test_xios.py (modified) (3 diffs)

codes/icosagcm/devel/Python/test/py/test_xios.py

@@ -4 +4 @@
 print '%d/%d starting'%(mpi_rank,mpi_size)
 
-import sys
 import numpy as np
-#import dynamico.wrap as wrap
 print 'import dynamico.unstructured'
 import dynamico.unstructured as unst
@@ -12 +10 @@
 import dynamico.xios as xios
 print 'Done.'
+from dynamico.meshes import MPAS_Mesh as Mesh, radian
 
 def boundaries(degree,points,lon,lat):
@@ -29 +28 @@
 #----------------------------- read MPAS mesh --------------------------------
 
-radian=unst.radian
 grid, llm, nqdyn, nqtot = 10242, 1,1,1 # 2562, 10242, 40962  
 Omega, radius, g, gh0 = 2.*np.pi/86400., 6.4e6, 1., 2.94e4
 N, T, courant = 40, 10800., 1.2 # simulation length = N*T
 print 'Omega, planetary PV', Omega, 2*Omega/gh0
-sys.stdout.flush()
 
 def f(lon,lat): return 2*Omega*np.sin(lat) # Coriolis parameter                                                                                                                
 print 'Reading MPAS mesh ...'
-sys.stdout.flush()
-mesh = unst.MPAS_Mesh('grids/x1.%d.grid.nc'%grid, llm, nqdyn, radius, f)
+mesh = Mesh('grids/x1.%d.grid.nc'%grid, llm, nqdyn, radius, f)
 print '...Done'
-sys.stdout.flush()
 
 lon_i, lat_i, lon_v, lat_v = [x*radian for x in mesh.lon_i, mesh.lat_i, mesh.lon_v, mesh.lat_v]