Changeset 238 for trunk/MESOSCALE_DEV/PLOT

Timestamp:

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

Author:

aslmd

Message:

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

Location:

trunk/MESOSCALE_DEV/PLOT/PYTHON

Files:

• mylib/myplot.py (modified) (3 diffs)
• scripts/winds.py (modified) (6 diffs)

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

@@ -400 +400 @@
     if char in ["cyl","lcc","merc","nsper","laea"]:   step = findstep(wlon)
     else:                                             step = 10.
+    steplon = step*2.
+    #if back in ["geolocal"]:                          
+    #    step = np.min([5.,step])
+    #    steplon = step
     print step
-    m.drawmeridians(np.r_[-180.:180.:step*2.], labels=[0,0,0,1], color='grey', fontsize=fontsizemer)
+    m.drawmeridians(np.r_[-180.:180.:steplon], 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)
     if back: m.warpimage(marsmap(back),scale=0.75)
@@ -494 +498 @@
              "HFX":          "RdYlBu",\
              "ICETOT":       "YlGnBu",\
-             "TAU_ICE":      "YlGnBu",\
+             "TAU_ICE":      "Blues",\
                      }
     if whichone not in whichcolorb:
@@ -535 +539 @@
                 "mola":        domain+"mars-mola-2k.jpg",\
                 "molabw":      domain+"MarsElevation_2500x1250.jpg",\
+                "clouds":      "http://www.johnstonsarchive.net/spaceart/marswcloudmap.jpg",\
+                "jupiter":     "http://www.mmedia.is/~bjj/data/jupiter_css/jupiter_css.jpg",\
+                "jupiter_voy": "http://www.mmedia.is/~bjj/data/jupiter/jupiter_vgr2.jpg",\
                         }
+        ### see http://www.mmedia.is/~bjj/planetary_maps.html
         if whichone not in whichlink: 
                 print "marsmap: choice not defined... you'll get the default one... "

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