Changeset 717

Timestamp:

Jul 4, 2012, 12:46:05 PM (12 years ago)

Author:

acolaitis

Message:

Python. Bug correction for 1D plots. This bug was leading to completely wrong profiles when asking for a unidim plot from a testphys1d output file.

Location:

trunk/UTIL/PYTHON

Files:

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

trunk/UTIL/PYTHON/myplot.py

@@ -204 +204 @@
 
 ## Author: AS + TN + AC
-def reducefield (input,d4=None,d3=None,d2=None,d1=None,yint=False,alt=None,anomaly=False,redope=None,mesharea=None):
+def reducefield (input,d4=None,d3=None,d2=None,d1=None,yint=False,alt=None,anomaly=False,redope=None,mesharea=None,unidim=999):
     ### 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
@@ -234 +234 @@
     ### now the main part
     if dimension == 2:
+        #### this is needed for 1d-type files (where dim=2 but axes are time-vert and not lat-lon)
+        if unidim == 1: d2=d4 ; d1=d3 ; d4=None ; d3=None
         if mesharea is None: mesharea=np.ones(shape)
         if   max(d2) >= shape[0]: error = True

trunk/UTIL/PYTHON/planetoplot.py

@@ -418 +418 @@
        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, redope=redope, mesharea=area )
+                                             yint=yintegral, alt=vert, anomaly=anomaly, redope=redope, mesharea=area, unidim=is1d)
            if mult != 2718.:  what_I_plot = what_I_plot*mult 
            else:              what_I_plot = np.log10(what_I_plot) ; print "log plot"