From: Rudy <ru...@ba...> - 2012-12-03 07:41:15
|
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 |