Changeset 9392

Timestamp:

Jan 13, 2015, 2:06:42 PM (10 years ago)

Author:

davies

Message:

More informative error message in quantity_setting_functions

File:

• trunk/anuga_core/source/anuga/utilities/quantity_setting_functions.py (modified) (1 diff)

trunk/anuga_core/source/anuga/utilities/quantity_setting_functions.py

@@ -369 +369 @@
                     quantityVal[fInds], upper_bound)
 
+        # End of loop
+
         if( min(isSet) != 1):
-            raise Exception('Some points were not inside any polygon')
+            print 'Some points were not inside any polygon, ',\
+                  'or evalute to nan over all datasets'
+            unset_inds = (isSet!=1).nonzero()[0]
+            lui = min(5, length(unset_inds)) 
+            print 'There are ', length(unset_inds), ' such points'
+            print 'Here are a few:'
+            for i in range(lui):
+                print x[i]+xll, y[i]+yll
+            raise Exception('It seems the input data needs to be fixed')
 
         return quantityVal