Ignore:
Timestamp:
Jul 30, 2008, 4:06:55 PM (14 years ago)
Author:
ole
Message:

Multibarrel culvert functionality.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/culvert_flows/culvert_class.py

    r5584 r5585  
    44from anuga.utilities.polygon import inside_polygon
    55from anuga.utilities.polygon import is_inside_polygon
     6from anuga.utilities.polygon import plot_polygons
     7
    68
    79
     
    5355                 blockage_bottup=None,
    5456                 culvert_routine=None,
     57                 number_of_barrels=1,
    5558                 verbose=False):
    5659       
     
    8891        assert self.culvert_type in ['circle', 'square', 'rectangle']
    8992       
     93        assert number_of_barrels >= 1
     94        self.number_of_barrels = number_of_barrels
     95       
     96       
    9097        # Set defaults
    9198        if manning is None: manning = 0.012   # Set a Default Mannings Roughness for Pipe
     
    118125                                    end_point1,
    119126                                    width=width,   
    120                                     height=height)
     127                                    height=height,
     128                                    number_of_barrels=number_of_barrels)
    121129       
    122130        if verbose is True:
     
    128136            #               P['enquiry_polygon1']],
    129137            #              figname='culvert_polygon_output')
     138            #import sys; sys.exit()                           
    130139
    131140
     
    362371        log_to_file(log_filename, s)
    363372
    364        
     373
     374        # Calculate discharge for one barrel
    365375        Q, barrel_velocity, culvert_outlet_depth = self.culvert_routine(self, inlet, outlet, delta_Et, g)
    366         #####################################################
     376       
     377        # Adjust discharge for multiple barrels
     378        Q *= self.number_of_barrels
     379
     380        # Compute barrel momentum
    367381        barrel_momentum = barrel_velocity*culvert_outlet_depth
    368382               
Note: See TracChangeset for help on using the changeset viewer.