Ignore:
Timestamp:
Sep 1, 2010, 6:11:25 PM (12 years ago)
Author:
steve
Message:

Added enquiry points to culverts

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/source/anuga/structures/boyd_box_culvert_routine.py

    r7982 r7984  
    7070        local_debug ='false'
    7171       
    72         if self.inflow.get_average_height() > 0.01: #this value was 0.01:
     72        if self.inflow.get_enquiry_height() > 0.01: #this value was 0.01:
    7373            if local_debug =='true':
    7474                log.critical('Specific E & Deltat Tot E = %s, %s'
    75                              % (str(self.inflow.get_average_specific_energy()),
     75                             % (str(self.inflow.get_enquiry_specific_energy()),
    7676                                str(self.delta_total_energy)))
    7777                log.critical('culvert type = %s' % str(culvert_type))
     
    7979
    8080            if self.log_filename is not None:
    81                 s = 'Specific energy  = %f m' % self.inflow.get_average_specific_energy()
     81                s = 'Specific energy  = %f m' % self.inflow.get_enquiry_specific_energy()
    8282                log_to_file(self.log_filename, s)
    8383
    8484            msg = 'Specific energy at inlet is negative'
    85             assert self.inflow.get_average_specific_energy() >= 0.0, msg
     85            assert self.inflow.get_enquiry_specific_energy() >= 0.0, msg
    8686
    8787            height = self.culvert_height
     
    8989            flow_width = self.culvert_width
    9090
    91             Q_inlet_unsubmerged = 0.540*g**0.5*width*self.inflow.get_average_specific_energy()**1.50 # Flow based on Inlet Ctrl Inlet Unsubmerged
    92             Q_inlet_submerged = 0.702*g**0.5*width*height**0.89*self.inflow.get_average_specific_energy()**0.61  # Flow based on Inlet Ctrl Inlet Submerged
     91            Q_inlet_unsubmerged = 0.540*g**0.5*width*self.inflow.get_enquiry_specific_energy()**1.50 # Flow based on Inlet Ctrl Inlet Unsubmerged
     92            Q_inlet_submerged = 0.702*g**0.5*width*height**0.89*self.inflow.get_enquiry_specific_energy()**0.61  # Flow based on Inlet Ctrl Inlet Submerged
    9393
    9494            # FIXME(Ole): Are these functions really for inlet control?
     
    120120                case = 'INLET CTRL Culvert is open channel flow we will for now assume critical depth'
    121121
    122             if self.delta_total_energy < self.inflow.get_average_specific_energy():
     122            if self.delta_total_energy < self.inflow.get_enquiry_specific_energy():
    123123                # Calculate flows for outlet control
    124124
    125125                # Determine the depth at the outlet relative to the depth of flow in the Culvert
    126                 if self.outflow.get_average_height() > height:        # The Outlet is Submerged
     126                if self.outflow.get_enquiry_height() > height:        # The Outlet is Submerged
    127127                    outlet_culvert_depth=height
    128128                    flow_area=width*height       # Cross sectional area of flow in the culvert
     
    172172        # END CODE BLOCK for DEPTH  > Required depth for CULVERT Flow
    173173
    174         else: # self.inflow.get_average_height() < 0.01:
     174        else: # self.inflow.get_enquiry_height() < 0.01:
    175175            Q = barrel_velocity = outlet_culvert_depth = 0.0
    176176
Note: See TracChangeset for help on using the changeset viewer.