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

Last change on this file since 3668 was 3668, checked in by duncan, 18 years ago

path fix - also, no extent and different grid size

File size: 6.8 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
9               
10#Making assumptions about the location of scenario data
11state = 'tasmania'
12scenario_dir_name = 'hobart_tsunami_scenario_2006'
13
14# data provided by Tas SES and checked by NM&I
15#onshore_name = 'hob3_topo' # original
16onshore_name = 'hob5_topo' # 12.5m grid and clipped to 100m elevation or 3000m from coast
17#onshore_name_25 = 'hob5_topo_25m' # 25m grid and clipped to 100m elevation or 3000m from coast
18onshore_name_25 = 'hob6_topo_25m' # 25m grid NOT clipped
19offshore_name_tas1 = 'derwent_2m'
20offshore_name_tas2 = 'derwent_5m'
21offshore_name_tas3 = 'south_east_tas' #actually this is AHO
22offshore_name_tas4 = 'hobart_1m'
23
24# AHO data and checked by NM&I
25offshore_name1 = 'xy100003760'
26offshore_name2 = 'xy100003761'
27offshore_name3 = 'xy100003762'
28offshore_name4 = 'xy100003907'
29offshore_name5 = 'xy100003908'
30offshore_name6 = 'xy100003909'
31offshore_name7 = 'xy100003910'
32offshore_name8 = 'xy100003932'
33offshore_name9 = 'xy100003933'
34offshore_name10 = 'xy100003934'
35offshore_name11 = 'xy100003935'
36offshore_name12 = 'xy100003936'
37offshore_name13 = 'xy100003964'
38offshore_name14 = 'xy100014250'
39offshore_name15 = 'xy100014253'
40offshore_name16 = 'xy100016142'
41
42# developed by NM&I
43coast_name = 'coastline_points'
44
45boundary_basename = 'SU-AU_clip'
46
47#swollen/ all data output
48basename = 'source'
49
50codename = 'project.py'
51
52if sys.platform == 'win32':
53    home = getenv('INUNDATIONHOME')
54    user = getenv('USERPROFILE')
55
56else:   
57    home = getenv('INUNDATIONHOME', sep+'d'+sep+'xrd'+sep+'gem'+sep+'2'+sep+'ramp'+sep+'risk_assessment_methods_project'+sep+'inundation')     
58    user = getenv('LOGNAME')
59    print 'USER:', user
60
61# INUNDATIONHOME is the inundation directory, not the data directory.
62home += sep +'data'
63
64#Derive subdirectories and filenames
65#time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
66local_time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
67
68meshdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'meshes'+sep
69datadir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'topographies'+sep
70gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep
71polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
72boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
73#output dir without time
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'
81
82gaugetimeseries = gaugedir + 'hobart'
83
84# boundary source data
85#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
86
87codedir = getcwd()+sep
88                               
89codedirname = codedir + 'project.py'
90
91meshname = outputtimedir + 'mesh_' + basename
92
93onshore_dem_name = datadir + onshore_name
94onshore_dem_name_25 = datadir + onshore_name_25
95all_onshore_dem_name = datadir + 'combined_onshore'
96offshore_dem_name_local1 = datadir + offshore_name_tas1
97offshore_dem_name_local2 = datadir + offshore_name_tas2
98offshore_dem_name_local3 = datadir + offshore_name_tas3
99offshore_dem_name_local4 = datadir + offshore_name_tas4
100offshore_dem_name_aho1 = datadir + offshore_name1
101offshore_dem_name_aho2 = datadir + offshore_name2
102offshore_dem_name_aho3 = datadir + offshore_name3
103offshore_dem_name_aho4 = datadir + offshore_name4
104offshore_dem_name_aho5 = datadir + offshore_name5
105offshore_dem_name_aho6 = datadir + offshore_name6
106offshore_dem_name_aho7 = datadir + offshore_name7
107offshore_dem_name_aho8 = datadir + offshore_name8
108offshore_dem_name_aho9 = datadir + offshore_name9
109offshore_dem_name_aho10 = datadir + offshore_name10
110offshore_dem_name_aho11 = datadir + offshore_name11
111offshore_dem_name_aho12 = datadir + offshore_name12
112offshore_dem_name_aho13 = datadir + offshore_name13
113offshore_dem_name_aho14 = datadir + offshore_name14
114offshore_dem_name_aho15 = datadir + offshore_name15
115offshore_dem_name_aho16 = datadir + offshore_name16
116coast_dem_name = datadir + coast_name
117combined_dem_name = datadir + 'hobart_combined_elevation'
118
119outputname = outputtimedir + basename  #Used by post processing
120
121# clipping 12.5m onshore data set
122eastingmin = 520000
123eastingmax = 536000
124northingmin = 5245000
125northingmax = 5265000
126# clipping 25m onshore data set
127eastingmin25 = 520000
128eastingmax25 = 580000
129northingmin25 = 5170000
130northingmax25 = 5260000
131
132# region to export for inundation map
133e_min_area = 521000#490000
134e_max_area = 522000#580000
135n_min_area = 5190000#5160000
136n_max_area = 5200000#5275000
137
138from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
139                       
140# bounding box
141south = degminsec2decimal_degrees(-43,45,0)
142north = degminsec2decimal_degrees(-42,30,0)
143west = degminsec2decimal_degrees(146,45,0)
144east = degminsec2decimal_degrees(148,0,0)#degminsec2decimal_degrees(148,15,0)
145
146#Main Domain of Hobart: first run JS 18/9/06
147d0 = [south, west]
148d1 = [south, east]
149d2 = [north, east]
150d3 = [north, west]
151
152polyAll, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
153refzone = zone
154#polyAll = [[500000, 5200000],[580000, 5200000],[580000, 5270000],[500000,5270000]] #073937
155#polyAll = [[500000, 5160000],[580000, 5160000],[580000, 5270000],[500000,5270000]] # try this is the above works
156#polyAll = [[510000, 5170000],[570000, 5170000],[570000, 5250000],[510000,5250000]] # try this is the above works
157#polyAll = [[520000, 5160000],[580000, 5170000],[580000, 5200000],[600000,5260000],[510000,5280000]] # try this is the above works
158polyAll = [[520000, 5170000],[580000, 5170000],[580000, 5200000],[590000,5240000],[520000,5260000]] # try this is the above works
159# bottom bright, topr, top, left
160
161# region to export for Alex to make bathymetry map
162e_min_area = 560000#500000 #490000
163e_max_area = 570000#560000 #580000
164n_min_area = 5200000#5240000#5160000
165n_max_area = 5230000#5260000#5270000
166
167#Interior region - Hobart city area + Glenorchy, Kingston
168i0 = [517000, 5267000]
169i1 = [517000, 5255000]
170i2 = [520000, 5250000]
171i3 = [522000, 5239000]
172i4 = [524000, 5238000]
173i5 = [526000, 5236000]
174i6 = [530000, 5244000]
175i7 = [530000, 5250000]
176i8 = [534000, 5254000]
177i9 = [520000, 5270000]
178
179poly_hobart = [i0, i1, i2, i3, i4, i5, i6, i7, i8, i9]
180
181# Tasman Peninsula
182l0 = [550000, 5247000]
183l1 = [550000, 5211000]
184l2 = [583000, 5211000]
185l3 = [583000, 5247000]
186
187poly_tasman_peninsula = [l0, l1, l2, l3]
188
189# Bruny Island
190from anuga.utilities.polygon import read_polygon
191poly_bruny = read_polygon(polygondir+'bruny.csv')
Note: See TracBrowser for help on using the repository browser.