Anuga-user Digest, Vol 41, Issue 1

28 views
Skip to first unread message

anuga-use...@lists.sourceforge.net

unread,
Jan 9, 2012, 8:39:03 AM1/9/12
to anuga...@lists.sourceforge.net
Send Anuga-user mailing list submissions to
anuga...@lists.sourceforge.net

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
anuga-use...@lists.sourceforge.net

You can reach the person managing the list at
anuga-us...@lists.sourceforge.net

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Anuga-user digest..."


Today's Topics:

1. queries related to simulation and sww file (rajesh tailor)
2. Catchment Finder (Peter Wells)
3. Hello (parul mogra)
4. Anuga-validation (Aurelien Hazart)
5. mesh file (C.Preciado)


----------------------------------------------------------------------

Message: 1
Date: Thu, 20 Oct 2011 16:49:52 +0530
From: rajesh tailor <tailor....@gmail.com>
Subject: [Anuga-user] queries related to simulation and sww file
To: anuga...@lists.sourceforge.net
Message-ID:
<CAA1D9o65hYV7m51RZm2CDA_O...@mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"

Hi All,

1) Previously I queried about anuga. Now I want to know that what
is the minimum system requirement to generate a fine mesh upto resolution of
50m. My domain area is about 1000km^2. I am using the machine which have
16GB memory, even thought it is not generating the mesh with 2000
resolution.

2) How to convert sww file to avi or other video format.

3) I wish to simulate dam break phenomena in the following scenario:
I have the volume/quantum of water as my initial input. I am unable to
determine the flow rate with which water will flow from dam.
I am using Inflow function for the same which require rate, polygon, domain
as required inputs. How to set the rate of water outflow when the volume of
lake is given.

regards
Rajesh Tailor
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 2
Date: Thu, 27 Oct 2011 17:51:29 +0100
From: Peter Wells <wel...@faunalia.co.uk>
Subject: [Anuga-user] Catchment Finder
To: anuga...@lists.sourceforge.net
Message-ID: <1319734289.3142.8.camel@klaus>
Content-Type: text/plain; charset="utf-8"

Hello all,

You may be interested in a free hydrology/water related WebGIS
application that Faunalia have just released called Catchment Finder.

It allows you to calculate (and download) the hydrological catchment
area draining to a point anywhere in the UK through the web browser.

Feel free to give it a go (requires internet explorer 7 or above):

http://www.faunalia.co.uk/catchmentfinder

Kind regards,

Pete

Peter Wells

Faunalia UK
23 Chestnut Close, Burgess Hill, RH15 8HN
www.faunalia.co.uk
Tel: +44 (0)1444 848012
Mob: +44 (0)7704 781683
SIP: 150...@sipgate.co.uk
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 3
Date: Fri, 25 Nov 2011 16:43:43 +0530
From: parul mogra <scori...@gmail.com>
Subject: [Anuga-user] Hello
To: anuga...@lists.sourceforge.net
Message-ID:
<CAPkZ3MTovrXYBFirAnC_xk97...@mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"


-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 4
Date: Wed, 28 Dec 2011 13:43:23 +0900
From: "Aurelien Hazart" <auha...@getc.co.jp>
Subject: [Anuga-user] Anuga-validation
To: <anuga...@lists.sourceforge.net>
Message-ID: <000301ccc51b$3b51d7b0$b1f58710$@co.jp>
Content-Type: text/plain; charset="iso-2022-jp"

Dear Ole (and all)

Thank you for the great work on Anuga.

I am now trying to run the validation cases but some errors messages
happened, mostly regarding the import process.

I am running Anuga 1.2.1. Maybe I didn?t get the last version of the
validation scripts. Is there any updated version of the validations scripts
?

Thanks a lot,

Aurelien Hazart

Dear all

The Patong Bay validation scenario which is the a cornerstone in the

validation of ANUGA has now been rolled into the validation regressions test

suite validate_all.py (as of revision 7669) along with previous validation

results. It takes a few days to run on our server (Intel I7) but it serves

to verify that any future release of ANUGA still produces the results report

in our publication to within less than 1% error. Timeseries at selected

gauges are produced by the validation code for visual inspection as well.

To run all components of this you'll need something like 8GB of RAM but the

script will not attempt to run Patong unless this is the case.

Validate all runs successfully on our 64 bit cluster but please let me know

if anyone encounter problems. And let me know if it works for you too.

All the best

Ole

--------------------------

????????? ?HAZART Aurelien

??????????????

?101-0052?????????????????????

?? : 03-5283-5825 FAX : 03-5283-5826?

E-mail : <mailto:sa...@getc.co.jp> auha...@getc.co.jp

URL : <http://www.getc.co.jp> http://www.getc.co.jp

-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 5
Date: Mon, 9 Jan 2012 14:14:24 +0100
From: "C.Preciado" <p...@ces-info.de>
Subject: [Anuga-user] mesh file
To: <anuga...@lists.sourceforge.net>
Message-ID: <0A03B22AEA7A476EB31C879C3624E083@BCPCIBM18>
Content-Type: text/plain; charset="iso-8859-1"

Dear ANUGA users,

First of all I would like to introduce a little my work within my company
and current situation as a user of the ANUGA.

We would like to simulate a flow in a river bank in current situation and
with new scenarios with polder construction in order to prevent flood
events. From a topodata XYZ file (.csv) we have created a mesh file (.msh)
using Anuga.

The output mesh file can not be opened apparently because it is based on the
NetCDF format and the software I am using (Global Mapper) only supports
raster/grid data from NetCDF files and can?t load generically formatted
vector data like these meshes.

Since I am not a programmer and thus not an expert on scripting, I was
wondering if anyone could assist me with ideas on improving the edition of
my script and also in advising which might be the best software for the mesh
file to be managed (opened, visualised).

So far, my scrip looks as it follows:

""" Converts ASC Points Data (csv) to NETCDF DEM

This Script is used to Reformat a XYZ Data file from Ascii to NetCDF

NetCDF if the Format ANUGA Uses

"""

#---------------------------------------------------------------------------
---

# Preparation of topographic data

# Convert ASC 2 DEM 2 PTS using source data and store result in source data

#---------------------------------------------------------------------------
---

#---------------------------------------------------------------------------
---

# Import necessary modules

#---------------------------------------------------------------------------
---

# Standard modules

import anuga

from pylab import figure, plot, axis, quiver, quiverkey, show, title,
axhline

from pylab import cos, sin, pi

import numpy

import csv

import time

import os

from os.path import join

import time

import sys

import os.path

import glob

from anuga.abstract_2d_finite_volumes.util import file_function

from anuga.geometry.polygon_function import Polygon_function

from anuga.geometry.polygon import read_polygon, plot_polygons, polygon_area

from anuga.abstract_2d_finite_volumes.generic_domain import Quantity

from anuga.pmesh.mesh_interface import create_mesh_from_regions

from anuga.shallow_water.shallow_water_domain import Domain

from anuga.shallow_water.forcing import Rainfall

from anuga.shallow_water.forcing import Inflow

# from model_build_tools import get_polygon_value_list

select = """

Available Scenarios for this PROJECT

1. xyzfromRH.csv

"""

answer = int(raw_input(select))

if answer == 1:

basename = 'xyzfromRH'

if answer == 1:

elevation_file = 'xyzfromRH'

# Filenames

meshname = basename+'.msh'

meshname =join('xyzfromRH.msh') # mesh name, used to store mesh, if not
changed will speed up run

print 'Basename ',basename

#+++++END_SCENARIOS_BLOCK

#+++++START_PREPARE TERRAIN_BLOCK

#---------------------------------------------------------------------------
---

# Preparation of topographic data

# Convert ASC 2 DEM 2 PTS using source data and store result in source data

#---------------------------------------------------------------------------
---

print ' Preparing Terrain'

# Create pts file

print 'Reading and Converting', elevation_file

print '... This might take a minute or so depending on file size.....'

W=440747.9226

N=251067.1197

E=458455.4778

S=227167.9114

RangeX=E-W

RangeY=N-S

domain_Area= RangeX*RangeY

print' Domain is',RangeX,' x ',RangeY,' = ',domain_Area,' m2'

bounding_polygon = [[W, S], [E, S], [E, N], [W, N]] # should read this
polygon directly from a file be carefull of tags !!

G = Geospatial_data(elevation_file+'.csv')

G.export_points_file(basename+'.pts')

print 'Creating', meshname

print ' Creating Mesh'

maximum_triangle_area =1000.0 # basic mesh size for the catchment boundary

create_mesh_from_regions(bounding_polygon,

boundary_tags={'south': [0],

'east': [1],

'north': [2],

'west': [3]},

maximum_triangle_area=maximum_triangle_area,

filename=meshname,

use_cache=True,

verbose=True)

Thanks in advance.

Best regards,

C.Preciado
___________________________________________________________

CES Consulting Engineers Salzgitter GmbH
Hamburger Stra?e 277, D-38114 Braunschweig, Germany
Tel.: (49) +531 / 31704-164 / Fax: (49) +531 / 31704-199 / Skype:
ces-phone-pre
_______________________________________________________________________
Registered Office (Sitz der Gesellschaft): D-38114 Braunschweig
District Court of Braunschweig (Amtsgericht Braunschweig):
Register of Companies HRB 101260
Managing Directors (Gesch?ftsf?hrung):
Dr.-Ing. Detlef H?gin, Dr.-Ing. Ralf Meyerhoff
_______________________________________________________________________


-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

------------------------------------------------------------------------------
Ridiculously easy VDI. With Citrix VDI-in-a-Box, you don't need a complex
infrastructure or vast IT resources to deliver seamless, secure access to
virtual desktops. With this all-in-one solution, easily deploy virtual
desktops for less than the cost of PCs and save 60% on VDI infrastructure
costs. Try it free! http://p.sf.net/sfu/Citrix-VDIinabox

------------------------------

_______________________________________________
Anuga-user mailing list
Anuga...@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/anuga-user


End of Anuga-user Digest, Vol 41, Issue 1
*****************************************

Reply all
Reply to author
Forward
0 new messages