Changeset 801 for codes/icosagcm/devel/Python/test

Timestamp:

01/17/19 10:42:38 (5 years ago)

Author:

dubos

Message:

devel/unstructured : implement reference vs physical mesh for spherical meshes

File:

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

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

@@ -11 +11 @@
 from dynamico.meshes import MPAS_Format, Unstructured_PMesh as PMesh, Local_Mesh as Mesh
 from dynamico import time_step
+from dynamico import maps
 print '...Done'
 
@@ -24 +25 @@
 
 print 'Omega, planetary PV', Omega, 2*Omega/gh0
+planet = maps.SphereMap(radius, Omega)
 
-def f(lon,lat): return 2*Omega*np.sin(lat) # Coriolis parameter
 print 'Reading MPAS mesh ...'
 meshfile = MPAS_Format('grids/x1.%d.grid.nc'%grid)
 pmesh = PMesh(comm,meshfile)
 pmesh.partition_metis()
-mesh = Mesh(pmesh, llm, nqdyn, radius, f)
+mesh = Mesh(pmesh, llm, nqdyn, planet)
 print '...Done'
 lon, lat = mesh.lon_i, mesh.lat_i