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

[Stephen R. <ste...@an...>]

Hi Stefan,

Strange. I guess we need to find out if the problem is with the sww file or the res = get_output(self.sww_file, timeSlices=max) code.

I assume you actually use the code

res = get_output(self.sww_file, timeSlices="max")

You could use an application  like ncdump to look at the sww file directly.

Is your problem only with elev, or all the quantities?

Could you provide the actual code that is producing the error.

Cheers
Steve




==============================
Stephen Roberts
Undergraduate Convenor
Mathematical Sciences Institute
Room 4.74 Hanna Neumann Building #145
The Australian National University
Canberra, ACT 2600 AUSTRALIA
Ph: +61 2 61254445
CRICOS: 00120C

________________________________
From: S.Kummer <ste...@we...>
Sent: Friday, 29 March 2019 5:04 AM
To: anu...@li...
Subject: [Anuga-user] Vertice values stored wrong in sww-file?

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


_______________________________________________
Anuga-user mailing list
Anu...@li...
https://lists.sourceforge.net/lists/listinfo/anuga-user