Changeset 9390

Timestamp:

Jan 8, 2015, 3:48:15 PM (10 years ago)

Author:

steve

Message:

Updated region.py to produce a warning if region doesn't intersect domain

Location:

trunk/anuga_core

Files:

• source/anuga/abstract_2d_finite_volumes/region.py (modified) (1 diff)
• source/anuga/fit_interpolate/interpolate.py (modified) (2 diffs)
• validation_tests/case_studies/merewether/plot_hydrograph.py (modified) (1 diff)

trunk/anuga_core/source/anuga/abstract_2d_finite_volumes/region.py

@@ -241 +241 @@
         
         if not self.domain.parallel:
-            msg = 'No centroids found for polygon %s '% str(self.polygon)
-            if len(indices) is 0: raise Exception(msg)
+            # only warn if not parallel as we should get lots of subdomains without indices
+            if len(indices) is 0:
+                msg = 'No centroids found for polygon %s '% str(self.polygon)
+                import warnings
+                warnings.warn(msg)
 
             

trunk/anuga_core/source/anuga/fit_interpolate/interpolate.py

@@ -914 +914 @@
             for i, t in enumerate(self.time):
                 # Interpolate quantities at this timestep
-                if verbose and i%((p+10)/10) == 0:
+                #if verbose and i%((p+10)/10) == 0:
+                if verbose:
                     log.critical('  time step %d of %d' % (i, p))
 
@@ -923 +924 @@
                         Q = quantities[name][:]   # No time dependency
 
-                    if verbose and i%((p+10)/10) == 0:
+                    #if verbose and i%((p+10)/10) == 0:
+                    if verbose:
                         log.critical('    quantity %s, size=%d' % (name, len(Q)))
 

trunk/anuga_core/validation_tests/case_studies/merewether/plot_hydrograph.py

@@ -100 +100 @@
 
 
-def get_interpolated_quantities_at_polyline_midpoints(filename,
-                                                      quantity_names=None,
-                                                      polyline=None,
-                                                      verbose=False):
-
-
-
-
-    return segments, I
-
-
 #===================================================================================