Index: trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/makefig.sh
===================================================================
--- trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/makefig.sh	(revision 233)
+++ trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/makefig.sh	(revision 237)
@@ -4,6 +4,6 @@
  file1="$fold/POLAR_THOMAS_ls50/wrfout_d03_2024-02-45_06:00:00"
  file2="$fold/POLAR_THOMAS_ls20/wrfout_d03_2024-01-42_06:00:00"
- file3="$fold/POLAR_THOMAS_ls0/wrfout_d03_2024-01-06_06:00:00"
- file4="$fold/POLAR_THOMAS_ls10/wrfout_d03_2024-01-19_06:00:00"
+ #file3="$fold/POLAR_THOMAS_ls0/wrfout_d03_2024-01-06_06:00:00"
+ file3="$fold/POLAR_THOMAS_ls10/wrfout_d03_2024-01-19_06:00:00"
  dest="/u/aslmd/WWW/antichambre/thomas/"
 
@@ -13,4 +13,4 @@
            -v USTM -m   0. -M 0.8 \
            -v HGT  -m   0  -M 0 \
-           -f $file1 -f $file2 -f $file3 -f $file4 \
+           -f $file1 -f $file2 -f $file3 \
            -d
Index: trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py
===================================================================
--- trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py	(revision 233)
+++ trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/winds.py	(revision 237)
@@ -33,5 +33,5 @@
     ### Load librairies and functions
     from netCDF4 import Dataset
-    from myplot import getcoord2d,define_proj,makeplotpng,simplinterv,vectorfield,ptitle,latinterv,getproj,wrfinterv,dumpbdy,\
+    from myplot import getcoord2d,define_proj,makeplotpng,makeplotpngres,simplinterv,vectorfield,ptitle,latinterv,getproj,wrfinterv,dumpbdy,\
                        fmtvar,definecolorvec,defcolorb,getprefix,putpoints,calculate_bounds,errormess,definesubplot,\
                        zoomset,getcoorddef,getwinddef,whatkindfile,reducefield,bounds,getstralt,getfield
@@ -87,5 +87,5 @@
 
        ### Which local time ?
-       print interv[0], interv[1], itime
+       #print interv[0], interv[1], itime
        ltst = ( interv[0] + 0.5*(wlon[0]+wlon[1])/15.) + itime*interv[1]
        ltst = int (ltst * 10) / 10.
@@ -122,6 +122,5 @@
                palette = get_cmap(name=colorb)
                if not tile:
-                   hole = True
-                   hole = None
+                   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)
                    zelevels = np.linspace(zevmin,zevmax)
@@ -178,5 +177,7 @@
     else:               zeplot = target + "/" + zeplot  
     ###
-    if found_lct:     makeplotpng(zeplot,pad_inches_value=0.35,disp=display)   
+    #if found_lct:     makeplotpng(zeplot,pad_inches_value=0.35,disp=display)   
+    if found_lct:     makeplotpng(zeplot,disp=display)
+    #if found_lct:     makeplotpngres(zeplot,200,disp=display)
     else:             print "Local time not found"