Changeset 1054 for trunk/UTIL

Timestamp:

Sep 27, 2013, 3:48:25 AM (11 years ago)

Author:

aslmd

Message:

planetoplot. added possibility to retrieve attribute from netCDF file. added an option to set value formatting through ppclass or pp.py. added a nice example used in a press release.

Location:

trunk/UTIL/PYTHON/planetoplot_v2

Files:

• examples/ppclass_additional/archean.sh (added)
• pp.py (modified) (1 diff)
• ppclass.py (modified) (6 diffs)
• ppplot.py (modified) (1 diff)

trunk/UTIL/PYTHON/planetoplot_v2/pp.py

@@ -71 +71 @@
 parser.add_option('-V','--void',action='store_true',dest='void',default=False,help="no colorbar, no title, no labels")
 parser.add_option('-U','--units',action='append',dest='units',type="string",default=None,help="units for the field")
+parser.add_option('-F','--fmt',action='append',dest='fmt',type="string",default=None,help="values formatting. ex: '%.0f' '%3.1e'")
 # -- 1D plot
 parser.add_option('-L','--linestyle',action='append',dest='linestyle',type="string",default=None,help="[1D] linestyle: '-' '--' '.' '..'")

trunk/UTIL/PYTHON/planetoplot_v2/ppclass.py

@@ -56 +56 @@
         except: pass
         print output ; output = ""
+
+# request a given attribute (and test if it is here)
+def ncattr(filename,char):
+    nc = netCDF4.Dataset(filename)
+    if hasattr(nc,char): ncattr=getattr(nc,char)
+    return ncattr
 
 # check a tab and exit if wrong. if just one string make it a list.
@@ -164 +170 @@
                       savtxt=False,\
                       modx=None,\
+                      fmt=None,\
                       xp=16,yp=8,\
                       missing=1.e25,\
@@ -223 +230 @@
         self.xp = xp ; self.yp = yp
         self.nxticks = nxticks ; self.nyticks = nyticks
+        self.fmt = fmt
 
     # print status
@@ -282 +290 @@
             self.missing = other.missing
             self.nxticks = other.nxticks ; self.nyticks = other.nyticks
+            self.fmt = other.fmt
         else:
             print "!! ERROR !! argument must be a pp object." ; exit()
@@ -825 +834 @@
                     if self.nxticks is not None: plobj.nxticks = self.nxticks
                     if self.nyticks is not None: plobj.nyticks = self.nyticks
+                    if self.fmt is not None: plobj.fmt = self.fmt
                     # -- 1D specific
                     if dp == 1:
@@ -1210 +1220 @@
             except:
                 try: self.p[iii].modx = opt.modx[0]
+                except: pass
+            ###
+            try: self.p[iii].fmt = opt.fmt[iii]
+            except:
+                try: self.p[iii].fmt = opt.fmt[0]
                 except: pass
 

trunk/UTIL/PYTHON/planetoplot_v2/ppplot.py

@@ -815 +815 @@
                 objC2 = m.contour(x, y, what_I_contour, \
                             zelevelsc, colors = ccol, linewidths = cline)
-                #mpl.clabel(objC2, inline=1, fontsize=10)
+                #mpl.clabel(objC2, inline=1, fontsize=10,manual=True,fmt='-%2.0f$^{\circ}$C',colors='r')
             m.contourf(x, y, what_I_plot, zelevels, cmap = palette, alpha = self.trans)
         ############################################################################################