Changeset 2477 in lmdz_wrf

Timestamp:

Apr 29, 2019, 4:27:12 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `from360d_reg': Function to transform from 360d = 12 * 30d calendar to a regular one

File:

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

trunk/tools/generic_tools.py

@@ -551 +551 @@
 ###
 #
+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_CFYmd('0001-01-01')
+    '0001-01-01'
+    >>> fromanydate_CFYmd('1850-1-1')
+    '1850-01-01'
+    >>> fromanydate_CFYmd('1850/1/1')
+    '1850-01-01'
+    >>> fromanydate_CFYmd('1976-02-17_08:32:27')
+    '1976-02-17'
+    """
+    fname = 'fromanydate_CFYmd'
+
+    if date.count('-') + date.count('/') == 0 + date.count('.') == 0:
+        print errormsg
+        print '  ' + fname + ": date without either '-' or '/' or '.' not ready !!"
+        print "    passed date '" + date + "' not manegeable"
+        quit(-1)
+
+    # Replacing '/', '.' by '-'
+    if date.count('/') != 0: date = date.replace('/','-')
+    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
+
+class CFtimeU_inf(object):
+    """ Classs to provide information of CFtime units
+      cftimeu: CF time units to provide information for [Timeu] since [RefDate]
+      * attributes:
+        refdate: reference date in Y-m-d H:M:S format
+        refdateYMHdms = reference date in YmdHMS format
+        refdatemat = reference date in matricial format
+        refdatejuliand = reference date as julian date (days since -4712 January 1st 
+          00 UTC)
+        refdateDT = reference date as datetime object
+        Tunits = origtunits
+        utoweeks = Factor to convert values from original temporal units to weeks
+        utodays = Factor to convert values from original temporal units to days
+        utohours = Factor to convert values from original temporal units to hours
+        utominutes = Factor to convert values from original temporal units to minutes
+        utoseconds = Factor to convert values from original temporal units to seconds
+    """
+
+    def __init__(self, cftimeu):
+        import datetime as dt
+        fname = 'CFtimeU_inf'
+
+        availtu = ['weeks', 'days', 'hours', 'minutes', 'seconds']
+
+        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 = cftimeu.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
+
+        horrefS = str(horref).zfill(2)
+        minrefS = str(minref).zfill(2)
+        secrefS = str(secref).zfill(2)
+
+        # Transforming factors
+        if origtunits == 'weeks':
+            self.utoweeks = 1. 
+            self.utodays = 7.
+            self.utohours = 7. * 24.
+            self.utominutes = 7. * 24. * 60.
+            self.utoseconds = 7. * 24. * 3600.
+        elif origtunits == 'days':
+            self.utoweeks = 1. / 7.
+            self.utodays = 1.
+            self.utohours = 24.
+            self.utominutes = 24. * 60.
+            self.utoseconds = 24. * 3600.
+        elif origtunits == 'hours':
+            self.utoweeks = 1. / ( 7. * 24.)
+            self.utodays = 1. / 24.
+            self.utohours = 1.
+            self.utominutes = 60.
+            self.utoseconds = 3600.
+        elif origtunits == 'minutes':
+            self.utoweeks = 1. / ( 7. * 24. * 60.)
+            self.utodays = 1. / ( 24. * 60.)
+            self.utohours = 1. / 60.
+            self.utominutes = 1.
+            self.utoseconds = 60.
+        elif origtunits == 'seconds':
+            self.utoweeks = 1. / ( 7. * 24. * 3600.)
+            self.utodays = 1. / ( 24. * 3600.)
+            self.utohours = 1. / 3600.
+            self.utominutes = 1. / 60.
+            self.utoseconds = 1.
+        else:
+            print errormsg
+            print '  ' + fname + ": new time-units '" + origtunits + "' not ready !!"
+            print '    available ones: ', availtu
+            quit(-1)
+
+        self.refdate = origSrefdate + ' ' + horrefS + ':' + minrefS + ':' + secrefS
+        self.refdateYMHdms = origSrefdate[0:4] + origSrefdate[5:7] +                 \
+          origSrefdate[8:10] + horrefS + ':' + minrefS + ':' + secrefS
+        self.refdatemat = [yrref, monref, dayref, horref, minref, secref]
+        self.refdatejuliand = juliandate(yrref,monref,dayref,horref,minref,secref)
+        if yrref >= 1900: self.refdateDT = dt.datetime(yrref,monref,dayref,horref,   \
+          minref,secref)
+        else: self.refdateDT = None
+
+        self.Tunits = origtunits
+
+
 def days_month(year,month):
     """ Function to give the number of days of a month of a given year
@@ -16231 +16379 @@
 #  '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_CFYmd('0001-01-01')
-    '0001-01-01'
-    >>> fromanydate_CFYmd('1850-1-1')
-    '1850-01-01'
-    >>> fromanydate_CFYmd('1850/1/1')
-    '1850-01-01'
-    >>> fromanydate_CFYmd('1976-02-17_08:32:27')
-    '1976-02-17'
-    """
-    fname = 'fromanydate_CFYmd'
-
-    if date.count('-') + date.count('/') == 0 + date.count('.') == 0:
-        print errormsg
-        print '  ' + fname + ": date without either '-' or '/' or '.' not ready !!"
-        print "    passed date '" + date + "' not manegeable"
-        quit(-1)
-
-    # Replacing '/', '.' by '-'
-    if date.count('/') != 0: date = date.replace('/','-')
-    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'):
@@ -16291 +16397 @@
 
     # original
-    origtv = origcftimeu.split(' ')
-    origtunits = origtv[0]
-    origSrefdate0 = origtv[2]
-
-    origSrefdate = fromanydate_CFYmd(origSrefdate)
-
-    yrref = int(origSrefdate[0:4])
-    monref = int(origSrefdate[5:7])
-    dayref = int(origSrefdate[8:10])
+    inftunits = CFtimeU_inf(origcftimeu)
+    origtunits = inftunits.Tunits
+    origSrefdate = inftunits.refdate
+
+    [yrref, monref, dayref, horref, minref, secref]= inftunits.refdatemat
 
     if int(yrref) < 1970: origsing = -1.
     else: origsing = 1.
-
-# 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
 
     # Transforming to julian days (reference -4712 January 1st 00 UTC)
@@ -16402 +16486 @@
 
     # original
-    origtv = origcftimeu.split(' ')
-    origtunits = origtv[0]
-    if len(origtv) == 4:
-        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_CFYmd(origtv[2]) + ' 00:00:00'
+    inftunits = CFtimeU_inf(origcftimeu)
+    origtunits = inftunits.Tunits
+    origSrefdate = inftunits.refdate
 
     # From original to days
@@ -16458 +16536 @@
 
     return newtimevals, newcftimeu
-
-class CFtimeU_inf(object):
-    """ Classs to provide information of CFtime units
-      cftimeu: CF time units to provide information for [Timeu] since [RefDate]
-      * attributes:
-        refdate: reference date in Y-m-d H:M:S format
-        refdateYMHdms = reference date in YmdHMS format
-        refdatemat = reference date in matricial format
-        refdatejuliand = reference date as julian date (days since -4712 January 1st 
-          00 UTC)
-        refdateDT = reference date as datetime object
-        Tunits = origtunits
-        utoweeks = Factor to convert values from original temporal units to weeks
-        utodays = Factor to convert values from original temporal units to days
-        utohours = Factor to convert values from original temporal units to hours
-        utominutes = Factor to convert values from original temporal units to minutes
-        utoseconds = Factor to convert values from original temporal units to seconds
-    """
-
-    def __init__(self, cftimeu):
-        import datetime as dt
-        fname = 'CFtimeU_inf'
-
-        availtu = ['weeks', 'days', 'hours', 'minutes', 'seconds']
-
-        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 = cftimeu.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
-
-        horrefS = str(horref).zfill(2)
-        minrefS = str(minref).zfill(2)
-        secrefS = str(secref).zfill(2)
-
-        # Transforming factors
-        if origtunits == 'weeks':
-            self.utoweeks = 1. 
-            self.utodays = 7.
-            self.utohours = 7. * 24.
-            self.utominutes = 7. * 24. * 60.
-            self.utoseconds = 7. * 24. * 3600.
-        elif origtunits == 'days':
-            self.utoweeks = 1. / 7.
-            self.utodays = 1.
-            self.utohours = 24.
-            self.utominutes = 24. * 60.
-            self.utoseconds = 24. * 3600.
-        elif origtunits == 'hours':
-            self.utoweeks = 1. / ( 7. * 24.)
-            self.utodays = 1. / 24.
-            self.utohours = 1.
-            self.utominutes = 60.
-            self.utoseconds = 3600.
-        elif origtunits == 'minutes':
-            self.utoweeks = 1. / ( 7. * 24. * 60.)
-            self.utodays = 1. / ( 24. * 60.)
-            self.utohours = 1. / 60.
-            self.utominutes = 1.
-            self.utoseconds = 60.
-        elif origtunits == 'seconds':
-            self.utoweeks = 1. / ( 7. * 24. * 3600.)
-            self.utodays = 1. / ( 24. * 3600.)
-            self.utohours = 1. / 3600.
-            self.utominutes = 1. / 60.
-            self.utoseconds = 1.
-        else:
-            print errormsg
-            print '  ' + fname + ": new time-units '" + origtunits + "' not ready !!"
-            print '    available ones: ', availtu
-            quit(-1)
-
-        self.refdate = origSrefdate + ' ' + horrefS + ':' + minrefS + ':' + secrefS
-        self.refdateYMHdms = origSrefdate[0:4] + origSrefdate[5:7] +                 \
-          origSrefdate[8:10] + horrefS + ':' + minrefS + ':' + secrefS
-        self.refdatemat = [yrref, monref, dayref, horref, minref, secref]
-        self.refdatejuliand = juliandate(yrref,monref,dayref,horref,minref,secref)
-        if yrref >= 1900: self.refdateDT = dt.datetime(yrref,monref,dayref,horref,   \
-          minref,secref)
-        else: self.refdateDT = None
-
-        self.Tunits = origtunits
 
 def from360d_reg(tvals, tunits, kind='right'):
@@ -16712 +16686 @@
 
     elif origcal == '3650d':
-        # We only can fix the 12x30 days calendar for monthly values...
+        print '  ' + fname + ": found non-standard calendar '" + origcal + "' !!"
+        # 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]
@@ -16718 +16693 @@
             dt = 30.
         # Transforming to regular months
+        if newcal == 'gregorian' or newcal == 'standard':
+            newtimevals = from360d_reg(tvals, tunits, kind='monpercen')
 
     else: