Changeset 251 for trunk/MESOSCALE_DEV/PLOT

Timestamp:

Aug 2, 2011, 1:45:34 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: python graphics: improved polar plots. added a detection for dust devils + stats.

Location:

trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts

Files:

• find_devils.py (added)
• winds.py (modified) (7 diffs)

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

@@ -40 +40 @@
     from myplot import getcoord2d,define_proj,makeplotres,simplinterv,vectorfield,ptitle,latinterv,getproj,wrfinterv,dumpbdy,\
                        fmtvar,definecolorvec,defcolorb,getprefix,putpoints,calculate_bounds,errormess,definesubplot,\
-                       zoomset,getcoorddef,getwinddef,whatkindfile,reducefield,bounds,getstralt,getfield,smooth,nolow
+                       zoomset,getcoorddef,getwinddef,whatkindfile,reducefield,bounds,getstralt,getfield,smooth,nolow,\
+                       getname,localtime,polarinterv
     from mymath import deg,max,min,mean
     from matplotlib.pyplot import contour,contourf, subplot, figure, rcParams, savefig, colorbar, pcolor, show
@@ -60 +61 @@
         [uchar,vchar,metwind] = getwinddef(nc)             
         if uchar == 'not found': winds = False
+    if not var and not winds: errormess("please set at least winds or var",printvar=nc.variables)
     [lon2d,lat2d] = getcoorddef(nc)                              ## COORDINATES, could be moved below
     if proj == None:   proj = getproj(nc)                        ## PROJECTION
@@ -65 +67 @@
     ##########################
     ### Define plot boundaries
-    if proj == "npstere":             [wlon,wlat] = latinterv("North_Pole")
-    elif proj == "spstere":           [wlon,wlat] = latinterv("Far_South_Pole")
+    ### todo: possible areas in latinterv in argument (ex: "Far_South_Pole")
+    if proj in ["npstere","spstere"]: [wlon,wlat] = polarinterv(lon2d,lat2d)
     elif proj in ["lcc","laea"]:      [wlon,wlat] = wrfinterv(lon2d,lat2d)
     else:                             [wlon,wlat] = simplinterv(lon2d,lat2d)
@@ -73 +75 @@
     #########################################
     ### Name for title and graphics save file
-    if var and winds:     basename = var + '_UV'
-    elif var:             basename = var
-    elif winds:           basename = 'UV'
-    else:
-        print nc.variables                 
-        errormess("please set at least winds or var")
-    if anomaly:           basename = 'd' + basename
+    basename = getname(var=var,winds=winds,anomaly=anomaly)
     basename = basename + getstralt(nc,nvert)  ## can be moved elsewhere for a more generic routine
 
@@ -98 +94 @@
 
        ### Which local time ?
-       #print interv[0], interv[1], itime
-       ltst = ( interv[0] + 0.5*(wlon[0]+wlon[1])/15.) + itime*interv[1]
-       ltst = int (ltst * 10) / 10.
-       ltst = ltst % 24
+       ltst = localtime ( interv[0]+itime*interv[1], 0.5*(wlon[0]+wlon[1]) )
 
        ### General plot settings
@@ -109 +102 @@
            if numplot > 1:     
                if typefile not in ['geo']:  subplot(sub+nplot-1)
-
            found_lct = True
        ### If only one local time is requested (numplot < 0)
@@ -157 +149 @@
                            vmin=zevmin, vmax=zevmax )
                if colorb != 'nobar' and var != 'HGT':              
-                         colorbar(fraction=0.05,pad=0.1,format=fmtvar(fvar),\
+                         colorbar(fraction=0.05,pad=0.03,format=fmtvar(fvar),\
                                            ticks=np.linspace(zevmin,zevmax,ndiv+1),\
                                            extend='neither',spacing='proportional')