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

Timestamp:

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

Author:

dubos

Message:

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

File:

• codes/icosagcm/devel/Python/test/py/RSW_MPAS_W02.py (modified) (2 diffs)

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

@@ -1 @@
-import sys
 print 'Loading modules ...'
-sys.stdout.flush()
-
 import math as math
-# select non-interactive backend, cf http://stackoverflow.com/questions/4931376/generating-matplotlib-graphs-without-a-running-x-server
-import matplotlib
-matplotlib.use('Agg') 
 import matplotlib.pyplot as plt
 import numpy as np
+print '...Done'
 
 print 'Loading DYNAMICO modules ...'
-sys.stdout.flush()
 from dynamico import unstructured as unst
+from dynamico.meshes import MPAS_Mesh as Mesh
 from dynamico import time_step
 print '...Done'
-sys.stdout.flush()
 
 grid, llm, nqdyn = 10242, 1,1 # 2562, 10242, 40962
@@ -22 +16 @@
 
 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, lat = mesh.lon_i, mesh.lat_i
 x,y,z = np.cos(lat)*np.cos(lon), np.cos(lat)*np.sin(lon), np.sin(lat)