Changeset 5585 for anuga_core/source/anuga/culvert_flows/culvert_class.py

Timestamp:

Jul 30, 2008, 4:06:55 PM (16 years ago)

Author:

ole

Message:

Multibarrel culvert functionality.

File:

• anuga_core/source/anuga/culvert_flows/culvert_class.py (modified) (6 diffs)

anuga_core/source/anuga/culvert_flows/culvert_class.py

@@ -4 +4 @@
 from anuga.utilities.polygon import inside_polygon
 from anuga.utilities.polygon import is_inside_polygon
+from anuga.utilities.polygon import plot_polygons
+
 
 
@@ -53 +55 @@
                  blockage_bottup=None,
                  culvert_routine=None,
+                 number_of_barrels=1,
                  verbose=False):
         
@@ -88 +91 @@
         assert self.culvert_type in ['circle', 'square', 'rectangle']
         
+        assert number_of_barrels >= 1
+        self.number_of_barrels = number_of_barrels
+        
+        
         # Set defaults
         if manning is None: manning = 0.012   # Set a Default Mannings Roughness for Pipe
@@ -118 +125 @@
                                     end_point1,
                                     width=width,   
-                                    height=height)
+                                    height=height,
+                                    number_of_barrels=number_of_barrels)
         
         if verbose is True:
@@ -128 +136 @@
             #               P['enquiry_polygon1']],
             #              figname='culvert_polygon_output')
+            #import sys; sys.exit()                           
 
 
@@ -362 +371 @@
         log_to_file(log_filename, s)
 
-        
+
+        # Calculate discharge for one barrel
         Q, barrel_velocity, culvert_outlet_depth = self.culvert_routine(self, inlet, outlet, delta_Et, g)
-        #####################################################
+        
+        # Adjust discharge for multiple barrels
+        Q *= self.number_of_barrels
+
+        # Compute barrel momentum
         barrel_momentum = barrel_velocity*culvert_outlet_depth