Changeset 8050 for trunk/anuga_core/source/anuga/structures/testing_culvert.py

Timestamp:

Oct 22, 2010, 5:10:19 PM (14 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:

• trunk/anuga_core/source/anuga/structures/testing_culvert.py (modified) (3 diffs)

trunk/anuga_core/source/anuga/structures/testing_culvert.py

@@ -11 +11 @@
 
 from anuga.structures.boyd_pipe_operator import Boyd_pipe_operator
+from anuga.structures.inlet_operator import Inlet_operator
                             
 #from anuga.culvert_flows.culvert_routines import boyd_generalised_culvert_model
@@ -101 +102 @@
 
 
-#culvert2 = Generic_box_culvert(domain,
-#                              end_point0=[19.0, 2.5],
-#                              end_point1=[25.0, 2.5],
-#                              width=1.00,
-#                              verbose=False)
+line = [[0.0, 5.0], [0.0, 10.0]]
+Q = 5.0
+#Inlet_operator(domain, line, Q)
 
 
 
-
-#print domain.fractional_step_operators
-
-#domain.apply_fractional_steps()
 
 ##-----------------------------------------------------------------------
@@ -143 +138 @@
     domain.write_time()
 
-    if domain.get_time() > 150.5 and domain.get_time() < 151.5 :
-        Bi = anuga.Dirichlet_boundary([0.0, 0.0, 0.0])
-        domain.set_boundary({'left': Bi, 'right': Br, 'top': Br, 'bottom': Br})
+    #if domain.get_time() > 150.5 and domain.get_time() < 151.5 :
+        #Bi = anuga.Dirichlet_boundary([0.0, 0.0, 0.0])
+        #domain.set_boundary({'left': Bi, 'right': Br, 'top': Br, 'bottom': Br})
 
     #delta_w = culvert.inlet.stage - culvert.outlet.stage