Changeset 1749 in lmdz_wrf

Timestamp:

Dec 17, 2017, 4:45:04 PM (8 years ago)

Author:

lfita

Message:

Final version:

• Adding year attribution as argument
• Adding time-units

File:

• trunk/tools/ORforcing_reconstruct.py (modified) (3 diffs)

trunk/tools/ORforcing_reconstruct.py

@@ -12 +12 @@
 import subprocess as sub
 import module_ForSci as Sci
-
-filen = 'cruncep_halfdeg_1950.nc'
+from optparse import OptionParser
+
+parser = OptionParser()
+parser.add_option("-y", "--year", dest="year", help="year to process",               \
+  metavar="VALUE")
+
+(opts, args) = parser.parse_args()
+
+####### ###### ##### #### ### ## #
+
+filen = 'cruncep_halfdeg_' + opts.year + '.nc'
 
 # Variable whcih provides the indices of a 1D vector from the dimy, dimx space
@@ -53 +62 @@
 availProj = ['latlon']
 
-ofilen = 'reconstruct_matrix.nc'
+ofilen = 'reconstruct_matrix_' + opts.year + '.nc'
 
 ncvars = onc.variables.keys()
@@ -260 +269 @@
 onewnc.sync()
 ncvar.add_global_PyNCplot(onewnc, main, fname, '0.1')
-
 onc.close()
+
+# Reconstructing times
+otime = onewnc.variables['time']
+ncvar.set_attribute(otime, 'units', 'seconds since ' + opts.year + '-01-01 00:00:00')
+
 onewnc.sync()
 onewnc.close()
 
-print fname + ": Successful writing of file 'reconstruct_matrix.nc' !!"
+print fname + ": Successful writing of file '" + ofilen + ".nc' !!"