Changeset 5450 for anuga_core/source/anuga/shallow_water/shallow_water_domain.py

Timestamp:

Jun 30, 2008, 10:41:27 AM (15 years ago)

Author:

ole

Message:

Added exception is inlet_region for General flow does not contain any triangles.
An input check was also added.

File:

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

anuga_core/source/anuga/shallow_water/shallow_water_domain.py

@@ -1491 +1491 @@
                  verbose=False):
                      
-
+        if center is None:
+            msg = 'I got radius but no center.'        
+            assert radius is None, msg
+            
+        if radius is None:
+            msg += 'I got center but no radius.'        
+            assert center is None, msg
+            
+            
+                     
         from math import pi
 
@@ -1529 +1538 @@
             # Inlet is circular
             
+            inlet_region = 'center=%s, radius=%s' %(self.center, self.radius) 
+            
             self.exchange_indices = []
             for k in range(N):
@@ -1537 +1548 @@
         if self.polygon is not None:                    
             # Inlet is polygon
+            
+            inlet_region = 'polygon=%s' %(self.polygon) 
+                        
             self.exchange_indices = inside_polygon(points, self.polygon)
             
-
-            
+            
+        if self.exchange_indices is not None:
+            #print inlet_region
+        
+            if len(self.exchange_indices) == 0:
+                msg = 'No triangles have been identified in specified region: %s' %inlet_region
+                raise Exception, msg