Ignore:
Timestamp:
Oct 22, 2010, 5:10:19 PM (12 years ago)
Author:
steve
Message:

Added in an inlet_operator class to implement inflow "boundary condition" by setting up a inlet defined by a line close to the boundary. Look at testing_culvert_inlet.py as an example.

File:
1 edited

Legend:

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

    r8007 r8050  
    1111
    1212from anuga.structures.boyd_pipe_operator import Boyd_pipe_operator
     13from anuga.structures.inlet_operator import Inlet_operator
    1314                           
    1415#from anuga.culvert_flows.culvert_routines import boyd_generalised_culvert_model
     
    101102
    102103
    103 #culvert2 = Generic_box_culvert(domain,
    104 #                              end_point0=[19.0, 2.5],
    105 #                              end_point1=[25.0, 2.5],
    106 #                              width=1.00,
    107 #                              verbose=False)
     104line = [[0.0, 5.0], [0.0, 10.0]]
     105Q = 5.0
     106#Inlet_operator(domain, line, Q)
    108107
    109108
    110109
    111 
    112 #print domain.fractional_step_operators
    113 
    114 #domain.apply_fractional_steps()
    115110
    116111##-----------------------------------------------------------------------
     
    143138    domain.write_time()
    144139
    145     if domain.get_time() > 150.5 and domain.get_time() < 151.5 :
    146         Bi = anuga.Dirichlet_boundary([0.0, 0.0, 0.0])
    147         domain.set_boundary({'left': Bi, 'right': Br, 'top': Br, 'bottom': Br})
     140    #if domain.get_time() > 150.5 and domain.get_time() < 151.5 :
     141        #Bi = anuga.Dirichlet_boundary([0.0, 0.0, 0.0])
     142        #domain.set_boundary({'left': Bi, 'right': Br, 'top': Br, 'bottom': Br})
    148143
    149144    #delta_w = culvert.inlet.stage - culvert.outlet.stage
Note: See TracChangeset for help on using the changeset viewer.