Changeset 2471 in lmdz_wrf for trunk/tools

Timestamp:

Apr 24, 2019, 8:43:23 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `fromanydate_CFYmd': Function to transform from any string format of date to Y-m-d used in CF-conventions

File:

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

trunk/tools/generic_tools.py

@@ -122 +122 @@
 # fill_Narray: Function to fill a n-dimensional array with an arrary of lesser rank
 # fix_CFdates: Fixing CF-time values with wrong setting
+# fromanydate_CFYmd: Function to transform from any string format of date to Y-m-d used in 
+#   CF-conventions
 # gen_impose_gregorian: Function to impose gregorian calendar to a series of times with a 
 #   non-standard calendar independently from datetime (imposes yr > 1900)
@@ -16239 +16241 @@
 #  'days since 1949-12-01 00:00:00')
 
+def fromanydate_CFYmd(date):
+    """ Function to transform from any string format of date to Y-m-d used in 
+        CF-conventions
+      date: date to transform
+    >>> fromanydate_Ymd('0001-01-01')
+    '0001-01-01'
+    >>> fromanydate_Ymd('1850-1-1')
+    '1850-01-01'
+    >>> fromanydate_Ymd('1850/1/1')
+    '1850-01-01'
+    >>> fromanydate_Ymd('1976-02-17_08:32:27')
+    '1976-02-17'
+    """
+    fname = 'fromanydate_CFYmd'
+
+    if date.count('-') + date.count('/') == 0:
+        print errormsg
+        print '  ' + fname + ": date without either '-' or '/' not ready !!"
+        print "    passed date '" + date + "' not manegeable"
+        quit(-1)
+
+    # Replacing '/' by '-'
+    if date.count('/') != 0: date = date.replace('/','-')
+
+    # Removing possible [datre]_[time]
+    if date.count('_') != 0: date = date.split('_')[0]
+
+    ldate = date.split('-')
+    if len(ldate) != 3:
+        print errormsg
+        print '  ' + fname + ": date '" + date + "' not ready !!"
+        print "     it requires 3 '-' sections '[yr]-[mo]-[dd]'"
+        quit(-1)
+
+    yr = str(int(ldate[0])).zfill(4)
+    mo = str(int(ldate[1])).zfill(2)
+    dd = str(int(ldate[2])).zfill(2)
+
+    newdate = yr + '-' + mo + '-' + dd
+
+    return newdate
+
 def change_CFRefdate(timevals, origcftimeu, newSrefdate='1949-12-01 00:00:00'):
     """ Change CF-time values with a new reference date
@@ -16258 +16302 @@
     origtv = origcftimeu.split(' ')
     origtunits = origtv[0]
-    origSrefdate = origtv[2]
+    origSrefdate0 = origtv[2]
+
+    origSrefdate = fromanydate_Ymd(origSrefdate)
 
     yrref = int(origSrefdate[0:4])
@@ -16368 +16414 @@
     origtunits = origtv[0]
     if len(origtv) == 4:
-        origSrefdate = origtv[2] + ' ' + origtv[3]
+        origSrefdate = fromanydate_Ymd(origtv[2]) + ' ' + origtv[3]
     else:
         print warnmsg
         print '  ' + fname + ": original time units wihout time !!"
         print "   adding it as '00:00:00'"
-        origSrefdate = origtv[2] + ' 00:00:00'
+        origSrefdate = fromanydate_Ymd(origtv[2]) + ' 00:00:00'
 
     # From original to days
@@ -16437 +16483 @@
     """
     fname = 'change_CFcalendar'
-    availorigcal = ['noleap', '365d']
+    availorigcal = ['noleap', '365d', '365_day']
     availnewcal = ['gregorian', 'standard']
 
     # If the calendar is not gregorian, it needs to be fixed from origSrefdate!
-    if origcal == 'noleap' or origcal == '365d':
+    if origcal == 'noleap' or origcal == '365d' or origcal == '365_day' :
         print '  ' + fname + ": found non-standard calendar '" + origcal + "' !!"
         if newcal == 'gregorian' or newcal == 'standard':