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

Last change on this file since 5498 was 4580, checked in by sexton, 17 years ago

minor updates for gauge names plus plotting files

File size: 10.5 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'
129beach_gauges = gaugedir + 'beach_gauges.csv'
130gauge_checking = gaugedir + 'gauge_checking_test.csv'
131buildings_filename = gaugedir + 'Port_Hedland_res_Project.csv'
132buildings_filename_out = 'Port_Hedland_res_Project_modified.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
196community_filename = gaugedir + 'CHINS_v2.csv'
197community_pt_hedland = gaugedir + 'community_pt_hedland.csv'
198
199# clipping region to make DEM (pts file) from onshore data
200eastingmin = 594000
201eastingmax = 715000
202northingmin = 7720000
203northingmax = 7880000
204
205# for ferret2sww
206south = degminsec2decimal_degrees(-20,30,0)
207north = degminsec2decimal_degrees(-17,10,0)
208west = degminsec2decimal_degrees(117,00,0)
209east = degminsec2decimal_degrees(120,00,0)
210
211# region to export (used from export_results.py)
212e_min_area = 648000
213e_max_area = 675000
214n_min_area = 7745000
215n_max_area = 7761000
216
217export_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]]
218                 
219refzone = 50 
220
221from anuga.coordinate_transforms.redfearn import redfearn
222# boundary up to 50 m contour
223lat1_50 = degminsec2decimal_degrees(-19,20,0)
224lat2_50 = degminsec2decimal_degrees(-19,30,0)
225lat3_50 = degminsec2decimal_degrees(-19,45,0)
226lon1_50 = degminsec2decimal_degrees(119,05,0)
227lon2_50 = degminsec2decimal_degrees(118,20,0)
228lon3_50 = degminsec2decimal_degrees(117,45,0)
229z, easting, northing = redfearn(lat1_50, lon1_50)
230d0_50 = [easting, northing]
231z, easting, northing = redfearn(lat2_50, lon2_50)
232d1_50 = [easting, northing]
233z, easting, northing= redfearn(lat3_50, lon3_50)
234d2_50 = [easting, northing]
235
236d4_50 = [285000, 7585000]
237d6_50 = [330000, 7605000]
238#bounding_poly50 = [p0_50, p1_50, p2_50, d6_50, d5, d4_50]
239
240d0 = [763852.0, 7934358.0]
241d1 = [710987.0, 7925797.0]
242d2 = [658264.0, 7926314.0]
243d3 = [552686.0, 7871580.0]
244#d4 = [604415.81, 7733013.56]
245d4 = [638000.0, 7733013.56]
246#d5 = [656561.15, 7732615.11]
247d5 = [662000.0, 7732615.11]
248#d6 = [708940.32, 7750510.33]
249d6 = [690000.0, 7740510.33]
250#polyAll = [d0, d1, d2, d3, d4, d5, d6]
251#polyAll = [d0_50, d1_50, d2_50, d4, d5, d6]
252# from Hamish
253h0=[629262.17, 7747205.47]
254h1=[552686.00, 7871579.99] #d3
255h2=[658264.00, 7926314.00] #d2
256h3=[710986.99, 7925796.99] #d1
257h4=[763851.99, 7934357.99] #d0
258h5=[701485.21, 7770656.86]
259h6=[698273.75, 7762227.38]
260h7=[698194.23, 7762018.65]
261h8=[691627.41, 7744781.98]
262h9=[679220.75, 7743604.59]
263h10=[653512.59, 7740528.56]
264h11=[634777.71, 7738247.17]
265h12=[629443.86, 7746910.37]
266h13=[629396.84, 7746986.75]
267h14=[629352.32, 7747059.06]
268h15=[629276.24, 7747182.63]
269h16=[629262.17, 7747205.47] #repeat of h0
270# using Hamish's new bounding polygon
271#polyAll = [d0_50, d1_50, d2_50, h16,h15,h14,h13,h12,h11,h10,h9,h8,h7,h6,h5]
272polyAll = [d0_50, d1_50, d2_50, h16,h11,h8,h6, h5]
273
274#Interior region - Pt Hedland town
275i0 = [668000, 7757000]
276i1 = [659000, 7755000]
277i2 = [660000, 7749000]
278i3 = [667000, 7746000]
279i4 = [678000, 7751000]
280
281poly_pt_hedland = [i0, i1, i2, i3, i4]
282
283#Are there other significant features?
284j0 = [670000, 7760000]
285j1 = [633000, 7745000]
286j2 = [665000, 7743000]
287j3 = [690000, 7755000]
288
289poly_region = [j0, j1, j2, j3]
Note: See TracBrowser for help on using the repository browser.