Changeset 181 for trunk/MESOSCALE/PLOT/PYTHON/scripts/winds.py

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.

File:

• trunk/MESOSCALE/PLOT/PYTHON/scripts/winds.py (modified) (3 diffs)

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