Changeset 207 for trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib

Timestamp:

Jul 12, 2011, 11:51:13 AM (15 years ago)

Author:

aslmd

Message:

MESOSCALE: A GENERAL CLEAN-UP FOLLOWING UPDATING THE USER MANUAL. EVERYTHING ESSENTIAL IS IN MESOSCALE (much lighter than before). EVERYTHING FOR DEVELOPPERS OR EXPERTS IS IN MESOSCALE_DEV.

Location:

trunk/MESOSCALE_DEV

Files:

• . (added)
• PLOT (copied) (copied from trunk/MESOSCALE/PLOT)
• PLOT/PYTHON/mylib/api_wrapper.py (modified) (2 diffs)
• PLOT/PYTHON/mylib/myplot.py (copied) (copied from trunk/MESOSCALE/PLOT/PYTHON/mylib/myplot.py) (3 diffs)

trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/api_wrapper.py

@@ -14 +14 @@
                     extrapolate     = 0, \
                     unstagger_grid  = False, \
-                    onelevel        = 0.020  ):
+                    onelevel        = 0.020, \
+                    nocall          = False ):
     import api
     import numpy as np
@@ -27 +28 @@
     print input_name, output_name
 
-    api.api_main ( path_to_input, input_name, path_to_output, 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), \
                    interp_method, extrapolate, unstagger_grid, onelevel )

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

@@ -68 +68 @@
     return prefix
 
-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')
@@ -240 +218 @@
     elif wlat[0] <= -80.:  blat = -40. 
     else:                  blat = wlat[0]
-    h = 2000.
+    h = 50.  ## en km
     radius = 3397200.
     if   char == "cyl":     m = Basemap(rsphere=radius,projection='cyl',\
@@ -250 +228 @@
     elif char == "npstere": m = Basemap(rsphere=radius,projection='npstere', boundinglat=blat, lon_0=0.)
     elif char == "spstere": m = Basemap(rsphere=radius,projection='spstere', boundinglat=blat, lon_0=0.)
+    elif char == "nplaea":  m = Basemap(rsphere=radius,projection='nplaea', boundinglat=wlat[0], lon_0=meanlon)
+    elif char == "laea":    m = Basemap(rsphere=radius,projection='laea',lon_0=meanlon,lat_0=meanlat,lat_ts=meanlat,\
+                              llcrnrlat=wlat[0],urcrnrlat=wlat[1],llcrnrlon=wlon[0],urcrnrlon=wlon[1])
     elif char == "nsper":   m = Basemap(rsphere=radius,projection='nsper',lon_0=meanlon,lat_0=meanlat,satellite_height=h*1000.)
     elif char == "merc":    m = Basemap(rsphere=radius,projection='merc',lat_ts=0.,\
                               llcrnrlat=wlat[0],urcrnrlat=wlat[1],llcrnrlon=wlon[0],urcrnrlon=wlon[1])
     fontsizemer = int(mpl.rcParams['font.size']*3./4.)
-    if char in ["cyl","lcc","merc"]:   step = findstep(wlon)
-    else:                              step = 10.
+    if char in ["cyl","lcc","merc","nsper","laea"]:   step = findstep(wlon)
+    else:                                             step = 10.
     m.drawmeridians(np.r_[-180.:180.:step*2.], labels=[0,0,0,1], color='grey', fontsize=fontsizemer)
     m.drawparallels(np.r_[-90.:90.:step], labels=[1,0,0,0], color='grey', fontsize=fontsizemer)