Changeset 610 for trunk/UTIL/PYTHON/pp.py

Timestamp:

Mar 31, 2012, 7:33:41 PM (13 years ago)

Author:

aslmd

Message:

UTIL PYTHON : for mesoscale plots 1. improved the case of pressure interpolation -i 1 or 2 ; 2. fixed the mean capabilities e.g. --lat 15.,35. now working fine the problem was the wrong dimension being considered in slon, my bad.

File:

• trunk/UTIL/PYTHON/pp.py (modified) (1 diff)

trunk/UTIL/PYTHON/pp.py

@@ -79 +79 @@
             zelevel = float(inputnvert[0])
             ze_interp_levels = [-9999.]
-        elif np.array(inputnvert).size == 2:
+        elif np.array(inputnvert).size > 2:
             zelevel = -99.
-            ze_interp_levels = np.linspace(float(inputnvert[0]),float(inputnvert[1]),20)
-        else:
-            zelevel = -99.
-            ze_interp_levels = np.linspace(float(inputnvert[0]),float(inputnvert[1]),float(inputnvert[2]))
+            start = float(inputnvert[0])
+            stop = float(inputnvert[1])
+            if np.array(inputnvert).size == 2:  numsample = 20
+            else:                               numsample = float(inputnvert[2])
+            if stop > start:   
+               ## altitude coordinates
+               ze_interp_levels = np.linspace(start,stop,numsample)
+            else:
+               ## pressure coordinates
+               ze_interp_levels = np.logspace(np.log10(start),np.log10(stop),numsample)
 
         #########################################################