Changeset 5413
- Timestamp:
- Jun 24, 2008, 8:53:59 AM (17 years ago)
- Location:
- anuga_work/development/Hinwood_2008
- Files:
-
- 1 added
- 2 edited
-
plot.py (modified) (2 diffs)
-
run_dam.py (modified) (1 diff)
-
slope.py (added)
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/Hinwood_2008/plot.py
r5410 r5413 82 82 """ 83 83 from scenarios import scenarios 84 outputdir_tag = "_good_tri_area_0.0 01_A"84 outputdir_tag = "_good_tri_area_0.01_A" 85 85 #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!! 86 86 # Supported formats: emf, eps, pdf, png, ps, raw, rgba, svg, svgz. … … 91 91 file_exp, file_sim, location_sims, location_exps, outputdir_name, \ 92 92 save_as_list = temp 93 print "run_data['scenario_id']", run_data['scenario_id'] 94 #location_sims = [location_sims[0]] 95 #location_exps = [location_exps[0]] 96 #save_as_list = [save_as_list[0]] 93 97 for loc_sim, loc_exp, save_as, gauge in map(None, location_sims, 94 98 location_exps, -
anuga_work/development/Hinwood_2008/run_dam.py
r5410 r5413 248 248 #------------------------------------------------------------- 249 249 if __name__ == "__main__": 250 #slopes = [[-4.5,0.0],[0.0,0.0],[1.285,0.090],[16.1,.960]] 251 # Note, gauge A has been removed, since it is used as the 252 # boundary. 253 run_type = 6 254 255 # T1R5 256 run_data = {'xleft':[-3.106,0.0], # Av' of ADV and Gauge A 257 'xtoe':[0.0,0.0], 258 'xbeach':[1.285,0.090], 259 'xright':[16.1,.960], 260 'offshore_water_depth':.4, 261 'scenario_id':'T1R5', 262 'gauge_names':['B','1','2','3','4','5','6','7','8', 263 '9','10','11','12','13','14'], 264 'gauge_x':[-0.68, 1.572, 2.572, 3.572, 4.572, 5.572, 265 6.572, 7.572, 8.572, 9.572, 10.572, 11.572, 266 12.572, 13.572, 14.572], 267 'gauge_bed_elevation':[-0.400000, -0.293158, 268 -0.234473, -0.175788, -0.117104, -0.058419, 0.000266, 269 0.058950, 0.117635, 0.176320, 0.235004, 0.293689, 270 0.352374, 0.411058, 0.469743] 271 } 272 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 273 run_type = run_type, 274 outputdir_name=run_data['scenario_id']) 275 276 # T1R3 277 run_data = {'xleft':[-3.106,0.0], # Av' of ADV and Gauge A 278 'xtoe':[0.0,0.0], 279 'xbeach':[1.285,0.090], 280 'xright':[16.1,.960], 281 'offshore_water_depth':.4, 282 'scenario_id':'T1R3', 283 'gauge_names':['B','1','2','3','4','5','6','7','8', 284 '9','10','11','12','13','14'], 285 'gauge_x':[-0.68, 1.572, 2.572, 3.572, 4.572, 5.572, 286 6.572, 7.572, 8.572, 9.572, 10.572, 11.572, 287 12.572, 13.572, 14.572], 288 'gauge_bed_elevation':[-0.400000, -0.293158, 289 -0.234473, -0.175788, -0.117104, -0.058419, 0.000266, 290 0.058950, 0.117635, 0.176320, 0.235004, 0.293689, 291 0.352374, 0.411058, 0.469743] 292 } 293 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 294 run_type = run_type, 295 outputdir_name=run_data['scenario_id']) 296 297 # #T2R7 298 # xleft is different 299 run_data = {'xleft':[-4.586,0.0], # Av' of ADV and Gauge A 300 'xtoe':[0.0,0.0], 301 'xbeach':[1.285,0.090], 302 'xright':[16.1,.960], 303 'offshore_water_depth':.4, 304 'scenario_id':'T2R7', 305 'gauge_names':['B','1','2','3','4','5','6','7','8', 306 '9','10','11','12','13','14'], 307 'gauge_x':[-0.68, 1.572, 2.572, 3.572, 4.572, 5.572, 308 6.572, 7.572, 8.572, 9.572, 10.572, 11.572, 309 12.572, 13.572, 14.572], 310 'gauge_bed_elevation':[-0.400000, -0.293158, 311 -0.234473, -0.175788, -0.117104, -0.058419, 0.000266, 312 0.058950, 0.117635, 0.176320, 0.235004, 0.293689, 313 0.352374, 0.411058, 0.469743] 314 } 315 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 316 run_type = run_type, 317 outputdir_name=run_data['scenario_id']) 318 319 # #T2R8 320 # xleft is different 321 run_data = {'xleft':[-4.586,0.0], # Av' of ADV and Gauge A 322 'xtoe':[0.0,0.0], 323 'xbeach':[1.285,0.090], 324 'xright':[16.1,.960], 325 'offshore_water_depth':.4, 326 'scenario_id':'T2R8', 327 'gauge_names':['B','1','2','3','4','5','6','7','8', 328 '9','10','11','12','13','14'], 329 'gauge_x':[-0.68, 1.572, 2.572, 3.572, 4.572, 5.572, 330 6.572, 7.572, 8.572, 9.572, 10.572, 11.572, 331 12.572, 13.572, 14.572], 332 'gauge_bed_elevation':[-0.400000, -0.293158, 333 -0.234473, -0.175788, -0.117104, -0.058419, 0.000266, 334 0.058950, 0.117635, 0.176320, 0.235004, 0.293689, 335 0.352374, 0.411058, 0.469743] 336 } 337 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 338 run_type = run_type, 339 outputdir_name=run_data['scenario_id']) 340 341 # #T3R29 342 # xleft is different 343 run_data = {'xleft':[-3.875,0.0], # Av' of ADV and Gauge A 344 'xtoe':[0.0,0.0], 345 'xbeach':[1.285,0.090], 346 'xright':[16.1,.440], 347 'offshore_water_depth':.336, 348 'scenario_id':'T3R29', 349 'gauge_names':['1','2','3','4','5','6','7','8', 350 '9','10','11','12','13','14','B'], 351 'gauge_x':[1.572, 2.572, 3.572, 4.572, 5.572, 352 6.572, 7.572, 8.572, 9.572, 10.572, 11.572, 353 12.572, 13.572, 14.572, -0.325], 354 'gauge_bed_elevation':[-0.237263, -0.213789, -0.190315, 355 -0.166841, -0.143368, -0.119894, 356 -0.096420, -0.072946, -0.049472, 357 -0.025998, -0.002524, 0.020949, 358 0.044423, 0.067897, -0.336000] 359 } 360 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 361 run_type = run_type, 362 outputdir_name=run_data['scenario_id']) 363 364 # #T3R28 365 run_data = {'xleft':[-3.875,0.0], # Av' of ADV and Gauge A 366 'xtoe':[0.0,0.0], 367 'xbeach':[1.285,0.090], 368 'xright':[16.1,.440], 369 'offshore_water_depth':.336, 370 'scenario_id':'T3R28', 371 'gauge_names':['1','2','3','4','5','6','7','8', 372 '9','10','11','12','13','14','B'], 373 'gauge_x':[1.572, 2.572, 3.572, 4.572, 5.572, 374 6.572, 7.572, 8.572, 9.572, 10.572, 11.572, 375 12.572, 13.572, 14.572, -0.325], 376 'gauge_bed_elevation':[-0.237263, -0.213789, -0.190315, 377 -0.166841, -0.143368, -0.119894, 378 -0.096420, -0.072946, -0.049472, 379 -0.025998, -0.002524, 0.020949, 380 0.044423, 0.067897, -0.336000] 381 } 382 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 383 run_type = run_type, 384 outputdir_name=run_data['scenario_id']) 385 386 # #T4R31 387 # xleft is different 388 run_data = {'xleft':[-2.43,0.0], # Av' of ADV and Gauge A 389 'xtoe':[0.0,0.0], 390 'xbeach':[1.285,0.090], 391 'xright':[16.1,.440], 392 'offshore_water_depth':.336, 393 'scenario_id':'T4R31', 394 'gauge_names':['1','2','3','4','5','6','7','8', 395 '9','10','11','12','13','14','B'], 396 'gauge_x':[1.572, 2.572, 3.572, 4.572, 5.572, 397 6.572, 7.572, 8.572, 9.572, 10.572, 11.572, 398 12.572, 13.572, 14.572, -0.325], 399 'gauge_bed_elevation':[-0.237263, -0.213789, -0.190315, 400 -0.166841, -0.143368, -0.119894, 401 -0.096420, -0.072946, -0.049472, 402 -0.025998, -0.002524, 0.020949, 403 0.044423, 0.067897, -0.336000] 404 } 405 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 406 run_type = run_type, 407 outputdir_name=run_data['scenario_id']) 408 409 # #T4R32 410 run_data = {'xleft':[-2.43,0.0], # Av' of ADV and Gauge A 411 'xtoe':[0.0,0.0], 412 'xbeach':[1.285,0.090], 413 'xright':[16.1,.440], 414 'offshore_water_depth':.336, 415 'scenario_id':'T4R32', 416 'gauge_names':['1','2','3','4','5','6','7','8', 417 '9','10','11','12','13','14','B'], 418 'gauge_x':[1.572, 2.572, 3.572, 4.572, 5.572, 419 6.572, 7.572, 8.572, 9.572, 10.572, 11.572, 420 12.572, 13.572, 14.572, -0.325], 421 'gauge_bed_elevation':[-0.237263, -0.213789, -0.190315, 422 -0.166841, -0.143368, -0.119894, 423 -0.096420, -0.072946, -0.049472, 424 -0.025998, -0.002524, 0.020949, 425 0.044423, 0.067897, -0.336000] 426 } 427 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 428 run_type = run_type, 429 outputdir_name=run_data['scenario_id']) 250 251 from scenarios import scenarios 252 253 run_type = 5 254 for run_data in [scenarios[0]]: 255 main( run_data['scenario_id'] + '_boundary.tsm' , run_data, 256 run_type = run_type, 257 outputdir_name=run_data['scenario_id']) 258
Note: See TracChangeset
for help on using the changeset viewer.
