Changeset 238 for trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py

Timestamp:

Jul 20, 2011, 5:49:17 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: small adjustments to graphical plots routines, including additional maps.

File:

• trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py (modified) (6 diffs)

trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py

@@ -25 +25 @@
            zoom=None,\
            display=True,\
-           itstep=None):
+           itstep=None,\
+           hole=False):
 
     ####################################################################################################################
@@ -50 +51 @@
     typefile = whatkindfile(nc)                                  ## TYPEFILE
     if var not in nc.variables: var = False                      ## VAR
-    if winds: [uchar,vchar,metwind] = getwinddef(nc)             ## WINDS
-    if uchar == 'not found': winds = False
+    if winds:                                                    ## WINDS
+        [uchar,vchar,metwind] = getwinddef(nc)             
+        if uchar == 'not found': winds = False
     [lon2d,lat2d] = getcoorddef(nc)                              ## COORDINATES, could be moved below
     if proj == None:   proj = getproj(nc)                        ## PROJECTION
@@ -122 +124 @@
                palette = get_cmap(name=colorb)
                if not tile:
-                   if var in ["TAU_ICE","ICETOT"]: hole = True #nice plots with vis img
-                   if not hole:  what_I_plot = bounds(what_I_plot,zevmin,zevmax)
+                   if not hole: what_I_plot = bounds(what_I_plot,zevmin,zevmax)
                    zelevels = np.linspace(zevmin,zevmax)
                    contourf( x, y, what_I_plot, zelevels, cmap = palette )
@@ -177 +178 @@
     else:               zeplot = target + "/" + zeplot  
     ###
-    #if found_lct:     makeplotpng(zeplot,pad_inches_value=0.35,disp=display)   
-    if found_lct:     makeplotpng(zeplot,disp=display)
+    if found_lct:     makeplotpng(zeplot,pad_inches_value=0.35,disp=display)   
     #if found_lct:     makeplotpngres(zeplot,200,disp=display)
     else:             print "Local time not found"
@@ -232 +232 @@
     parser.add_option('-d', action='store_false',dest='display',                default=True,  help='do not pop up created images')
     parser.add_option('-e', action='store',dest='itstep',       type="int",     default=None,  help='stride time (def=4)')
+    parser.add_option('-H', action='store_true',dest='hole',                    default=False, help='holes above max and below min')
     #parser.add_option('-V', action='store', dest='comb',        type="float",   default=None,  help='a defined combination of variables')
     (opt,args) = parser.parse_args()
@@ -302 +303 @@
                 addchar=lschar,interv=[zehour,zehourin],vmin=argvmin,vmax=argvmax,\
                 tile=opt.tile,zoom=opt.zoom,display=opt.display,\
-                itstep=opt.itstep)
+                itstep=opt.itstep,hole=opt.hole)
             print 'Done: '+name