Changeset 181 for trunk/MESOSCALE/PLOT/PYTHON/scripts

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/scripts

Files:

• domain.py (modified) (2 diffs)
• winds.py (modified) (3 diffs)

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]))