Changeset 310 in lmdz_wrf

Timestamp:

Feb 26, 2015, 11:22:37 AM (10 years ago)

Author:

lfita

Message:

Working 1-var version of the 'ncstediffp'

File:

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

trunk/tools/nc_var_tools.py

@@ -2651 +2651 @@
                   '" does not have attribute: "units"'
                 units = None
+                self.unitsval = None
+                self.units = None
+                self.calendar = None
+                self.Urefdate = None
+                self.Srefdate = None
+                self.refdate = None
+                if len(times.shape) == 1:
+                    self.dimt = len(times[:])
+                    self.dt = times[1] - times[0]
+                    self.firstTu = times[0]
+                    self.firstTt = None
+                    self.firstTS = None
+                    self.firstTm = None
+                    self.lastTu = times[self.dimt]
+                    self.lastTt = None
+                    self.lastTS = None
+                    self.lastTm = None
+                else:
+                    self.dimt = None
+                    self.dt = None
+                    self.firstTu = None
+                    self.firstTt = None
+                    self.firstTS = None
+                    self.firstTm = None
+                    self.lastTu = None
+                    self.lastTt = None
+                    self.lastTS = None
+                    self.lastTm = None
             else:
                 units = times.getncattr('units')
@@ -13498 +13526 @@
             else:
                 print warnmsg
-                print '  ' + fname ": variable '"vn + "de' already exist!!"
+                print '  ' + fname + ": variable '" + vn + "de' already exist!!"
                 print '    replacing values'
                 newvar = onc.variables[vn + 'de']
@@ -13506 +13534 @@
             newattr = set_attribute(newvar, 'deaccumulated', 'True')
             if timeinf.dt is not None:
+                print timeinf.dt
                 newattr = set_attributek(newvar, 'time_step_diff', timeinf.dt,       \
                   'npfloat')