source: anuga_work/production/karratha_2006/project.py @ 3669

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

reflecting PATH change in project scripts

File size: 6.2 KB
Line 
1"""Common filenames and locations for topographic data, meshes and outputs.
2"""
3
4from os import sep, environ, getenv, getcwd
5from os.path import expanduser, basename
6#from anuga.utilities.polygon import read_polygon
7import sys
8from anuga.coordinate_transforms.redfearn import\
9     degminsec2decimal_degrees,\
10     convert_points_from_latlon_to_utm
11
12from time import localtime, strftime
13from anuga.geospatial_data.geospatial_data import *
14
15# Location and naming of scenario data
16state = 'western_australia'
17scenario_dir_name = 'karratha_tsunami_scenario_2005'
18
19basename = 'elevation50m'
20boundary_basename = 'SU-AU'
21
22
23
24# Data output
25codename = 'project.py'
26
27home = getenv('INUNDATIONHOME') #Sandpit's parent dir   
28
29if sys.platform == 'win32':
30    user = getenv('USERPROFILE')
31else:
32    user = getenv('LOGNAME')
33
34# INUNDATIONHOME is the inundation directory, not the data directory.
35home += sep +'data'
36
37#Derive subdirectories and filenames
38time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
39outputtimedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'output'+sep+time+sep
40
41meshdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'meshes'+sep
42datadir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'topographies'+sep
43gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep
44polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
45boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
46outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'output'+sep
47tidedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'tide_data'+sep
48
49#gauge_filename = gaugedir + 'gauge_location_broome.csv'
50
51meshname = meshdir + basename
52outputname = outputtimedir + basename  #Used by post processing
53#boundaryname = outputtimedir + basename  #Used by post processing
54
55
56
57# Regions
58
59refzone = 50 
60south = degminsec2decimal_degrees(-20,55,0)
61north = degminsec2decimal_degrees(-20,15,0)
62west = degminsec2decimal_degrees(116,17,0)
63east = degminsec2decimal_degrees(117,10,0)
64
65p0 = [south, degminsec2decimal_degrees(116,32,0)]
66p1 = [south, west]
67p2 = [degminsec2decimal_degrees(-20,23,0), west]
68p3 = [north, degminsec2decimal_degrees(116,45,0)]
69p4 = [north, degminsec2decimal_degrees(117,0,0)]
70p5 = [p2[0], degminsec2decimal_degrees(117,8,0)]
71p6 = [degminsec2decimal_degrees(-20,30,0), east]
72p7 = [degminsec2decimal_degrees(-20,38,0), east]
73p8 = [south, east]
74
75
76bounding_polygon, zone = convert_points_from_latlon_to_utm([p0, p1, p2, p3, p4, p5, p6, p7, p8])
77refzone = zone
78
79
80#Interior regions -  NOT CONVERTED YET TO UTM
81karratha_south = degminsec2decimal_degrees(-20,44,0)
82karratha_north = degminsec2decimal_degrees(-20,42,0)
83karratha_west = degminsec2decimal_degrees(116,48,0)
84karratha_east = degminsec2decimal_degrees(116,53,30)
85
86k0 = [karratha_south, karratha_west]
87k1 = [karratha_south, karratha_east]
88k2 = [karratha_north, karratha_east]
89k3 = [karratha_north, karratha_west]   
90
91karratha_polygon = [k0, k1, k2, k3]
92
93
94#Interior regions
95dampier_south = degminsec2decimal_degrees(-20,40,0)
96dampier_north = degminsec2decimal_degrees(-20,38,10)
97dampier_west = degminsec2decimal_degrees(116,43,0)
98dampier_east = degminsec2decimal_degrees(116,45,0)
99
100d0 = [dampier_south, dampier_west]
101d1 = [dampier_south, dampier_east]
102d2 = [dampier_north, dampier_east]
103d3 = [dampier_north, dampier_west]   
104
105dampier_polygon = [d0, d1, d2, d3]
106
107
108#Interior regions
109refinery_south = degminsec2decimal_degrees(-20,37,50)
110refinery_north = degminsec2decimal_degrees(-20,36,0)
111refinery_west = degminsec2decimal_degrees(116,44,0)
112refinery_east = degminsec2decimal_degrees(116,46,10)
113
114d0 = [refinery_south, refinery_west]
115d1 = [refinery_south, refinery_east]
116d2 = [refinery_north, refinery_east]
117d3 = [refinery_north, refinery_west]   
118
119refinery_polygon = [d0, d1, d2, d3]
120
121
122#Interior region around 468899, 7715177:
123#lat (-20, 39, 44.93753), lon (116, 42, 5.09106)
124
125point_south = degminsec2decimal_degrees(-20,39,46)
126point_north = degminsec2decimal_degrees(-20,39,42)
127point_west = degminsec2decimal_degrees(116,42,0)
128point_east = degminsec2decimal_degrees(116,42,10)
129
130d0 = [point_south, point_west]
131d1 = [point_south, point_east]
132d2 = [point_north, point_east]
133d3 = [point_north, point_west]   
134
135point_polygon = [d0, d1, d2, d3]
136
137
138#Neils areas around interesting points
139neil1_point1 = [degminsec2decimal_degrees(-20,35,34),
140                degminsec2decimal_degrees(116,45,18)]
141neil1_point2 = [degminsec2decimal_degrees(-20,36,15),
142                degminsec2decimal_degrees(116,46,18)]
143neil1_point3 = [degminsec2decimal_degrees(-20,35,9),
144                degminsec2decimal_degrees(116,47,17)]
145neil1_point4 = [degminsec2decimal_degrees(-20,34,26),
146                degminsec2decimal_degrees(116,46,17)]
147
148neil1_polygon = [neil1_point1, neil1_point2, neil1_point3, neil1_point4]
149
150
151
152neil2_point1 = [degminsec2decimal_degrees(-20,39,36),
153                degminsec2decimal_degrees(116,41,33)]
154neil2_point2 = [degminsec2decimal_degrees(-20,40,10),
155                degminsec2decimal_degrees(116,42,13)]
156neil2_point3 = [degminsec2decimal_degrees(-20,38,39),
157                degminsec2decimal_degrees(116,43,49)]
158neil2_point4 = [degminsec2decimal_degrees(-20,38,5),
159                degminsec2decimal_degrees(116,43,9)]
160
161neil2_polygon = [neil2_point1, neil2_point2, neil2_point3, neil2_point4]
162
163
164
165
166#Withnell bay
167wb_point1 = [degminsec2decimal_degrees(-20,35,34),
168                degminsec2decimal_degrees(116,45,18)]
169wb_point2 = [degminsec2decimal_degrees(-20,36,15),
170                degminsec2decimal_degrees(116,46,18)]
171wb_point3 = [degminsec2decimal_degrees(-20,35,9),
172                degminsec2decimal_degrees(116,47,17)]
173wb_point4 = [degminsec2decimal_degrees(-20,34,26),
174                degminsec2decimal_degrees(116,46,17)]
175
176wb_polygon = [wb_point1, wb_point2, wb_point3, wb_point4]
177
178
179
180
181#Larger Withnell bay
182lwb_point1 = [degminsec2decimal_degrees(-20,35,59),
183                degminsec2decimal_degrees(116,42,00)]
184lwb_point2 = [degminsec2decimal_degrees(-20,36,50),
185                degminsec2decimal_degrees(116,46,50)]
186lwb_point3 = [degminsec2decimal_degrees(-20,34,00),
187                degminsec2decimal_degrees(116,47,39)]
188lwb_point4 = [degminsec2decimal_degrees(-20,33,00),
189                degminsec2decimal_degrees(116,42,50)]
190
191lwb_polygon = [lwb_point1, lwb_point2, lwb_point3, lwb_point4]
192
193
194
Note: See TracBrowser for help on using the repository browser.