Re: [Anuga-user] set stage value for interior regions

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

Hi Tejaswi,

You can pass a function f(x,y) to domain.set_quantity, like:

 > domain.set_quantity('stage', f)

So you just need to make function f(x,y) which takes the x/y coordinates 
and returns the desired stage at x/y (presumably with a different 
approach depending on whether you are in the 'river basin' or not).

One routine that might help you build f(x,y) is 
anuga.utilities.quantity_setting_functions.composite_quantity_setting_function 
-- this is designed to apply different functions in polygonal regions -- 
it is fairly well documented.

If you use another approach, beware that f(x,y) assumes x/y are in 
ANUGA's **internal** coordinate system -- which has a lower-left of 
(0,0) -- so you will generally need to translate the coordinates inside 
your function.

Cheers,
Gareth.


On 27/07/17 14:41, Tejaswi Jadhav wrote:
> Hi
>     I am software developer and new to anuga I want set stage value 
> for interior regions with stage= elevation+depth.
> I set stage for whole domain
>       h = 0.05 # Constant depth
> domain.set_quantity('stage', expression='elevation + %f' % h)
>
> But I want set stage for river basin ( interior region ) how to do that?
>
> Thank you.
>
>
> Regards,
> Tejaswi
>
>
> ------------------------------------------------------------------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, Slashdot.org! http://sdm.link/slashdot
>
>
> _______________________________________________
> Anuga-user mailing list
> Anu...@li...
> https://lists.sourceforge.net/lists/listinfo/anuga-user