Changeset 181

Timestamp:

Jun 29, 2011, 7:00:20 PM (14 years ago)

Author:

aslmd

Message:

MESOSCALE:PYTHON: petites ameliorations pour rendre les scripts applicables si on n'a pas les droits dans le dossier ou se trouvent les donnees.

Location:

trunk/MESOSCALE/PLOT/PYTHON

Files:

• map_ecmwf.py (modified) (1 diff)
• mylib/myplot.py (modified) (3 diffs)
• scripts/domain.py (modified) (2 diffs)
• scripts/winds.py (modified) (3 diffs)

trunk/MESOSCALE/PLOT/PYTHON/map_ecmwf.py

@@ -42 +42 @@
                         m = myp.define_proj(proj,wlon,wlat)
                         x, y = m(lon2d, lat2d)
-                        m.drawmeridians(np.r_[wlon[0]:wlon[1]:step*2], labels=[0,0,0,1], color='grey')
-                        m.drawparallels(np.r_[wlat[0]:wlat[1]:step  ], labels=[1,0,0,0], color='grey')
-                        #m.bluemarble()
-                        m.warpimage(myp.earthmap("nice"))
+                        #m.drawmeridians(np.r_[wlon[0]:wlon[1]:step*2], labels=[0,0,0,1], color='grey')
+                        #m.drawparallels(np.r_[wlat[0]:wlat[1]:step  ], labels=[1,0,0,0], color='grey')
+                        m.bluemarble()
+                        #m.warpimage(myp.earthmap("nice"))
                         zeplot = m.contour(x, y, field, 20)
                         plt.clabel(zeplot, inline=1, inline_spacing=1, fontsize=7, fmt='%0i')

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

@@ -47 +47 @@
     import  matplotlib.pyplot as plt
     res = int(res)
-    if folder != '':    name = folder+'/'+filename+str(res)+".png"
-    else:               name = filename+str(res)+".png"
+    if folder != '':      name = folder+'/'+filename+str(res)+".png"
+    else:                 name = filename+str(res)+".png"
     plt.savefig(name,dpi=res,bbox_inches='tight',pad_inches=pad_inches_value)
-    if disp:            display(name)           
+    if disp:              display(name)         
     return
 
@@ -118 +118 @@
 def display (name):
         from    os                              import system
-        system("eog "+name+" > /dev/null 2> /dev/null &")
+        system("display "+name+" > /dev/null 2> /dev/null &")
         return name
 
@@ -138 +138 @@
     meanlat = 0.5*(wlat[0]+wlat[1])
     h = 2000.
-    if   char == "cyl":     m = Basemap(projection='cyl',llcrnrlat=wlat[0],urcrnrlat=wlat[1],llcrnrlon=wlon[0],urcrnrlon=wlon[1])#,suppress_ticks=False)
+    if   char == "cyl":     m = Basemap(projection='cyl',llcrnrlat=wlat[0],urcrnrlat=wlat[1],llcrnrlon=wlon[0],urcrnrlon=wlon[1])
     elif char == "moll":    m = Basemap(projection='moll',lon_0=meanlon)
     elif char == "ortho":   m = Basemap(projection='ortho',lon_0=meanlon,lat_0=meanlat)

trunk/MESOSCALE/PLOT/PYTHON/scripts/domain.py

@@ -4 +4 @@
 
 def usage():
-    print 'USAGE:  plot.py file '
+    print 'USAGE:  plot.py file (target)'
     print 'file  : name of input netcdf file.'
-    print 'Example:  domain.py /tmp7/aslmd/HELLAS/wrfout_d02_2024-10-03_06:00:00'
+    print '(target) : a directory with write rights.'
+    print 'Example:  domain.py /d5/aslmd/LMD_MM_MARS_SIMUS/OM/OM6_TI85/wrfout_d01_2024-06-43_06:00:00 ~/'
 
-def domain (namefile,proj="ortho",back="molabw"): #vishires
+def domain (namefile,proj="ortho",back="molabw",target=None): #vishires
     from netCDF4 import Dataset
     from myplot import getcoord2d,define_proj,makeplotpng,simplinterv
@@ -18 +19 @@
     x, y = m(lon2d, lat2d)
     contourf(x, y, nc.variables['HGT'][0,:,:] / 1000., 50)
-    makeplotpng(namefile+".domain",pad_inches_value=0.05)
+    if not target:   zeplot = namefile+".domain"
+    else:            zeplot = target+"domain"
+    makeplotpng(zeplot,pad_inches_value=0.35)
 
 if __name__ == "__main__":
     import sys
-    domain(str(sys.argv[1]))
-
+    if (len(sys.argv)) == 2:     domain( str(sys.argv[1]) )
+    elif (len(sys.argv)) == 3:   domain( str(sys.argv[1]) , target = str(sys.argv[2]) )

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

@@ -4 +4 @@
 
 def usage():
-    print 'USAGE:  winds.py ntime nvert file'
+    print 'USAGE:  winds.py nvert file (target)'
     print 'file  : name of input netcdf file.'
     print 'nvert : vertical level.'
-    print 'Example: winds.py 0 /d5/aslmd/LMD_MM_MARS_SIMUS/OM/OM6_TI85/wrfout_d01_2024-06-43_06:00:00_zabg'
+    print '(target) : a directory with write rights.'
+    print 'Example: winds.py 0 /d5/aslmd/LMD_MM_MARS_SIMUS/OM/OM6_TI85/wrfout_d01_2024-06-43_06:00:00_zabg ~/'
 
-def winds (namefile,nvert,proj="cyl",back=None): 
+def winds (namefile,nvert,proj="cyl",back=None,target=None): 
     from netCDF4 import Dataset
     from myplot import getcoord2d,define_proj,makeplotpng,simplinterv,vectorfield,ptitle
@@ -31 +32 @@
                       scale=20., factor=300., color='k')
        sub += 1
-    makeplotpng(namefile+".winds"+str(nvert),pad_inches_value=0.35)   
+    if not target:   zeplot = namefile+".winds"+str(nvert)
+    else:            zeplot = target+"winds"+str(nvert)
+    makeplotpng(zeplot,pad_inches_value=0.35)   
 
 def getwinds (nc,charu='Um',charv='Vm'):
@@ -40 +43 @@
 if __name__ == "__main__":
     import sys
-    winds(str(sys.argv[2]),int(sys.argv[1]))
+    if (len(sys.argv)) == 3:     winds( str(sys.argv[2]) , int(sys.argv[1]) )
+    elif (len(sys.argv)) == 4:   winds( str(sys.argv[2]) , int(sys.argv[1]) , target=str(sys.argv[3]))