| 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 | # run_boundary_points |
|---|
| 27 | #----------------------------------------------------------------------------- |
|---|
| 28 | |
|---|
| 29 | LL_LAT = -50.0 #ll_lat - lower left latitude, in decimal degrees |
|---|
| 30 | LL_LONG = 80.0 #ll-long - lower left longitude, in decimal degrees |
|---|
| 31 | GRID_SPACING = 2.0/60.0 #grid_spacing - in deciamal degrees |
|---|
| 32 | LAT_AMOUNT = 4800 #lat_amount - number of latitudes |
|---|
| 33 | LONG_AMOUNT = 3600 #long_amount- number of longs |
|---|
| 34 | ZONE = 51 |
|---|
| 35 | |
|---|
| 36 | #------------------------------------------------------------------------------ |
|---|
| 37 | # Copy scripts to time stamped output directory and capture screen |
|---|
| 38 | # output to file |
|---|
| 39 | #------------------------------------------------------------------------------ |
|---|
| 40 | poly = project_urs.poly_all |
|---|
| 41 | print "poly", poly |
|---|
| 42 | print "project_urs.boundaries_in_dir_name",project_urs.boundaries_in_dir_name |
|---|
| 43 | URS_points_needed_to_file(project_urs.boundaries_in_dir_name,poly, ZONE, |
|---|
| 44 | LL_LAT, LL_LONG,GRID_SPACING,LAT_AMOUNT, LONG_AMOUNT,export_csv=True) |
|---|
| 45 | #URS_points_needed_to_file(project.boundaries_dir_name5,poly, export_csv=True, zone=50) |
|---|
| 46 | #URS_points_needed_to_file(project.boundaries_dir_name5,poly, export_csv=True) |
|---|
| 47 | #urs_ungridded2sww(project.boundaries_in_dir_name2, project.boundaries_dir_name5, mean_stage= project.tide, verbose=True) |
|---|
