source: anuga_work/development/Hinwood_2008/scenarios.py @ 6823

Last change on this file since 6823 was 5710, checked in by duncan, 16 years ago

comments

File size: 9.1 KB
Line 
1"""
2Scenario information, used to run simulations and create graphs.
3"""
4
5# A list of scenario dictionaries
6scenarios = []
7
8# 0 - t1r3
9# 1 - t1r5
10# 2 - t2r7
11# 3 - t2r8
12# 4 - t3r28
13# 5 - t3r29
14# 6 - t4r31
15# 7 - t4r32
16
17# T1R3
18data = {'xleft':[-3.106,0.0],  # Av' of ADV and Gauge A
19        'xtoe':[0.0,0.0],
20        'xbeach':[1.285,0.090],
21        'xright':[8,.4843334],
22        'offshore_water_depth':.4,
23        'scenario_id':'T1R3',
24        'gauge_names':['B','1','2','3','4','5','7'],
25        'gauge_x':[-0.68, 1.569, 2.568, 3.566, 4.564, 5.562, 7.559],
26        'gauge_bed_elevation':[-0.400000, -0.293158,
27                               -0.234473, -0.175788, -0.117104,
28                               -0.058419, 0.058950],
29        'start_slope_x':4,
30        'finish_slope_x':6,
31        'break_times':[34.15, 38.35, 42.95, 47.85, 53.15, 58.15],
32        'break_xs':[5.45, 5.45, 5.28, 5.28, 5.29, 5.29],
33        'break_type':['none', 'plunge', 'plunge', 'plunge', 'plunge', \
34                      'plunge'],
35        'axis':[0,80,-0.040,0.040],
36        'axis_maximum_x':6.0,
37        'ANUGA_start_time':1.83,
38        'band_offset':-1.0,
39        'wave_times':[23.0,59.0]  # this is in Anuga time
40                 }
41scenarios.append(data)
42
43# T1R5
44data = {'xleft':[-3.106,0.0],  # Av' of ADV and Gauge A
45        'xtoe':[0.0,0.0],
46        'xbeach':[1.285,0.090],
47        'xright':[8,.4843334],
48        'offshore_water_depth':.4,
49        'scenario_id':'T1R5',
50        'gauge_names':['B','1','2','3','4','5','7'],
51        'gauge_x':[-0.68, 1.569, 2.568, 3.566, 4.564, 5.562, 7.559],
52        'gauge_bed_elevation':[-0.400000, -0.293158,
53                               -0.234473, -0.175788, -0.117104,
54                               -0.058419, 0.058950],
55        'start_slope_x':4,
56        'finish_slope_x':6,
57        'break_times':[28.68, 33.18, 37.68, 42.18, 47.38, 52.18],
58        'break_xs':[5.45, 5.58, 5.29, 5.245, 5.21, 5.21],
59        'break_type':['none', 'plunge', 'plunge', 'plunge', 'plunge', \
60                      'plunge'],
61        'axis':[0,80,-0.04,0.04],
62        'axis_maximum_x':6.0,
63        'ANUGA_start_time':0.38,
64        'band_offset':-0.5,
65        'wave_times':[20.0,59.0] # this is in Anuga time
66        }
67scenarios.append(data)
68
69
70
71# #T2R7
72# xleft is different
73data = {'xleft':[-4.586,0.0],  # Av' of ADV and Gauge A
74        'xtoe':[0.0,0.0],
75        'xbeach':[1.285,0.090],
76        'xright':[8,.4843334],
77        'offshore_water_depth':.4,
78        'scenario_id':'T2R7',
79        'gauge_names':['B','1','2','3','4','5'],
80        'gauge_x':[-0.68, 1.569, 2.568, 3.566, 4.564, 5.562],
81        'gauge_bed_elevation':[-0.400000, -0.293158,
82                               -0.234473, -0.175788, -0.117104,
83                               -0.058419],
84        'start_slope_x':5,
85        'finish_slope_x':6.5,
86        'break_times':[33.95, 41.25, 49.55, 56.75], # last 2 removed
87        'break_xs':[5.93,5.93,5.907,5.78], # since depth finishes at 65 sec
88        'break_type':['none', 'front steepened', 'front steepened', 'spill'],
89        'axis':[0,80,-0.040,0.040],
90        'axis_maximum_x':6.0,
91        'ANUGA_start_time':0.05,
92        'band_offset':2.0,
93        'wave_times':[15.0,64.0] # this is in Anuga time
94                 }
95scenarios.append(data)
96
97# #T2R8
98# xleft is different
99data = {'xleft':[-4.586,0.0],  # Av' of ADV and Gauge A
100        'xtoe':[0.0,0.0],
101        'xbeach':[1.285,0.090],
102        'xright':[8,.4843334],
103        'offshore_water_depth':.4,
104        'scenario_id':'T2R8',
105        'gauge_names':['B','1','2','3','4','5'],
106        'gauge_x':[-0.68, 1.569, 2.568, 3.566, 4.564, 5.562],
107        'gauge_bed_elevation':[-0.400000, -0.293158,
108                               -0.234473, -0.175788, -0.117104,
109                               -0.058419],
110        'start_slope_x':5,
111        'finish_slope_x':6.5,
112        'break_times':[47.00,55.00,62.00,70.00], # last 2 removed
113        'break_xs':[6.32,6.32,6.033,5.935],  # since depth finishes at 75 sec
114        'break_type':['none', 'front steepened', 'front steepened',
115                      'weak break'],
116        'axis':[0,80,-0.040,0.040],
117        'axis_maximum_x':6.0,
118        'ANUGA_start_time':0.04,
119        'band_offset':-2.0,
120        'wave_times':[34.0,74.0] # this is in Anuga time
121        }
122scenarios.append(data)
123
124# #T3R28
125# xleft is different
126data = {'xleft':[-3.875,0.0],  # Av' of ADV and Gauge A
127        'xtoe':[0.0,0.0],
128        'xbeach':[1.285,0.090],
129        'xright':[14.,.3903881],
130        'offshore_water_depth':.336,
131        'scenario_id':'T3R28',
132        'gauge_names':['B','1','2','3','5','7',
133                       '9','10','11','12'],
134
135        'gauge_x':[-0.325, 1.572, 2.571, 3.571, 5.57,
136                   7.57, 9.569, 10.569, 11.569,
137                   12.569],
138
139        # remove
140        'gauge_bed_elevation':[-0.336000, -0.237263, -0.213789, -0.190315,
141                               -0.143368,
142                               -0.096420, -0.049472,
143                               -0.025998, -0.002524, 0.020949],
144        'start_slope_x':8.,
145        'finish_slope_x':9.5,
146        'break_times':[61.21,68.51,76.11,84.11],
147        'break_xs':[9.17,9.135,9.135,9.105],
148        'break_type':['collapse', 'collapse', 'collapse',
149                      'collapse'],
150        'axis':[0,80,-0.040,0.040],
151        'axis_maximum_x':12.0,
152        'ANUGA_start_time':12.18,
153        'band_offset':-0.5,
154        'wave_times':[30.0,85.0] # this is in Anuga time
155        }
156scenarios.append(data)
157
158# #T3R29
159data = {'xleft':[-3.875,0.0],  # Av' of ADV and Gauge A
160        'xtoe':[0.0,0.0],
161        'xbeach':[1.285,0.090],
162        'xright':[14.,.3903881],
163        'offshore_water_depth':.336,
164        'scenario_id':'T3R29',
165        'gauge_names':['B','1','2','3','5','7',
166                       '9','10','11','12'],
167
168        'gauge_x':[-0.325, 1.572, 2.571, 3.571, 5.57,
169                   7.57, 9.569, 10.569, 11.569,
170                   12.569],
171
172        # remove
173        'gauge_bed_elevation':[-0.336000, -0.237263, -0.213789, -0.190315,
174                               -0.143368,
175                               -0.096420, -0.049472,
176                               -0.025998, -0.002524, 0.020949],
177           
178        'start_slope_x':8.,
179        'finish_slope_x':9.5,
180        'break_times':[61.61,66.61,74.31,81.71],
181        'break_xs':[9.063,9.063,9.043,9.043],
182        'break_type':['none','collapse', 'collapse', 'collapse'],
183        'axis':[0,80,-0.040,0.040],
184        'axis_maximum_x':12.0,
185        'ANUGA_start_time':10.48,
186        'band_offset':-2.0,
187        'wave_times':[30.0,74.0] # this is in Anuga time
188                 }
189scenarios.append(data)
190
191# #T4R31
192# xleft is different
193data = {'xleft':[-2.43,0.0],  # Av' of ADV and Gauge A
194        'xtoe':[0.0,0.0],
195        'xbeach':[1.285,0.090],
196        'xright':[14.,.3903881],
197        'offshore_water_depth':.336,
198        'scenario_id':'T4R31',     
199        'gauge_names':['B','1','2','3','5','7',
200                       '9','10','11','12'],
201
202        'gauge_x':[-0.325, 1.572, 2.571, 3.571, 5.57,
203                   7.57, 9.569, 10.569, 11.569,
204                   12.569],
205
206        # remove
207        'gauge_bed_elevation':[-0.336000, -0.237263, -0.213789, -0.190315,
208                               -0.143368,
209                               -0.096420, -0.049472,
210                               -0.025998, -0.002524, 0.020949],
211       
212        'start_slope_x':6,
213        'finish_slope_x':9,
214        'break_times':[65.99,69.89,75.49,81.49,86.19],
215        'break_xs':[8.53,7.46,7.492,7.492,7.444],
216        'break_type':['collapse', 'spill', 'spill', 'spill', 'spill'],
217        'axis':[0,80,-0.040,0.040],
218        'axis_maximum_x':12.0,
219        'ANUGA_start_time':11.63,
220        'band_offset':-0.5,
221        'wave_times':[34.0,75.0] # this is in Anuga time
222        }
223scenarios.append(data)
224
225# #T4R32
226data = {'xleft':[-2.43,0.0],  # Av' of ADV and Gauge A
227        'xtoe':[0.0,0.0],
228        'xbeach':[1.285,0.090],
229        'xright':[14.,.3903881],
230        'offshore_water_depth':.336,
231        'scenario_id':'T4R32',   
232        'gauge_names':['B','1','2','3','5','7',
233                       '9','10','11','12'],
234
235        'gauge_x':[-0.325, 1.572, 2.571, 3.571, 5.57,
236                   7.57, 9.569, 10.569, 11.569,
237                   12.569],
238
239        # remove
240        'gauge_bed_elevation':[-0.336000, -0.237263, -0.213789, -0.190315,
241                               -0.143368,
242                               -0.096420, -0.049472,
243                               -0.025998, -0.002524, 0.020949],
244       
245        'start_slope_x':6,
246        'finish_slope_x':9,
247        # The ANUGA_start_time has not been subtracted from these values
248        'break_times':[64.29,68.69,73.49,79.39,84.69], # last wave removed
249        'break_xs':[7.673,7.673,7.43,7.448,7.448], # Pressure info ends early
250        'break_type':['none','spill', 'spill', 'spill', 'spill'],
251        'axis':[0,80,-0.040,0.040],
252        'axis_maximum_x':12.0,
253        'ANUGA_start_time':12.68,
254        'band_offset':-1.0,
255        'wave_times':[34.0,75.0] # this is in Anuga time
256        }
257scenarios.append(data)
258
259# These are the indexes from 'gauge_x' of the gauges for the journal article.
260to_publish_indexes = {'T1R5':[0, 5, 6],
261              'T2R7':[0, 3, 5],
262              'T3R29':[0, 6, 8],
263              'T4R32':[0, 5, 8]
264              }
265
266
267
Note: See TracBrowser for help on using the repository browser.