Index: /trunk/MESOSCALE/LMD_MM_MARS/makemeso
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/makemeso	(revision 246)
+++ /trunk/MESOSCALE/LMD_MM_MARS/makemeso	(revision 247)
@@ -188,5 +188,5 @@
     if [[ "${phys}" == "newphys_" ]]
     then
-     \rm copy_model_tmp > /dev/null
+     \rm copy_model_tmp > /dev/null 2> /dev/null
         if [[ "${scenario}" == "" ]]
         then
Index: /trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/myplot.py
===================================================================
--- /trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/myplot.py	(revision 246)
+++ /trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/myplot.py	(revision 247)
@@ -362,6 +362,7 @@
     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,8 +442,11 @@
 
 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
 
Index: /trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py
===================================================================
--- /trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py	(revision 246)
+++ /trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py	(revision 247)
@@ -150,5 +150,5 @@
                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,5 +177,6 @@
                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