Changeset 8095


Ignore:
Timestamp:
Dec 7, 2010, 10:56:18 PM (8 years ago)
Author:
habili
Message:

changing of water height to depth
Fixes for instability

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

Legend:

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

    r8073 r8095  
    3636        anuga.Structure_operator.__init__(self,
    3737                                          domain,
    38                                           width,
    39                                           height,
    4038                                          end_points,
    4139                                          exchange_lines,
    4240                                          enquiry_points,
     41                                          width,
     42                                          height,
    4343                                          apron,
    4444                                          manning,
     
    4949                                          logging,
    5050                                          verbose)     
    51        
    5251       
    5352        if type(losses) == types.DictType:
  • trunk/anuga_core/source/anuga/structures/structure_operator.py

    r8073 r8095  
    2424    def __init__(self,
    2525                 domain,
    26                  width,
    27                  height,
    2826                 end_points,
    2927                 exchange_lines,
    3028                 enquiry_points,
     29                 width,
     30                 height,
    3131                 apron,
    3232                 manning,
     
    246246            self.culvert_vector = self.end_points[1] - end_points[0]
    247247       
    248         self.culvert_length = math.sqrt(num.sum(self.culvert_vector**2))   
     248        self.culvert_length = math.sqrt(num.sum(self.culvert_vector**2))
    249249        assert self.culvert_length > 0.0, 'The length of culvert is less than 0'
    250250       
  • trunk/anuga_core/source/anuga/structures/testing_outlet_control.py

    r8073 r8095  
    186186
    187187    Boyd_box_operator(domain,
    188                             end_points=[ep0,ep1],
    189                             losses=losses,
    190                             width=culvert_width,
    191                             depth=10.0,
    192                             apron=6.0,
    193                             use_momentum_jet=True,
    194                             use_velocity_head=True,
    195                             manning=0.013,
    196                             logging=True,
    197                             label='box_culvert',
    198                             verbose=False)
    199 
    200                        
     188                      losses=losses,
     189                      width=culvert_width,
     190                      height=10.0,
     191                      end_points=[ep0,ep1],
     192                      apron=6.0,
     193                      manning=0.013,
     194                      label='box_culvert',
     195                      logging=True,
     196                      verbose=False)
    201197
    202198#losses = {'inlet':1, 'outlet':1, 'bend':1, 'grate':1, 'pier': 1, 'other': 1}
  • trunk/anuga_core/source/anuga/structures/testing_wide_bridge.py

    r8073 r8095  
    184184    culvert_width = 50.0/number_of_culverts
    185185    y = 100-i*culvert_width - culvert_width/2.0
    186     ep0 = [40.0, y]
    187     ep1 = [50.0, y]
     186    ep0 = num.array([40.0, y])
     187    ep1 = num.array([50.0, y])
    188188    culverts.append(anuga.Boyd_box_operator(domain,
    189                             end_point0=ep0,
    190                             end_point1=ep1,
    191                             losses=1.5,
    192                             width=3.658,
    193                             depth=3.658,
    194                             apron=5.0,
    195                             use_momentum_jet=True,
    196                             use_velocity_head=True,
    197                             manning=0.013,
    198                             description='bridge culvert',
    199                             verbose=False))
    200 
    201 
    202 
     189                                            losses=1.5,
     190                                            width=3.658,
     191                                            height=3.658,
     192                                            end_points=[ep0, ep1],
     193                                            apron=0.5,
     194                                            manning=0.013,
     195                                            enquiry_gap=1.0,
     196                                            description='bridge culvert',
     197                                            verbose=False))
    203198
    204199#culvert2 = Culvert_operator(domain,
     
    260255    for i, culvert in enumerate(culverts):
    261256        print 'culvert: ', i
    262         print culvert.structure_statistics()
     257        print culvert.statistics()
    263258   
    264259
Note: See TracChangeset for help on using the changeset viewer.