Changeset 827 for trunk/UTIL/PYTHON/mcd/mcd.py

Timestamp:

Nov 2, 2012, 5:15:38 PM (12 years ago)

Author:

aslmd

Message:

PYTHON UTIL mcd interface. few improvements: set bounds, set better resolution for figs, address Steve remarks.

File:

• trunk/UTIL/PYTHON/mcd/mcd.py (modified) (7 diffs)

trunk/UTIL/PYTHON/mcd/mcd.py

@@ -80 +80 @@
         self.fixedlt = False
         self.zonmean = False
+        self.min2d = None
+        self.max2d = None
+        self.dpi = 80.
 
     def viking1(self): self.name = "Viking 1 site. MCD v4.3 output" ; self.lat = 22.48 ; self.lon = -49.97 ; self.xdate = 97.
@@ -101 +104 @@
         if not oneline: self.title = self.title + "\n"
         if self.lats is None:  self.title = self.title + " Latitude " + str(self.lat) + "N"
-        if self.zonmean: self.title = self.title + "Zonal mean over all longitudes."
-        elif self.lons is None: self.title = self.title + " Longitude " + str(self.lon) + "E"
+        if self.zonmean and self.lats is not None and self.xzs is not None: 
+            self.title = self.title + "Zonal mean over all longitudes."
+        elif self.lons is None: 
+            self.title = self.title + " Longitude " + str(self.lon) + "E"
         if self.xzs is None:   
             self.vertunits()
@@ -513 +518 @@
       # The size * the dpi gives the final image size
       #   a4"x4" image * 80 dpi ==> 320x320 pixel image
-      canvas.print_figure(figname, dpi=80)
+      canvas.print_figure(figname, dpi=self.dpi)
 
 ###################
@@ -704 +709 @@
         ## define field. bound field.
         what_I_plot = np.transpose(field)
-        zevmin, zevmax = myplot.calculate_bounds(what_I_plot)  ## vmin=min(what_I_plot_frame), vmax=max(what_I_plot_frame))
+        zevmin, zevmax = myplot.calculate_bounds(what_I_plot,vmin=self.min2d,vmax=self.max2d)  
         what_I_plot = myplot.bounds(what_I_plot,zevmin,zevmax)
         ## define contour field levels. define color palette
@@ -734 +739 @@
       # The size * the dpi gives the final image size
       #   a4"x4" image * 80 dpi ==> 320x320 pixel image
-      canvas.print_figure(figname, dpi=80)
+      canvas.print_figure(figname, dpi=self.dpi)
 
     def htmlplot2d(self,tabtodo,figname="temp.png"):
@@ -775 +780 @@
         ## define field. bound field.
         what_I_plot = np.transpose(field)
-        zevmin, zevmax = myplot.calculate_bounds(what_I_plot)  ## vmin=min(what_I_plot_frame), vmax=max(what_I_plot_frame))
+        zevmin, zevmax = myplot.calculate_bounds(what_I_plot,vmin=self.min2d,vmax=self.max2d)  
         what_I_plot = myplot.bounds(what_I_plot,zevmin,zevmax)
         ## define contour field levels. define color palette
@@ -806 +811 @@
       # The size * the dpi gives the final image size
       #   a4"x4" image * 80 dpi ==> 320x320 pixel image
-      canvas.print_figure(figname, dpi=80)
+      canvas.print_figure(figname, dpi=self.dpi)