Changeset 8097


Ignore:
Timestamp:
Dec 8, 2010, 1:54:41 PM (9 years ago)
Author:
habili
Message:

Updating of code and test case

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

Legend:

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

    r8073 r8097  
    3838                                          enquiry_points,
    3939                                          width=diameter,
    40                                           depth=None,
     40                                          height=None,
    4141                                          apron=apron,
    4242                                          manning=manning,
     
    7373        self.case = 'N/A'
    7474       
    75 
    76     def __determine_inflow_outflow(self):
    77         # Determine flow direction based on total energy difference
    78 
     75    def discharge_routine(self):
     76
     77        #self.__determine_inflow_outflow()
     78
     79        local_debug ='false'
     80       
    7981        if self.use_velocity_head:
    8082            self.delta_total_energy = self.inlets[0].get_enquiry_total_energy() - self.inlets[1].get_enquiry_total_energy()
     
    9092            self.delta_total_energy = -self.delta_total_energy
    9193
    92     def discharge_routine(self):
    93 
    94         self.__determine_inflow_outflow()
    95 
    96         local_debug ='false'
    97        
    98         #import pdb
    99         #pdb.set_trace()
    100        
    10194        if self.inflow.get_enquiry_depth() > 0.01: #this value was 0.01: Remember this needs to be compared to the Invert Lvl
    10295            if local_debug =='true':
  • trunk/anuga_core/source/anuga/structures/testing_culvert.py

    r8050 r8097  
    4141                                               len2=width)
    4242domain = anuga.Domain(points, vertices, boundary)   
     43domain.set_starttime(10)
    4344domain.set_name('Test_culvert')                 # Output name
    4445domain.set_default_order(2)
     
    8889filename=os.path.join(path, 'example_rating_curve.csv')
    8990
    90 
     91end_point0 = num.array([9.0, 2.5])
     92end_point1 = num.array([13.0, 2.5])
    9193
    9294Boyd_pipe_operator(domain,
    93                             end_point0=[9.0, 2.5],
    94                             end_point1=[13.0, 2.5],
    95                             losses=1.5,
    96                             diameter=1.5,
    97                             apron=5.0,
    98                             use_momentum_jet=True,
    99                             use_velocity_head=False,
    100                             manning=0.013,
    101                             verbose=False)
     95                    #end_point0=[9.0, 2.5],
     96                    #end_point1=[13.0, 2.5],
     97                    #exchange_line0=[[9.0, 1.75],[9.0, 3.25]],
     98                    #exchange_line1=[[13.0, 1.75],[13.0, 3.25]],
     99                    losses=1.5,
     100                    end_points=[end_point0, end_point1],
     101                    diameter=1.5,
     102                    apron=0.5,
     103                    use_momentum_jet=True,
     104                    use_velocity_head=False,
     105                    manning=0.013,
     106                    verbose=False)
    102107
    103108
     
    135140#min_delta_w = sys.maxint
    136141#max_delta_w = -min_delta_w
    137 for t in domain.evolve(yieldstep = 1.0, finaltime = 200):
     142for t in domain.evolve(yieldstep=1.0, finaltime=50.0):
    138143    domain.write_time()
    139144
    140     #if domain.get_time() > 150.5 and domain.get_time() < 151.5 :
     145    #if domain.get_time() > 150.5 and domain.ge t_time() < 151.5 :
    141146        #Bi = anuga.Dirichlet_boundary([0.0, 0.0, 0.0])
    142147        #domain.set_boundary({'left': Bi, 'right': Br, 'top': Br, 'bottom': Br})
Note: See TracChangeset for help on using the changeset viewer.