Changeset 516 for trunk/UTIL/PYTHON

Timestamp:

Feb 4, 2012, 1:16:59 PM (13 years ago)

Author:

aslmd

Message:

UTIL PYTHON : a few stupid tricks added. if you ask --mult 2718 you'll have a log(field). if you ask moll projection the colorbar will be horizontal and title under the colorbar. if you use -T (tile) with a 1D plot you will have individual points instead of a line.

Location:

trunk/UTIL/PYTHON

Files:

• myplot.py (modified) (2 diffs)
• myscript.py (modified) (1 diff)
• planetoplot.py (modified) (3 diffs)

trunk/UTIL/PYTHON/myplot.py

@@ -623 +623 @@
     #steplon = 5.625
     #zelatmax = 89.9
-    m.drawmeridians(np.r_[-180.:180.:steplon], labels=[0,0,0,1], color=zecolor, linewidth=zelinewidth, fontsize=fontsizemer, latmax=zelatmax)
-    m.drawparallels(np.r_[-90.:90.:step], labels=[1,0,0,0], color=zecolor, linewidth=zelinewidth, fontsize=fontsizemer, latmax=zelatmax)
+    if char not in ["moll"]:
+        m.drawmeridians(np.r_[-180.:180.:steplon], labels=[0,0,0,1], color=zecolor, linewidth=zelinewidth, fontsize=fontsizemer, latmax=zelatmax)
+        m.drawparallels(np.r_[-90.:90.:step], labels=[1,0,0,0], color=zecolor, linewidth=zelinewidth, fontsize=fontsizemer, latmax=zelatmax)
     if back: m.warpimage(marsmap(back),scale=0.75)
             #if not back:
@@ -726 +727 @@
              "DOWNDRAFT":    "%.0f",\
              "TK":           "%.0f",\
+             #"ZMAX_TH":      "%.0f",\
+             #"WSTAR":        "%.0f",\
              # Variables from TES ncdf format
              "T_NADIR_DAY":  "%.0f",\

trunk/UTIL/PYTHON/myscript.py

@@ -15 +15 @@
     parser.add_option('-w', '--with',   action='store',dest='var2',      type="string",  default=None,  help='variable contoured')
     parser.add_option('-a', '--anomaly',action='store_true',dest='anomaly',              default=False, help='compute and plot relative anomaly in %')
-    parser.add_option('--mult',         action='store',dest='mult',      type="float",   default=1.,    help='a multiplicative factor to plotted field')
+    parser.add_option('--mult',         action='store',dest='mult',      type="float",   default=1.,    help='multiplicative factor to plotted field (2718=log)')
     parser.add_option('-m', '--min',    action='append',dest='vmin',     type="float",   default=None,  help='bounding minimum value [min]')    
     parser.add_option('-M', '--max',    action='append',dest='vmax',     type="float",   default=None,  help='bounding maximum value [max]') 

trunk/UTIL/PYTHON/planetoplot.py

@@ -342 +342 @@
            what_I_plot, error = reducefield( all_var[index_f], d4=indextime, d1=indexlon, d2=indexlat, d3=indexvert, \
                                              yint=yintegral, alt=vert, anomaly=anomaly, redope=redope )
-           what_I_plot = what_I_plot*mult
+           if mult != 2718.:  what_I_plot = what_I_plot*mult 
+           else:              what_I_plot = np.log10(what_I_plot) ; print "log plot"
        if var2:      ### what is contoured.
            what_I_plot_contour, error = reducefield( all_var2[index_f], d4=indextime, d1=indexlon, d2=indexlat , d3=indexvert, \
@@ -446 +447 @@
                         if mrate is not None: x = y  ## because swapaxes...
                         #what_I_plot_frame = np.diff(what_I_plot_frame, n=1) ; x = x[1:]
-                        if indexvert is not None or indextime is None:    plot(x,what_I_plot_frame,label=lbl)  ## regular plot
-                        else:                                             plot(what_I_plot_frame,x,label=lbl)  ## vertical profile
+                        
+                        if not tile:  zeline='-'
+                        else:         zeline=','
+                        if indexvert is not None or indextime is None:    plot(x,what_I_plot_frame,zeline,label=lbl)  ## regular plot
+                        else:                                             plot(what_I_plot_frame,x,zeline,label=lbl)  ## vertical profile
                         if nplot > 1: legend(loc='best')
                         if indextime is None and axtime is not None and xlab is None:    xlabel(axtime.upper()) ## define the right label
@@ -469 +473 @@
                             elif mult != 1:                                        daformat = "%.1f"
                             else:                                                  daformat = fmtvar(fvar.upper())
-                            colorbar( fraction=0.05,pad=0.03,format=daformat,\
+                            if proj in ['moll']:  zeorientation="horizontal"
+                            else:                 zeorientation="vertical"
+                            zecb = colorbar( fraction=0.05,pad=0.03,format=daformat,orientation=zeorientation,\
                                       ticks=np.linspace(zevmin,zevmax,num=min([ticks/2+1,21])),extend='neither',spacing='proportional' ) 
+                            if zeorientation == "horizontal":
+                                zecb.ax.set_xlabel(zetitle)
+                                zetitle=""
                         if winds:
                             if typefile in ['mesoapi','meso']: