- Timestamp:
- Jun 15, 2010, 9:11:59 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_core/documentation/user_manual/demos/cairns/project.py
r7816 r7838 3 3 """ 4 4 5 from anuga.geometry.polygon import read_polygon, plot_polygons, \ 6 polygon_area, is_inside_polygon 5 import anuga 7 6 8 7 #------------------------------------------------------------------------------ 9 8 # Define scenario as either slide or fixed_wave. Choose one. 10 9 #------------------------------------------------------------------------------ 11 #scenario = 'fixed_wave' # Huge wave applied at the boundary12 scenario = 'slide' # Slide wave form applied inside the domain10 scenario = 'fixed_wave' # Huge wave applied at the boundary 11 #scenario = 'slide' # Slide wave form applied inside the domain 13 12 14 13 #------------------------------------------------------------------------------ 15 14 # Filenames 16 15 #------------------------------------------------------------------------------ 17 demname_stem = 'cairns'18 meshname = demname+ '.msh'16 name_stem = 'cairns' 17 meshname = name_stem + '.msh' 19 18 20 19 # Filename for locations where timeseries are to be produced … … 25 24 #------------------------------------------------------------------------------ 26 25 # bounding polygon for study area 27 bounding_polygon = read_polygon('extent.csv')26 bounding_polygon = anuga.read_polygon('extent.csv') 28 27 29 A = polygon_area(bounding_polygon) / 1000000.028 A = anuga.polygon_area(bounding_polygon) / 1000000.0 30 29 print 'Area of bounding polygon = %.2f km^2' % A 31 30 … … 34 33 #------------------------------------------------------------------------------ 35 34 # Read interior polygons 36 poly_cairns = read_polygon('cairns.csv')37 poly_island0 = read_polygon('islands.csv')38 poly_island1 = read_polygon('islands1.csv')39 poly_island2 = read_polygon('islands2.csv')40 poly_island3 = read_polygon('islands3.csv')41 poly_shallow = read_polygon('shallow.csv')35 poly_cairns = anuga.read_polygon('cairns.csv') 36 poly_island0 = anuga.read_polygon('islands.csv') 37 poly_island1 = anuga.read_polygon('islands1.csv') 38 poly_island2 = anuga.read_polygon('islands2.csv') 39 poly_island3 = anuga.read_polygon('islands3.csv') 40 poly_shallow = anuga.read_polygon('shallow.csv') 42 41 43 42 # Optionally plot points making up these polygons
Note: See TracChangeset
for help on using the changeset viewer.