1 | """Common filenames and locations for topographic data, meshes and outputs. |
---|
2 | Also includes origin for slump scenario. |
---|
3 | """ |
---|
4 | |
---|
5 | from os import sep, environ, getenv, getcwd |
---|
6 | from os.path import expanduser, basename |
---|
7 | #from anuga.utilities.polygon import read_polygon |
---|
8 | import sys |
---|
9 | from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees |
---|
10 | from time import localtime, strftime |
---|
11 | from anuga.geospatial_data.geospatial_data import * |
---|
12 | |
---|
13 | #Making assumptions about the location of scenario data |
---|
14 | state = 'western_australia' |
---|
15 | scenario_dir_name = 'pt_hedland_tsunami_scenario_2006' |
---|
16 | |
---|
17 | # onshore data from 30m DTED level 2 |
---|
18 | onshore_name_dted = 'pt_hedland_onshore_30m_dted' # get from Neil/Ingo (DEM or topo data) |
---|
19 | onshore_name_dli = 'pt_hedland_onshore_20m_dli' # get from Neil/Ingo (DEM or topo data) |
---|
20 | |
---|
21 | # offshore data from GA digitised charts |
---|
22 | offshore_name1 = 'pt_hedland_offshore_points_orig' |
---|
23 | |
---|
24 | # offshore data from AHO fairsheets |
---|
25 | offshore_name2 = 'pt_hedland_offshore_points_fairsheet_orig' |
---|
26 | |
---|
27 | # reissued offshore data from PMD and NM&I Nov 2006 |
---|
28 | offshore_name0 = 'hydrographic_charts' |
---|
29 | offshore_name1 = 'xy100003817' |
---|
30 | offshore_name2 = 'xy100003819' |
---|
31 | offshore_name3 = 'xy100003822' |
---|
32 | offshore_name4 = 'xy100003824' |
---|
33 | offshore_name5 = 'xy100003825' |
---|
34 | offshore_name6 = 'xy100003826' |
---|
35 | offshore_name7 = 'xy100003827' |
---|
36 | offshore_name8 = 'xy100003828' |
---|
37 | offshore_name9 = 'xy100003829' |
---|
38 | offshore_name10 = 'xy100003830' |
---|
39 | offshore_name11 = 'xy100003831' |
---|
40 | offshore_name12 = 'xy100003832' |
---|
41 | offshore_name13 = 'xy100003837' |
---|
42 | offshore_name14 = 'xy100003838' |
---|
43 | offshore_name15 = 'xy100003840' |
---|
44 | offshore_name16 = 'xy100003841' |
---|
45 | offshore_name17 = 'xy100003843' |
---|
46 | offshore_name18 = 'xy100003845' |
---|
47 | offshore_name19 = 'xy100003846' |
---|
48 | offshore_name20 = 'xy100003847' |
---|
49 | offshore_name21 = 'xy100003945' |
---|
50 | offshore_name22 = 'xy100003946' |
---|
51 | offshore_name23 = 'xy100003947' |
---|
52 | offshore_name24 = 'xy100003948' |
---|
53 | offshore_name25 = 'xy100004071' |
---|
54 | offshore_name26 = 'xy100004072' |
---|
55 | offshore_name27 = 'xy100004074' |
---|
56 | offshore_name28 = 'xy100004078' |
---|
57 | offshore_name29 = 'xy100004079' |
---|
58 | offshore_name30 = 'xy100004080' |
---|
59 | offshore_name31 = 'xy100004081' |
---|
60 | offshore_name32 = 'xy100004082' |
---|
61 | offshore_name33 = 'xy100004083' |
---|
62 | offshore_name34 = 'xy100004084' |
---|
63 | offshore_name35 = 'xy100004085' |
---|
64 | offshore_name36 = 'xy100004086' |
---|
65 | offshore_name37 = 'xy100004087' |
---|
66 | offshore_name38 = 'xy100004088' |
---|
67 | offshore_name39 = 'xy100004089' |
---|
68 | offshore_name40 = 'xy100004090' |
---|
69 | offshore_name41 = 'xy100011753' |
---|
70 | offshore_name42 = 'xy100014967' |
---|
71 | offshore_name43 = 'xy100014968' |
---|
72 | offshore_name44 = 'xy100014969' |
---|
73 | offshore_name45 = 'xy100021030' |
---|
74 | offshore_name46 = 'xy100021031' |
---|
75 | offshore_name47 = 'xy100021032' |
---|
76 | offshore_name48 = 'xy100021033' |
---|
77 | offshore_name49 = 'xy100021035' |
---|
78 | offshore_name50 = 'xy100022767' |
---|
79 | |
---|
80 | # coastline developed from aerial photography and 1.5m DLI contour |
---|
81 | #coast_name = 'pt_hedland_coastline_points_dted2_new' |
---|
82 | coast_name = 'coast_with_extracted_z' |
---|
83 | |
---|
84 | |
---|
85 | boundary_basename = 'SU-AU_clip' |
---|
86 | |
---|
87 | #swollen/ all data output |
---|
88 | basename = 'source' |
---|
89 | |
---|
90 | codename = 'project.py' |
---|
91 | |
---|
92 | if sys.platform == 'win32': |
---|
93 | home = getenv('INUNDATIONHOME') #Sandpit's parent dir |
---|
94 | # python_home = getenv('PWD') |
---|
95 | user = getenv('USERPROFILE') |
---|
96 | #print 'USER:', user |
---|
97 | else: |
---|
98 | # original |
---|
99 | #home = getenv('INUNDATIONHOME', sep+'d'+sep+'cit'+sep+'1'+sep+'cit'+sep+'risk_assessment_methods_project'+sep+'inundation') |
---|
100 | # update to perlite 2 |
---|
101 | home = getenv('INUNDATIONHOME', sep+'d'+sep+'cit'+sep+'2'+sep+'cit'+sep+'inundation'+sep+'data') |
---|
102 | user = getenv('LOGNAME') |
---|
103 | |
---|
104 | # INUNDATIONHOME is the inundation directory, not the data directory. |
---|
105 | home += sep +'data' |
---|
106 | |
---|
107 | #Derive subdirectories and filenames |
---|
108 | time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir |
---|
109 | outputtimedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep+time+sep |
---|
110 | |
---|
111 | #print 'outputtimedir', outputtimedir |
---|
112 | #meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep |
---|
113 | #datadir = home+sep+scenario_dir_name+sep+'topographies'+sep |
---|
114 | #gaugedir = home+sep+scenario_dir_name+sep+'gauges'+sep |
---|
115 | #polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep |
---|
116 | #boundarydir = home+sep+scenario_dir_name+sep+'boundaries'+sep |
---|
117 | #outputdir = home+sep+scenario_dir_name+sep+'outputs'+sep |
---|
118 | #tidedir = home+sep+scenario_dir_name+sep+'tide_data'+sep |
---|
119 | |
---|
120 | meshdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'meshes'+sep |
---|
121 | datadir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'topographies'+sep |
---|
122 | gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep |
---|
123 | polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep |
---|
124 | boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep |
---|
125 | outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep |
---|
126 | tidedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'tide_data'+sep |
---|
127 | |
---|
128 | gauge_filename = gaugedir + 'gauge_location_port_hedland.csv' |
---|
129 | gauge_checking = gaugedir + 'gauge_checking_test.csv' |
---|
130 | buildings_filename = gaugedir + 'pt_hedland_res.csv' |
---|
131 | buildings_filename_out = 'pt_hedland_res_modified.csv' |
---|
132 | buildings_filename_damage_out = 'pt_hedland_res_modified_damage.csv' |
---|
133 | community_filename = gaugedir + 'CHINS_v2.csv' |
---|
134 | community_scenario = gaugedir + 'community_pt_hedland.csv' |
---|
135 | tidal_filename = tidedir + 'pt_hedland_tide.txt' |
---|
136 | |
---|
137 | meshname = meshdir + basename |
---|
138 | #onshore_dem_name = datadir + onshore_name_dted |
---|
139 | onshore_dem_name = datadir + onshore_name_dli |
---|
140 | offshore_dem_name0 = datadir + offshore_name0 |
---|
141 | offshore_dem_name1 = datadir + offshore_name1 |
---|
142 | offshore_dem_name2 = datadir + offshore_name2 |
---|
143 | offshore_dem_name3 = datadir + offshore_name3 |
---|
144 | offshore_dem_name4 = datadir + offshore_name4 |
---|
145 | offshore_dem_name5 = datadir + offshore_name5 |
---|
146 | offshore_dem_name6 = datadir + offshore_name6 |
---|
147 | offshore_dem_name7 = datadir + offshore_name7 |
---|
148 | offshore_dem_name8 = datadir + offshore_name8 |
---|
149 | offshore_dem_name9 = datadir + offshore_name9 |
---|
150 | offshore_dem_name10 = datadir + offshore_name10 |
---|
151 | offshore_dem_name11 = datadir + offshore_name11 |
---|
152 | offshore_dem_name12 = datadir + offshore_name12 |
---|
153 | offshore_dem_name13 = datadir + offshore_name13 |
---|
154 | offshore_dem_name14 = datadir + offshore_name14 |
---|
155 | offshore_dem_name15 = datadir + offshore_name15 |
---|
156 | offshore_dem_name16 = datadir + offshore_name16 |
---|
157 | offshore_dem_name17 = datadir + offshore_name17 |
---|
158 | offshore_dem_name18 = datadir + offshore_name18 |
---|
159 | offshore_dem_name19 = datadir + offshore_name19 |
---|
160 | offshore_dem_name20 = datadir + offshore_name20 |
---|
161 | offshore_dem_name21 = datadir + offshore_name21 |
---|
162 | offshore_dem_name22 = datadir + offshore_name22 |
---|
163 | offshore_dem_name23 = datadir + offshore_name23 |
---|
164 | offshore_dem_name24 = datadir + offshore_name24 |
---|
165 | offshore_dem_name25 = datadir + offshore_name25 |
---|
166 | offshore_dem_name26 = datadir + offshore_name26 |
---|
167 | offshore_dem_name27 = datadir + offshore_name27 |
---|
168 | offshore_dem_name28 = datadir + offshore_name28 |
---|
169 | offshore_dem_name29 = datadir + offshore_name29 |
---|
170 | offshore_dem_name30 = datadir + offshore_name30 |
---|
171 | offshore_dem_name31 = datadir + offshore_name31 |
---|
172 | offshore_dem_name32 = datadir + offshore_name32 |
---|
173 | offshore_dem_name33 = datadir + offshore_name33 |
---|
174 | offshore_dem_name34 = datadir + offshore_name34 |
---|
175 | offshore_dem_name35 = datadir + offshore_name35 |
---|
176 | offshore_dem_name36 = datadir + offshore_name36 |
---|
177 | offshore_dem_name37 = datadir + offshore_name37 |
---|
178 | offshore_dem_name38 = datadir + offshore_name38 |
---|
179 | offshore_dem_name39 = datadir + offshore_name39 |
---|
180 | offshore_dem_name40 = datadir + offshore_name40 |
---|
181 | offshore_dem_name41 = datadir + offshore_name41 |
---|
182 | offshore_dem_name42 = datadir + offshore_name42 |
---|
183 | offshore_dem_name43 = datadir + offshore_name43 |
---|
184 | offshore_dem_name44 = datadir + offshore_name44 |
---|
185 | offshore_dem_name45 = datadir + offshore_name45 |
---|
186 | offshore_dem_name46 = datadir + offshore_name46 |
---|
187 | offshore_dem_name47 = datadir + offshore_name47 |
---|
188 | offshore_dem_name48 = datadir + offshore_name48 |
---|
189 | offshore_dem_name49 = datadir + offshore_name49 |
---|
190 | offshore_dem_name50 = datadir + offshore_name50 |
---|
191 | |
---|
192 | coast_dem_name = datadir + coast_name |
---|
193 | combined_dem_name = datadir + 'pt_hedland_combined_elevation' |
---|
194 | outputname = outputtimedir + basename #Used by post processing |
---|
195 | |
---|
196 | # clipping region to make DEM (pts file) from onshore data |
---|
197 | eastingmin = 594000 |
---|
198 | eastingmax = 715000 |
---|
199 | northingmin = 7720000 |
---|
200 | northingmax = 7880000 |
---|
201 | |
---|
202 | # for ferret2sww |
---|
203 | south = degminsec2decimal_degrees(-20,30,0) |
---|
204 | north = degminsec2decimal_degrees(-17,10,0) |
---|
205 | west = degminsec2decimal_degrees(117,00,0) |
---|
206 | east = degminsec2decimal_degrees(120,00,0) |
---|
207 | |
---|
208 | # region to export (used from export_results.py) |
---|
209 | e_min_area = 648000 |
---|
210 | e_max_area = 675000 |
---|
211 | n_min_area = 7745000 |
---|
212 | n_max_area = 7761000 |
---|
213 | |
---|
214 | export_region = [[e_min_area,n_min_area],[e_min_area,n_max_area],[e_max_area,n_max_area],[e_max_area,n_min_area]] |
---|
215 | |
---|
216 | refzone = 50 |
---|
217 | |
---|
218 | from anuga.coordinate_transforms.redfearn import redfearn |
---|
219 | # boundary up to 50 m contour |
---|
220 | lat1_50 = degminsec2decimal_degrees(-19,20,0) |
---|
221 | lat2_50 = degminsec2decimal_degrees(-19,30,0) |
---|
222 | lat3_50 = degminsec2decimal_degrees(-19,45,0) |
---|
223 | lon1_50 = degminsec2decimal_degrees(119,05,0) |
---|
224 | lon2_50 = degminsec2decimal_degrees(118,20,0) |
---|
225 | lon3_50 = degminsec2decimal_degrees(117,45,0) |
---|
226 | z, easting, northing = redfearn(lat1_50, lon1_50) |
---|
227 | d0_50 = [easting, northing] |
---|
228 | z, easting, northing = redfearn(lat2_50, lon2_50) |
---|
229 | d1_50 = [easting, northing] |
---|
230 | z, easting, northing= redfearn(lat3_50, lon3_50) |
---|
231 | d2_50 = [easting, northing] |
---|
232 | |
---|
233 | d4_50 = [285000, 7585000] |
---|
234 | d6_50 = [330000, 7605000] |
---|
235 | #bounding_poly50 = [p0_50, p1_50, p2_50, d6_50, d5, d4_50] |
---|
236 | |
---|
237 | d0 = [763852.0, 7934358.0] |
---|
238 | d1 = [710987.0, 7925797.0] |
---|
239 | d2 = [658264.0, 7926314.0] |
---|
240 | d3 = [552686.0, 7871580.0] |
---|
241 | #d4 = [604415.81, 7733013.56] |
---|
242 | d4 = [638000.0, 7733013.56] |
---|
243 | #d5 = [656561.15, 7732615.11] |
---|
244 | d5 = [662000.0, 7732615.11] |
---|
245 | #d6 = [708940.32, 7750510.33] |
---|
246 | d6 = [690000.0, 7740510.33] |
---|
247 | #polyAll = [d0, d1, d2, d3, d4, d5, d6] |
---|
248 | #polyAll = [d0_50, d1_50, d2_50, d4, d5, d6] |
---|
249 | # from Hamish |
---|
250 | h0=[629262.17, 7747205.47] |
---|
251 | h1=[552686.00, 7871579.99] #d3 |
---|
252 | h2=[658264.00, 7926314.00] #d2 |
---|
253 | h3=[710986.99, 7925796.99] #d1 |
---|
254 | h4=[763851.99, 7934357.99] #d0 |
---|
255 | h5=[701485.21, 7770656.86] |
---|
256 | h6=[698273.75, 7762227.38] |
---|
257 | h7=[698194.23, 7762018.65] |
---|
258 | h8=[691627.41, 7744781.98] |
---|
259 | h9=[679220.75, 7743604.59] |
---|
260 | h10=[653512.59, 7740528.56] |
---|
261 | h11=[634777.71, 7738247.17] |
---|
262 | h12=[629443.86, 7746910.37] |
---|
263 | h13=[629396.84, 7746986.75] |
---|
264 | h14=[629352.32, 7747059.06] |
---|
265 | h15=[629276.24, 7747182.63] |
---|
266 | h16=[629262.17, 7747205.47] #repeat of h0 |
---|
267 | # using Hamish's new bounding polygon |
---|
268 | #polyAll = [d0_50, d1_50, d2_50, h16,h15,h14,h13,h12,h11,h10,h9,h8,h7,h6,h5] |
---|
269 | polyAll = [d0_50, d1_50, d2_50, h16,h11,h8,h6, h5] |
---|
270 | |
---|
271 | #Interior region - Pt Hedland town |
---|
272 | i0 = [668000, 7757000] |
---|
273 | i1 = [659000, 7755000] |
---|
274 | i2 = [660000, 7749000] |
---|
275 | i3 = [667000, 7746000] |
---|
276 | i4 = [678000, 7751000] |
---|
277 | |
---|
278 | poly_pt_hedland = [i0, i1, i2, i3, i4] |
---|
279 | |
---|
280 | #Are there other significant features? |
---|
281 | j0 = [670000, 7760000] |
---|
282 | j1 = [633000, 7745000] |
---|
283 | j2 = [665000, 7743000] |
---|
284 | j3 = [690000, 7755000] |
---|
285 | |
---|
286 | poly_region = [j0, j1, j2, j3] |
---|