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

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

updates

File size: 7.9 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.csv'
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
173##eastingmin25  = 522000
174##eastingmax25  = 555000
175##northingmin25 = 5234000
176##northingmax25 = 5260000
177
178# clipping 25m data set - Hobart
179eastingmin25  = 524208.387
180eastingmax25  = 554867.24
181northingmin25 = 5229154.555
182northingmax25 = 5258511.857
183
184
185# clipping 25m data set - Bruny
186eastingmin25_2 = 523672.502
187eastingmax25_2 = 536020.057
188northingmin25_2 = 5203120.035
189northingmax25_2 = 5212052.309
190
191
192
193###############################
194# Interior region definitions
195###############################
196
197###Interior region - Hobart city area + Glenorchy, Kingston
198##i0 = [517000, 5267000]
199##i1 = [517000, 5255000]
200##i2 = [520000, 5250000]
201##i3 = [522000, 5239000]
202##i4 = [524000, 5238000]
203##i5 = [526000, 5236000]
204##i6 = [530000, 5244000]
205##i7 = [530000, 5250000]
206##i8 = [534000, 5254000]
207##i9 = [520000, 5270000]
208##
209##poly_hobart = [i0, i1, i2, i3, i4, i5, i6, i7, i8, i9]
210##
211### Tasman Peninsula
212##l0 = [550000, 5247000]
213##l1 = [550000, 5211000]
214##l2 = [583000, 5211000]
215##l3 = [583000, 5247000]
216##
217##poly_tasman_peninsula = [l0, l1, l2, l3]
218##
219### Bruny Island
220##poly_bruny = read_polygon(polygondir+'bruny.csv')
221
222
223# Hobart digitized polygons
224poly_hobart1 = read_polygon(polygondir+'Hob_poly1.csv')
225poly_hobart2 = read_polygon(polygondir+'Hob_poly2.csv')
226poly_hobart3 = read_polygon(polygondir+'Hob_poly3.csv')
227poly_hobart4 = read_polygon(polygondir+'Hob_poly4.csv')
228
229plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3, poly_hobart4],'boundingpoly2',verbose=False)
Note: See TracBrowser for help on using the repository browser.