Changeset 3309

Timestamp:

Jul 11, 2006, 3:57:26 PM (18 years ago)

Author:

sexton

Message:

handling INF's

File:

• inundation/damage/inundation_damage.py (modified) (3 diffs)

inundation/damage/inundation_damage.py

@@ -10 +10 @@
 from pyvolution.util import file_function
 from geospatial_data.geospatial_data import ensure_absolute
-
+from utilities.numerical_tools import INF
 try:  
     import kinds  
@@ -84 +84 @@
         for time in callable_sww.get_time():
             quantities = callable_sww(time,point_i)
+            print "quantities", quantities
+            
             w = quantities[0]
             z = quantities[1]
@@ -89 +91 @@
             vh = quantities[3]
             depth = w - z - ground_floor_height
-            momentum = sqrt(uh*uh + vh*vh)  
+              
             if depth > max_depth[point_i]:
                 max_depth[point_i] = depth
+            if w == INF or z == INF or uh == INF or vh == INF:
+                continue
+            momentum = sqrt(uh*uh + vh*vh)
             if momentum > max_momentum[point_i]:
                 max_momentum[point_i] = momentum