Changeset 2178 in lmdz_wrf for trunk

Timestamp:

Oct 11, 2018, 1:58:56 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `generate_CFtimes': Function to generate CFtimes for a given period, frequency and units

File:

• trunk/tools/generic_tools.py (modified) (4 diffs)

trunk/tools/generic_tools.py

@@ -113 +113 @@
 # get_specdictionary_HMT: Function to get specific values from a dictionary by selcting that keys with H*M*T
 # getting_fixedline: Function to get the values from a line of text with fixed lenght of different values
+# generate_CFtimes: Function to generate CFtimes for a given period, frequency and units
 # grib_CFequiv: Function to provide the CF name of a GRIB variable code number
 # inf_operSlist: Function to provide information from a string as a list separated by 
@@ -14539 +14540 @@
       datei: initial date of the period ([YYYY][MM][DD][HH][MI][SS] format)
       datef: final date of the period ([YYYY][MM][DD][HH][MI][SS] format)
-      units: units of time (one from 'years', 'months', 'days', 'minutes', 'hours', 
-        'seconds')
+      units: units of time (one from 'days', 'hours', 'minutes', 'seconds')
       Tfreq: frequency of time-steps in number of [units] time
       dateref: reference date ([YYYY][MM][DD][HH][MI][SS] format)
+    >>> generate_CFtimes('20021230000000', '20030103000000', 'days', 1)
+    [19387.0, 19388.0, 19389.0, 19390.0, 19391.0, 19392.0], 'days since 1949-12-01 00:00:00'
     """
     import datetime as dt
     fname = 'generate_CFtimes'
 
-    availtunits = ['years', 'months', 'days', 'minutes', 'hours', 'seconds']
+    availtunits = ['days', 'minutes', 'hours', 'seconds']
+
+    if units == 'days': ddt = 3600.*24.
+    elif units == 'hours': ddt = 3600.
+    elif units == 'minutes': ddt = 60.
+    elif units == 'seconds': ddt = 1.
+    else:
+        print errormsg
+        print '  ' + fname + ": time units '" + uints + "' not ready !!"
+        print '    available ones:', availtunits
+        quit(-1)
 
     idate = dt.datetime.strptime(datei, '%Y%m%d%H%M%S')
@@ -14561 +14573 @@
         diffT = delta.days*24*3600. + delta.seconds
 
-    cftimes = [diffT]
+    cftimes = [np.float64(diffT/ddt)]
 
     date = idate
-    while date < edate:
+    while date <= edate:
         date = advance_date(date, units, freq)
 
@@ -14573 +14585 @@
         else:
             diffT = delta.days*24*3600. + delta.seconds
-        cftimes.append(diffT)
-        print 'Lluis:', date, ':', diffT, 'frq:', freq
+        cftimes.append(np.float64(diffT/ddt))
     
-    delta = edate - refdate
-    # Python version depending
-    if searchInlist(dir(delta), 'total_seconds'):
-        diffT = delta.total_seconds()
-    else:
-        diffT = delta.days*24*3600. + delta.seconds
-
-    cftimes.append(diffT)
     timeu = units + ' since ' + refdateS
 
     return cftimes, timeu
 
-print generate_CFtimes('20021201000000', '20030228000000', 'days', 1)
-
 #quit()