Changeset 8859

Timestamp:

May 9, 2013, 12:17:00 PM (12 years ago)

Author:

steve

Message:

Adding facility to change culvert_height

Location:

trunk/anuga_core/source/anuga/structures

Files:

• run_gate_operator.py (modified) (7 diffs)
• structure_operator.py (modified) (1 diff)

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

@@ -10 +10 @@
 import anuga
 
-from anuga.structures.boyd_pipe_operator import Boyd_pipe_operator
+from anuga.structures.boyd_box_operator import Boyd_box_operator
 from anuga.structures.inlet_operator import Inlet_operator
                             
@@ -41 +41 @@
                                                len2=width)
 domain = anuga.Domain(points, vertices, boundary)   
-domain.set_name('run_culvert_inlet')                 # Output name
+domain.set_name('run_gate_operator')                 # Output name
 domain.set_default_order(2)
 #domain.set_beta(1.5)
@@ -89 +89 @@
 
 
-
-Boyd_pipe_operator(domain,
+gate = Boyd_box_operator(domain,
                             end_points=[[9.0, 2.5],[13.0, 2.5]],
                             losses=1.5,
-                            diameter=1.5,
+                            width=1.5,
                             apron=5.0,
                             use_momentum_jet=True,
@@ -101 +100 @@
 
 
+gate.set_culvert_height(10.0)
+
 line = [[0.0, 5.0], [0.0, 10.0]]
-Q = 5.0
+Q = 1.0
 Inlet_operator(domain, line, Q)
 
@@ -113 +114 @@
 
 ## Inflow based on Flow Depth and Approaching Momentum
-Bi = anuga.Dirichlet_boundary([2.0, 0.0, 0.0])
 Br = anuga.Reflective_boundary(domain)              # Solid reflective wall
-#Bo = anuga.Dirichlet_boundary([-5, 0, 0])           # Outflow
 
-## Upstream and downstream conditions that will exceed the rating curve
-## I.e produce delta_h outside the range [0, 10] specified in the the 
-## file example_rating_curve.csv
-#Btus = anuga.Time_boundary(domain, \
-            #lambda t: [100*num.sin(2*pi*(t-4)/10), 0.0, 0.0])
-#Btds = anuga.Time_boundary(domain, \
-            #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': Br, 'right': Br, 'top': Br, 'bottom': Br})
 
@@ -134 +125 @@
 #min_delta_w = sys.maxint 
 #max_delta_w = -min_delta_w
-for t in domain.evolve(yieldstep = 1.0, finaltime = 200):
+for t in domain.evolve(yieldstep = 1.0, finaltime = 50):
     domain.write_time()
+
+
+    if num.allclose(t, 10.0):
+        gate.set_culvert_height(0.000001)
+
+    Q, velocity, depth = gate.discharge_routine()
+
+    print gate.culvert_height
+    print Q
+    print velocity
+    print depth
+        
 
     #if domain.get_time() > 150.5 and domain.get_time() < 151.5 :
@@ -146 +149 @@
     #if delta_w < min_delta_w: min_delta_w = delta_w
 
-    print domain.volumetric_balance_statistics()
+    #print domain.volumetric_balance_statistics()
     
     pass

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

@@ -210 +210 @@
 
 
-    def set_height(self, height):
-
-        self.height = height
-
-    def set_width(self, width):
-
-        self.width = width
+    def set_culvert_height(self, height):
+
+        self.culvert_height = height
+
+    def set_culvert_width(self, width):
+
+        self.culvert_width = width