Changeset 1027 for trunk/UTIL/PYTHON

Timestamp:

Sep 2, 2013, 4:21:01 PM (11 years ago)

Author:

aslmd

Message:

PYTHON planetoplot_v2. possibility to invert xaxis if xmin>xmax. added modulo and more ticks for xlabel non-projected 2d plots

Location:

trunk/UTIL/PYTHON/planetoplot_v2

Files:

• ppplot.py (modified) (4 diffs)
• set_ppclass.txt (modified) (1 diff)

trunk/UTIL/PYTHON/planetoplot_v2/ppplot.py

@@ -262 +262 @@
     strtab = []
     for tick in ax.get_xaxis().get_ticklabels():
-        num = float(tick.get_text())
-        strtab.append(num % mod)
+        onetick = tick.get_text()
+        if len(onetick) > 0:
+          num = float(onetick)
+          strtab.append(num % mod)
     ax.get_xaxis().set_ticklabels(strtab)
     return ax
@@ -490 +492 @@
         else:
             mpl.plot(x,y,marker=self.marker,label=self.label)
+        # AXES
+        ax = mpl.gca()
         # make log axes and/or invert ordinate
         # ... this must be after plot so that axis bounds are well-defined
@@ -495 +499 @@
         if self.logx: mpl.xscale("log") # not mpl.semilogx() because excludes log on y
         if self.logy: mpl.yscale("log") # not mpl.semilogy() because excludes log on x
-        if self.invert: ax = mpl.gca() ; ax.set_ylim(ax.get_ylim()[::-1])
+        if self.invert: ax.set_ylim(ax.get_ylim()[::-1])
+        if self.xmin is not None and self.xmax is not None:
+          if self.xmin > self.xmax:
+            ax.set_xlim(ax.get_xlim()[::-1])
+            self.xmin,self.xmax = self.xmax,self.xmin
         # add a label for line(s)
         if self.label is not None:
             if self.label != "":
                 mpl.legend(loc="best",fancybox=True)
-        # AXES
-        ax = mpl.gca()
         # format labels
         if self.swap: ax.xaxis.set_major_formatter(FormatStrFormatter(self.fmt))
@@ -670 +676 @@
             if self.logy: mpl.semilogy()
             if self.invert: ax.set_ylim(ax.get_ylim()[::-1])
+            if self.xmin is not None and self.xmax is not None:
+              if self.xmin > self.xmax: 
+                ax.set_xlim(ax.get_xlim()[::-1])
+                self.xmin,self.xmax = self.xmax,self.xmin
             if self.xmin is not None: ax.set_xbound(lower=self.xmin)
             if self.xmax is not None: ax.set_xbound(upper=self.xmax)
             if self.ymin is not None: ax.set_ybound(lower=self.ymin)
             if self.ymax is not None: ax.set_ybound(upper=self.ymax)
+            # set the number of ticks (hardcoded)
+            ax.xaxis.set_major_locator(MaxNLocator(10))
+            ## specific modulo labels
+            if self.modx is not None:
+                ax = labelmodulo(ax,self.modx)
         else:
             ## A 2D MAP USING PROJECTIONS (basemap)

trunk/UTIL/PYTHON/planetoplot_v2/set_ppclass.txt

@@ -10 +10 @@
   # z-axis: possible names to search for. add in the line with a separating ;
   #
-altitude;Alt;presnivs;bottom_top;altitude_abg;lev;PTOT;vert
+altitude;Alt;presnivs;bottom_top;altitude_abg;lev;PTOT;vert;p
   # t-axis: possible names to search for. add in the line with a separating ;
   #