Changeset 822

Timestamp:

Oct 29, 2012, 12:19:22 PM (12 years ago)

Author:

aslmd

Message:

UTIL PYTHON. corrected a problem introduced in a previous commit. added --axtime ind to ask for time subscript instead of values.

Location:

trunk/UTIL/PYTHON

Files:

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

trunk/UTIL/PYTHON/myplot.py

@@ -1041 +1041 @@
              "TK":           "%.0f",\
              "T":            "%.0f",\
+             "MARS_TI":      "%.0f",\
+          "THERMAL_INERTIA": "%.0f",\
              #"ZMAX_TH":      "%.0f",\
              #"WSTAR":        "%.0f",\

trunk/UTIL/PYTHON/planetoplot.py

@@ -248 +248 @@
               if initime == -1: initime=input("Please type initial local time:")
               time = (initime+time*24)%24 ; print "LOCAL TIMES.... ", time
+          ### (simply ask for subscript)
+          if axtime in ["ind"]:
+              dadim = getdimfromvar(nc)
+              if   len(dadim) == 4:  time = np.arange(dadim[-4])
+              elif len(dadim) == 3:  time = np.arange(dadim[-3])
+              elif len(dadim) == 2:  time = np.arange(dadim[-2]) 
 
     ### TYPE 2 : MESOSCALE FILES
@@ -630 +636 @@
                         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":
+                        if axtime == "lt" and indextime is None:
                             ax = mpl.pyplot.gca()
                             # set ticks where your images will be

trunk/UTIL/PYTHON/pp.py

@@ -83 +83 @@
             zelevel = float(inputnvert[0])
             ze_interp_levels = [-9999.]
-            zesvert = readslices([str(zelevel)])
         elif np.array(inputnvert).size > 2:
             zelevel = -99.
@@ -105 +104 @@
           if typefile in ["meso"]:
             if zelevel == 0. and opt.itp == 4:  zelevel = 0.010
+            if np.array(inputnvert).size == 1:  zesvert = readslices([str(zelevel)])
             ### winds or no winds
             if opt.winds            :  zefields = 'uvmet'