Changeset 8541 for trunk/anuga_core/source/anuga/file_conversion/sww2dem.py

Timestamp:

Aug 30, 2012, 10:35:14 AM (13 years ago)

Author:

steve

Message:

Dealing with holes

File:

• trunk/anuga_core/source/anuga/file_conversion/sww2dem.py (modified) (2 diffs)

trunk/anuga_core/source/anuga/file_conversion/sww2dem.py

@@ -330 +330 @@
     if verbose: log.critical('Interpolating')
     grid_values = interp.interpolate(result, grid_points).flatten()
+    outside_indices = interp.get_outside_poly_indices()
 
     if verbose:
@@ -336 +337 @@
 
     # Assign NODATA_value to all points outside bounding polygon (from interpolation mesh)
-    P = interp.mesh.get_boundary_polygon()
-    outside_indices = outside_polygon(grid_points, P, closed=True)
+    
+#    P = interp.mesh.get_boundary_polygon()
+#    outside_indices = outside_polygon(grid_points, P, closed=True)
 
     for i in outside_indices: