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

Timestamp:

Oct 25, 2012, 11:24:38 PM (12 years ago)

Author:

aslmd

Message:

UTIL PYTHON : improved projections. improved local time ticks. small bug fix with vertical coordinates. added a keyword nocolorb.

File:

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

trunk/UTIL/PYTHON/planetoplot.py

@@ -76 +76 @@
            trycol=False,\
            streamflag=False,\
+           nocolorb=False,\
            analysis=None):
 
@@ -429 +430 @@
     ### Open a figure and set subplots
     fig = figure()
-    subv,subh = definesubplot( numplot, fig ) 
+    subv,subh = definesubplot( numplot, fig, ipreferline = (mapmode == 1) ) 
     if ope in ['-','-%','-_histo'] and len(namefiles) ==1 : subv,subh = 2,2
     elif ope in ['-','-%'] and len(namefiles)>1 : subv, subh = len(namefiles), 3
@@ -464 +465 @@
        ltst = None 
        if typefile in ['meso'] and indextime is not None and len(indextime) < 2: ltst = localtime ( indextime, slon , all_namefile[index_f])
+
+
+       #if mapmode == 0:
+       # if xlab is None:
+       #  if indexlon is None:  xlab = "Longitude"
+       #  elif indexlat is None: xlab = "Latitude"
+       # if ylab is None:
+       #  if indexvert is None: ylab = "Altitude"
+
+
        print "********** INDEX LON:",indexlon," LAT:",indexlat," VERT:",indexvert," TIME:",indextime
        ##var = nc.variables["phisinit"][:,:]
@@ -547 +558 @@
                else:                                                   palette = get_cmap(name=colorb)
                #palette = cm.GMT_split
+               #palette = cm.GMT_globe
                ##### 1. ELIMINATE 0D or >3D CASES
                if len(what_I_plot.shape) == 0:   
@@ -618 +630 @@
                         if indextime is None and axtime is not None and xlab is None:    xlabel(axtime.upper()) ## define the right label
                         if save == 'txt':  writeascii(np.transpose([x,np.transpose(what_I_plot)]),'profile'+str(nplot*1000+imov)+'.txt')
-
+                        if axtime == "lt":
+                            ax = mpl.pyplot.gca()
+                            # set ticks where your images will be
+                            ax.get_xaxis().set_ticks(np.arange(0,48,2))
+                            # rename tick labels
+                            ax.get_xaxis().set_ticklabels(["0","2","4","6","8","10","12","14","16","18","20","22",\
+                                                       "0","2","4","6","8","10","12","14","16","18","20","22"])
+                            ## rebound everyone
+                            ax.set_xbound(lower=min(x), upper=max(x))
                     elif which == "xy":
                         if lbls is not None: lbl=lbls[index_f]
@@ -707 +727 @@
                             mpl.pyplot.plot([idx],[idy],'+k',mew=0.5,ms=18)
 
-                        if colorb not in ['nobar','onebar']:        
+                        if not nocolorb:
+                          if colorb not in ['nobar','onebar']:
                             if (fileref is not None) and ((index_f+1)%3 == 0):   daformat = "%.3f" 
                             elif mult != 1:                                        daformat = "%.1f"