Changeset 8049 for trunk/anuga_core/source/anuga/structures/structure_operator.py

Timestamp:

Oct 22, 2010, 10:36:57 AM (13 years ago)

Author:

habili

Message:

Changed from using polygons to polylines. The inlets are now represented by triangles that intersect or contain a polyline - The enquiry gap is now computed as: gap = (self.apron + self.enquiry_gap)*self.culvert_vector

File:

• trunk/anuga_core/source/anuga/structures/structure_operator.py (modified) (4 diffs)

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

@@ -83 +83 @@
         self.driving_energy = 0.0
         
-        self.__create_exchange_polygons()
+        self.__create_exchange_polylines()
 
         self.inlets = []
@@ -194 +194 @@
 
 
-    def __create_exchange_polygons(self):
+    def __create_exchange_polylines(self):
 
         """Create polylines at the end of a culvert inlet and outlet.
@@ -222 +222 @@
 
         #gap = 1.5*h + self.enquiry_gap
-        gap = 1.5*self.culvert_vector + self.enquiry_gap
+        gap = (self.apron+ self.enquiry_gap)*self.culvert_vector
 
         self.inlet_polylines = []
         self.inlet_equiry_points = []
 
-        # Build exchange polygon and enquiry point
+        # Build exchange polyline and enquiry point
         for i in [0, 1]:
             i0 = (2*i-1) #i0 determines the sign of the points
@@ -267 +267 @@
             message += '\n'
 
-            message += 'polygon\n'
+            message += 'polyline\n'
             message += '%s' % inlet.polyline
             message += '\n'