Changeset 1843 in lmdz_wrf

Timestamp:

Mar 22, 2018, 8:28:48 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `CFtime_freq': Function to provide the frequency of CF times (in seconds)

File:

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

trunk/tools/generic_tools.py

@@ -57 +57 @@
 # CFcorValues: Function to provide CF information about a given coordinate
 # CFmonthU_daysU: Function to transform from a CF date series with units as 'months since [DATE]' to 'days since [DATE]'
+# CFtime_freq: Function to provide the frequency of CF times (in seconds)
 # CFvar_DIAGvar: Function to provide which model diagnostic values can provide a CF-variable from ASCII file
 # CFvar_MODvar: Function to provide which model values can provide a CF-variable from ASCII file
@@ -664 +665 @@
           str(da).zfill(2) + ' ' + str(ho).zfill(2) + '-' + str(mi).zfill(2) + '-' + \
           str(se).zfill(2) 
-    elif typeSo = 'Y-m-dTH:M:SZ':
+    elif typeSo == 'Y-m-dTH:M:SZ':
         dateYmdHMS = str(yr).zfill(4) + '-' + str(mo).zfill(2) + '-' +               \
           str(da).zfill(2) + 'T' + str(ho).zfill(2) + ':' + str(mi).zfill(2) + ':' + \
@@ -12700 +12701 @@
     # Dictionary for each dimension
     CFlon= {'dimn': 'lon', 'stdn': 'longitude', 'longname': 'Longitude', 'units':    \
-      'degrees_east', 'axis': 'X', '_CoordinateAxisType': 'Lon', 'maxrank:', 2,      \
+      'degrees_east', 'axis': 'X', '_CoordinateAxisType': 'Lon', 'maxrank': 2,       \
       'length': 0}
     CFlat= {'dimn': 'lat', 'stdn': 'latitude', 'longname': 'Latitude', 'units':      \
-      'degrees_north', 'axis': 'Y', '_CoordinateAxisType': 'Lat', 'maxrank:', 2,     \
+      'degrees_north', 'axis': 'Y', '_CoordinateAxisType': 'Lat', 'maxrank':  2,     \
       'length': 0}
     CFpress= {'dimn': 'pressure', 'stdn': 'pressure', 'longname': 'Pressure',        \
       'units': 'Pa', 'axis': 'Z', '_CoordinateAxisType': 'Pres', 'positive': 'down', \
-      'maxrank:', 1, 'length': 0}
+      'maxrank': 1, 'length': 0}
     CFheight= {'dimn': 'height', 'stdn': 'height', 'longname': 'Height',             \
       'units': 'm', 'axis': 'Z', '_CoordinateAxisType': 'Height', 'positive': 'up',  \
-      'maxrank:', 1, 'length': 0}
+      'maxrank': 1, 'length': 0}
     CFtime= {'dimn': 'time', 'stdn': 'time', 'longname': 'Time',                     \
       'units': '[t] since [date]', 'axis': 'T', '_CoordinateAxisType': 'Time',       \
-      'calendar': '[calendar]', 'maxrank:', 1, 'length': -1}
+      'calendar': '[calendar]', 'maxrank': 1, 'length': -1}
 
     # names to search for
@@ -12744 +12745 @@
     return cdfim
 
+def CFtime_freq(time1, time2, tunits):
+    """ Function to provide the frequency of CF times (in seconds)
+      time1: time-step value
+      time2: following time-step value
+      tunits: time untis following CF conventions
+    >>> CFtime_freq(521438., 521439., 'hours since 1950-01-01_00:00:00')
+    3600.0
+    """
+    fname = 'CFtime_freq'
+
+    timediff = time2 - time1
+
+    tunitsval = tunits.split(' ')[0]
+
+    # Units for the variable 'time' to seconds
+    availtunitvals = ['centuries', 'years', 'days', 'hours', 'minutes', 'seconds']
+    if tunitsval == 'centuries':
+        tsecsuv = 100 * 365 * 24 * 3600.
+    elif tunitsval == 'years':
+        tsecsuv = 365 * 24 * 3600.
+    elif tunitsval == 'weeks':
+        tsecsuv = 7 * 24 * 3600.
+    elif tunitsval == 'days':
+        tsecsuv = 24 * 3600.
+    elif tunitsval == 'hours':
+        tsecsuv = 3600.
+    elif tunitsval == 'minutes':
+        tsecsuv = 60.
+    elif tunitsval == 'seconds':
+        tsecsuv = 1.
+    else:
+        print errormsg
+        print '  ' + fname + ":' tunitsvals= '" + tunitsval + "' not ready!!"
+        print '    available ones: ', availtunitvals
+        quit(-1)
+
+    freq = timediff*tsecsuv
+
+    return freq
+#print CFtime_freq(521438., 521439., 'hours since 1950-01-01_00:00:00')
+
 #quit()