Changeset 602 in lmdz_wrf for trunk/tools/validation_sim.py

Timestamp:

Jul 16, 2015, 3:25:23 PM (9 years ago)

Author:

lfita

Message:

Changing reference units on `coincident_CFtimes' to seconds

File:

• trunk/tools/validation_sim.py (modified) (2 diffs)

trunk/tools/validation_sim.py

@@ -795 +795 @@
 
         difft = trefTB - trefTA
-        diffv = difft.days*24.*3600.*10.e6 + difft.seconds*10.e6 + difft.microseconds
+        diffv = difft.days*24.*3600. + difft.seconds + difft.microseconds/1.e6
         print '  ' + fname + ': different reference refA:',trefTA,'refB',trefTB
-        print '    difference:',difft,':',diffv,'microseconds'
+        print '    difference:',difft,':',diffv,'(in seconds)'
 
         if tuA == 'microseconds':
+            tB = tB + diffv*1.e6
+        elif tuA == 'seconds':
             tB = tB + diffv
-        elif tuA == 'seconds':
-            tB = tB + diffv/10.e6
         elif tuA == 'minutes':
-            tB = tB + diffv/(60.*10.e6)
+            tB = tB + diffv/(60.)
         elif tuA == 'hours':
-            tB = tB + diffv/(3600.*10.e6)
-        elif tuA == 'dayss':
-            tB = tB + diffv/(24.*3600.*10.e6)
+            tB = tB + diffv/(3600.)
+        elif tuA == 'days':
+            tB = tB + diffv/(24.*3600.)
         else:
             print errormsg
@@ -1937 +1937 @@
 iobst = CFtimes_datetime_NOfile([valdimobs['T'][0]], obstunits, 'standard')
 fobst = CFtimes_datetime_NOfile([valdimobs['T'][dimtobs-1]], obstunits, 'standard')
+print 'Lluis:', simobstimes[0],'tobs:',obstunits
 isimt = CFtimes_datetime_NOfile([simobstimes[0]], obstunits, 'standard')
 fsimt = CFtimes_datetime_NOfile([simobstimes[dimtsim-1]], obstunits, 'standard')