[Anuga-user] Obtaining maximum stage values from internal timesteps

[Ole N. <ole...@gm...>]

Hi Miguel (and everyone else)

Thanks for your posting and your interest in ANUGA

What you are asking is not currently implemented in ANUGA (alas there'll
always be more things we could do),
but it seems very doable (I always say that) - by copying and modifying  the
code behind get_maximum_value().
We understand that you need the maximum across all internal timesteps (hence
using the sww file won't be adequate for your need) and you need it for all
triangles.

So, lets try to hone in on your requirements:


   - You want the maximum value of a quantity for each triangle
   (centroid) accumulated over time?
      - Does this mean since the simulation started and with an update
      at every timestep?
   - Suppose your function was called
   quantity.get_maximum_values_over_time()
      - How would you like to call it?
      - Would you need to restrict the time interval and/or the mesh
      over which the calculations are done?
      - And what would you ideally like it to return (e.g. an array of
      values - one for each triangle)?

The only issue I can see is that of storage and overhead as we'd need to ask
each quantity to record maxima for each triangle and update it at each
timestep. But if that it what you need, it is certainly doable.

Best regards
Ole







Message: 2
Date: Tue, 6 May 2008 18:52:31 +0100
From: "Miguel Fernandes" <mig...@gm...>
Subject: [Anuga-user] Obtaining maximum stage values from internal
       timesteps
To: anu...@li...
Message-ID:
       <5c0...@ma...>
Content-Type: text/plain; charset="iso-8859-1"

Hello all,

I would like to obtain the the maximum stage value in each of the triangles.
But from a comparison between the stages calculated in the internal
timesteps and not from the each slice in the output sww. Something like
get_maximum_value() but for each of the triangles.

I have been testing with get_maximum_value() but it only output the highest
and lowest value and where it occured, also tried during the evolve loop
with get_quantity('stage').get_values() and compare from the last iteration
in loop with a yieldstep of 100s and a finaltime of 4500s but seems that i'm
doing the same of comparing each stage slice in the sww and get the maximum
values there.

Hope to have clear in what i tried to explain, any help would be very
important.