Changeset 2188 in lmdz_wrf for trunk/tools/generic_tools.py

Timestamp:

Oct 17, 2018, 7:49:15 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `get_right_CFtimeunits': Function to get the right CFtime units from any given format

File:

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

trunk/tools/generic_tools.py

@@ -111 +111 @@
 # fill_Narray: Function to fill a n-dimensional array with an arrary of lesser rank
 # get_configuration: Function to get the configuration from an ASCII external file
+# get_right_CFtimeunits: Function to get the right CFtime units from any given format
 # 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
@@ -14631 +14632 @@
     return [ix, iy], mindist
 
+def get_right_CFtimeunits(CFtimeu):
+    """ Function to get the right CFtime units from any given format
+       CF-time units must be like:
+         [time_unit] since [YYYY]-[MM]-[DD] [HH]:[MI]:[SS]
+      CFtimeu: CF-time units to look
+    >>> get_right_CFtimeunits('minutes since 1949-12-01 00:00:00')
+    minutes simnce 1949-12-01 00:00:00
+    >>> get_right_CFtimeunits('minutes since 1949-12-01')
+    minutes simnce 1949-12-01 00:00:00
+    >>> get_right_CFtimeunits('minutes since 19491201 000000')
+    minutes simnce 1949-12-01 00:00:00
+    >>> get_right_CFtimeunits('days since 1850-1-1')
+    days simnce 1850-01-01 00:00:00
+    """
+    fname = 'get_right_CFtimeunits'
+
+    if CFtimeu.find('since') == -1:
+        print errormsg
+        print '  ' + fname + ": provided time-units '" + CFtimeu + "' does not " +   \
+          "have even 'since', impossible to proceed !!"
+        quit(-1)
+
+    uvalsecs = CFtimeu.split(' ')
+
+    # Getting four sections
+    if len(uvalsecs) == 4:
+        date=uvalsecs[2]
+        time=uvalsecs[3]
+    else:
+        if uvalsecs[2].find('_') == -1:
+            print infmsg
+            print '  ' + fname + ": CF-time units without time !!"
+            print "    adding '00:00:00'"
+            date=uvalsecs[2]
+            time='00:00:00'
+        else:
+            date=uvalsecs[2].split('_')[0]
+            time=uvalsecs[2].split('_')[1]
+
+    # Getting four digits
+    if date.find('-') == -1: 
+        if len(date) == 8: date = date[0:4] + '-' + date[4:6] + '-' + date[6:8]
+        else: date = date.replace('/','-')
+    datesecs=date.split('-')
+    date = datesecs[0].zfill(4)+'-'+datesecs[1].zfill(2)+'-'+datesecs[2].zfill(2)
+    #print '  ' + fname + ': final date:', date
+
+    if time.find(':') == -1:
+        if len(time) == 6: time = time[0:2] + ':' + time[2:4] + ':' + time[4:6]
+        else: date = date.replace('_',':')
+    timesecs=time.split(':')
+    time = timesecs[0].zfill(2)+':'+timesecs[1].zfill(2)+':'+timesecs[2].zfill(2)
+    #print '  ' + fname + ': final time:', time
+
+    refdateS = date + ' ' + time
+
+    if int(date.split('-')[0]) > 1900:
+        refdate = datetimeStr_datetime(refdateS)
+        refdateS = refdate.strftime("%Y-%m-%d %H:%M:%S")
+    else:
+        print warnmsg
+        print '  ' + fname + ": reference date: '" + date + "' below 1900 !!"
+
+    time_unit = uvalsecs[0]
+
+    rightCFtimeu = time_unit + ' simnce ' + refdateS
+
+    return rightCFtimeu
 
 #quit()