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

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.

File:

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

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)