- Timestamp:
- 10/24/19 19:00:48 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/dynamico/meshes.py
r984 r987 20 20 21 21 radian=180/math.pi # convert from radians to degrees 22 def Triangulation(lon,lat): return tri.Triangulation((lon*radian+360.)%360., radian*lat) 22 23 23 24 #------------------- Hybrid mass-based coordinate ------------- … … 451 452 plt.colorbar() 452 453 plt.ylim((-90,90)) 453 plt.xlim(( -180,180))454 plt.xlim((0,360)) 454 455 def plot_i(self,data, **kwargs): 455 456 self.plot(self.primal,data, **kwargs) … … 512 513 gridfile.normalize(self) 513 514 self.to_dynamico() 514 self.primal = tri.Triangulation(lon_i*radian, lat_i*radian)515 self.dual = tri.Triangulation(lon_v*radian, lat_v*radian)516 self.triedge = tri.Triangulation(lon_e*radian, lat_e*radian)515 self.primal = Triangulation(lon_i, lat_i) 516 self.dual = Triangulation(lon_v, lat_v) 517 self.triedge = Triangulation(lon_e, lat_e) 517 518 518 519 class Unstructured_PMesh(Abstract_Mesh): # Mesh data distributed across MPI processes … … 712 713 self.com_edges.set_dynamico_transfer('edge') 713 714 714 self.primal = tri.Triangulation(self.lon_i*radian, self.lat_i*radian)715 self.dual = tri.Triangulation(self.lon_v*radian, self.lat_v*radian)716 self.triedge = tri.Triangulation(self.lon_e*radian, self.lat_e*radian)715 self.primal = Triangulation(lon_i, self.lat_i) 716 self.dual = Triangulation(lon_v, self.lat_v) 717 self.triedge = Triangulation(lon_e, self.lat_e) 717 718 718 719 def plot_patches(self, ax, clim, degree, bounds, lon, lat, data):
Note: See TracChangeset
for help on using the changeset viewer.