Changeset 949 in lmdz_wrf for trunk/tools

Timestamp:

Jun 24, 2016, 5:13:03 PM (8 years ago)

Author:

lfita

Message:

Fixing fractional months for 'CFUmonhts_Udays'

File:

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

trunk/tools/generic_tools.py

@@ -7968 +7968 @@
         newmm = mm + datev[itv] - refm
         if newmm/int(newmm) - 1. != 0.:
-            print errormsg
-            print '  ' + fname + ': fractional months values', newmm, 'not ready!!'
-            quit(-1)
-        newmm = np.int(newmm)
-
-        idate = dt.datetime(yr, newmm, dd, hh, mi, ss)
+            print warnmsg
+            print '  ' + fname + ': fractional months values', newmm, '!!'
+            mmdec = newmm - int(newmm)*1.
+            newmm = np.int(newmm)
+            daysnextmon = days_month(yr,newmm+1)
+            newdd = daysnextmon*mmdec 
+            dddec = newdd - int(newdd)*1.
+            newdd = int(newdd)
+            newhh = dddec*24.
+            hhdec = newhh - int(newhh)*1.
+            newhh = int(newhh)
+            newmi = hhdec*60.
+            midec = newmi - int(newmi)*1.
+            newmi = int(newmi)
+            newss = midec*60.
+            ssdec = newss - int(newss)*1.
+            newss = int(newss)
+            if ssdec != 0.: print '    residual seconds after adjustmen:', ssdec
+        else:
+            newmm = np.int(newmm)
+            newdd = np.int(newdd)
+            newhh = np.int(newhh)
+            newmi = np.int(newmi)
+            newss = np.int(newss)
+
+        idate = dt.datetime(yr, newmm, newdd, newhh, newmi, newss)
         
         diffT = idate - refT