Changeset 1673 in lmdz_wrf for trunk/tools/nc_var_tools.py

Timestamp:

Dec 4, 2017, 12:02:14 AM (7 years ago)

Author:

lfita

Message:

Adding:

• `compute_WRFtime': function to copmute CFtimes from WRFtime variable

File:

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

trunk/tools/nc_var_tools.py

@@ -46 +46 @@
 # compute_tevolboxtraj: Function to compute tevolboxtraj: temporal evolution at a given point along a box following a trajectory
 # compute_tevolboxtraj_radialsec: Function to compute tevolboxtraj_radialsec: temporal evolution at a given point along a number of radii at a given frequency following a trajectory
+# compute_WRFtime: unction to copmute CFtimes from WRFtime variable
 # curve_section: Function to provide a section of a file following a given curve
 # DataSetSection: Function to get a section (values along a dimension) of a given data-set
@@ -20662 +20663 @@
     return
 
+def compute_WRFtime(timewrfv, refdate='19491201000000', tunitsval='minutes'):
+    """ Function to copmute CFtimes from WRFtime variable
+      refdate= [YYYYMMDDMIHHSS] format of reference date
+      tunitsval= CF time units
+      timewrfv= matrix string values of WRF 'Times' variable
+    """
+    fname = 'compute_WRFtime'
+
+    yrref=refdate[0:4]
+    monref=refdate[4:6]
+    dayref=refdate[6:8]
+    horref=refdate[8:10]
+    minref=refdate[10:12]
+    secref=refdate[12:14]
+
+    refdateS = yrref + '-' + monref + '-' + dayref + ' ' + horref + ':' + minref +   \
+      ':' + secref
+
+    dt = timewrfv.shape[0]
+    WRFtime = np.zeros((dt), dtype=np.float)
+
+    for it in range(dt):
+        wrfdates = gen.datetimeStr_conversion(timewrfv[it,:],'WRFdatetime', 'matYmdHMS')
+        WRFtime[it] = gen.realdatetime1_CFcompilant(wrfdates, refdate, tunitsval)
+
+    tunits = tunitsval + ' since ' + refdateS
+
+    return WRFtime, tunits
+
 #quit()