Changeset 2137 in lmdz_wrf

Timestamp:

Sep 13, 2018, 12:00:30 AM (6 years ago)

Author:

lfita

Message:

Fixing netcdf_fold_concatenation_HMT by `WRFtime'?

File:

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

trunk/tools/nc_var_tools.py

@@ -11538 +11538 @@
         print infmsg
         print '    ' + fname + ": creation of variable 'time' from WRF 'Times' !!"
-        odimvar = ncobjs[0].variables['Times']
-        timewrfv = odimvar[:]
-        tvals, urefvals = compute_WRFtime(timewrfv, refdate='19491201000000', tunitsval='minutes')
+        tvals = []
+        for ifile in range(Nfiles):
+            odimvar = ncobjs[ifile].variables['Times']
+            timewrfv = odimvar[:]
+            tvals0, urefvals = compute_WRFtime(timewrfv, refdate='19491201000000',   \
+              tunitsval='minutes')
+            tvals = tvals + list(tvals0)
+        tvals = np.array(tvals)
         dimtwrf = len(tvals)
         if not gen.searchInlist(newnc.dimensions,'time'):
             newnc.createDimension('time',None)
         newvar = newnc.createVariable('time', 'f8', 'time')
+        newvar[:] = np.zeros((totcondim), dtype=np.float64)
         basicvardef(newvar, 'time', 'Time', urefvals)
         newvar.setncattr('axis', 'T')
@@ -11601 +11607 @@
             timewrfvf = odimvarf[:]
             tvalsf, utimef = compute_WRFtime(timewrfvf, refdate='19491201000000', tunitsval='minutes')
-            ldimfile = tvals.shape[0] 
+            ldimfile = tvalsf.shape[0] 
             dimn = 'time'
             condimn = 'time'
-            if dimn == condim:
+            if dimn == condim or (condim == 'Time' and vardimname == 'WRFtime'):
                 slicevartot.append(slice(begslicetot,begslicetot+ldimfile))
                 begslicetot = begslicetot + ldimfile