Changeset 2344 in lmdz_wrf for trunk/tools

Timestamp:

Feb 15, 2019, 9:21:39 PM (6 years ago)

Author:

lfita

Message:

Fixing `except_fillValue' 'exact' filtering

File:

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

trunk/tools/nc_var_tools.py

@@ -29474 +29474 @@
     if not onc.variables.has_key(varname):
         print errormsg
-        print '  ' + fname + ": file '" + + "' does not have variable '" + varname + \
-          "' !!"
+        print '  ' + fname + ": file '" + ncfile + "' does not have variable '" +    \
+          varname + "' !!"
         Varns = list(onc.variables.keys())
         Varns.sort()
@@ -29488 +29488 @@
             if not gen.searchInlist(onc.dimensions,dn):
                 print errormsg    
-                print '  ' + fname + ": file '" + + "' does not have dimension '" +  \
-                  dn + "' !!"
+                print '  ' + fname + ": file '" + ncfile + "' does not have " +      \
+                  "dimension '" + dn + "' !!"
                 dimns = list(onc.dimensions)
                 dimns.sort()
@@ -29530 +29530 @@
         if not onc.variables.has_key(varn):
             print errormsg
-            print '  ' + fname + ": file '" + + "' does not have variable '" + varn+ \
-              "' !!"
+            print '  ' +fname + ": file '" + ncfile + "' does not have variable '" + \
+              varn + "' !!"
             Varns = list(onc.variables.keys())
             Varns.sort()
@@ -29547 +29547 @@
     varns = newvarns + []
 
+    fvals = []
     if kind == 'exact':
         varvalue = gen.typemod(varvalue, varv.dtype)
         mavarv = ma.masked_not_equal(varv, varvalue)
+        fvals.append(varvalue)
     elif kind == 'range':
         ivarvalue = gen.typemod(ivarvalue, varv.dtype)
         evarvalue = gen.typemod(evarvalue, varv.dtype)
         mavarv = ma.masked_outside(varv, ivarvalue, evarvalue)
+        fvals.append(ivarvalue)
+        fvals.append(evarvalue)
     mamask = mavarv.mask
+    Nmask = np.sum(mamask)
+    Tshape = np.prod(varv.shape)
+    if Nmask == Tshape:
+        print errormsg
+        print '  ' + fname + ": filtering '" + varname + "' by '" + kind + "' =",    \
+          fvals, ' less without values !!'
+        quit(-1)
+
+    print '  ' + fname + ': filtering', Tshape-Nmask, 'from:', Tshape
 
     if ovar.dtype == type(1):