Changeset 1088

Timestamp:

Mar 16, 2005, 5:46:16 PM (20 years ago)

Author:

prow

Message:

sww2domain now throws an exception if NaNs? are in the file. Option to set NaNs? to a default value.

File:

• inundation/ga/storm_surge/pyvolution/data_manager.py (modified) (3 diffs)

inundation/ga/storm_surge/pyvolution/data_manager.py

@@ -1533 +1533 @@
 
 
-def sww2domain(filename,t=None):
+def sww2domain(filename,t=None,fail_if_NaN=True,NaN_filler=0):
     """Read sww Net CDF file containing Shallow Water Wave simulation
 
@@ -1541 +1541 @@
 
     """
+    NaN=9.969209968386869e+036
     from Scientific.IO.NetCDF import NetCDFFile
     from domain import Domain
@@ -1598 +1599 @@
     for quantity in other_quantities:
         X = fid.variables[quantity][:]
+        print quantity
+        print 'max(X)'
+        print max(X)
+        print 'max(X)==NaN'
+        print max(X)==NaN
+        if (max(X)==NaN) or (min(X)==NaN):
+            if fail_if_NaN:
+                msg = 'quantity %s contains no_data entry'%quantity
+                raise msg
+            else:
+                data = (X<>NaN)
+                X = (X*data)+(data==0)*NaN_filler
         domain.set_quantity(quantity,X)
 #
     for quantity in conserved_quantities:
         X = interpolated_quantities[quantity]
+        print quantity
+        print 'max(X)'
+        print max(X)
+        print 'max(X)==NaN'
+        print max(X)==NaN
+        if (max(X)==NaN) or (min(X)==NaN):
+            if fail_if_NaN:
+                msg = 'quantity %s contains no_data entry'%quantity
+                raise msg
+            else:
+                data = (X<>NaN)
+                X = (X*data)+(data==0)*NaN_filler
         domain.set_quantity(quantity,X)
 #