source: anuga_work/production/busselton/project.py @ 5444

Last change on this file since 5444 was 5444, checked in by kristy, 16 years ago
File size: 7.0 KB
Line 
1# -*- coding: cp1252 -*-
2"""Common filenames and locations for topographic data, meshes and outputs.
3"""
4
5from os import sep, environ, getenv, getcwd ,umask
6from os.path import expanduser
7import sys
8from time import localtime, strftime, gmtime
9from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles
10#from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
11from anuga.utilities.system_tools import get_user_name, get_host_name
12
13# file and system info
14#---------------------------------
15#codename = 'project.py'
16
17home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir   
18user = get_user_name()
19host = get_host_name()
20
21# INUNDATIONHOME is the inundation directory, not the data directory.
22
23#time stuff
24time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
25gtime = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
26build_time = time+'_build'
27run_time = time+'_run'
28print 'gtime: ', gtime
29
30#Making assumptions about the location of scenario data
31state = 'western_australia'
32scenario_name = 'busselton'
33scenario = 'busselton_tsunami_scenario'
34
35tide = 0.6
36
37alpha = 0.1 #changed for testing purposes
38friction=0.01
39starttime=0
40finaltime=80000
41export_cellsize=25
42setup='final'
43source='exmouth'
44
45
46if setup =='trial':
47    print'trial'
48    res_factor=10
49    time_thinning=48
50    yieldstep=240
51if setup =='basic': 
52    print'basic'
53    res_factor=4
54    time_thinning=12
55    yieldstep=120
56if setup =='final': 
57    print'final'
58    res_factor=1
59    time_thinning=4
60    yieldstep=60
61
62rev_num = 'current'
63##rev_num = '4695' # 2nd Sept 2007
64##rev_num = '4743' # 3nd Oct 2007
65##rev_num = '4777' # 1st Nov 2007
66##rev_num = '4874' # 3rd Dec 2007
67##rev_num = '4901' # 3rd Jan 2007
68##rev_num = '4990' # 5th Feb 2007
69##rev_num = '5103' # 3rd March 2007
70##rev_num = '5185' # 1st April 2007
71##rev_num = '5273' # 2nd May 2007
72
73
74dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+str(rev_num)+'_'+str(user)
75
76
77# onshore data provided by WA DLI - provided by Hamish on the 17th June 2008
78
79onshore_name = 'busselton_v2_gda94_mga50' # original
80
81# AHO + DPI data
82coast_name = 'Contour_0' # provided by hamish, represent better coastline than the 100km as compared to charts
83coast_name1 = 'beach_survey_final'
84offshore_name = 'Bathymetry_final'
85offshore_name1 = 'Busselton_Chart'
86offshore_name2 = 'Busselton_Digitised'
87
88#final topo name
89combined_name ='busselton_combined_elevation'
90combined_name_small = 'busselton_combined_elevation_smaller'
91
92anuga_dir = home+state+sep+scenario+sep+'anuga'+sep
93
94topographies_in_dir = home+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep
95topographies_dir = anuga_dir+'topographies'+sep
96
97# input topo file location
98onshore_in_dir_name = topographies_in_dir + onshore_name #topo
99
100coast_in_dir_name = topographies_in_dir + coast_name #coastline
101coast_in_dir_name1 = topographies_in_dir + coast_name1 #beach survey
102
103offshore_in_dir_name = topographies_in_dir + offshore_name #bathymetry
104offshore_in_dir_name1 = topographies_in_dir + offshore_name1 #bathymetry Charts
105offshore_in_dir_name2 = topographies_in_dir + offshore_name2 #Digitised Fairsheet
106
107#output to anuga from build file
108onshore_dir_name = topographies_dir + onshore_name
109
110coast_dir_name = topographies_dir + coast_name
111coast_dir_name1 = topographies_dir + coast_name1
112
113offshore_dir_name = topographies_dir + offshore_name
114offshore_dir_name1 = topographies_dir + offshore_name1
115offshore_dir_name2 = topographies_dir + offshore_name2
116
117#final topo files
118combined_dir_name = topographies_dir + combined_name
119combined_dir_name_small = topographies_dir + combined_name_small
120
121meshes_dir = anuga_dir+'meshes'+sep
122meshes_dir_name = meshes_dir + scenario_name
123
124polygons_dir = anuga_dir+'polygons'+sep
125tide_dir = anuga_dir+'tide_data'+sep
126
127#boundaries_source = '1'
128
129if source=='exmouth':
130    boundaries_name = 'busselton_3103_30052008' #exmouth gun
131    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
132
133if source=='test':
134    boundaries_name = 'other' #exmouth gun
135    boundaries_in_dir = anuga_dir+'boundaries'+sep
136
137
138#boundaries locations
139boundaries_in_dir_name = boundaries_in_dir + boundaries_name
140boundaries_dir = anuga_dir+'boundaries'+sep
141boundaries_dir_name = boundaries_dir + scenario_name
142
143#output locations
144output_dir = anuga_dir+'outputs'+sep
145output_build_time_dir = output_dir +'test'+sep + build_time + dir_comment + sep
146#output_run_time_dir = output_dir +run_time+dir_comment+sep
147output_run_time_dir = output_dir +'test'+sep + run_time + dir_comment +sep # take test out
148output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
149
150#gauges
151gauge_name = 'busselton.csv'
152gauges_dir = home+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
153gauges_dir_name = gauges_dir + gauge_name
154
155buildings_filename = gauges_dir + 'Busselton_res_Project.csv'
156buildings_filename_out = 'Busselton_res_Project_modified.csv'
157
158community_filename = gauges_dir +''
159community_broome = gauges_dir + ''
160
161
162###############################
163# Domain definitions
164###############################
165
166from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
167
168# bounding polygon for study area
169poly_all = read_polygon(polygons_dir+'poly_all.csv')
170res_poly_all = 100000*res_factor
171
172
173###############################
174# Interior region definitions
175###############################
176
177#digitized polygons
178poly_large = read_polygon(polygons_dir+'coast_5km_d20m.csv')
179res_large = 40000*res_factor
180
181poly_busselton = read_polygon(polygons_dir+'busselton_1km.csv')
182res_busselton = 500*res_factor
183
184poly_bunbury = read_polygon(polygons_dir+'bunbury_1km.csv')
185res_bunbury = 500*res_factor
186
187poly_dunsborough = read_polygon(polygons_dir+'dunsborough_1km.csv')
188res_dunsborough = 500*res_factor
189
190poly_busselton2 = read_polygon(polygons_dir+'busselton_2km.csv')
191res_busselton2 = 10000*res_factor
192
193poly_bunbury2 = read_polygon(polygons_dir+'bunbury_2km.csv')
194res_bunbury2 = 10000*res_factor
195
196
197interior_regions = [[poly_large,res_large],[poly_busselton,res_busselton],[poly_bunbury,res_bunbury]
198                    ,[poly_dunsborough,res_dunsborough],[poly_busselton2,res_busselton2]
199                    ,[poly_bunbury2,res_bunbury2]]
200
201
202boundary_tags={'back': [0, 1, 2, 3, 4, 9],
203               'side': [5], 'ocean': [6, 7, 8]}
204
205trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
206
207poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')
208
209print 'min number triangles', trigs_min
210
211
212###################################################################
213# Clipping regions for export to asc and regions for clipping data
214###################################################################
215
216# exporting asc grid for Busselton
217xminBusselton = 340000
218xmaxBusselton = 352000
219yminBusselton = 6271500
220ymaxBusselton = 6280000
221
222# exporting asc grid for Bunbury
223xminBunbury = 369000
224xmaxBunbury = 381000
225yminBunbury = 6308000
226ymaxBunbury = 6308000
227
228# exporting asc grid for Dunsborough
229xminDunsborough = 321000
230xmaxDunsborough = 327500
231yminDunsborough = 6277000
232ymaxDunsborough = 6282000
233
Note: See TracBrowser for help on using the repository browser.