Changeset 6006

Timestamp:

Nov 26, 2008, 6:15:27 PM (16 years ago)

Author:

ole

Message:

More error messages

File:

• anuga_core/source/anuga/shallow_water/shallow_water_domain.py (modified) (3 diffs)

anuga_core/source/anuga/shallow_water/shallow_water_domain.py

@@ -1691 +1691 @@
                          # previous timestep in order to obtain rate
 
-        bounding_polygon = domain.get_boundary_polygon() # Returns absolute coordinates
+        # Get boundary (in absolute coordinates)
+        bounding_polygon = domain.get_boundary_polygon()
 
 
@@ -1728 +1729 @@
         
         if polygon is not None:
-            self.exchange_area = polygon_area(self.polygon)
 
             # Check that polygon lies within the mesh.
             for point in self.polygon:
-                msg = 'Point %s in polygon for forcing term' %(str(point))
+                msg = 'Point %s in polygon for forcing term' %(point)
                 msg += ' did not fall within the domain boundary.'
                 assert is_inside_polygon(point, bounding_polygon), msg
-
+                
+            # Compute area and check that it is greater than 0    
+            self.exchange_area = polygon_area(self.polygon)
+            
+            msg = 'Polygon %s in forcing term' %(self.polygon)
+            msg += ' has area = %f' %self.exchange_area
+            assert self.exchange_area > 0.0            
+            
+                
+                
 
         # Pointer to update vector
@@ -1762 +1771 @@
             # Inlet is polygon
             
-            inlet_region = 'polygon=%s' %(self.polygon) 
+            inlet_region = 'polygon=%s, area=%f m^2' %(self.polygon, 
+                                                       self.exchange_area) 
                         
             self.exchange_indices = inside_polygon(points, self.polygon)
+