Changeset 874

Timestamp:

Feb 11, 2013, 9:04:23 AM (12 years ago)

Author:

aslmd

Message:

UTIL PYTHON. planetoplot, better-sized colorbars

Location:

trunk/UTIL/PYTHON

Files:

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

trunk/UTIL/PYTHON/myplot.py

@@ -1138 +1138 @@
              "W":            "jet",\
              "WMAX_TH":      "spectral",\
-             "ANOMALY":      "RdBu_r",\
+             "ANOMALY":      "gist_ncar",\
+             #"ANOMALY":      "RdBu_r",\
              "QSURFICE":     "hot_r",\
              "ALBBARE":      "spectral",\

trunk/UTIL/PYTHON/planetoplot.py

@@ -757 +757 @@
                             elif mult != 1:                                        daformat = "%.1f"
                             else:                                                  daformat = fmtvar(fvar.upper())
-                            #if proj in ['moll','cyl']:  zeorientation="horizontal" ; zepad = 0.05
-                            if proj in ['moll']:        zeorientation="horizontal" ; zepad = 0.05
-                            else:                       zeorientation="vertical" ; zepad = 0.03
-                            zecb = colorbar( fraction=0.05,pad=zepad,format=daformat,orientation=zeorientation,\
+                            if proj in ['moll']:  zeorientation="horizontal" ; zepad = 0.07 ; zefrac = 0.1 #zepad=0.05
+                            elif proj in ['cyl']: zeorientation="vertical" ; zepad = 0.03 ; zefrac = 0.023
+                            else:                 zeorientation="vertical" ; zepad = 0.03 ; zefrac = 0.05
+                            if mapmode == 0:      zefrac = 0.1
+                            zecb = colorbar( fraction=zefrac,pad=zepad,format=daformat,orientation=zeorientation,\
                                       ticks=np.linspace(zevmin,zevmax,num=min([ticks/2+1,21])),extend='neither',spacing='proportional' ) 
+                            #if zeorientation == "horizontal": 
+                            #    daticks = zecb.ax.get_xticks()
+                            #    zecb.ax.set_xticklabels(daticks,rotation='vertical')
                             if zeorientation == "horizontal" and zetitle[0] != "fill": zecb.ax.set_xlabel(zetitle[index_f]) ; zetitle[index_f]=""
                         if winds:
@@ -807 +811 @@
                         if var2 == 'HGT':        zelevels = np.arange(-10000.,30000.,250.) #1000.)
                         elif var2 == 'tpot':     zelevels = np.arange(270,370,5)
-                        elif var2 == 'tk':       zelevels = np.arange(150,250,5)
+                        elif var2 in ['tk','temp']:     zelevels = np.arange(0.,1000.,5.)#(150,250,5)
                         elif var2 == 'wstar':    zelevels = np.arange(0,10,1.0)
                         elif var2 == 'zmax_th':  zelevels = np.arange(0,10,2.0) ; what_I_plot_frame = what_I_plot_frame / 1000.
@@ -816 +820 @@
                             ## is this needed? only if len(all_var2[index_f].shape) != len(all_var[index_f].shape)
                             cs = contour( x,y,what_I_plot_frame, zelevels, colors='k', linewidths = 0.33)#, alpha=0.5, linestyles='solid')
-                            #cs = contour( x,y,what_I_plot_frame, zelevels, colors='w', linewidths = 0.5)#, alpha=0.5, linestyles='solid')
-                            clabel(cs, inline=1, fontsize = 4.*rcParams['font.size']/5., fmt=fmtvar(var2.upper()))
+                            ##cs = contour( x,y,what_I_plot_frame, zelevels, colors='w', linewidths = 0.5)#, alpha=0.5, linestyles='solid')
+                            #clabel(cs, inline=1, fontsize = 4.*rcParams['font.size']/5., fmt=fmtvar(var2.upper()))
                         elif mapmode == 1:  
                             cs = m.contour( x,y,what_I_plot_frame, zelevels, colors='k', linewidths = 0.33)#, alpha=0.5, linestyles='solid')