Changeset 6905 for branches/numpy_anuga_validation/automated_validation_tests/flow_tests/test_inflow_using_flowline.py

Timestamp:

Apr 26, 2009, 8:36:31 AM (15 years ago)

Author:

rwilson

Message:

Back-merge from main trunk and ongoing Patong Beach changes.

File:

• branches/numpy_anuga_validation/automated_validation_tests/flow_tests/test_inflow_using_flowline.py (modified) (5 diffs)

branches/numpy_anuga_validation/automated_validation_tests/flow_tests/test_inflow_using_flowline.py

@@ -18 +18 @@
 from anuga.shallow_water.shallow_water_domain import Reflective_boundary
 from anuga.shallow_water.shallow_water_domain import Dirichlet_boundary
-from anuga.shallow_water.shallow_water_domain import Transmissive_Momentum_Set_Stage_boundary
+from anuga.shallow_water.shallow_water_domain import Transmissive_momentum_set_stage_boundary
 from anuga.shallow_water.shallow_water_domain import Inflow
 from anuga.shallow_water.data_manager import get_flow_through_cross_section
@@ -46 +46 @@
 
         domain = Domain(points, vertices, boundary)   
-        domain.set_name('inflow_flowline_test')              # Output name
+        domain.set_name('inflow_flowline_test')
 
 
@@ -57 +57 @@
             return z
 
-        domain.set_quantity('elevation', topography)  # Use function for elevation
-        domain.set_quantity('friction', mannings_n)   # Constant friction of conc surface   
+        domain.set_quantity('elevation', topography)
+        domain.set_quantity('friction', mannings_n)
         domain.set_quantity('stage',
-                            expression='elevation')   # Dry initial condition
+                            expression='elevation') # Dry initial condition
 
 
@@ -93 +93 @@
         #----------------------------------------------------------------------
 
-        Br = Reflective_boundary(domain)              # Solid reflective wall
-        #Bo = Dirichlet_boundary([-(slope*length)+normal_depth, 0, 0])  # Outflow stage at normal depth
+        Br = Reflective_boundary(domain) # Solid reflective wall
+
+        # Define downstream boundary based on predicted depth
         def normal_depth_stage_downstream(t):
-            return (-slope*length) + normal_depth 
-        Bt = Transmissive_Momentum_Set_Stage_boundary(domain=domain,function=normal_depth_stage_downstream)
+            return (-slope*length) + normal_depth
+        
+        Bt = Transmissive_momentum_set_stage_boundary(domain=domain,
+                                                      function=normal_depth_stage_downstream)
+        
         domain.set_boundary({'left': Br, 'right': Bt, 'top': Br, 'bottom': Br})
 
@@ -110 +114 @@
             if verbose :
                 print domain.timestepping_statistics()
+                print domain.volumetric_balance_statistics()