Changeset 237 for trunk/MESOSCALE_DEV/PLOT/PYTHON

Timestamp:

Jul 20, 2011, 4:11:27 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: small adjustments to graphical Python routines (to make nice cloud plots)

Location:

trunk/MESOSCALE_DEV/PLOT/PYTHON

Files:

• mylib/myplot.py (modified) (3 diffs)
• scripts/makefig.sh (modified) (2 diffs)
• scripts/winds.py (modified) (4 diffs)

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

@@ -165 +165 @@
     from numpy import array
     nc  = Dataset(namefile)
+    zetime = None
     if 'Times' in nc.variables: 
         zetime = nc.variables['Times'][0]
@@ -476 +477 @@
              "HFX":          "%.0f",\
              "ICETOT":       "%.1e",\
+             "TAU_ICE":      "%.2f",\
                     }
     if whichvar not in fmtvar:
@@ -492 +494 @@
              "HFX":          "RdYlBu",\
              "ICETOT":       "YlGnBu",\
+             "TAU_ICE":      "YlGnBu",\
                      }
     if whichone not in whichcolorb:

trunk/MESOSCALE_DEV/PLOT/PYTHON/scripts/makefig.sh

@@ -4 +4 @@
  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 +13 @@
            -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

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

@@ -33 +33 @@
     ### 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 +87 @@
 
        ### 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 +122 @@
                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 +177 @@
     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"