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

Timestamp:

Aug 20, 2010, 5:27:03 PM (14 years ago)

Author:

steve

Message:

Got a flow from one region to another.

File:

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

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

@@ -33 +33 @@
 width = 5.
 
-dx = dy = 1           # Resolution: Length of subdivisions on both axes
+dx = dy = 0.5          # Resolution: Length of subdivisions on both axes
 
 points, vertices, boundary = rectangular_cross(int(length/dx),
@@ -85 +85 @@
 
 filename=os.path.join(path, 'example_rating_curve.csv')
-culvert = Generic_box_culvert(domain,
+culvert1 = Generic_box_culvert(domain,
                               end_point0=[9.0, 2.5], 
                               end_point1=[13.0, 2.5],
@@ -91 +91 @@
                               verbose=False)
 
-#domain.forcing_terms.append(culvert)
 
+#culvert2 = Generic_box_culvert(domain,
+#                              end_point0=[19.0, 2.5],
+#                              end_point1=[25.0, 2.5],
+#                              width=1.00,
+#                              verbose=False)
+
+
+
+
+#print domain.fractional_step_operators
+
+#domain.apply_fractional_steps()
 
 ##-----------------------------------------------------------------------
@@ -99 +110 @@
 
 ## Inflow based on Flow Depth and Approaching Momentum
-#Bi = anuga.Dirichlet_boundary([0.0, 0.0, 0.0])
-#Br = anuga.Reflective_boundary(domain)              # Solid reflective wall
+Bi = anuga.Dirichlet_boundary([1.0, 0.0, 0.0])
+Br = anuga.Reflective_boundary(domain)              # Solid reflective wall
 #Bo = anuga.Dirichlet_boundary([-5, 0, 0])           # Outflow
 
@@ -111 +122 @@
             #lambda t: [-5*(num.cos(2*pi*(t-4)/20)), 0.0, 0.0])
 #domain.set_boundary({'left': Btus, 'right': Btds, 'top': Br, 'bottom': Br})
+domain.set_boundary({'left': Bi, 'right': Br, 'top': Br, 'bottom': Br})
 
 
@@ -119 +131 @@
 #min_delta_w = sys.maxint 
 #max_delta_w = -min_delta_w
-#for t in domain.evolve(yieldstep = 1, finaltime = 25):
+for t in domain.evolve(yieldstep = 1.0, finaltime = 300):
+    domain.write_time()
     #delta_w = culvert.inlet.stage - culvert.outlet.stage
     
@@ -125 +138 @@
     #if delta_w < min_delta_w: min_delta_w = delta_w            
     
-    #pass
+    pass
 
 ## Check that extreme values in rating curve have been exceeded