Changeset 1551 in lmdz_wrf

Timestamp:

May 1, 2017, 7:36:14 PM (8 years ago)

Author:

lfita

Message:

Fixing `compute_opersvarsfiles':

• Adding complement variables for the non-sliced dimension-variable
• Adding a counter increase which was missed

File:

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

trunk/tools/nc_var_tools.py

@@ -7391 +7391 @@
                 if dimv != 'WRFtime':
                     odimv = objnc.variables[dimv]
-                    print fname + ': Lluis dictslice:', dictslice
-                    print fname + ': Lluis odimv:', odimv.dimensions
                     varslice, Dddv[dimn] = SliceVarDict(odimv, dictslice)
-                    print fname + ': Lluis varslice:', varslice
-                    print fname + ': Lluis Dddv:', Dddv
                     Ddv[dimn] = odimv[tuple(varslice)]
                     Dtdv[dimn] = odimv.dtype
@@ -7447 +7443 @@
                 found = False
                 if not gen.searchInlist(Dddv.keys(), dimn):
+                    ovar = objnc.variables[dimn]
                     print warnmsg
                     print '  ' + fname + ": sliced variable has dimension '" + dimn+ \
@@ -7452 +7449 @@
                     print "    add range for the dimension as '[dim]|[dimv]|[range]'"
                     print '    trying now tacking all values of the dimension-variable'
-                    print "    ", varn + "'s dimensions:", dimvarvals
+                    print "    " + dimn + "'s dimensions:", ovar.dimensions
                     print '    provided slicing dimension-variables:', Dddv.keys()
-                    Dddv[dimn] = dimn
+                    Dddv[dimn] = [dimn]
+                    Ddv[dimn] = ovar[:]
+                    Dtdv[dimn] = ovar.dtype
+                    dicattrs = {}
+                    for attrn in odimv.ncattrs():
+                        attrval = odimv.getncattr(attrn)
+                        dicattrs[attrn] = attrval
+                    Dadv[dimn] = dicattrs
                     #quit(-1)
         else:
@@ -7557 +7561 @@
     #   dimension name: dim + '_' + operS + '_' + var
     #   dimension size: value from the shape of the additional values
+
 
     for operSvarn in Dadditionalv.keys():
@@ -7583 +7588 @@
                     newdim = objofile.createDimension(newdimn,dsize)
                     newvdims.append(newdimn)
+                    idim = idim + 1
 
                 for Hvarn in Hvarns: