Changeset 516 for trunk/UTIL/PYTHON/planetoplot.py

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.

File:

• trunk/UTIL/PYTHON/planetoplot.py (modified) (3 diffs)

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']: