source: anuga_work/production/dampier_2006/project.py @ 3870

Last change on this file since 3870 was 3870, checked in by nick, 17 years ago

dampier updates

File size: 10.5 KB
Line 
1"""Common filenames and locations for topographic data, meshes and outputs.
2"""
3
4import sys
5from os import sep, environ, getenv, getcwd
6from os.path import expanduser, basename
7
8from anuga.coordinate_transforms.redfearn import\
9     degminsec2decimal_degrees,\
10     convert_from_latlon_to_utm
11
12from time import localtime, strftime
13from anuga.geospatial_data.geospatial_data import *
14
15# file and system info
16#---------------------------------
17codename = 'project.py'
18
19home = getenv('INUNDATIONHOME') #Sandpit's parent dir   
20
21if sys.platform == 'win32':
22    user = getenv('USERPROFILE')
23else:
24    user = getenv('LOGNAME')
25
26# INUNDATIONHOME is the inundation directory, not the data directory.
27home += sep +'data'
28#----------------------------------
29# Location and naming of scenario data
30#----------------------------------
31state = 'western_australia'
32scenario_name = 'dampier_tsunami'
33scenario_datas_name = 'dampier_tsunami_scenario_2006'  #name of the directory where the data is stored
34#scenario_datas_name = 'karratha_tsunami_scenario_2005' # Tmp location
35
36#mesh_name = 'elevation50m'
37boundaries_name = 'dampier'
38boundaries_source = 'test'
39
40# topography file names
41onshore_name = 'dli_no_islands'
42coast_name = 'DTED_05_Contour'
43islands_name = 'dted_islands'
44offshore_name = 'XY100003902'
45offshore_name1 = 'XY100003903'
46offshore_name2 = 'XY100003951'
47offshore_name3 = 'XY100006321'
48offshore_name4 = 'XY100011756'
49offshore_name5 = 'XY100014243'
50offshore_name6 = 'XY100014244'
51offshore_name7 = 'XY100021081'
52offshore_name8 = 'XY100021082'
53offshore_name9 = 'XY100021083'
54offshore_name10 = 'XY100021085'
55offshore_name11 = 'XY100021086'
56offshore_name12 = 'XY100026309'
57offshore_name13 = 'XY100026338'
58offshore_name14 = 'XYDM83'
59
60combined_name ='dampier_combined_elevation'
61
62
63
64#Derive subdirectories and filenames
65
66meshes_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'meshes'+sep
67topographies_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'topographies'+sep
68gauges_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'gauges'+sep
69polygons_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'polygons'+sep
70boundaries_in_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'boundaries'+sep+'urs'+sep+boundaries_source+sep
71#outputdir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'output'+sep
72tide_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'tide_data'+sep
73
74time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
75build_time = time+'_build'
76run_time = time+'_run'
77
78output_build_time_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'outputs'+sep+build_time+sep
79output_run_time_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'outputs'+sep+run_time+sep
80topographies_time_dir = topographies_dir+build_time+sep
81boundaries_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'boundaries'+sep
82boundaries_time_dir = home+sep+state+sep+scenario_datas_name+sep+'anuga'+sep+'boundaries'+sep+build_time+sep
83
84#ideas
85#boundaries_time_dir = boundaries_in_dir+'urs'+sep+boundaries_source+sep
86
87
88#gauge_filename = gaugedir + 'gauge_location_broome.csv'
89
90onshore_dir_name = topographies_dir + onshore_name
91coast_dir_name = topographies_dir + coast_name
92islands_dir_name = topographies_dir + islands_name
93offshore_dir_name = topographies_dir + offshore_name
94offshore_dir_name1 = topographies_dir + offshore_name1
95offshore_dir_name2 = topographies_dir + offshore_name2
96offshore_dir_name3 = topographies_dir + offshore_name3
97offshore_dir_name4 = topographies_dir + offshore_name4
98offshore_dir_name5 = topographies_dir + offshore_name5
99offshore_dir_name6 = topographies_dir + offshore_name6
100offshore_dir_name7 = topographies_dir + offshore_name7
101offshore_dir_name8 = topographies_dir + offshore_name8
102offshore_dir_name9 = topographies_dir + offshore_name9
103offshore_dir_name10 = topographies_dir + offshore_name10
104offshore_dir_name11 = topographies_dir + offshore_name11
105offshore_dir_name12 = topographies_dir + offshore_name12
106offshore_dir_name13 = topographies_dir + offshore_name13
107offshore_dir_name14 = topographies_dir + offshore_name14
108
109#output dir
110combined_dir_name = topographies_dir + combined_name
111combined_time_dir_name = topographies_time_dir + combined_name
112
113meshes_dir_name = meshes_dir + scenario_name
114
115#output_build_time_dir_name = output_build_time_dir + scenario_name  #Used by post processing
116output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
117boundaries_in_dir_name = boundaries_in_dir + boundaries_name
118boundaries_time_dir_name = boundaries_time_dir + boundaries_name  #Used by post processing
119
120
121
122
123# Regions
124
125refzone = 50 
126south = degminsec2decimal_degrees(-20,55,0)
127north = degminsec2decimal_degrees(-20,15,0)
128#north = degminsec2decimal_degrees(-19,15,0)
129west = degminsec2decimal_degrees(116,17,0)
130east = degminsec2decimal_degrees(117,10,0)
131
132#only used to clip boundary condition
133north_boundary = north + 0.04
134south_boundary = south - 0.04
135west_boundary = west - 0.04
136east_boundary = east + 0.04
137
138
139p0 = [south, degminsec2decimal_degrees(116,32,0)]
140p1 = [south, west]
141p2 = [degminsec2decimal_degrees(-20,23,0), west]
142p3 = [north, degminsec2decimal_degrees(116,45,0)]
143p4 = [north, degminsec2decimal_degrees(117,0,0)]
144p5 = [p2[0], degminsec2decimal_degrees(117,8,0)]
145p6 = [degminsec2decimal_degrees(-20,30,0), east]
146p7 = [degminsec2decimal_degrees(-20,38,0), east]
147p8 = [south, east]
148
149bounding_polygon, zone =\
150                  convert_from_latlon_to_utm([p0, p1, p2, p3, p4,
151                                              p5, p6, p7, p8])
152refzone = zone
153
154
155#Interior regions
156
157# CIPMA point of interest
158cipma_latitude = -20.588456
159cipma_longitude = 116.771527
160
161
162k0 = [cipma_latitude-0.02, cipma_longitude-0.02]
163k1 = [cipma_latitude-0.02, cipma_longitude+0.02]
164k2 = [cipma_latitude+0.02, cipma_longitude+0.02]
165k3 = [cipma_latitude+0.02, cipma_longitude-0.02]
166
167cipma_polygon, zone = convert_from_latlon_to_utm([k0, k1, k2, k3])
168assert zone == refzone
169
170e_min_area = 474000
171e_max_area = 480000
172n_min_area = 7719000
173n_max_area = 7725000
174
175#Interior regions
176karratha_south = degminsec2decimal_degrees(-20,44,0)
177karratha_north = degminsec2decimal_degrees(-20,42,0)
178karratha_west = degminsec2decimal_degrees(116,48,0)
179karratha_east = degminsec2decimal_degrees(116,53,30)
180
181k0 = [karratha_south, karratha_west]
182k1 = [karratha_south, karratha_east]
183k2 = [karratha_north, karratha_east]
184k3 = [karratha_north, karratha_west]   
185
186karratha_polygon, zone = convert_from_latlon_to_utm([k0, k1, k2, k3])
187assert zone == refzone
188
189
190#Interior regions
191dampier_south = degminsec2decimal_degrees(-20,40,0)
192dampier_north = degminsec2decimal_degrees(-20,38,10)
193dampier_west = degminsec2decimal_degrees(116,43,0)
194dampier_east = degminsec2decimal_degrees(116,45,0)
195
196d0 = [dampier_south, dampier_west]
197d1 = [dampier_south, dampier_east]
198d2 = [dampier_north, dampier_east]
199d3 = [dampier_north, dampier_west]   
200
201dampier_polygon, zone = convert_from_latlon_to_utm([d0, d1, d2, d3])
202assert zone == refzone
203
204
205#Interior regions
206refinery_south = degminsec2decimal_degrees(-20,37,50)
207refinery_north = degminsec2decimal_degrees(-20,36,0)
208refinery_west = degminsec2decimal_degrees(116,44,0)
209refinery_east = degminsec2decimal_degrees(116,46,10)
210
211d0 = [refinery_south, refinery_west]
212d1 = [refinery_south, refinery_east]
213d2 = [refinery_north, refinery_east]
214d3 = [refinery_north, refinery_west]   
215
216refinery_polygon, zone = convert_from_latlon_to_utm([d0, d1, d2, d3])
217assert zone == refzone
218
219
220#Interior region around 468899, 7715177:
221#lat (-20, 39, 44.93753), lon (116, 42, 5.09106)
222
223point_south = degminsec2decimal_degrees(-20,39,46)
224point_north = degminsec2decimal_degrees(-20,39,42)
225point_west = degminsec2decimal_degrees(116,42,0)
226point_east = degminsec2decimal_degrees(116,42,10)
227
228d0 = [point_south, point_west]
229d1 = [point_south, point_east]
230d2 = [point_north, point_east]
231d3 = [point_north, point_west]   
232
233point_polygon, zone = convert_from_latlon_to_utm([d0, d1, d2, d3])
234assert zone == refzone
235
236
237#Neils areas around interesting points
238neil1_point1 = [degminsec2decimal_degrees(-20,35,34),
239                degminsec2decimal_degrees(116,45,18)]
240neil1_point2 = [degminsec2decimal_degrees(-20,36,15),
241                degminsec2decimal_degrees(116,46,18)]
242neil1_point3 = [degminsec2decimal_degrees(-20,35,9),
243                degminsec2decimal_degrees(116,47,17)]
244neil1_point4 = [degminsec2decimal_degrees(-20,34,26),
245                degminsec2decimal_degrees(116,46,17)]
246
247neil1_polygon, zone = convert_from_latlon_to_utm([neil1_point1,
248                                                  neil1_point2,
249                                                  neil1_point3,
250                                                  neil1_point4])
251assert zone == refzone
252
253
254
255
256neil2_point1 = [degminsec2decimal_degrees(-20,39,36),
257                degminsec2decimal_degrees(116,41,33)]
258neil2_point2 = [degminsec2decimal_degrees(-20,40,10),
259                degminsec2decimal_degrees(116,42,13)]
260neil2_point3 = [degminsec2decimal_degrees(-20,38,39),
261                degminsec2decimal_degrees(116,43,49)]
262neil2_point4 = [degminsec2decimal_degrees(-20,38,5),
263                degminsec2decimal_degrees(116,43,9)]
264
265neil2_polygon, zone = convert_from_latlon_to_utm([neil2_point1,
266                                                  neil2_point2,
267                                                  neil2_point3,
268                                                  neil2_point4])
269assert zone == refzone
270
271
272
273
274
275#Withnell bay
276wb_point1 = [degminsec2decimal_degrees(-20,35,34),
277                degminsec2decimal_degrees(116,45,18)]
278wb_point2 = [degminsec2decimal_degrees(-20,36,15),
279                degminsec2decimal_degrees(116,46,18)]
280wb_point3 = [degminsec2decimal_degrees(-20,35,9),
281                degminsec2decimal_degrees(116,47,17)]
282wb_point4 = [degminsec2decimal_degrees(-20,34,26),
283                degminsec2decimal_degrees(116,46,17)]
284
285wb_polygon, zone = convert_from_latlon_to_utm([wb_point1, wb_point2,
286                                               wb_point3, wb_point4])
287assert zone == refzone
288
289
290
291
292
293#Larger Withnell bay
294lwb_point1 = [degminsec2decimal_degrees(-20,35,59),
295                degminsec2decimal_degrees(116,42,00)]
296lwb_point2 = [degminsec2decimal_degrees(-20,36,50),
297                degminsec2decimal_degrees(116,46,50)]
298lwb_point3 = [degminsec2decimal_degrees(-20,34,00),
299                degminsec2decimal_degrees(116,47,39)]
300lwb_point4 = [degminsec2decimal_degrees(-20,33,00),
301                degminsec2decimal_degrees(116,42,50)]
302
303lwb_polygon, zone = convert_from_latlon_to_utm([lwb_point1, lwb_point2,
304                                                lwb_point3, lwb_point4])
305                                                     
306assert zone == refzone
307
308
309
310
Note: See TracBrowser for help on using the repository browser.