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

Timestamp:

Feb 21, 2018, 4:28:22 PM (7 years ago)

Author:

lfita

Message:

Adding `WRFtime' on 'netcdf_fold_concatenation_HMT'

File:

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

trunk/tools/nc_var_tools.py

@@ -11210 +11210 @@
       giving Header, Middle, Tail for the name files
     netcdf_fold_concateation(values, ncfile, varn)
-      [values]= [fold],[dimname]
+      [values]= [fold],[dimname],[vardimname]
         [fold]: folder with the location of the netCDF files
         [dimname]: dimension along which files should be concatenated
-        [vardimname]: name of the dimension variable with the values
+        [vardimname]: name of the dimension variable with the values ('WRFtime' for 
+          CF time conversion from WRF 'Times')
       [ncfile]= [header],[middle],[tail]
         [header] = header of the name of the files to concatenate [ncfile]*[middle]*[tail]
@@ -11306 +11307 @@
 # Dimension variable for the concatenation
     
-    odimvar = ncobjs[0].variables[vardimname]
-    odimvardims = odimvar.dimensions
-    odimattrs = odimvar.ncattrs()
-    for dimn in odimvardims:
-        if not gen.searchInlist(newnc.dimensions,dimn):
-            Ldim = len(odimvar.dimensions[dimn])
-    newvar = newnc.createVariable(vardimname, odimvar.dtype, odimvardims)
-    for attrn in odimattrs:
-        attrv = odimvar.getncattr(attrn)
-        newattr = set_attribute(newvar, attrn, attrv)
-        if attrn == 'units': urefvals = attrv
+    if vardimname != 'WRFtime': 
+        odimvar = ncobjs[0].variables[vardimname]
+        odimvardims = odimvar.dimensions
+        odimattrs = odimvar.ncattrs()
+        for dimn in odimvardims:
+            if not gen.searchInlist(newnc.dimensions,dimn):
+                Ldim = len(odimvar.dimensions[dimn])
+        newvar = newnc.createVariable(vardimname, odimvar.dtype, odimvardims)
+        for attrn in odimattrs:
+            attrv = odimvar.getncattr(attrn)
+            newattr = set_attribute(newvar, attrn, attrv)
+            if attrn == 'units': urefvals = attrv
+    else:
+        print infmsg
+        print '    ' + fname + ": creation of variable 'time' from WRF 'Times' !!"
+        odimvar = ncobjs[0].variables[vardimname]
+        timewrfv = odimvar[:]
+        tvals, tunits = compute_WRFtime(timewrfv, refdate='19491201000000', tunitsval='minutes')
+        dimtwrf = len(tvals)
+        if not gen.searchInlist(newnc.dimensions,'time'):
+            newnc.createDimension('time',None)
+        newvar = newnc.createVariable('time', 'f8', 'time')
+        basicvardef(newvar, 'time', 'Time', tunits)
+        newvar.setncattr('axis', 'T')
+        newvar.setncattr('_CoordinateAxisType', 'Time')
+        set_attribute(newvar, 'calendar', 'standard')
 
     # Checking units of condimvar