Changeset 2465 in lmdz_wrf

Timestamp:

Apr 24, 2019, 2:38:00 PM (6 years ago)

Author:

lfita

Message:

Update `varinfile'

File:

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

trunk/tools/nc_var_tools.py

@@ -299 +299 @@
     varn= variable name
     """
+    fname = 'varinfile'
     if not ncf.variables.has_key(varn):
         print emsg
-        print '  file "' + filen + '" does not have ' + vmsg + ' variable "' + varn + '" !!'
-        print emsg
+        print '  file "'+filen+'" does not have '+vmsg+ ' variable "' + varn + '" !!'
         ncf.close()
-        quit(-1)    
+        Varns = list(ncf.variables.keys())
+        Varns.sort()
+        print '    avaialble ones:', Varns
+        quit(-1)
 
     return
@@ -31284 +31287 @@
 #compress_data(values, '/home/lluis/PY/wrfout_d01_1995-01-01_00:00:00', 'T2,LU_INDEX')
 
-
 def CFfile_creation(values, ncfile, variables):
     """ Operation to create a file folowing CF-conventions 
@@ -31500 +31502 @@
 #CFfile_creation(values, 'CFtest.nc', varsS)
 
+def CFfile_fixTime(values, ncfile, variable):
+    """ Operation to fix a netCDF file with a time axis folowing CF-conventions, but
+        with a wrong set-up
+      values= [operations]
+        [operations]: ';' list of names of operations to perform
+          'setRefDate',[YYYYMMDDHHMMSS]: re-set the reference date to [YYYYMMDDHHMMSS]
+          'setCalendar': re-set the calendar to gregorian or standard
+      ncfile= name of the file to fix its time-axis
+      variable= name of the variable time 
+    """
+    fname = 'CFfile_fixTime'
+
+    if values == 'h':
+        print fname + '_____________________________________________________________'
+        print CFfile_fixTime.__doc__
+        quit()
+
+    expectargs = '[operations]'
+ 
+    gen.check_arguments(fname,values,expectargs,':')
+
+    operations = values.split(':')[0]
+
+    onc = NetCDFFile(ncfile,'a')
+    if not onc.variables.has_key(variable):
+
+
+    otime = onc.variables[variable]
+
+    for op in operations:
+
+