Changeset 6223 for anuga_core/source/anuga/utilities/polygon.py

Timestamp:

Jan 21, 2009, 4:37:45 PM (15 years ago)

Author:

ole

Message:

Various cleanup

File:

• anuga_core/source/anuga/utilities/polygon.py (modified) (4 diffs)

anuga_core/source/anuga/utilities/polygon.py

@@ -661 +661 @@
     """
 
-    def __init__(self, regions, default=0.0, geo_reference=None):
+    def __init__(self,
+                 regions,
+                 default=0.0,
+                 geo_reference=None):
 
         try:
@@ -702 +705 @@
         for polygon, value in regions:
             P = geo_reference.change_points_geo_ref(polygon)
-            self.regions.append( (P, value) )
+            self.regions.append((P, value))
 
 
@@ -714 +717 @@
         assert len(y) == N
 
-        points = num.concatenate( (num.reshape(x, (N, 1)),
-                                   num.reshape(y, (N, 1))), axis=1 )
+        points = num.concatenate((num.reshape(x, (N, 1)),
+                                  num.reshape(y, (N, 1))), axis=1)
 
         if callable(self.default):
@@ -735 +738 @@
                     z[i] = value
 
+        if len(z) == 0:
+            msg = 'Warning: points provided to Polygon function did not fall within'
+            msg += 'its regions'
+            msg += 'x in [%.2f, %.2f], y in [%.2f, %.2f]' % (min(x), max(x),
+                                                             min(y), max(y))
+            print msg
+
+            
         return z