source: production/onslow_2006/project.py @ 3031

Last change on this file since 3031 was 2958, checked in by sexton, 19 years ago

new gauges for Onslow

File size: 4.8 KB
Line 
1"""Common filenames and locations for topographic data, meshes and outputs.
2Also includes origin for slump scenario.
3"""
4
5from os import sep, environ, getenv, getcwd
6from os.path import expanduser
7from utilities.polygon import read_polygon
8import sys
9
10from pmesh.create_mesh import convert_points_from_latlon_to_utm
11
12from coordinate_transforms.redfearn import degminsec2decimal_degrees
13
14from time import localtime, strftime
15
16               
17#Making assumptions about the location of scenario data
18scenario_dir_name = 'onslow_tsunami_scenario_2006'
19
20# 250m data to be provided
21coarsename = 'onsl_bathydem250' # get from Neil/Ingo (DEM or topo data)
22
23# 30m data to be provided
24onshore_name = 'onslow_onshore_30m_dted' # get from Neil/Ingo (DEM or topo data)
25
26offshore_name = 'onslow_offshore_points'
27
28boundary_basename = 'SU-AU_clip'
29
30#swollen/ all data output
31basename = 'source'
32
33codename = 'project.py'
34
35if sys.platform == 'win32':
36    home = environ['INUNDATIONHOME']     #Sandpit's parent dir
37#    home = environ['INUNDATIONHOME']     #Sandpit's parent dir
38else:   
39    home = getenv('INUNDATIONHOME', sep+'d'+sep+'cit'+sep+'1'+sep+'cit'+sep+'risk_assessment_methods_project'+sep+'inundation')     
40
41#Derive subdirectories and filenames
42time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
43outputtimedir = home+sep+scenario_dir_name+sep+'output'+sep+time+sep
44meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
45datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
46gaugedir = home+sep+scenario_dir_name+sep+'gauges'+sep
47polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep
48boundarydir = home+sep+scenario_dir_name+sep+'boundaries'+sep
49#output dir without time
50outputdir = home+sep+scenario_dir_name+sep+'output'+sep
51tidedir = home+sep+scenario_dir_name+sep+'tide_data'+sep
52
53print'bound', boundarydir
54
55#gauge_filename = gaugedir + 'onslow_gauges.xya'
56#for MOST
57#gauge_filename = gaugedir + 'onslow_gauges1.xya'
58gauge_filename = gaugedir + 'gauge_location_onslow.csv'
59buildings_filename = gaugedir + 'onslow_res.csv'
60community_filename = gaugedir + 'CHINS_v2.csv'
61community_scenario = gaugedir + 'community_onslow.csv'
62gaugetimeseries = gaugedir + 'onslow'
63
64tidal_filename = tidedir + 'onsl.txt'
65tidal_outname = tidedir + 'max_min.txt'
66
67# boundary source data
68#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
69
70codedir = getcwd()+sep
71                               
72codedirname = codedir + 'project.py'
73
74meshname = meshdir + basename
75
76coarsedemname = datadir + coarsename
77
78onshore_dem_name = datadir + onshore_name
79
80offshore_dem_name = datadir + offshore_name
81
82combined_dem_name = datadir + 'onslow_combined_elevation'
83
84outputname = outputtimedir + basename  #Used by post processing
85
86#!gauge_filename = outputdir + 'onslow_gauges.xya'
87#!gauge_outname = outputdir + 'gauges_max_output.xya'
88
89# clipping region to make DEM (pts file) from fine elevation data
90eastingmin = 240000
91eastingmax = 340000
92northingmin = 7580000
93northingmax = 7700000
94
95south = degminsec2decimal_degrees(-22,15,0)
96north = degminsec2decimal_degrees(-20,30,0)
97west = degminsec2decimal_degrees(114,15,0)
98east = degminsec2decimal_degrees(115,50,0)
99'''
100# region for visualisation
101eminviz = 260000
102emaxviz = 320000
103nminviz = 7590000
104nmaxviz = 7630000
105'''
106# region to export
107
108e_min_area = 300000
109e_max_area = 310000
110n_min_area = 7600000
111n_max_area = 7610000
112
113#Georeferencing
114from coordinate_transforms.redfearn import degminsec2decimal_degrees
115
116refzone = 50
117
118#Updated Main Domain of Onslow: first run NB 6/4/06
119d0 = [310000, 7690000]
120d1 = [280000, 7690000]
121d2 = [270000, 7645000]
122d3 = [240000, 7625000]
123d4 = [270000, 7580000]
124d5 = [300000, 7590000]
125d6 = [340000, 7610000]
126
127polyAll = [d0, d1, d2, d3, d4, d5, d6]
128
129#Interior region - Onslow town
130
131#i0 = [304000, 7608000]
132i0 = [304000, 7607000]
133i1 = [302000, 7605000]
134#i2 = [303000, 7602000]
135i2 = [304000, 7603000]
136#i3 = [305000, 7601000]
137i3 = [307000, 7602000]
138i4 = [309000, 7603000]
139#i5 = [307000, 7606500]
140i5 = [307000, 7606000]
141
142poly_onslow = [i0, i1, i2, i3, i4, i5]
143
144#Thevenard Island
145j0 = [294000, 7629000]
146j1 = [285000, 7625000]
147j2 = [294000, 7621000]
148j3 = [299000, 7625000]
149
150poly_thevenard = [j0, j1, j2, j3]
151
152'''
153# Direction Is
154k0 = [309000, 7619000]
155k1 = [304000, 7619000]
156k2 = [304000, 7616500]
157k3 = [309000, 7616500]
158
159poly_direction = [k0, k1, k2, k3]
160'''
161
162#med res around onslow
163l0 = [300000, 7610000]
164l1 = [285000, 7600000]
165l2 = [300000, 7597500]
166l3 = [310000, 7600000]
167l4 = [315000, 7610000]
168#l5 = [310000, 7612500]
169
170#poly_coast = [l0, l1, l2, l3, l4, l5]
171poly_coast = [l0, l1, l2, l3, l4]
172
173#general coast and local area to onslow region
174m0 = [270000, 7581000]
175m1 = [300000, 7591000]
176m2 = [339000, 7610000]
177m3 = [330000, 7630000]
178m4 = [290000, 7640000]
179m5 = [260000, 7600000]
180
181poly_region = [m0, m1, m2, m3, m4, m5]
Note: See TracBrowser for help on using the repository browser.