Changeset 483

Timestamp:

Dec 19, 2011, 9:57:26 PM (13 years ago)

Author:

aslmd

Message:

PYTHON UTIL: added a new keyword --redope which allow to plot the minimum/maximum value over times stored in the file. MESOSCALE: provided a small fix in runmeso for case mars=12.

Location:

trunk

Files:

• MESOSCALE/LMD_MM_MARS/SIMU/runmeso (modified) (1 diff)
• UTIL/PYTHON/README.PP (modified) (1 diff)
• UTIL/PYTHON/myplot.py (modified) (3 diffs)
• UTIL/PYTHON/myscript.py (modified) (3 diffs)
• UTIL/PYTHON/planetoplot.py (modified) (4 diffs)
• UTIL/PYTHON/pp.py (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

@@ -195 +195 @@
        3) scat=1 ;;
       11) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
+      12) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
       20) scat=1 ;;
        *) scat=1 ;;

trunk/UTIL/PYTHON/README.PP

@@ -71 +71 @@
 Goal: I want to plot averaged results in the file from one time to another time
 pp.py -f diagfi.nc -v tsurf --time 4,7
+
+Goal: I want to plot the minimum/maximum value over times stored in the file
+pp.py -f diagfi.nc -v tsurf --redope mint
+pp.py -f diagfi.nc -v tsurf --redope maxt
 
 Goal: I want to plot a globally-averaged 1D temperature profile

trunk/UTIL/PYTHON/myplot.py

@@ -44 +44 @@
       elif 'vert' in nc.variables:             typefile = 'mesoapi'
       elif 'U' in nc.variables:                typefile = 'meso'
+      #elif 'HGT_M' in nc.variables:            typefile = 'geo'
+      elif 'HGT' in nc.variables:              typefile = 'meso'
     elif 'HGT_M' in nc.variables:              typefile = 'geo'
-    #else:                            errormess("whatkindfile: typefile not supported.")
     else:                                      typefile = 'gcm' # for lslin-ed files from gcm
     return typefile
@@ -75 +76 @@
 
 ## Author: AS + TN + AC
-def reducefield (input,d4=None,d3=None,d2=None,d1=None,yint=False,alt=None,anomaly=False):
+def reducefield (input,d4=None,d3=None,d2=None,d1=None,yint=False,alt=None,anomaly=False,redope=None):
     ### we do it the reverse way to be compliant with netcdf "t z y x" or "t y x" or "y x"
     ### it would be actually better to name d4 d3 d2 d1 as t z y x
     ### ... note, anomaly is only computed over d1 and d2 for the moment
     import numpy as np
-    from mymath import max,mean
+    from mymath import max,mean,min
     csmooth = 12 ## a fair amount of grid points (too high results in high computation time)
+    if redope is not None:
+       if   redope == "mint":     input = min(input,axis=0) ; d1 = None
+       elif redope == "maxt":     input = max(input,axis=0) ; d1 = None
+       else:                      errormess("not supported. but try lines in reducefield beforehand.")
+       #elif redope == "minz":     input = min(input,axis=1) ; d2 = None
+       #elif redope == "maxz":     input = max(input,axis=1) ; d2 = None
+       #elif redope == "miny":     input = min(input,axis=2) ; d3 = None
+       #elif redope == "maxy":     input = max(input,axis=2) ; d3 = None
+       #elif redope == "minx":     input = min(input,axis=3) ; d4 = None
+       #elif redope == "maxx":     input = max(input,axis=3) ; d4 = None
     dimension = np.array(input).ndim
     shape = np.array(input).shape
@@ -201 +212 @@
         subh = 99999
     elif numplot == 2:
-        subv = 2
-        subh = 1
+        subv = 1 #2
+        subh = 2 #1
         fig.subplots_adjust(wspace = 0.35)
         rcParams['font.size'] = int( rcParams['font.size'] * 3. / 4. )

trunk/UTIL/PYTHON/myscript.py

@@ -7 +7 @@
     parser.add_option('-d', '--display',action='store_false',dest='display',             default=True,  help='do not pop up created images')
     parser.add_option('-O','--output',  action='store',dest='out',       type="string",  default=None,  help='output file name')
+    parser.add_option('--rate',         action='store'     ,dest='rate', type="int",   default=None, help='output is a movie along Time dimension [None]')
+    parser.add_option('--quality',      action='store_true',dest='quality',            default=False,help='For movie mode: improves movie quality.(slower)')
 
     ### WHAT I PLOT
@@ -16 +18 @@
     parser.add_option('-M', '--max',    action='append',dest='vmax',     type="float",   default=None,  help='bounding maximum value [max]') 
     parser.add_option('-H', '--hole',   action='store_true',dest='hole',                 default=False, help='holes above max and below min')
-    parser.add_option('--nolow',        action='store_true',dest='nolow',                default=False, help='plot only minima and maxima [false]')
+    parser.add_option('--nolow',        action='store_true',dest='nolow',                default=False, help='do not plot low |values| [False]')
+    parser.add_option('--redope',       action='store',dest='redope',    type="string",  default=None,  help='REDuce OPErators: mint,maxt for the moment [None]')
 
     ### VERTICAL INTERPOLATION
@@ -65 +68 @@
     ### SPECIAL
     parser.add_option('--tsat',         action='store_true',dest='tsat',               default=False,help='convert temperature field T in Tsat-T using pressure')
-    parser.add_option('--rate',         action='store'     ,dest='rate', type="int",   default=None, help='output is a movie along Time dimension [None]')
-    parser.add_option('--quality',      action='store_true',dest='quality',            default=False,help='For movie mode: improves movie quality.(slower)')
 
     return parser

trunk/UTIL/PYTHON/planetoplot.py

@@ -58 +58 @@
            trans=1,\
            zarea=None,\
-           axtime=None):
+           axtime=None,\
+           redope=None):
 
 
@@ -181 +182 @@
               if "Times" in nc.variables:   time = count + np.arange(0,len(nc.variables["Times"]),1)
               elif "Time" in nc.variables:  time = count + np.arange(0,len(nc.variables["Time"]),1)
+              else:                         time = count + np.arange(0,1,1)
               count = time[-1] + 1  ## so that a cat is possible with simple subscripts
           ###
@@ -303 +305 @@
        if varname:   ### what is shaded.
            what_I_plot, error = reducefield( all_var[index_f], d4=indextime, d1=indexlon, d2=indexlat, d3=indexvert, \
-                                             yint=yintegral, alt=vert, anomaly=anomaly )
+                                             yint=yintegral, alt=vert, anomaly=anomaly, redope=redope )
            what_I_plot = what_I_plot*mult
        if var2:      ### what is contoured.
@@ -473 +475 @@
 
        ### Next subplot
-       basename = getname(var=varname,var2=var2,winds=winds,anomaly=anomaly)
+       zevarname = varname
+       if redope is not None: zevarname = zevarname + "_" + redope 
+       basename = getname(var=zevarname,var2=var2,winds=winds,anomaly=anomaly)
        if len(what_I_plot.shape) > 3:
            basename = basename + getstralt(nc,level) 

trunk/UTIL/PYTHON/pp.py

@@ -152 +152 @@
                 invert_y=opt.inverty,xaxis=zexaxis,yaxis=zeyaxis,ylog=opt.logy,yintegral=opt.column,\
                 blat=opt.blat,blon=opt.blon,tsat=opt.tsat,flagnolow=opt.nolow,\
-                mrate=opt.rate,mquality=opt.quality,trans=opt.trans,zarea=opt.area,axtime=opt.axtime)
+                mrate=opt.rate,mquality=opt.quality,trans=opt.trans,zarea=opt.area,axtime=opt.axtime,\
+                redope=opt.redope)
         print 'DONE: '+name
         system("rm -f to_be_erased")