Changeset 600 in lmdz_wrf

Timestamp:

Jul 15, 2015, 11:11:44 AM (9 years ago)

Author:

lfita

Message:

Adding on `RAINTOT', what should happen with only 1 time-step

File:

• trunk/tools/diagnostics.py (modified) (1 diff)

trunk/tools/diagnostics.py

@@ -981 +981 @@
 
 # Transforming to a flux
-        if depvars[2] != 'WRFtime':
-            dtimeunits = var2.getncattr('units')
-            tunits = dtimeunits.split(' ')[0]
-
-            dtime = (var2[1] - var2[0])*timeunits_seconds(tunits)
+        if var2.shape[0] > 0:
+            if depvars[2] != 'WRFtime':
+                dtimeunits = var2.getncattr('units')
+                tunits = dtimeunits.split(' ')[0]
+   
+                dtime = (var2[1] - var2[0])*timeunits_seconds(tunits)
+            else:
+                var2 = ncobj.variables['Times']
+                time1 = var2[0,:]
+                time2 = var2[1,:]
+                tmf1 = ''
+                tmf2 = ''
+                for ic in range(len(time1)):
+                    tmf1 = tmf1 + time1[ic]
+                    tmf2 = tmf2 + time2[ic]
+                dtdate1 = dt.datetime.strptime(tmf1,"%Y-%m-%d_%H:%M:%S")
+                dtdate2 = dt.datetime.strptime(tmf2,"%Y-%m-%d_%H:%M:%S")
+                diffdate12 = dtdate2 - dtdate1
+                dtime = diffdate12.total_seconds()
+                print 'dtime:',dtime
         else:
-            var2 = ncobj.variables['Times']
-            time1 = var2[0,:]
-            time2 = var2[1,:]
-            tmf1 = ''
-            tmf2 = ''
-            for ic in range(len(time1)):
-                tmf1 = tmf1 + time1[ic]
-                tmf2 = tmf2 + time2[ic]
-            dtdate1 = dt.datetime.strptime(tmf1,"%Y-%m-%d_%H:%M:%S")
-            dtdate2 = dt.datetime.strptime(tmf2,"%Y-%m-%d_%H:%M:%S")
-            diffdate12 = dtdate2 - dtdate1
-            dtime = diffdate12.total_seconds()
-            print 'dtime:',dtime
+            print warnmsg
+            print '  ' + fname + ": only 1 time-step for '" + diag + "' !!"
+            print '    leaving a zero value!'
+            diagout = var0*0.
+            dtime=1.
 
         ncvar.insert_variable(ncobj, 'pr', diagout/dtime, diagoutd, diagoutvd, newnc)