Changeset 5453

Timestamp:

Jun 30, 2008, 4:44:52 PM (16 years ago)

Author:

ole

Message:

Fixed georeferencing problem in culvert class

File:

• anuga_core/source/anuga/culvert_flows/culvert_class.py (modified) (1 diff)

anuga_core/source/anuga/culvert_flows/culvert_class.py

@@ -215 +215 @@
 
             # Indices corresponding to energy enquiry field for this opening
-            coordinates = domain.get_centroid_coordinates() # Get all centroid points (x,y)
+            coordinates = domain.get_centroid_coordinates(absolute=True) # Get all centroid points (x,y)
             enquiry_indices = inside_polygon(coordinates, self.enquiry_polygons[i]) 
-            
+
+            if len(enquiry_indices) == 0:
+                msg = 'No triangles have been identified in specified region: %s' %str(self.enquiry_polygons[i])
+                raise Exception, msg                
             
             # Get model values for points in enquiry polygon for this opening