[[Image(Screenshot.png)]] = Buildings, Bridges and Dams = A new feature in ANUGA 1.2 lets you create create arbitrary geometry to represent solid structures. These take the form of holes within the mesh, and can be given arbitrary boundary properties, just like the edges of the domain. The demo code can be found within demos/buildings.py, part of which is reproduced here: {{{ length = 50 width = 10 resolution = 0.15 # make this number smaller to make the simulation more accurate # Create the "world" as a long, skinny channel boundary_poly = poly_from_box(0, length, 0, width) # Place 3 buildings downstream building_polys = [ poly_from_box(10, 15, 2.5, 7.5), # upstream box poly_from_box(22.5, 27.5, 1.5, 6.5), # middle box poly_from_box(35, 40, 3.5, 8.5)] # downstream box # create a domain mesh, with 3 building holes in it domain = anuga.create_domain_from_regions(boundary_poly, boundary_tags={'left': [0], 'bottom': [1], 'right': [2], 'top': [3]}, maximum_triangle_area = resolution, mesh_filename = 'building.msh', interior_holes = building_polys, use_cache=True, # to speed it up verbose=True) # log output on }}} Note the new parameter to create_domain_from_regions: {{{ interior_holes = building_polys }}} This lets you specify a list of polygons to be cut out of the mesh. The internal boundaries have their own boundary tag: {{{exterior}}} which can be set just like any other boundary.