1 | """Common filenames and locations for topographic data, meshes and outputs.""" |
---|
2 | |
---|
3 | from anuga.utilities.polygon import read_polygon, plot_polygons, \ |
---|
4 | polygon_area, is_inside_polygon |
---|
5 | |
---|
6 | #------------------------------------------------------------------------------ |
---|
7 | # Define scenario as either slide or fixed_wave. Choose one. |
---|
8 | #------------------------------------------------------------------------------ |
---|
9 | #scenario = 'fixed_wave' # Huge wave applied at the boundary |
---|
10 | scenario = 'slide' # Slide wave form applied inside the domain |
---|
11 | |
---|
12 | #------------------------------------------------------------------------------ |
---|
13 | # Filenames |
---|
14 | #------------------------------------------------------------------------------ |
---|
15 | demname = 'cairns' |
---|
16 | meshname = demname + '.msh' |
---|
17 | |
---|
18 | # Filename for locations where timeseries are to be produced |
---|
19 | gauge_filename = 'gauges.csv' |
---|
20 | |
---|
21 | #------------------------------------------------------------------------------ |
---|
22 | # Domain definitions |
---|
23 | #------------------------------------------------------------------------------ |
---|
24 | # bounding polygon for study area |
---|
25 | bounding_polygon = read_polygon('extent.csv') |
---|
26 | |
---|
27 | A = polygon_area(bounding_polygon) / 1000000.0 |
---|
28 | print 'Area of bounding polygon = %.2f km^2' % A |
---|
29 | |
---|
30 | #------------------------------------------------------------------------------ |
---|
31 | # Interior region definitions |
---|
32 | #------------------------------------------------------------------------------ |
---|
33 | # Read interior polygons |
---|
34 | poly_cairns = read_polygon('cairns.csv') |
---|
35 | poly_island0 = read_polygon('islands.csv') |
---|
36 | poly_island1 = read_polygon('islands1.csv') |
---|
37 | poly_island2 = read_polygon('islands2.csv') |
---|
38 | poly_island3 = read_polygon('islands3.csv') |
---|
39 | poly_shallow = read_polygon('shallow.csv') |
---|
40 | |
---|
41 | # Optionally plot points making up these polygons |
---|
42 | #plot_polygons([bounding_polygon, poly_cairns, poly_island0, poly_island1, |
---|
43 | # poly_island2, poly_island3, poly_shallow], |
---|
44 | # style='boundingpoly', verbose=False) |
---|
45 | |
---|
46 | # Define resolutions (max area per triangle) for each polygon |
---|
47 | default_res = 10000000 # Background resolution |
---|
48 | islands_res = 100000 |
---|
49 | cairns_res = 100000 |
---|
50 | shallow_res = 500000 |
---|
51 | |
---|
52 | # Define list of interior regions with associated resolutions |
---|
53 | interior_regions = [[poly_cairns, cairns_res], |
---|
54 | [poly_island0, islands_res], |
---|
55 | [poly_island1, islands_res], |
---|
56 | [poly_island2, islands_res], |
---|
57 | [poly_island3, islands_res], |
---|
58 | [poly_shallow, shallow_res]] |
---|
59 | |
---|
60 | #------------------------------------------------------------------------------ |
---|
61 | # Data for exporting ascii grid |
---|
62 | #------------------------------------------------------------------------------ |
---|
63 | eastingmin = 363000 |
---|
64 | eastingmax = 418000 |
---|
65 | northingmin = 8026600 |
---|
66 | northingmax = 8145700 |
---|
67 | |
---|
68 | #------------------------------------------------------------------------------ |
---|
69 | # Data for landslide |
---|
70 | #------------------------------------------------------------------------------ |
---|
71 | slide_origin = [451871, 8128376] # Assume to be on continental shelf |
---|
72 | slide_depth = 500. |
---|