| 1 | """Script for running tsunami inundation scenario for Dampier, WA, Australia. |
|---|
| 2 | |
|---|
| 3 | Source data such as elevation and boundary data is assumed to be available in |
|---|
| 4 | directories specified by project.py |
|---|
| 5 | The output sww file is stored in project.output_run_time_dir |
|---|
| 6 | |
|---|
| 7 | The scenario is defined by a triangular mesh created from project.polygon, |
|---|
| 8 | the elevation data and a simulated tsunami generated with URS code. |
|---|
| 9 | |
|---|
| 10 | Ole Nielsen and Duncan Gray, GA - 2005 and Jane Sexton, Nick Bartzis, GA - 2006 |
|---|
| 11 | """ |
|---|
| 12 | |
|---|
| 13 | #------------------------------------------------------------------------------ |
|---|
| 14 | # Import necessary modules |
|---|
| 15 | #------------------------------------------------------------------------------ |
|---|
| 16 | |
|---|
| 17 | # Standard modules |
|---|
| 18 | |
|---|
| 19 | # # Related major packages |
|---|
| 20 | # from anuga.shallow_water import Domain |
|---|
| 21 | from anuga.shallow_water.data_manager import URS_points_needed_to_file, \ |
|---|
| 22 | urs_ungridded2sww |
|---|
| 23 | import project_urs # Definition of file names and polygons |
|---|
| 24 | |
|---|
| 25 | #------------------------------------------------------------------------------ |
|---|
| 26 | # Copy scripts to time stamped output directory and capture screen |
|---|
| 27 | # output to file |
|---|
| 28 | #------------------------------------------------------------------------------ |
|---|
| 29 | poly = project_urs.poly_all |
|---|
| 30 | print "poly", poly |
|---|
| 31 | print "project_urs.boundaries_in_dir_name",project_urs.boundaries_in_dir_name |
|---|
| 32 | URS_points_needed_to_file(project_urs.boundaries_in_dir_name,poly, export_csv=True, zone=51) |
|---|
| 33 | #URS_points_needed_to_file(project.boundaries_dir_name5,poly, export_csv=True, zone=50) |
|---|
| 34 | #URS_points_needed_to_file(project.boundaries_dir_name5,poly, export_csv=True) |
|---|
| 35 | #urs_ungridded2sww(project.boundaries_in_dir_name2, project.boundaries_dir_name5, mean_stage= project.tide, verbose=True) |
|---|
