Changeset 186 for trunk/MESOSCALE/PLOT/PYTHON/mylib

Timestamp:

Jul 3, 2011, 4:27:24 AM (14 years ago)

Author:

aslmd

Message:

MESOSCALE: python post-processing. wrapper with my Fortran interpolator (small modifications done to api.F90). added the corresponding options to winds.py which is now a pretty complete script

Location:

trunk/MESOSCALE/PLOT/PYTHON/mylib

Files:

• api_g95.sh (added)
• api_wrapper.py (added)
• myplot.py (modified) (3 diffs)

trunk/MESOSCALE/PLOT/PYTHON/mylib/myplot.py

@@ -35 +35 @@
         return wlon,wlat
 
+def api_onelevel (  path_to_input   = None, \
+                    input_name      = 'wrfout_d0?_????-??-??_??:00:00', \
+                    path_to_output  = None, \
+                    output_name     = 'output.nc', \
+                    process         = 'list', \
+                    fields          = 'tk,W,uvmet,HGT', \
+                    debug           = False, \
+                    bit64           = False, \
+                    oldvar          = True, \
+                    interp_method   = 4, \
+                    extrapolate     = 0, \
+                    unstagger_grid  = False, \
+                    onelevel        = 0.020  ):
+    import api
+    import numpy as np
+    if not path_to_input:   path_to_input  = './'
+    if not path_to_output:  path_to_output = path_to_input
+    api.api_main ( path_to_input, input_name, path_to_output, output_name, \
+                   process, fields, debug, bit64, oldvar, np.arange (299), \
+                   interp_method, extrapolate, unstagger_grid, onelevel )
+    return
+
 def getproj (nc):
     map_proj = getattr(nc, 'MAP_PROJ')
@@ -51 +73 @@
         print "mercator projection"
         proj="merc"
+    else:
+        proj="merc"
     return proj    
 
@@ -70 +94 @@
     import  matplotlib.pyplot as plt
     res = int(res)
-    if folder != '':      name = folder+'/'+filename+str(res)+".png"
-    else:                 name = filename+str(res)+".png"
+    name = filename+"_"+str(res)+".png"
+    if folder != '':      name = folder+'/'+name
     plt.savefig(name,dpi=res,bbox_inches='tight',pad_inches=pad_inches_value)
     if disp:              display(name)