Changeset 1088 in lmdz_wrf for trunk

Timestamp:

Sep 9, 2016, 5:31:47 PM (9 years ago)

Author:

lfita

Message:

Adding CFRefTime in `LMDZ_toCF'

File:

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

trunk/tools/nc_var_tools.py

@@ -16737 +16737 @@
     return
 
-def LMDZ_toCF(ncfile):
+def LMDZ_toCF(values,ncfile):
     """ Function to pass a LMDZ original file to CF-conventions
+      values= [CFtimeUnits]
+        [CFtimeUnits]: [tunits] since [YYYYMMDDHHMISS], CF time-units to which variable time should be used ('!' for spaces)
       ncfile= file to transform 
-      variables= variables to transform
     """
     fname = 'LMDZ_toCF'
@@ -16748 +16749 @@
         print LMDZ_toCF.__doc__
         quit()
+
+    CFT = values.replace('!',' ')
 
     ofile = 'CF_LMDZ.nc'
@@ -16760 +16763 @@
         print fname + ": renaming variable 'pres' as 'p' !!"
         newname = ncf.renameVariable('pres','p')
+    ncf.sync()
+
+    # Re-aranging time-axis
+    tvaro = ncf.variables['time_counter']
+    tvals = tvaro[:]
+    tunits = tvaro.getncattr('units')
+
+    RefDate = CFT.split(' ')[2]
+    RefDateS = RefDate[0:4] + '-' + RefDate[4:6] + '-' + RefDate[6:8] + ' ' +        \
+      RefDate[8:10] + ':' + RefDate[10:12] + ':' + RefDate[12:14]
+    CFtimeUnits = ' '.join(CFT.split(' ')[0:2]) + ' '+ RefDateS
+
+    newtvals = gen.coincident_CFtimes(tvals, CFtimeUnits, tunits)
+    tvaro[:] = newtvals
+    newattr = set_attribute(tvaro, 'units', CFtimeUnits)
+    newattr = set_attribute(tvaro, 'time_origin', CFtimeUnits)
+
     ncf.sync()