source: anuga_work/production/pt_hedland_2006/project.py @ 4347

Last change on this file since 4347 was 4045, checked in by sexton, 18 years ago

updates to pt hedland script (new data provided) and sydney slide scenario

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