[Anuga-user] Fwd: Anuga-user Digest, Vol 51, Issue 4

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 454-5900

---------- Forwarded message ----------
From: Ole Nielsen <ole...@gm...>
Date: 7 December 2012 08:45
Subject: Re: Anuga-user Digest, Vol 51, Issue 4
To: rajesh tailor <tai...@gm...>

Hi Rajesh

Yes the number is relative to mean sea level (=0 in the model).
If you want to *add* water depth to the existing surface (or add heigh to
elevation data) you can use the method add_quantity instead of set_quantity.

In your case, I was just wondering if the topography is higher than 167.20?
If that is the case, you won't see any water because it is an absolute
number.

Cheers
Ole

On 6 December 2012 13:04, rajesh tailor <tai...@gm...> wrote:

> Hi Ole,
>
> I have done the same and triangles are falling in the specified polygon
> only... Can you tell me what height it is taking for reference to provide
> specified height in the following function..
>  domain.set_quantity('stage',Polygon_function([(poly,167.20)]))
>
> If I use this function this should assign 167.20 m height of water surface
> in specified polygon (according to you) but what reference it is taking to
> assign this height... Mean Sea Level...??
>
>
> On Wed, Dec 5, 2012 at 8:37 AM, Ole Nielsen <ole...@gm...>wrote:
>
>> Hi Rajesh
>>
>> What you are doing looks correct.
>> The call
>> domain.set_quantity('stage',Polygon_function([(poly,167.20)]))
>>
>> will set the water level 'stage' to 167.20 for every triangle that falls
>> inside the polygon.
>> I don't know why it doesn't work for you, but maybe your mesh is too
>> coarse so that no triangles fall inside the polygon?
>>
>> just a thought
>>
>> Cheers
>> Ole
>>
>>
>>
>> On 4 December 2012 19:33, <anu...@li...>wrote:
>>
>>> Send Anuga-user mailing list submissions to
>>>         anu...@li...
>>>
>>> To subscribe or unsubscribe via the World Wide Web, visit
>>>         https://lists.sourceforge.net/lists/listinfo/anuga-user
>>> or, via email, send a message with subject or body 'help' to
>>>         anu...@li...
>>>
>>> You can reach the person managing the list at
>>>         anu...@li...
>>>
>>> When replying, please edit your Subject line so it is more specific
>>> than "Re: Contents of Anuga-user digest..."
>>>
>>>
>>> Today's Topics:
>>>
>>>    1. Re: Anuga-user Digest, Vol 51, Issue 3 (rajesh tailor)
>>>
>>>
>>> ----------------------------------------------------------------------
>>>
>>> Message: 1
>>> Date: Tue, 4 Dec 2012 10:56:59 +0530
>>> From: rajesh tailor <tai...@gm...>
>>> Subject: Re: [Anuga-user] Anuga-user Digest, Vol 51, Issue 3
>>> To: anu...@li...
>>> Message-ID:
>>>         <
>>> CAA...@ma...>
>>> Content-Type: text/plain; charset="iso-8859-1"
>>>
>>> Hi Rudy,
>>> I tried this...
>>>
>>> domain.set_quantity('stage',Polygon_function([('Water_Storage_DAM_TWL_167_2.csv',167.20)]))
>>> But this is giving errror...
>>> ValueError: setting an array element with a sequence
>>>
>>> so I have tried like this...
>>> poly = anuga.read_polygon('Water_Storage_DAM_TWL_167_2.csv')
>>> domain.set_quantity('stage',Polygon_function([(poly,167.20)]))
>>>
>>> this is working without giving any error... but it is not adding anything
>>> to the polygon specified and nothing is changed in whole domain....
>>>
>>>
>>> In previous case whenever I was trying to set stage values within a
>>> polygon, that stage value was taking exact shape of terrain in that
>>> polygon
>>> but Ideally water surface in that polygon should be horizontal
>>> (constant),
>>> not varying according to terrain..
>>>
>>> thanks,
>>>
>>> regards,
>>> Rajesh Tailor
>>>
>>>
>>> On Mon, Dec 3, 2012 at 5:53 PM, <
>>> anu...@li...>wrote:
>>>
>>> > Send Anuga-user mailing list submissions to
>>> >         anu...@li...
>>> >
>>> > To subscribe or unsubscribe via the World Wide Web, visit
>>> >         https://lists.sourceforge.net/lists/listinfo/anuga-user
>>> > or, via email, send a message with subject or body 'help' to
>>> >         anu...@li...
>>> >
>>> > You can reach the person managing the list at
>>> >         anu...@li...
>>> >
>>> > When replying, please edit your Subject line so it is more specific
>>> > than "Re: Contents of Anuga-user digest..."
>>> >
>>> >
>>> > Today's Topics:
>>> >
>>> >    1. Re: Anuga-user Digest, Vol 51, Issue 2 (Rudy)
>>> >
>>> >
>>> > ----------------------------------------------------------------------
>>> >
>>> > Message: 1
>>> > Date: Mon, 03 Dec 2012 18:21:34 +1100
>>> > From: Rudy <ru...@ba...>
>>> > Subject: Re: [Anuga-user] Anuga-user Digest, Vol 51, Issue 2
>>> > To: anu...@li...
>>> > Message-ID: <50B...@ba...>
>>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>>> >
>>> > Ok
>>> > A. to Q.
>>> >
>>> > I am not getting it.. could you please elaborate this...
>>> >
>>> > The first example with Tide Sets the Stage over the entire domain to a
>>> > level of 3.75m
>>> > So in any of your terrain where the ground level is less than 3.75m the
>>> > depth will be greater than 0.0m
>>> > tide = 3.75
>>> > domain.set_quantity('stage', tide)
>>> >
>>> > The 2nd example does the same but restricts the area of coverage to
>>> just
>>> > inside the Polygon....
>>> >
>>> > Your Q. what is this...
>>> > (join('02_POLYS','04_STAGE','16720','Water_Storage_DAM_TWL_167_2.
>>> > csv',167.20)
>>> >
>>> > Using :
>>> >
>>> > from os.path import join
>>> >
>>> > ensures that you can use a file in a directory tree structure so that
>>> it
>>> > is platform independent....
>>> >
>>> > if your csv file is in the same directory as your ANUGA script then you
>>> > simply use...
>>> >
>>> >
>>> >
>>> domain.set_quantity('stage',Polygon_function([('Water_Storage_DAM_TWL_167_2.csv',167.20)]))
>>> >
>>> > You Q.   and what the above .csv file will contains and in what
>>> format...
>>> >
>>> > ASCII   X,Y on each line describing your polyline
>>> > X1,Y1
>>> > X2,Y2
>>> > X3,Y3
>>> > etc...
>>> >
>>> > Regards Rudy
>>> >
>>> >
>>> >
>>> > On 3/12/2012 5:31 PM, anu...@li... wrote:
>>> > > hi Rudy,
>>> > >
>>> > > I am not getting it.. could you please elaborate this...
>>> > >
>>> > > There are a couple of ways to achieve this depending on your
>>> particular
>>> > > circumstance...
>>> > >
>>> > >
>>> > > For example if the river is tidal and connected to a boundary, can
>>> set a
>>> > > boundary condition
>>> > >
>>> > > Or better set the stage at the start of the run as follows:
>>> > > #-----------------------------------------------------------
>>> > > -------------------
>>> > > # SET INTIAL WATER LEVEL OVER DOMAIN
>>> > > #-----------------------------------------------------------
>>> > > -------------------
>>> > > tide = 3.75
>>> > > domain.set_quantity('stage', tide)
>>> > >
>>> > >
>>> > > If you really want to restrict the coverage, as in a lake then you
>>> can
>>> > > define a polygon and fill the stage in that polygon
>>> > >
>>> > > #-----------------------------------------------------------
>>> > > -------------------
>>> > > # Set a Initial Water Level over the Domain or over an area defined
>>> by a
>>> > > polyline
>>> > > #  eg if u have a basin that has water in it already u can nominate a
>>> > > polyline and set the stage level inside that polyline
>>> > > #-----------------------------------------------------------
>>> > > -------------------
>>> > >
>>> > > domain.set_quantity('stage',Polygon_function([(join('02_
>>> > >
>>> POLYS','04_STAGE','16720','Water_Storage_DAM_TWL_167_2.csv',167.20)]))
>>> > > #  Sets the Water LEvel over the Polygon to 167.20m AHD ???
>>> > > domain.set_quantity('stage',Polygon_function([(join('02_
>>> > > POLYS','04_STAGE','19165','SED_DAM_TWL_191_65.csv',191.65)]))
>>> > > # Sets the Water LEvel over the Polygon to 191.65m AHD ???
>>> > >
>>> > >
>>> > >
>>> > >
>>> > > what is this...
>>> > > (join('02_POLYS','04_STAGE','16720','Water_Storage_DAM_TWL_167_2.
>>> > > csv',167.20)
>>> > >
>>> > > and what the above .csv file will contains and in what format...
>>> > >
>>> > > thanks of your kind suggestions..
>>> > >
>>> > > regards,
>>> > > Rajesh Tailo
>>> >
>>> >
>>> >
>>> >
>>> > ------------------------------
>>> >
>>> >
>>> >
>>> ------------------------------------------------------------------------------
>>> > Keep yourself connected to Go Parallel:
>>> > BUILD Helping you discover the best ways to construct your parallel
>>> > projects.
>>> > http://goparallel.sourceforge.net
>>> >
>>> > ------------------------------
>>> >
>>> > _______________________________________________
>>> > Anuga-user mailing list
>>> > Anu...@li...
>>> > https://lists.sourceforge.net/lists/listinfo/anuga-user
>>> >
>>> >
>>> > End of Anuga-user Digest, Vol 51, Issue 3
>>> > *****************************************
>>> >
>>> -------------- next part --------------
>>> An HTML attachment was scrubbed...
>>>
>>> ------------------------------
>>>
>>>
>>> ------------------------------------------------------------------------------
>>> LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial
>>> Remotely access PCs and mobile devices and provide instant support
>>> Improve your efficiency, and focus on delivering more value-add services
>>> Discover what IT Professionals Know. Rescue delivers
>>> http://p.sf.net/sfu/logmein_12329d2d
>>>
>>> ------------------------------
>>>
>>> _______________________________________________
>>> Anuga-user mailing list
>>> Anu...@li...
>>> https://lists.sourceforge.net/lists/listinfo/anuga-user
>>>
>>>
>>> End of Anuga-user Digest, Vol 51, Issue 4
>>> *****************************************
>>>
>>
>>
>