Index: /trunk/UTIL/PYTHON/myplot.py
===================================================================
--- /trunk/UTIL/PYTHON/myplot.py	(revision 821)
+++ /trunk/UTIL/PYTHON/myplot.py	(revision 822)
@@ -1041,4 +1041,6 @@
              "TK":           "%.0f",\
              "T":            "%.0f",\
+             "MARS_TI":      "%.0f",\
+          "THERMAL_INERTIA": "%.0f",\
              #"ZMAX_TH":      "%.0f",\
              #"WSTAR":        "%.0f",\
Index: /trunk/UTIL/PYTHON/planetoplot.py
===================================================================
--- /trunk/UTIL/PYTHON/planetoplot.py	(revision 821)
+++ /trunk/UTIL/PYTHON/planetoplot.py	(revision 822)
@@ -248,4 +248,10 @@
               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,5 +636,5 @@
                         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
Index: /trunk/UTIL/PYTHON/pp.py
===================================================================
--- /trunk/UTIL/PYTHON/pp.py	(revision 821)
+++ /trunk/UTIL/PYTHON/pp.py	(revision 822)
@@ -83,5 +83,4 @@
             zelevel = float(inputnvert[0])
             ze_interp_levels = [-9999.]
-            zesvert = readslices([str(zelevel)])
         elif np.array(inputnvert).size > 2:
             zelevel = -99.
@@ -105,4 +104,5 @@
           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'