Changeset 1091

Timestamp:

Nov 7, 2013, 10:18:48 AM (11 years ago)

Author:

aslmd

Message:

PLANETOPLOT. better handling of map projections.

File:

• trunk/UTIL/PYTHON/planetoplot_v2/ppplot.py (modified) (2 diffs)

trunk/UTIL/PYTHON/planetoplot_v2/ppplot.py

@@ -757 +757 @@
             #mertab = np.r_[wlon[0]:wlon[1]:steplon] ; merlab = [0,0,0,1]
             #partab = np.r_[wlat[0]:wlat[1]:steplat] ; parlab = [1,0,0,0]
-            mertab = np.r_[-180.:180.:steplon] ; merlab = [0,0,0,1] #-360:360.
+            if np.abs(wlon[0]) < 180.1 and np.abs(wlon[1]) < 180.1:
+                mertab = np.r_[-180.:180.:steplon]
+            else:
+                mertab = np.r_[0.:360.:steplon]
+            merlab = [0,0,0,1]
             partab = np.r_[-90.:90.:steplat] ; parlab = [1,0,0,0]
             format = '%.1f'
@@ -770 +774 @@
             elif self.proj in ["ortho","moll","robin"]:
                 wlat[0] = None ; wlat[1] = None ; wlon[0] = None ; wlon[1] = None
+                lon_0 = np.ceil(lon_0) # reverse map if lon_0 is slightly below 180 with [0,360]
                 steplon = 30. ; steplat = 30.
                 if self.proj in ["moll"]: steplon = 60.