Changeset 311 in lmdz_wrf

Timestamp:

Feb 26, 2015, 11:34:12 AM (10 years ago)

Author:

lfita

Message:

Final version of the operation 'ncstepdiff'

Location:

trunk/tools

Files:

• nc_var.py (modified) (2 diffs)
• nc_var_tools.py (modified) (6 diffs)

trunk/tools/nc_var.py

@@ -24 +24 @@
   'flipdim', 'fvaradd', 'gaddattrk', 'gaddattr', 'get_namelist_vars', 'grattr',      \
   'grmattr', 'igattrs', 'isgattrs', 'isvattrs', 'ivars', 'ivattrs', 'maskvar',       \
-  'ncreplace', 'netcdf_concatenation', 'netcdf_fold_concatenation',                  \
+  'ncreplace', 'ncstepdiff', 'netcdf_concatenation', 'netcdf_fold_concatenation',    \
   'compute_opersvarsfiles',                                                          \
   'seasmean', 'sellonlatbox', 'sellonlatlevbox', 'setvar_asciivalues',               \
@@ -168 +168 @@
 elif oper == 'ncreplace':
     ncvar.ncreplace(opts.values, opts.ncfile, opts.varname)
+elif oper == 'ncstepdiff':
+    ncvar.ncstepdiff(opts.values, opts.ncfile, opts.varname)
 elif oper == 'netcdf_concatenation':
     ncvar.netcdf_concatenation(opts.ncfile)

trunk/tools/nc_var_tools.py

@@ -2648 +2648 @@
             if not searchInlist(attvar, 'units'):
                 print warnmsg
-                print '    cls_time_information: time variable "', tname,            \
+                print '    cls_time_information: time variable "' + tname +          \
                   '" does not have attribute: "units"'
                 units = None
@@ -13452 +13452 @@
     return
 
-def ncstediffp(values, ncfile, varn):
+def ncstepdiff(values, ncfile, varn):
     """ Function to compute differencies between time-steps (deacumulate) a netCDF file
       new variables labelled [varn] + 'de'will be added to the file
@@ -13460 +13460 @@
       ncfile= netCDF file to replace something
       varn= ',' list of variables to deaccumulate ('all' for all)
-      ncstediffp('Time,Times', '/home/lluis/PY/wrfout_d01_2001-11-11_00:00:00.tests', 'Q2')
+      ncstepdiff('Time,Times', '/home/lluis/PY/wrfout_d01_2001-11-11_00:00:00.tests', 'Q2')
     """
-    fname = 'ncstediffp'
+    fname = 'ncstepdiff'
 
     if values == 'h':
         print fname + '_____________________________________________________________'
-        print ncstediffp.__doc__
+        print ncstepdiff.__doc__
         quit()
 
@@ -13480 +13480 @@
     if varn.find(',') != -1:
         if varn != 'all':
-            varns = varn.slipt(',')
+            varns = varn.split(',')
         else:
             varns = onc.variables
@@ -13526 +13526 @@
             else:
                 print warnmsg
-                print '  ' + fname + ": variable '" + vn + "de' already exist!!"
-                print '    replacing values'
+                print '  ' + fname + ": variable '" + vn + "de' already exist " +    \
+                  'replacing values !!'
                 newvar = onc.variables[vn + 'de']
 
@@ -13547 +13547 @@
 
     return
-ncstediffp('Time,Times', '/home/lluis/PY/wrfout_d01_2001-11-11_00:00:00.tests', 'Q2')
-quit()
+#ncstepdiff('Time,Times', '/home/lluis/PY/wrfout_d01_2001-11-11_00:00:00.tests', 'Q2,T2')
 
 #quit()