source: anuga_work/production/hobart_2006/project.py @ 3683

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

updates (clipping 25m data into 50m grid, putting MOST boundary on relevant boundary segments, change saving unique vertices to False)

File size: 7.7 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
10from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
11
12if sys.platform == 'win32':
13    home = getenv('INUNDATIONHOME')
14    user = getenv('USERPROFILE')
15
16else:   
17    home = getenv('INUNDATIONHOME', sep+'d'+sep+'xrd'+sep+'gem'+sep+'2'+sep+'ramp'+sep+'risk_assessment_methods_project'+sep+'inundation')     
18    user = getenv('LOGNAME')
19    print 'USER:', user
20
21# INUNDATIONHOME is the inundation directory, not the data directory.
22home += sep +'data'
23
24#Making assumptions about the location of scenario data
25state = 'tasmania'
26scenario_dir_name = 'hobart_tsunami_scenario_2006'
27
28# data provided by Tas SES and checked by NM&I
29#onshore_name = 'hob3_topo' # original
30onshore_name = 'hob5_topo' # 12.5m grid and clipped to 100m elevation or 3000m from coast
31#onshore_name_25 = 'hob5_topo_25m' # 25m grid and clipped to 100m elevation or 3000m from coast
32onshore_name_25 = 'hob6_topo_25m' # 25m grid NOT clipped
33offshore_name_tas1 = 'derwent_2m'
34offshore_name_tas2 = 'derwent_5m'
35offshore_name_tas3 = 'south_east_tas' #actually this is AHO
36offshore_name_tas4 = 'hobart_1m'
37
38# AHO data and checked by NM&I
39offshore_name1 = 'xy100003760'
40offshore_name2 = 'xy100003761'
41offshore_name3 = 'xy100003762'
42offshore_name4 = 'xy100003907'
43offshore_name5 = 'xy100003908'
44offshore_name6 = 'xy100003909'
45offshore_name7 = 'xy100003910'
46offshore_name8 = 'xy100003932'
47offshore_name9 = 'xy100003933'
48offshore_name10 = 'xy100003934'
49offshore_name11 = 'xy100003935'
50offshore_name12 = 'xy100003936'
51offshore_name13 = 'xy100003964'
52offshore_name14 = 'xy100014250'
53offshore_name15 = 'xy100014253'
54offshore_name16 = 'xy100016142'
55
56# developed by NM&I
57coast_name = 'coastline_points'
58
59boundary_basename = 'puysegur'
60
61#swollen/ all data output
62basename = 'source'
63
64codename = 'project.py'
65
66#Derive subdirectories and filenames
67#time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
68local_time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
69meshdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'meshes'+sep
70datadir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'topographies'+sep
71gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep
72polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
73boundarydir = home+sep+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
74outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep
75outputtimedir = outputdir + local_time + sep
76polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
77
78gauge_filename = gaugedir + 'hobart_gauges.xya'
79buildings_filename = gaugedir + 'hobart_res.csv'
80buildings_filename_damage_out = 'hobart_res_modified.csv'
81gaugetimeseries = gaugedir + 'hobart'
82
83# boundary source data
84#MOST_dir = 'f:'+sep+'3'+sep+'ehn'+sep+'users'+sep+'davidb'+sep+'tsunami'+sep+'WA_project'+sep+'SU-AU_90'+sep+'most_2'+sep+'detailed'+sep
85
86codedir = getcwd()+sep                           
87codedirname = codedir + 'project.py'
88meshname = outputtimedir + 'mesh_' + basename
89
90onshore_dem_name = datadir + onshore_name
91onshore_dem_name_25 = datadir + onshore_name_25
92all_onshore_dem_name = datadir + 'combined_onshore'
93offshore_dem_name_local1 = datadir + offshore_name_tas1
94offshore_dem_name_local2 = datadir + offshore_name_tas2
95offshore_dem_name_local3 = datadir + offshore_name_tas3
96offshore_dem_name_local4 = datadir + offshore_name_tas4
97offshore_dem_name_aho1 = datadir + offshore_name1
98offshore_dem_name_aho2 = datadir + offshore_name2
99offshore_dem_name_aho3 = datadir + offshore_name3
100offshore_dem_name_aho4 = datadir + offshore_name4
101offshore_dem_name_aho5 = datadir + offshore_name5
102offshore_dem_name_aho6 = datadir + offshore_name6
103offshore_dem_name_aho7 = datadir + offshore_name7
104offshore_dem_name_aho8 = datadir + offshore_name8
105offshore_dem_name_aho9 = datadir + offshore_name9
106offshore_dem_name_aho10 = datadir + offshore_name10
107offshore_dem_name_aho11 = datadir + offshore_name11
108offshore_dem_name_aho12 = datadir + offshore_name12
109offshore_dem_name_aho13 = datadir + offshore_name13
110offshore_dem_name_aho14 = datadir + offshore_name14
111offshore_dem_name_aho15 = datadir + offshore_name15
112offshore_dem_name_aho16 = datadir + offshore_name16
113coast_dem_name = datadir + coast_name
114# addition once total grid delivered
115onshore_offshore_dem_name_25 = datadir + '25m_se_tas' #25m grid
116onshore_offshore_dem_name = datadir + '50m_se_tas' #50m grid
117# output names
118bruny_dem_name_25 = datadir + 'bruny_25_dem'
119hobart_dem_name_25 = datadir + 'hobart_25_dem'
120combined_dem_name = datadir + 'hobart_combined_elevation'
121combined_dem_name_2 = datadir + 'hobart_combined_elevation_2'
122
123outputname = outputtimedir + basename  #Used by post processing
124
125###############################
126# Domain definitions
127###############################
128
129# bounding box
130south = degminsec2decimal_degrees(-44,45,0)
131north = degminsec2decimal_degrees(-42,0,0)
132west = degminsec2decimal_degrees(146,45,0)
133east = degminsec2decimal_degrees(148,25,0) #degminsec2decimal_degrees(148,15,0)
134
135#Main Domain of Hobart: first run JS 18/9/06
136d0 = [south, west]
137d1 = [south, east]
138d2 = [north, east]
139d3 = [north, west]
140polyAll2, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
141print "polyAll2", polyAll2
142refzone = zone
143
144# Second run - bottom bright, topr, top, left
145#polyAll = [[520000, 5170000],[580000, 5170000],[580000, 5200000],[590000,5240000],[520000,5260000]]
146
147# Third run - afternoon Wed 27 Sep; surrounds -100m and 20mish elevation
148#polyAll = read_polygon(polygondir+'new_extent.csv')
149#plot_polygons([polyAll],'boundingpoly',verbose=False)
150
151# Mark run - morning Fri 29 Sep; surrounds -100m and 20mish elevation
152polyAll = read_polygon(polygondir+'new_extent.csv')
153plot_polygons([polyAll, polyAll2],'boundingpoly',verbose=False)
154
155
156###################################################################
157# Clipping regions for export to asc and regions for clipping data
158###################################################################
159
160# region to export for inundation map
161e_min_area = 520000#490000
162e_max_area = 580000#580000
163n_min_area = 5180000#5160000
164n_max_area = 5260000#5275000
165
166# clipping 12.5m onshore data set
167eastingmin = 520000
168eastingmax = 536000
169northingmin = 5245000
170northingmax = 5265000
171
172# clipping 25m data set - Hobart
173eastingmin25 = 524208.387
174eastingmax25 = 554867.24
175northingmin25 = 5229154.555
176northingmax25 = 5258511.857
177
178# clipping 25m data set - Bruny
179eastingmin25_2 = 523672.502
180eastingmax25_2 = 536020.057
181northingmin25_2 = 5203120.035
182northingmax25_2 = 5212052.309
183
184
185
186###############################
187# Interior region definitions
188###############################
189
190###Interior region - Hobart city area + Glenorchy, Kingston
191##i0 = [517000, 5267000]
192##i1 = [517000, 5255000]
193##i2 = [520000, 5250000]
194##i3 = [522000, 5239000]
195##i4 = [524000, 5238000]
196##i5 = [526000, 5236000]
197##i6 = [530000, 5244000]
198##i7 = [530000, 5250000]
199##i8 = [534000, 5254000]
200##i9 = [520000, 5270000]
201##
202##poly_hobart = [i0, i1, i2, i3, i4, i5, i6, i7, i8, i9]
203##
204### Tasman Peninsula
205##l0 = [550000, 5247000]
206##l1 = [550000, 5211000]
207##l2 = [583000, 5211000]
208##l3 = [583000, 5247000]
209##
210##poly_tasman_peninsula = [l0, l1, l2, l3]
211##
212### Bruny Island
213##poly_bruny = read_polygon(polygondir+'bruny.csv')
214
215
216# Hobart digitized polygons
217poly_hobart1 = read_polygon(polygondir+'Hob_poly1.csv')
218poly_hobart2 = read_polygon(polygondir+'Hob_poly2.csv')
219poly_hobart3 = read_polygon(polygondir+'Hob_poly3.csv')
220
221plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3],'boundingpoly2',verbose=False)
Note: See TracBrowser for help on using the repository browser.