Re: [Anuga-user] Problem with water surface elevations artificially high at wet-dry transitions

[Gareth D. <gar...@gm...>]

One more thing in addition to the previous points,

If you're making a raster from the sww outputs, you also need to make 
sure that it will not create "interpolation artefacts" near wet-dry 
boundaries. For instance, by interpolating between wet and dry cells.

I'm not sure whether the raster creation methods you're using will do that.

If you still have problems even after storing vertices uniquely, then 
consider using a different raster-maker. In anuga.utilities.plot_utils 
there is a function Make_Geotiff that can make rasters using 
nearest-neighbour interpolation of centroids. You can suppress 
interpolation completely by setting the "k_nearest_neighbours" argument 
to 1.

See the code in this file: 
https://github.com/GeoscienceAustralia/anuga_core/blob/master/anuga/utilities/plot_utils.py