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

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

Hi Stehphen,

yes I actually use

res = get_output(self.sww_file, timeSlices="max")
(It's a typo in the first email).

The problem is with all quantities.

It works when I use the code:
domain.set_store_vertices_uniquely(True)

But I try to avoid this due to the larger file-size.

As a next step I will try to examine the sww-file directly.


Cheers
Stefan




Am Thu, 28 Mar 2019 22:23:32 +0000
schrieb Stephen Roberts <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