[Anuga-user] Vertice values stored wrong in sww-file?

[S.Kummer <ste...@we...>]

Hi,

something seems to be wrong with storing vertice values in the sww-file.

When I setup a model I use
"domain.set_store_vertices_uniquely(False)".
I also use ["domain.set_quantities_to_be_stored({'elevation':1,...]

Before and during simulation the vertice values (e.g. elevation) are
correct ["domain.get_quantity('elevation',
location="vertices").get_values(location="vertices")].

But when I read the vertice values (eg. elevation, friction, stage) from
the sww-file
["res = get_output(self.sww_file, timeSlices=max)" -> res.elev]
nearly all of them are zero. Only a few vertice values on the
boundary-polygon have values different from zero.

However all the centroid values are right.

The number of vertice values in this case is about 57000.

I have done this many times before in different projects but this is the
first time that the vertice values are stored wrong in the sww-file.

Any ideas?


Thanks
Stefan