- Timestamp:
- Jan 17, 2009, 7:36:16 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/abstract_2d_finite_volumes/util.py
r6178 r6189 393 393 if boundary_polygon is not None: 394 394 #removes sts points that do not lie on boundary 395 # FIXME(Ole): Why don't we just remove such points from the list of points and associated data? 396 # I am actually convinced we can get rid of neighbour_gauge_id altogether as the sts file is produced using the ordering file. 397 # All sts points are therefore always present in the boundary. In fact, they *define* parts of the boundary. 395 398 boundary_polygon=ensure_numeric(boundary_polygon) 396 399 boundary_polygon[:,0] -= xllcorner … … 434 437 gauge_neighbour_id=None 435 438 439 #print gauge_neighbour_id 440 436 441 if interpolation_points is not None: 437 442 # Adjust for georef
Note: See TracChangeset
for help on using the changeset viewer.