Changeset 7998


Ignore:
Timestamp:
Sep 3, 2010, 12:20:38 PM (8 years ago)
Author:
habili
Message:

deal with boyd pipes

Location:
trunk/anuga_core/source/anuga/structures
Files:
1 added
2 edited

Legend:

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

    r7996 r7998  
    199199       
    200200       
     201    def get_culvert_diameter(self):
     202   
     203            return self.width
     204       
     205       
    201206    def get_culvert_height(self):
    202207   
  • trunk/anuga_core/source/anuga/structures/test_Outlet_Ctrl.py

    r7997 r7998  
    2525
    2626#from anuga.culvert_flows.culvert_class import Culvert_flow
    27 from anuga.structures.boyd_box_operator import Boyd_box_operator
     27from anuga.structures.boyd_pipe_operator import Boyd_pipe_operator
    2828#from anuga.culvert_flows.culvert_routines import weir_orifice_channel_culvert_model
    2929from math import pi,pow,sqrt
     
    165165    ep0 = [40.0, y]
    166166    ep1 = [50.0, y]
    167     losses = {'inlet':1, 'outlet':1, 'bend':1, 'grate':1, 'pier': 1, 'other': 1}
    168     culverts.append(Boyd_box_operator(domain,
     167    losses = {'inlet':0.5, 'outlet':1, 'bend':0, 'grate':0, 'pier': 0, 'other': 0}
     168    culverts.append(Boyd_pipe_operator(domain,
    169169                            end_point0=ep0,
    170170                            end_point1=ep1,
    171                             losses,
    172                             width=3.658, #culvert_width, #3.658,
    173                             height=3.658,
     171                            losses=losses,
     172                            diameter=1.5, #culvert_width, #3.658,
    174173                            apron=6.0,
    175174                            use_momentum_jet=True,
     
    178177                            verbose=False))
    179178
    180 
    181 
    182 losses = {'inlet':1, 'outlet':1, 'bend':1, 'grate':1, 'pier': 1, 'other': 1}
    183 culvert2 = Culvert_operator(domain,
    184                             end_point0=[40.0, 62.5],
    185                             end_point1=[50.0, 62.5],
    186                             losses,
    187                             width=25.0,
    188                             height=10.0,
    189                             apron=5.0,
    190                             manning=0.013,
    191                             verbose=False)
     179                       
     180
     181#losses = {'inlet':1, 'outlet':1, 'bend':1, 'grate':1, 'pier': 1, 'other': 1}
     182#culvert2 = Culvert_operator(domain,
     183                            #end_point0=[40.0, 62.5],
     184                            #end_point1=[50.0, 62.5],
     185                            #losses,
     186                            #width=25.0,
     187                            #height=10.0,
     188                            #apron=5.0,
     189                            #manning=0.013,
     190                            #verbose=False)
    192191
    193192
Note: See TracChangeset for help on using the changeset viewer.