Index: /DVD_images/extra_files/Hobart/analysis.html =================================================================== --- /DVD_images/extra_files/Hobart/analysis.html (revision 7330) +++ /DVD_images/extra_files/Hobart/analysis.html (revision 7331) @@ -16,5 +16,5 @@
- This page describes the scripts and shows how to use them. Note that due to the ad-hoc nature of the scripts,
+ This page describes the scripts and shows how to use them. Note that due to the ad-hoc nature of the scripts,
you must modify these scripts before they can be used to extract information from any SWW files you
generate.
@@ -33,47 +33,52 @@
export_results_max.py
The export_results_max.py script can be used to create an ASC file that contains the maximum
- of an ANUGA variable or expression for each point in a user-defined raster on the simulation region.
+ of an ANUGA variable or expression for each point in a user-defined ArcGIS grid on the simulation region.
You must change certain values within the export_results_max.py file to get what you want:
time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt' time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'- as necessary. Then modify the line: + For information on output_folder_name. +
Then modify the line:
time_dirs = [time_dir1, time_dir2]to include each of the variables you defined above. -
-
+
-
times = [37860]+ We do this by specifying the start times of all the extra SWW files:
- Note we do not have to specify the first SWW file, only the extra file start times. If there are no extra SWW files, - just do: -
+
times = [16980, 33960, 50940]+
+ Note we do not have to specify the first SWW file, only the extra file start times. If there are no extra SWW files, + just do: +
times = []
-
cellsize = 20 # raster cell size in metres-
-
cellsize = 20 # ArcGIS grid cell size in metres+
+ Note: That this value should not go beyond the most refined part of the mesh + i.e. if mesh resolution is 500m2, cellsize = (500 x 2)1/2 = 32m +
+
timestep = None # over all timesteps #timestep = 0-
- Note that if you want the elevation use timestep = 0. Elevation does not change over time and you only need the - first SWW file.
-
+
area = ['Hobart', 'NW', 'South']- Note that the name strings you put into the area list must match the names used in the +
+ Note that the name strings you put into the area list must match the names used in the project.py file. That is, if you use a name such as 'NW', then the - project.py file must have lines like: + project.py file must have lines like:
xminNW = 548000 xmaxNW = 561000 @@ -85,9 +90,9 @@Doing the above in a Command Prompt window is temporary - if you open another window you will not have the ANUGADATA environment variable defined. To make @@ -84,17 +85,13 @@ tool (Advanced tab). See the Installation Guide if you don't know how to do this. -area = ['All'] # no region of interest-
- Finally you must decide which variable or expression values you want sampled on your raster. +
- Finally you must decide which variable or expression values you want sampled on your ArcGIS grid. Define a list var that contains strings defining the required variable or expression:
var = ['depth', 'speed']- Note that the strings you supply must be defined in the var_equations dictionary - prior to your definition of var: + Note that the strings you supply must be defined in the var_equations dictionary + prior to your definition of var:var_equations = {'stage': 'stage', 'momentum': '(xmomentum**2 + ymomentum**2)**0.5', @@ -100,9 +105,9 @@+ X:\ is the assumed to be the CD drive letter.get_timeseries.py
The get_timeseries.py script is used to get timeseries data for a selection of variable - data at one or more gauge points in a simulation. + data at one or more timeseries gauge points in a simulation.You must change certain values within the get_timeseries.py file to get what you want:
-
- You must determine how many SWW files your simulation produced and create as many lines like this: +
- Depending on how many simulations you have run will determine the number of output folder names you insert here:
time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt' time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'@@ -111,5 +116,5 @@ to include each of the variables you defined above.-
- Make sure that the gauges list in projects.py contains one or more CSV files +
- Make sure that the gauges list specified in projects.py contains one or more CSV files defining the gauges within the simulation that you want the timeseries data for. The gauge file must have this format: Index: D_images/extra_files/Hobart/auxiliary_scripts.html =================================================================== --- /DVD_images/extra_files/Hobart/auxiliary_scripts.html (revision 7330) +++ (revision ) @@ -1,132 +1,0 @@ - - - -
Geoscience Australia - - - - --
-- -- - -- - -- -Tsunami Inundation Models for the HOBART region
-
- -Auxiliary Scripts
- There are number of auxiliary scripts on this disk that you can use to extract information from the SWW files - produced by a simulation run. These scripts are not meant to be the definitive answer on how to get information - from SWW files, but are offered as an example of one way to do it. -- This page describes the scripts and shows how to use them. Note that due to the ad-hoc nature of the scripts, - you must modify these scripts before they can be used to extract information from any SWW files you - generate. - -
export_results_max.py
- The export_results_max.py script can be used to create an ASC file that contains the maximum - of an ANUGA variable or expression for each point in a user-defined raster on the simulation region. -- You must change certain values within the export_results_max.py file to get what you want: -
-
- -- You must determine how many SWW files your simulation produced and create as many lines like this: -
time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt' - time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt' - time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt' - time_dir4 = '20090505_150805_run_final_0.8_58280_None_kvanputt' - time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt' - time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt'- as necessary. Then modify the line: -time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]- to include each of the variables you defined above. - -- Modify the cellsize value to set the size of the raster you require. For example: -
cellsize = 20 # raster cell size in metres- -- Now set the timestep at which you want the raster generated. Either set the actual timestep required or use - None to indicate that you want the maximum values in the raster over all timesteps: -
timestep = 1000 - #timestep = None # over all timesteps-- Note that if you specify a particular timestep here you may not need to include all SWW files in the - time_dirs list above. For example, if you wanted data from timestep 0 only then you would need - to include only the first SWW file in the list. - -
- If you want to clip the raster to one or more small regions then put the names of the regions of interest - into the area list: -
area = ['Hobart', 'NW', 'South']- Note that the name strings you put into the area list must match the names used in the - project.py file. That is, if you use a name such as 'NW', then the - project.py file must have lines like: -xminNW = 548000 - xmaxNW = 561000 - yminNW = 5250000 - ymaxNW = 5258000- which set the maximum and minimum eastings and northings that define a rectangular region. -- If you don't want to clip to a region of interest, then do this: -
area = ['All'] # no region of interest- -- Now you must decide which variable or expression values you want sampled on your raster. - Define a list var that contains strings defining the required variable/expression: -
var = ['depth', 'speed']- Note that the strings you supply must be defined in the var_equations dictionary - prior to your definition of var: -var_equations = {'stage': 'stage', - 'momentum': '(xmomentum**2 + ymomentum**2)**0.5', - 'depth': 'stage-elevation', - 'speed': '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6)', - 'elevation': 'elevation' }- This dictionary maps your var strings to a variable or expression. -get_runup.py
- The get_runup.py script is used to get the maximum runup height in a specified region - or regions. -- You must change certain values within the get_runup.py file to get what you want: -
-
- -- You must determine how many SWW files your simulation produced and create as many lines like this: -
time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt' - time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt' - time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt' - time_dir4 = '20090505_150805_run_final_0.8_58280_None_kvanputt' - time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt' - time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt'- as necessary. Then modify the line: -time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]- to include each of the variables you defined above. - -- Make sure that the images list in projects.py contains one or more CSV files - defining the regions within the simulation that you want the maximum runup for. -
get_timeseries.py
- The get_timeseries.py script is used to get timeseries data for a selection of variable - data at one or more gauge points in a simulation. -- You must change certain values within the get_timeseries.py file to get what you want: -
-
- - Index: /DVD_images/extra_files/Hobart/index.html =================================================================== --- /DVD_images/extra_files/Hobart/index.html (revision 7330) +++ /DVD_images/extra_files/Hobart/index.html (revision 7331) @@ -16,5 +16,5 @@- You must determine how many SWW files your simulation produced and create as many lines like this: -
time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt' - time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt' - time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt' - time_dir4 = '20090505_150805_run_final_0.8_58280_None_kvanputt' - time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt' - time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt'- as necessary. Then modify the line: -time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]- to include each of the variables you defined above. - -- Make sure that the gauges list in projects.py contains one or more CSV files - defining the gauges within the simulation that you want the timeseries data for. -
- @@ -23,31 +23,29 @@Tsunami Inundation Models for the HOBART region
+Tsunami Inundation Models for south east Tasmania
Introduction
- - The information on this DVD and the associated report is intended to assist emergency managers in developing preparation and - response plans that can be used during a tsunami threat or emergency. -- The Attorney General's Department (AGD) has supported Geoscience Australia (GA) in developing a range of products to support - the understanding of tsunami hazard through the Australian Probabilistic Tsunami Warning System Project. The work reported here is intended - to further build the capacity of the Jurisdictions in developing inundation models for prioritised locations. + The information within this DVD is intended to assist emergency managers in tsunami planning and preparation + activities. The Australian Government Attorney Generals Department (AGD) has supported Geoscience Australia (GA) in developing a + range of products to support the understanding of tsunami hazard through the Australian Tsunami Warning System Project. The work + reported here is intended to further build the capacity of the Tasmanian State Government in developing inundation models for + prioritised locations.
While the associated report provides background, model results and interpretations, the DVD contains all data, scripts and - software necessary to reproduce and potentially augment the models underpinning the report. This will allow the Jurisdiction + instructions for installing software necessary to reproduce and potentially augment the models underpinning the report. This will allow Tasmanian State Government to rerun the models with minor modifications or new elevation data as needed and also potentially replicate the methodology - to other locations within the model extent, if desired. + for other locations within the model extent, if desired.
The data provided on this DVD is:
-
- The Professional Opinion Report 2009/XX entitled -
- Figures of maximum inundation and maximum speed -
- Generated Arc grids of maximum inundation and maximum speed for the areas of interest -
- Generated timeseries from the models +
- The Professional Opinion Report 2009/07 entitled +
- Figures of maximum inundation depth and maximum flow speed +
- Generated ArcGIS grids of maximum inundation depth and maximum flow speed for the areas of interest +
- Generated timeseries from the models
- Combined elevation data used by the simulation -
- The Arc grid of the resulting elevation generated by ANUGA -
- The Python scripts used to run the models +
- The ArcGIS grid of the resulting elevation + generated by the model software; ANUGA +
- The Python scripts used to run the models
- The tsunami wave data on this disk are all derived from events available in the Australian Tsunami - Hazard Map. As this dataset is very large, this DVD only contains the events + The input tsunami wave data on this disk are all derived from events available in the Probabilistic Tsunami Hazard Assessment (PTHA) for Australia. As this dataset is very large, this DVD only contains the events described in the report. Index: /DVD_images/extra_files/Hobart/installation.html =================================================================== --- /DVD_images/extra_files/Hobart/installation.html (revision 7330) +++ /DVD_images/extra_files/Hobart/installation.html (revision 7331) @@ -16,5 +16,5 @@
- @@ -31,7 +31,7 @@ The tsunami inundation simulations are based on the Open Source software package called ANUGA. - For more information about ANUGA visit https://datamining.anu.edu.au/anuga. + For more information about ANUGA visit https://datamining.anu.edu.au/anuga.Tsunami Inundation Models for the HOBART region
+Tsunami Inundation Models for south east Tasmania
- ANUGA has been developed and tested under the Windows XP and Linux (Ubuntu, Mint and Red Hat) operating systems. + ANUGA has been developed and tested under the Windows XP and Linux (Ubuntu, Mint and Red Hat) operating systems. ANUGA may also be installed on a Windows Vista system and Debian Linux, but these have not been extensively tested.
@@ -40,8 +40,10 @@
Installing ANUGA
- + Download Software: please note Internet Explorer has the habit of renaming the + .tgz file to .gz - the remedy is to rename them back or use another browser such as Firefox. +You must install ANUGA according to the Installation Guide. As the installation procedure comprises a number of steps we suggest that it is carried out by someone with a - basic understanding of how to install software packages on the operating system of choice. + basic understanding of how to install software packages on the operating system of choice.
How to run a simulation
@@ -52,5 +54,5 @@ can be used to get help.- +
@@ -61,22 +63,21 @@
- Next, you must copy the entire contents of the DVD to a place in your filesystem. Let's assume - that you want to put the data into C:\ANUGA. You would do this in a Command Prompt window: -
-
- C:
- cd \
- mkdir ANUGA
+ Next, you must copy the entire contents of the DVD to a place in your file system. Let's assume + that you want to put the data into C:\ANUGA. You could drag all files from the DVD into this folder with Explorer. + Or through the Command Prompt window: ++ C: + cd + mkdir ANUGA xcopy /e X:\ ANUGA -- Of course, you could just drag all files on the DVD (X:\ is the assumed drive letter) to the desired target directory with Explorer. +
Now you must create an environment variable ANUGADATA that points to the directory you just created: -
-
+ set ANUGADATA=C:\ANUGA -+
Next, you prepare and execute the model: -
- cd C:\ANUGA - cd project - python setup_model.py + + cd C:\ANUGA + cd project + python setup_model.py python run_model.py - |
+
This should all run to completion though it may take several days, depending on the hardware used. @@ -105,8 +102,8 @@ drilling down from that point until you find a directory called outputs. The latest directory in outputs should contain a file called screen_error.txt - with an explanation of the error at the end. + with an explanation of the error at the end.
Once the simulation has completed you may view the results using the ANUGA
- viewer as described in the Installation Guide
+ viewer as described in the Installation Guide
and User Manual. You may also
generate rasters suitable for GIS mapping or write dedicated scripts extracting timeseries at selected
@@ -120,5 +117,5 @@
After you have run a simulation you may wish to extract information from the generated SWW file.
On this disk there are a selection of example analysis scripts that you can use or modify.
- The analysis scripts page explains how to use the scripts.
+ The analysis page explains how to use the scripts.
How to modify a simulation
@@ -138,8 +135,6 @@
Frequently Asked Questions
- There is a Frequently Asked Questions (FAQ) page at
- https://datamining.anu.edu.au/anuga/wiki/FrequentlyAskedQuestions.
- Look in there and search the mailing list archives before asking questions on the ANUGA user's mailing list.
- Someone may already have asked your question!
+ Look in Frequently Asked Questions (FAQ) and search the mailing list archives before asking questions on the ANUGA user's mailing list.
+ Someone may already have asked and answered your question!
Note that you can edit the FAQ pages, so it is helpful if you can update the FAQ with your question and its answer if you feel Index: /DVD_images/extra_files/Hobart/modifications.html =================================================================== --- /DVD_images/extra_files/Hobart/modifications.html (revision 7330) +++ /DVD_images/extra_files/Hobart/modifications.html (revision 7331) @@ -16,5 +16,5 @@
- First we describe the uses of the script files you might change, then we walk through a few examples + First we describe the usage of the script files you might change, then we walk through a few examples of changes you might make.
@@ -37,5 +37,5 @@
project.py | Defines the input data used, where to place output, etc. |
build_elevation.py | Builds a PTS file of all elevation data specified by project.py |
build_elevation.py | Combines the elevation data specified by project.py into one file; with the extension .PTS |
setup_model.py | Prepares the simulation before actually running it |
run_model.py | Runs the simulation |
event_number = 51436 # 1 in 10000 yr event from New Hebrides |
+ The event numbers correspond to a quake ID from the Probabilistic Tsunami Hazard Map Assessment of Australia. + +
+ event_number = 58260 # 1 in 10000 yr event from Puysegur Trench+
Elevation
- Elevation data can be changed in the project.py elevation script.
+ Elevation data can be changed in the project.py script under ELEVATION DATA.
Elevation data can be read as either a point file, comma delimited, or as an ASCII grid file
- (ASC) with an accompanying projection file (PRJ). All elevation input must be projected in the correct UTM zone.
+ (ASC) with an accompanying projection file (PRJ). All elevation input should sit in topographies and must be projected in the correct UTM zone.
A header for a CSV file has the format: -
x,y,elevation |
+ x,y,elevation+
An ASC file header has the format: -
ncols 868 -nrows 856 -xllcorner 418933.86055096 -yllcorner 5151810.6668096 -cellsize 250 -NODATA_value -9999 |
+ ncols 868 + nrows 856 + xllcorner 418933.86055096 + yllcorner 5151810.6668096 + cellsize 250 + NODATA_value -9999+
The header of a PRJ file has the format: -
Projection UTM -Zone 55 -Datum D_GDA_1994 -Zunits NO -Units METERS -Spheroid GRS_1980 -Xshift 500000 -Yshift 10000000 -Parameters |
- - The elevation filenames must be listed in either point_filenames or ascii_grid_filenames + +
+ Projection UTM + Zone 55 + Datum D_GDA_1994 + Zunits NO + Units METERS + Spheroid GRS_1980 + Xshift 500000 + Yshift 10000000 + Parameters+ +
+ + The elevation filenames in project.py must be listed in either point_filenames or ascii_grid_filenames depending on their format. Point files need to have their extension shown however the ascii grid files have the .asc extension assumed: -
point_filenames = ['SD100031996_jgriffin_clip.csv', - 'tomaga_offshore_AHD_MGA_1997.csv', - 'hobart_BBHD_MGA_1995.csv', - 'hobart_AHD_MGA_2000.csv'] - -ascii_grid_filenames = ['sd100031996_p', - 'sd100031996_p2', - 'sd100031996_p3', - 'sd100031996_p4'] |
+ point_filenames = ['point1.csv', + 'point2.csv', + 'point3.csv'] + + ascii_grid_filenames = ['grid1', + 'grid2', + 'grid3']+
@@ -183,6 +188,39 @@
Interior regions
- The interior_regions parameter allows you to change the mesh of the model.
+ The user can specify a number of internal polygons within each of which the resolution of the mesh can be specified.
+ Mesh resolution is the maximum allowable area specified for each region, defining the largest area an indivdual
+ triangular element of the mesh can take (and therefore the minimum mesh resolution).
+ These polygons need to be nested within each other with no overlapping edges.
+
+ The interior regions can be changed in the project.py script under INTERIOR REGIONS. + Interior regions can be read as either seperate CSV files for each polygon displayed as a listed paired with its + resolution and/or one CSV file for all polygons, where its resolution is defined within the csv under 'id'. + All file inputs should sit in polygons and must be projected in the correct UTM zone. +
+ The format for a CSV file with ONE polygon has the format: +
+ easting,northing Note: NO Header+ + + The header for a CSV file with MANY polygons has the format: +
+ easting,northing,id,value+
+ Where id = polygon number and value = maximum allowable area. +
+ The project.py script for this section looks like this: + + +
+ interior_regions_list = [['aos1.csv', 1500], + ['aos2.csv', 1500], + ['sw.csv', 30000]] + interior_regions_multiple_csv = 'PriorityAreas.csv'+ +
+ + For further information on ANUGA file formats please see the ANUGA User Manual, section 5.1.