Changeset 802 in lmdz_wrf

Timestamp:

Jun 2, 2016, 12:42:21 PM (9 years ago)

Author:

lfita

Message:

Adding `secondsDate': Function to know how many units of time are from a given pair of dates

File:

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

trunk/tools/generic_tools.py

@@ -22 +22 @@
 # PolyArea: Function to compute the area of the polygon following 'Shoelace formula'
 # significant_decomposition: Function to decompose a given number by its signifcant potencies
+# secondsDate: Function to know how many units of time are from a given pair of dates
 
 def values_line(line, splitv, chars):
@@ -6357 +6358 @@
 
     return iipot, potencies, potval
+
+def secondsDate(refdate, date, units):
+    """ Function to know how many units of time are from a given pair of dates
+      refdate= date of reference (as [YYYY][MM][DD][HH][MI][SS])
+      date= date to look for (as [YYYY][MM][DD][HH][MI][SS])
+      units= units of time ('century', 'year', 'month', 'day', 'week', 'hour', 'minute', 'second')
+    >>> secondsDate('19490101000000', '19760217082932', 'second')
+    856081772.0
+    >>> secondsDate('19490101000000', '19760217082932', 'minute')
+    237800.492222
+    """
+    import datetime as dt
+    fname = 'SecondsDate'
+
+    yrref=int(refdate[0:4])
+    monref=int(refdate[4:6])
+    dayref=int(refdate[6:8])
+    horref=int(refdate[8:10])
+    minref=int(refdate[10:12])
+    secref=int(refdate[12:14])
+
+    yrd=int(date[0:4])
+    mond=int(date[4:6])
+    dayd=int(date[6:8])
+    hord=int(date[8:10])
+    mind=int(date[10:12])
+    secd=int(date[12:14])
+
+    refdateT = dt.datetime(yrref, monref, dayref, horref, minref, secref)
+    ddateT = dt.datetime(yrd, mond, dayd, hord, mind, secd)
+
+    diffT = ddateT - refdateT
+
+    if searchInlist(dir(diffT), 'total_seconds'):
+        totsecs = diffT.total_seconds()
+    else:
+        totsecs = diffT.days*24*3600. + diffT.seconds
+
+    if units == 'century':
+        times = totsecs / (100.*365.*24.*3600.)
+    elif units == 'year':
+        times = totsecs / (365.*24.*3600.)
+    elif units == 'month':
+        print errormsg
+        print '  ' + fname + "' units '" + units + "' not ready!!"
+        quit(-1)
+    elif units == 'week':
+        times = totsecs / (7.*24.*3600.)
+    elif units == 'day':
+        times = totsecs / (24.*3600.)
+    elif units == 'hour':
+        times = totsecs / (3600.)
+    elif units == 'minute':
+        times = totsecs / (3600.)
+    elif units == 'second':
+        times = totsecs / (1.)
+
+    return times
+
+print secondsDate('19490101000000', '19760217082932', 'second')
+print secondsDate('19490101000000', '19760217082932', 'minute')
+
+quit()