Changeset 8027


Ignore:
Timestamp:
Sep 30, 2010, 10:01:46 AM (14 years ago)
Author:
habili
Message:

Added control case text to print_timestepping_statistics

Location:
trunk/anuga_core/source/anuga/structures
Files:
3 edited

Legend:

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

    r8020 r8027  
    7373        self.discharge = 0.0
    7474        self.velocity = 0.0
     75       
     76        self.case = 'N/A'
    7577
    7678   
     
    117119                    perimeter= 2.0*dcrit+width
    118120                outlet_culvert_depth = dcrit
    119                 case = 'Inlet unsubmerged Box Acts as Weir'
     121                self.case = 'Inlet unsubmerged Box Acts as Weir'
    120122            else: # Inlet Submerged but check internal culvert flow depth
    121123                Q = Q_inlet_submerged
     
    129131                    perimeter= 2.0*dcrit+width
    130132                outlet_culvert_depth = dcrit
    131                 case = 'Inlet submerged Box Acts as Orifice'
     133                self.case = 'Inlet submerged Box Acts as Orifice'
    132134
    133135            dcrit = (Q**2/anuga.g/width**2)**0.333333
     
    138140                flow_area = width*height  # Cross sectional area of flow in the culvert
    139141                perimeter = 2*(width+height)
    140                 case = 'Inlet CTRL Outlet unsubmerged PIPE PART FULL'
     142                self.case = 'Inlet CTRL Outlet unsubmerged PIPE PART FULL'
    141143            else:
    142144                flow_area = width * outlet_culvert_depth
    143145                perimeter = width+2*outlet_culvert_depth
    144                 case = 'INLET CTRL Culvert is open channel flow we will for now assume critical depth'
     146                self.case = 'INLET CTRL Culvert is open channel flow we will for now assume critical depth'
    145147            # Initial Estimate of Flow for Outlet Control using energy slope
    146148            #( may need to include Culvert Bed Slope Comparison)
     
    158160                    flow_area=width*height       # Cross sectional area of flow in the culvert
    159161                    perimeter=2.0*(width+height)
    160                     case = 'Outlet submerged'
     162                    self.case = 'Outlet submerged'
    161163                else:   # Here really should use the Culvert Slope to calculate Actual Culvert Depth & Velocity
    162164                    dcrit = (Q**2/anuga.g/width**2)**0.333333
     
    166168                        flow_area=width*height
    167169                        perimeter=2.0*(width+height)
    168                         case = 'Outlet is Flowing Full'
     170                        self.case = 'Outlet is Flowing Full'
    169171                    else:
    170172                        flow_area=width*outlet_culvert_depth
    171173                        perimeter=(width+2.0*outlet_culvert_depth)
    172                         case = 'Outlet is open channel flow'
     174                        self.case = 'Outlet is open channel flow'
    173175
    174176                hyd_rad = flow_area/perimeter
  • trunk/anuga_core/source/anuga/structures/boyd_pipe_operator.py

    r8018 r8027  
    7171        self.velocity = 0.0
    7272       
     73        self.case = 'N/A'
     74       
    7375   
    7476    def discharge_routine(self):
     
    140142                perimeter = diameter * math.pi
    141143                flow_width= diameter
    142                 case = 'Inlet CTRL Outlet submerged Circular PIPE FULL'
     144                self.case = 'Inlet CTRL Outlet submerged Circular PIPE FULL'
    143145                if local_debug == 'true':
    144146                    anuga.log.critical('Inlet CTRL Outlet submerged Circular '
     
    151153                flow_width= diameter*math.sin(alpha/2.0)
    152154                perimeter = alpha*diameter/2.0
    153                 case = 'INLET CTRL Culvert is open channel flow we will for now assume critical depth'
     155                self.case = 'INLET CTRL Culvert is open channel flow we will for now assume critical depth'
    154156                if local_debug =='true':
    155157                    anuga.log.critical('INLET CTRL Culvert is open channel flow '
     
    167169                    perimeter = diameter * math.pi
    168170                    flow_width= diameter
    169                     case = 'Outlet submerged'
     171                    self.case = 'Outlet submerged'
    170172                    if local_debug =='true':
    171173                        anuga.log.critical('Outlet submerged')
     
    183185                        perimeter = diameter * math.pi
    184186                        flow_width= diameter
    185                         case = 'Outlet unsubmerged PIPE FULL'
     187                        self.case = 'Outlet unsubmerged PIPE FULL'
    186188                        if local_debug =='true':
    187189                            anuga.log.critical('Outlet unsubmerged PIPE FULL')
     
    191193                        flow_width= diameter*math.sin(alpha/2.0)
    192194                        perimeter = alpha*diameter/2.0
    193                         case = 'Outlet is open channel flow we will for now assume critical depth'
     195                        self.case = 'Outlet is open channel flow we will for now assume critical depth'
    194196                        if local_debug == 'true':
    195197                            anuga.log.critical('Q Outlet Depth and ALPHA = %s, %s, %s'
  • trunk/anuga_core/source/anuga/structures/structure_operator.py

    r8024 r8027  
    307307        message += 'Inlet Driving Energy %.2f\n' % self.driving_energy
    308308        message += 'Delta Total Energy %.2f\n' % self.delta_total_energy
     309        message += 'Control at this instant: %s\n' % self.case
    309310
    310311        print message
Note: See TracChangeset for help on using the changeset viewer.