Changeset 1781 in lmdz_wrf

Timestamp:

Feb 23, 2018, 2:14:03 PM (7 years ago)

Author:

lfita

Message:

Fixing use of 'WRFtime' in 'netcdf_fold_concatenation_HMT'

File:

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

trunk/tools/nc_var_tools.py

@@ -11322 +11322 @@
         print infmsg
         print '    ' + fname + ": creation of variable 'time' from WRF 'Times' !!"
-        odimvar = ncobjs[0].variables[vardimname]
+        odimvar = ncobjs[0].variables['Times']
         timewrfv = odimvar[:]
-        tvals, tunits = compute_WRFtime(timewrfv, refdate='19491201000000', tunitsval='minutes')
+        tvals, urefvals = 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)
+        basicvardef(newvar, 'time', 'Time', urefvals)
         newvar.setncattr('axis', 'T')
         newvar.setncattr('_CoordinateAxisType', 'Time')
@@ -11355 +11355 @@
     for ifile in range(Nfiles):
         slicevartot = []
-        if not ncobjs[ifile].variables.has_key(vardimname):
+        if vardimname != 'WRFtime' and not ncobjs[ifile].variables.has_key(vardimname):
             print errormsg
             print '  ' + fname + ": file '" + confiles[ifile] + "' does not have " + \
               " dimension-variable '" + vardimname + "' !!"
             quit(-1)
-        objvar = ncobjs[ifile].variables[vardimname]
-        for dimn in objvar.dimensions:
-            ldimfile = len(ncobjs[ifile].dimensions[dimn])
+        if vardimname != 'WRFtime':
+            objvar = ncobjs[ifile].variables[vardimname]
+            for dimn in objvar.dimensions:
+                ldimfile = len(ncobjs[ifile].dimensions[dimn])
+                if dimn == condim:
+                    slicevartot.append(slice(begslicetot,begslicetot+ldimfile))
+                    begslicetot = begslicetot + ldimfile
+                else:
+                    slicevartot.append(slice(0,ldimfile))
+                if tunits:
+                    uvals = objvar.getncattr('units')
+                    #urefvals, objvals = gen.check_timestep(odimvar)
+                    #uvals, objvals = gen.check_timestep(objvar)
+                    newtimes = gen.coincident_CFtimes(objvar[:], urefvals, uvals)
+                    newvar[tuple(slicevartot)] = newtimes[:]
+                else:
+                    newvar[tuple(slicevartot)] = objvar[:]
+                newnc.sync()
+        else:
+            odimvarf = ncobjs[0].variables['Times']
+            timewrfvf = odimvarf[:]
+            tvalsf, utimef = compute_WRFtime(timewrfvf, refdate='19491201000000', tunitsval='minutes')
+            ldimfile = tvals.shape[0] 
+            dimn = 'time'
+            condimn = 'time'
             if dimn == condim:
                 slicevartot.append(slice(begslicetot,begslicetot+ldimfile))
@@ -11369 +11391 @@
                 slicevartot.append(slice(0,ldimfile))
             if tunits:
-                uvals = objvar.getncattr('units')
+                uvals = utimef
                 #urefvals, objvals = gen.check_timestep(odimvar)
                 #uvals, objvals = gen.check_timestep(objvar)
-                newtimes = gen.coincident_CFtimes(objvar[:], urefvals, uvals)
+                newtimes = gen.coincident_CFtimes(tvalsf, urefvals, uvals)
                 newvar[tuple(slicevartot)] = newtimes[:]
-            else:
-                newvar[tuple(slicevartot)] = objvar[:]
             newnc.sync()