Changeset 692 for codes/icosagcm/devel
- Timestamp:
- 05/03/18 17:27:43 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/test/py/partition.py
r680 r692 80 80 meshfile = meshes.MPAS_Format('grids/%s.grid.nc'%grid) 81 81 pmesh = meshes.Unstructured_PMesh(comm, meshfile) 82 lmesh = meshes.Local_Mesh(pmesh) 82 83 def coriolis(lon,lat): return 0.*lat 84 llm, nqdyn, radius = 1,1,1. 85 lmesh = meshes.Local_Mesh(pmesh, llm, nqdyn, radius, coriolis) 83 86 84 87 (primal_deg, primal_vertex, dim_vertex, dim_cell, cell_owner, … … 103 106 buf.data = np.cos(10.*buf.data) 104 107 buf.update() # updates halo 105 lmesh.plot (ax,[-math.pi/2,math.pi/2], mydegree, mybounds, mylon, mylat, buf.data)108 lmesh.plot_patches(ax,[-math.pi/2,math.pi/2], mydegree, mybounds, mylon, mylat, buf.data) 106 109 plt.xlim(-190.,190.) 107 110 plt.ylim(-90.,90.) … … 111 114 buf.read_own(cell_owner) 112 115 buf.update() 113 lmesh.plot (ax,[0,mpi_rank+1], mydegree, mybounds, mylon, mylat, buf.data)116 lmesh.plot_patches(ax,[0,mpi_rank+1], mydegree, mybounds, mylon, mylat, buf.data) 114 117 plt.xlim(-190.,190.) 115 118 plt.ylim(-90.,90.)
Note: See TracChangeset
for help on using the changeset viewer.