Changeset 987

Timestamp:

10/24/19 19:00:48 (5 years ago)

Author:

dubos

Message:

devel/Python : ensure positive longitudes when plotting

File:

• codes/icosagcm/devel/Python/dynamico/meshes.py (modified) (4 diffs)

codes/icosagcm/devel/Python/dynamico/meshes.py

@@ -20 +20 @@
 
 radian=180/math.pi # convert from radians to degrees
+def Triangulation(lon,lat): return tri.Triangulation((lon*radian+360.)%360., radian*lat)
 
 #------------------- Hybrid mass-based coordinate -------------
@@ -451 +452 @@
         plt.colorbar()
         plt.ylim((-90,90))
-        plt.xlim((-180,180))
+        plt.xlim((0,360))
     def plot_i(self,data, **kwargs):
         self.plot(self.primal,data, **kwargs)
@@ -512 +513 @@
         gridfile.normalize(self)
         self.to_dynamico()
-        self.primal  = tri.Triangulation(lon_i*radian, lat_i*radian)
-        self.dual    = tri.Triangulation(lon_v*radian, lat_v*radian)
-        self.triedge = tri.Triangulation(lon_e*radian, lat_e*radian)        
+        self.primal  = Triangulation(lon_i, lat_i)
+        self.dual    = Triangulation(lon_v, lat_v)
+        self.triedge = Triangulation(lon_e, lat_e)
         
 class Unstructured_PMesh(Abstract_Mesh): # Mesh data distributed across MPI processes
@@ -712 +713 @@
         self.com_edges.set_dynamico_transfer('edge')
 
-        self.primal  = tri.Triangulation(self.lon_i*radian, self.lat_i*radian)
-        self.dual    = tri.Triangulation(self.lon_v*radian, self.lat_v*radian)
-        self.triedge = tri.Triangulation(self.lon_e*radian, self.lat_e*radian)        
+        self.primal  = Triangulation(lon_i, self.lat_i)
+        self.dual    = Triangulation(lon_v, self.lat_v)
+        self.triedge = Triangulation(lon_e, self.lat_e)        
 
     def plot_patches(self, ax, clim, degree, bounds, lon, lat, data):