Changeset 8034

Timestamp:

Oct 14, 2010, 3:20:51 PM (14 years ago)

Author:

habili

Message:

Added in Rudy's amended boyd_box_operator

Location:

trunk/anuga_core/source/anuga/structures

Files:

• boyd_box_operator_Amended3.py (added)
• boyd_pipe_operator.py (modified) (1 diff)
• inlet.py (modified) (2 diffs)
• structure_operator.py (modified) (1 diff)

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

@@ -77 +77 @@
 
         local_debug ='false'
+        
+        #import pdb
+        #pdb.set_trace()
         
         if self.inflow.get_enquiry_height() > 0.01: #this value was 0.01: Remember this needs to be compared to the Invert Lvl

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

@@ -30 +30 @@
 
         # Check that polygon lies within the mesh.
-        for point in self.polygon:
+        for point in self.polygon:  
                 msg = 'Point %s ' %  str(point)
                 msg += ' did not fall within the domain boundary.'
@@ -42 +42 @@
         self.triangle_indices = polygon_overlap(vertex_coordinates, self.polygon)
         #self.triangle_indices_cen = inside_polygon(domain_centroids, self.polygon, verbose=self.verbose)
+        #import pdb
+        #pdb.set_trace()
 
         if len(self.triangle_indices) == 0:

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

@@ -228 +228 @@
         # Build exchange polygon and enquiry point
         for i in [0, 1]:
-            i0 = (2*i-1)
+            i0 = (2*i-1) #i0 determines the sign of the points
             p0 = self.end_points[i] + w
             p1 = self.end_points[i] - w