Changeset 351 for trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/api_wrapper.py

Timestamp:

Nov 7, 2011, 12:53:40 AM (13 years ago)

Author:

aslmd

Message:

PYTHON GRAPHICS: for mesoscale files section with vertically interpolated fields with alitude or pressure (through api wrapper) is now supported. the same thing might be useful to be done with zrecast for GCM.

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/api_wrapper.py (modified) (2 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/api_wrapper.py

@@ -14 +14 @@
                     extrapolate     = 0, \
                     unstagger_grid  = False, \
+                    interp_level    = [-9999.], \
                     onelevel        = 0.020, \
                     nocall          = False ):
@@ -26 +27 @@
         if interp_method == 4:    output_name = input_name+'_zabg'
 
+    if interp_level[0] == -9999.:  
+        interp_level = np.arange(299)
+    else:                 
+        zelen = len(interp_level)
+        zetemp = np.zeros(299)
+        zetemp[0:zelen] = interp_level[0:zelen]
+        zetemp[zelen] = -99999.
+        interp_level = zetemp
+        onelevel = -99999.
+
     #print input_name, output_name
 
     if nocall:     pass
     else:          api.api_main ( path_to_input, input_name, path_to_output, output_name, \
-                   process, fields, debug, bit64, oldvar, np.arange(299), \
+                   process, fields, debug, bit64, oldvar, interp_level, \
                    interp_method, extrapolate, unstagger_grid, onelevel )