Changeset 247 for trunk/MESOSCALE_DEV/PLOT

Timestamp:

Jul 29, 2011, 3:47:26 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: minor bug fixes in execution script and graphical python routines.

Location:

trunk/MESOSCALE_DEV/PLOT/PYTHON

Files:

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

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

@@ -362 +362 @@
     meanlat = 0.5*(wlat[0]+wlat[1])
     if   wlat[0] >= 80.:   blat =  40. 
-    elif wlat[1] <= -80.:  blat = -40. 
-    else:                  blat = max([wlat[0],wlat[1]])
+    elif wlat[1] <= -80.:  blat = -40.
+    elif wlat[1] >= 0.:    blat = wlat[0] 
+    elif wlat[0] <= 0.:    blat = wlat[1]
     print blat
     h = 50.  ## en km
@@ -441 +442 @@
 
 def bounds(what_I_plot,zevmin,zevmax):
+    from mymath import max,min,mean
     ### might be convenient to add the missing value in arguments
     what_I_plot[ what_I_plot < zevmin ] = zevmin*(1. + 1.e-7)
+    print "new min ", min(what_I_plot)
     what_I_plot[ what_I_plot > 9e+35  ] = -9e+35
     what_I_plot[ what_I_plot > zevmax ] = zevmax*(1. - 1.e-7)
+    print "new max ", max(what_I_plot)
     return what_I_plot
 

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

@@ -150 +150 @@
                if not tile:
                    if not hole: what_I_plot = bounds(what_I_plot,zevmin,zevmax)
-                   zelevels = np.linspace(zevmin,zevmax) #,num=20)
+                   zelevels = np.linspace(zevmin*(1. + 1.e-7),zevmax*(1. - 1.e-7)) #,num=20)
                    contourf( x, y, what_I_plot, zelevels, cmap = palette )
                else:
@@ -177 +177 @@
                vectorfield(vecx, vecy,\
                           x, y, stride=stride, csmooth=2,\
-                          scale=15., factor=300., color=colorvec, key=key)
+                          #scale=15., factor=300., color=colorvec, key=key)
+                          scale=20., factor=250., color=colorvec, key=key)
                                             #200.         ## or csmooth=stride