source: anuga_work/production/broome_2006/project.py @ 4477

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

update project file to include newly acquired residential building data from EVP

File size: 6.1 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
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
12
13# file and system info
14#---------------------------------
15codename = 'project.py'
16
17home = getenv('INUNDATIONHOME') #Sandpit's parent dir   
18user = get_user_name()
19
20# INUNDATIONHOME is the inundation directory, not the data directory.
21home += sep +'data'
22
23#time stuff
24time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
25#time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
26gtime = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
27build_time = time+'_build'
28run_time = time+'_run'
29print 'gtime: ', gtime
30
31#tide = -5.3
32tide = 0
33#tide = 4.9
34
35#Making assumptions about the location of scenario data
36state = 'western_australia'
37scenario_name = 'broome'
38scenario = 'broome_tsunami_scenario_2006'
39
40# onshore data provided by WA DLI
41onshore_name = 'Broome_topography_DLI' # original
42
43#island
44#island_name = 'rott_dli_ext' # original
45
46# offshore
47coast_name = 'Broome_coastline'
48
49offshore_name = 'Broome_Bathymetry'
50offshore_name1 = 'inferred_north'
51offshore_name2 = 'inferred_south'
52
53#final topo name
54combined_name ='broome_combined_elevation'
55combined_name1 ='broome_combined_elevation1'
56combined_name_unclipped1 ='broome_combined_elevation_unclipped1'
57combined_small_name = 'broome_combined_elevation_small'
58
59topographies_in_dir = home+sep+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep+'Feb07'+sep
60topographies_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep
61topographies_time_dir = topographies_dir+build_time+sep
62
63#input topo file location
64onshore_in_dir_name = topographies_in_dir + onshore_name
65#island_in_dir_name = topographies_in_dir + island_name
66
67coast_in_dir_name = topographies_in_dir + coast_name
68
69offshore_in_dir_name = topographies_in_dir + offshore_name
70offshore_in_dir_name1 = topographies_in_dir + offshore_name1
71offshore_in_dir_name2 = topographies_in_dir + offshore_name2
72
73onshore_dir_name = topographies_dir + onshore_name
74
75#island_dir_name = topographies_dir + island_name
76
77coast_dir_name = topographies_dir + coast_name
78
79offshore_dir_name = topographies_dir + offshore_name
80offshore_dir_name1 = topographies_dir + offshore_name1
81offshore_dir_name2 = topographies_dir + offshore_name2
82
83#final topo files
84combined_dir_name = topographies_dir + combined_name
85combined_dir_name_unclipped1 = topographies_dir + combined_name_unclipped1
86combined_dir_name1 = topographies_dir + combined_name1
87combined_small_name_dir = topographies_dir + combined_small_name
88
89meshes_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'meshes'+sep
90meshes_dir_name = meshes_dir + scenario_name
91
92polygons_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'polygons'+sep
93tide_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'tide_data'+sep
94
95
96boundaries_source = 'broome_3854_17042007'
97#boundaries locations
98boundaries_in_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep+'urs'+sep
99#boundaries_in_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep+'urs'+sep+boundaries_source+sep
100boundaries_in_dir_name = boundaries_in_dir + boundaries_source
101boundaries_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep
102boundaries_dir_name = boundaries_dir + boundaries_source
103#boundaries_time_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep+build_time+sep
104#boundaries_time_dir_name = boundaries_time_dir + boundaries_name  #Used by post processing
105
106#output locations
107output_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'outputs'+sep
108output_build_time_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'outputs'+sep+build_time+sep
109output_run_time_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'outputs'+sep+run_time+sep
110output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
111
112#gauges
113gauge_name = 'broome_gauges.csv'
114gauges_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
115gauges_dir_name = gauges_dir + gauge_name
116
117community_filename = gauges_dir + 'CHINS_v2.csv'
118community_broome = gauges_dir + 'community_broome.csv'
119
120buildings_filename = gauges_dir + 'Broome_res_Project.csv'
121buildings_filename_out = 'Broome_res_Project_modified.csv'
122
123###############################
124# Domain definitions
125###############################
126from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
127
128poly_all = read_polygon(polygons_dir+'extent_small.csv')
129#poly_all = read_polygon(polygons_dir+'extent.csv')
130res_factor = 2
131#res_poly_all = 500000
132res_poly_all = 150000*res_factor
133
134###############################
135# Interior region definitions
136###############################
137
138poly_0 = read_polygon(polygons_dir+'neg20_coast_contour_pts.csv')
139#res_0 = 100000
140res_0 = 20000*res_factor
141
142poly_1 = read_polygon(polygons_dir+'broome_north_coast_inside_extent.csv')
143#res_1 = 50000
144res_1 = 5000*res_factor
145
146poly_2 = read_polygon(polygons_dir+'broome_south_coast_inside_extent.csv')
147#res_2 = 50000
148res_2 = 5000*res_factor
149
150poly_3 = read_polygon(polygons_dir+'Broome_town_pts.csv')
151#res_3 = 20000
152res_3 = 2000*res_factor
153
154poly_4 = read_polygon(polygons_dir+'Broome_inner_town_pts.csv')
155#res_4 = 5000
156res_4 = 500*res_factor
157#assert zone == refzone
158
159interior_regions = [[poly_0,res_0],[poly_1,res_1],[poly_2,res_2]
160                     ,[poly_3,res_3],[poly_4,res_4]]
161
162trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
163
164print 'min number triangles', trigs_min
165
166poly_mainland = read_polygon(polygons_dir+'Initial_Condition.csv')
167
168###################################################################
169# Clipping regions for export to asc and regions for clipping data
170###################################################################
171
172# exporting asc grid
173e_min_area = 412000.0
174e_max_area = 423000.0
175n_min_area = 8007000.0
176n_max_area = 8022000.0
177
178
179
Note: See TracBrowser for help on using the repository browser.