Changeset 425

Timestamp:

Nov 25, 2011, 5:06:34 PM (14 years ago)

Author:

aslmd

Message:

GRAPHICS: cosmetic changes. corrected a bug in treating local times with two input meso files.

Location:

trunk/UTIL/PYTHON

Files:

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

trunk/UTIL/PYTHON/myplot.py

@@ -678 +678 @@
     fmtvar    =     { \
              "TK":           "%.0f",\
-# Variables from TES ncdf format
+             # Variables from TES ncdf format
              "T_NADIR_DAY":  "%.0f",\
              "T_NADIR_NIT":  "%.0f",\
-# Variables from tes.py ncdf format
+             # Variables from tes.py ncdf format
              "TEMP_DAY":     "%.0f",\
              "TEMP_NIGHT":   "%.0f",\
-# Variables from MCS and mcs.py ncdf format
+             # Variables from MCS and mcs.py ncdf format
              "DTEMP":     "%.0f",\
              "NTEMP":   "%.0f",\
+             # other stuff
              "TPOT":         "%.0f",\
              "TSURF":        "%.0f",\
@@ -726 +727 @@
              "HGT":          "spectral",\
              "TK":           "gist_heat",\
+             "TPOT":         "Paired",\
              "TSURF":        "RdBu_r",\
              "QH2O":         "PuBu",\
@@ -927 +929 @@
 # output : all indexes to be taken into account for reducing field
   import numpy as np
+  ### added by AS to treat the case of stime = - LT
+  if saxis is not None:
+      if saxis[0][0] < 0: saxis = - saxis
+  ###
   if ( np.array(axis).ndim == 2):
       axis = axis[:,0]

trunk/UTIL/PYTHON/myscript.py

@@ -21 +21 @@
     parser.add_option('-l', '--level',  action='store',dest='lvl',       type="string",  default="0",   help='level / start,stop,step (-i 2: p,mb)(-i 3,4: z,km) [0]')
     parser.add_option('-i', '--interp', action='store',dest='itp',       type="int",     default=None,  help='interpolation (2: p, 3: z-amr, 4:z-als)')
-    parser.add_option('--intas',        action='store',dest='intas',     type="string",  default=None,  help='specify "mcs" or "tes" for specific gcm pressure interpolation grid')
+    parser.add_option('--intas',        action='store',dest='intas',     type="string",  default=None,  help='specify "mcs" or "tes" for gcm P interpolation grid')
     parser.add_option('-N', '--no-api', action='store_true',dest='nocall',               default=False, help='do not recreate api file')
 

trunk/UTIL/PYTHON/planetoplot.py

@@ -257 +257 @@
                if typefile in ['mesoapi','meso']:
                    for i in range(len(time)):  time[i] = localtime ( interv[0]+time[i]*interv[1], 0.5*(wlon[0]+wlon[1]) )
-                   stime = -stime ; print "OK... WORKING WITH LOCAL TIMES"
+                   print "OK... WORKING WITH LOCAL TIMES"
                else: errormess("local times not supported for GCM files. not too hard to modify the code though.")
        if mapmode == 1 and stime is None:   indextime = 1
        else:                                indextime = getsindex(stime,((nplot-1)//(nlon*nlat*nvert))%ntime,time)
        if typefile in ['mesoapi','meso'] and indextime is not None:  ltst = localtime ( interv[0]+indextime*interv[1], 0.5*(wlon[0]+wlon[1]) )  
+       print "********** index lon, lat, vert, time ",indexlon,indexlat,indexvert,indextime
        ####################################################################