Changeset 2473 in lmdz_wrf for trunk/tools

Timestamp:

Apr 25, 2019, 2:12:54 PM (6 years ago)

Author:

lfita

Message:

 

File:

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

trunk/tools/generic_tools.py

@@ -16245 +16245 @@
         CF-conventions
       date: date to transform
-    >>> fromanydate_Ymd('0001-01-01')
+    >>> fromanydate_CFYmd('0001-01-01')
     '0001-01-01'
-    >>> fromanydate_Ymd('1850-1-1')
+    >>> fromanydate_CFYmd('1850-1-1')
     '1850-01-01'
-    >>> fromanydate_Ymd('1850/1/1')
+    >>> fromanydate_CFYmd('1850/1/1')
     '1850-01-01'
-    >>> fromanydate_Ymd('1976-02-17_08:32:27')
+    >>> fromanydate_CFYmd('1976-02-17_08:32:27')
     '1976-02-17'
     """
     fname = 'fromanydate_CFYmd'
 
-    if date.count('-') + date.count('/') == 0:
+    if date.count('-') + date.count('/') == 0 + date.count('.') == 0:
         print errormsg
-        print '  ' + fname + ": date without either '-' or '/' not ready !!"
+        print '  ' + fname + ": date without either '-' or '/' or '.' not ready !!"
         print "    passed date '" + date + "' not manegeable"
         quit(-1)
 
-    # Replacing '/' by '-'
+    # Replacing '/', '.' by '-'
     if date.count('/') != 0: date = date.replace('/','-')
+    if date.count('.') != 0: date = date.replace('.','-')
 
     # Removing possible [datre]_[time]
@@ -16304 +16305 @@
     origSrefdate0 = origtv[2]
 
-    origSrefdate = fromanydate_Ymd(origSrefdate)
+    origSrefdate = fromanydate_CFYmd(origSrefdate)
 
     yrref = int(origSrefdate[0:4])
@@ -16414 +16415 @@
     origtunits = origtv[0]
     if len(origtv) == 4:
-        origSrefdate = fromanydate_Ymd(origtv[2]) + ' ' + origtv[3]
+        origSrefdate = fromanydate_CFYmd(origtv[2]) + ' ' + origtv[3]
     else:
         print warnmsg
         print '  ' + fname + ": original time units wihout time !!"
         print "   adding it as '00:00:00'"
-        origSrefdate = fromanydate_Ymd(origtv[2]) + ' 00:00:00'
+        origSrefdate = fromanydate_CFYmd(origtv[2]) + ' 00:00:00'
 
     # From original to days
@@ -16468 +16469 @@
     return newtimevals, newcftimeu
 
+print fromanydate_CFYmd('1979.12.07')
+
+class CFtimeU_inf(object):
+    """ Classs to provide information of CFtime units
+      cftimeu: CF time units to provide information for [Timeu] since [RefDate]
+    """
+
+    def __init__(self, cftimeu):
+        fname = 'CFtimeU_inf'
+        origtv = cftimeu.split(' ')
+        origtunits = origtv[0]
+        origSrefdate0 = origtv[2]
+
+        origSrefdate = fromanydate_CFYmd(origSrefdate0)
+
+        yrref = int(origSrefdate[0:4])
+        monref = int(origSrefdate[5:7])
+        dayref = int(origSrefdate[8:10])
+
+# Does original reference date contain a time value [YYYY]-[MM]-[DD] [HH]:[MI]:[SS]
+##
+        trefT = origcftimeu.find(':')
+        if not trefT == -1:
+            if len(origtv) == 3:
+                horref=int(origSrefdate[11:13])
+                minref=int(origSrefdate[14:16])
+                secref=int(origSrefdate[17:19])
+            else:
+                origSreftime = origtv[3]
+                horref=int(origSreftime[0:2])
+                minref=int(origSreftime[3:5])
+                secref=int(origSreftime[6:8])
+        else:
+            horref = 0
+            minref = 0
+            secref = 0
+
+        self.refdate = origSrefdate + ' ' + horref + ':' + minref + ':' + secref
+        self.refdateYMHdms = yrref + monref + dayref + horref + minref + secref
+        self.Tunits = origtunits
+
+cftimev = CFtimeU_inf('hours since 1979-12-1')
+printing_class(cftimev)
+
+quit()
+
+def from360d_reg(tvals, tunits):
+    """ Function to transform from 360d = 12 * 30d calendar to a regular one
+      tvals: temporal values
+      tunits: units of the values in CF format [Tunits] since [RefDate]
+    """
+    fname = 'from360d_reg'
+
+    
+
+    return newtvals
+
 def change_CFcalendar(timevals, origcftimeu, origcal, newcal='gregorian'):
     """ Changing CF-time values to a new calendar
@@ -16498 +16556 @@
             quit(-1)
 
+    elif origcal == '3650d':
+        # We only can fix the 12x30 days calendar for monthly values...
+        if type(timevals) == type(range(2)) or type(timevals) == type(np.arange(2)):
+            dt = timevals[1] - timevals[0]
+        else:
+            dt = 30.
+        # Transforming to regular months
+
     else:
         print errormsg