Changeset 6104

Timestamp:

Dec 24, 2008, 1:39:10 PM (14 years ago)

Author:

ole

Message:

Tested that input to culvert rating curve are guaranteed to
exercise the out of bounds condition

Location:

anuga_core/source/anuga/culvert_flows

Files:

• culvert_class.py (modified) (1 diff)
• test_culvert_class.py (modified) (2 diffs)

anuga_core/source/anuga/culvert_flows/culvert_class.py

@@ -231 +231 @@
                                                    indices=[self.enquiry_indices[i]])            
             opening.elevation = elevation
+            opening.stage = elevation # Simple assumption that culvert is dry initially
 
         # Assume two openings for now: Referred to as 0 and 1

anuga_core/source/anuga/culvert_flows/test_culvert_class.py

@@ -4 +4 @@
 import unittest
 import os.path
+import sys
 
 from anuga.utilities.system_tools import get_pathname_from_package
@@ -125 +126 @@
         #-----------------------------------------------------------------------
 
+        min_delta_w = sys.maxint 
+        max_delta_w = -min_delta_w
         for t in domain.evolve(yieldstep = 1, finaltime = 25):
+            delta_w = culvert.inlet.stage - culvert.outlet.stage
+            
+            if delta_w > max_delta_w: max_delta_w = delta_w
+            if delta_w < min_delta_w: min_delta_w = delta_w            
+            
             #print domain.timestepping_statistics()
             pass
-    
-
+
+        # Check that extreme values in rating curve have been exceeded
+        # so that we know that condition has been exercised
+        assert min_delta_w < 0
+        assert max_delta_w > 10        
+        
 
     def test_that_culvert_dry_bed_rating(self):