Changeset 7984 for trunk/anuga_core/source/anuga/structures/testing_wide_bridge.py

Timestamp:

Sep 1, 2010, 6:11:25 PM (13 years ago)

Author:

steve

Message:

Added enquiry points to culverts

File:

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

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

@@ -46 +46 @@
 
 points, vertices, boundary = rectangular_cross(int(length/dx), int(width/dy),
-                                               len1=length, len2=width)
+                                                    len1=length, len2=width)
 domain = Domain(points, vertices, boundary)   
 domain.set_name('Test_WIDE_BRIDGE')                 # Output name
@@ -168 +168 @@
 """
 from anuga.structures.culvert_operator import Culvert_operator
+#culvert0 = Culvert_operator(domain,
+#                            end_point0=[40.0, 75.0],
+#                            end_point1=[50.0, 75.0],
+#                            width=50.0,
+#                            height=10.0,
+#                            apron=5.0,
+#                            verbose=False)
+
+
 culvert1 = Culvert_operator(domain,
-                            end_point0=[40.0, 75.0],
-                            end_point1=[50.0, 75.0],
-                            width=50.0,
+                            end_point0=[40.0, 87.5],
+                            end_point1=[50.0, 87.5],
+                            width=25.0,
                             height=10.0,
+                            apron=5.0,
+                            verbose=False)
+
+
+culvert2 = Culvert_operator(domain,
+                            end_point0=[40.0, 62.5],
+                            end_point1=[50.0, 62.5],
+                            width=25.0,
+                            height=10.0,
+                            apron=5.0,
                             verbose=False)
 
@@ -220 +239 @@
 for t in domain.evolve(yieldstep = 1, finaltime = 100):
     print domain.timestepping_statistics()
+    print domain.volumetric_balance_statistics()