Index: /anuga_work/production/dampier_2006/make_report.py
===================================================================
--- /anuga_work/production/dampier_2006/make_report.py	(revision 4134)
+++ /anuga_work/production/dampier_2006/make_report.py	(revision 4134)
@@ -0,0 +1,328 @@
+"""
+Generate report for production run
+
+NOTE - this will only work on Windows system (as pylab not installed on unix).
+
+Inputs:
+
+report_title:    title to be included in tex file
+production dirs: dictionary of production directories with a
+                 association to that simulation run, eg high tide,
+                 low tide and MSL.
+                    
+
+Outputs:
+
+* Report generated to scenario_report.tex where the scenario relates
+the name of the production directory this script is being run from.
+* figures used for report stored in the report_figure directory
+NOTE, this directory will need to be committed, as well as
+the latex files.
+
+The report structure is
+
+* Executive Summary
+* Introduction
+* Modelling Methodology
+* Tsunami scenario
+* Data sources
+* Inundation modelling results
+* Impact modelling
+* Summary
+* Acknowledgements
+* References
+* Appendix: Metadata
+* Appendix: Inundation model metadata
+* Appendix: Time series outputs
+
+Other files included in document which require manual intervention:
+
+* an abstract must be written in abstract.tex
+* an introduction must be written in introduction.tex; a basic outline and
+  some of the core inputs are already in place
+* the tsunami-genic event should be discussed in tsunami_scenario.tex
+* a discussion of the ANUGA model is required in anuga.tex
+* a computational_setup.tex file needs to be written for the particular scenario
+* the interpretation of the results needs to be written to interpretation.tex
+* maximum inundation map names need to be included in HAT_map and LAT_map etc.
+* damage modelling map names need to be included in HAT_damage and LAT_damage etc.
+* a summary must be written into summary.tex
+* metadata for the scenario data to be included in metadata.tex
+
+May, June 2006                    
+"""
+
+from os import getcwd, sep, altsep, mkdir, access, F_OK
+import project
+from anuga.abstract_2d_finite_volumes.util import sww2timeseries, get_gauges_from_file
+
+# Derive scenario name
+p = getcwd().split(sep)
+scenario = p[-1] # Last element of absolute CWD path
+scenario_name = scenario.split('_2006')[0] # Strip any text past `_2006`
+test = scenario_name.split('_')
+if len(test) <> 1:
+    scenario_name = '%s %s' %(test[0], test[1])
+
+# Create report directory
+reportdir = getcwd()+sep+'report07'+sep
+if access(reportdir,F_OK) == 0:
+    mkdir (reportdir)
+    
+# User defined inputs
+report_title = 'Tsunami impact modelling for %s' %scenario_name.title()
+
+# WA DLI data
+production_dirs = {#'20061101_053132_run': 'big event'} # parallel
+                    #'20061101_054432_run': 'big event'} # sequential
+                    #'20061107_070805_run': 'Mw 9.0',
+                    '20061116_055135_run': 'Mw 9.0'} # multiple interior regions
+                    #'20061116_062924_run': 'Mw 9.0'} # one interior region
+
+is_parallel = True
+if is_parallel == True:
+    nodes = 8
+
+max_maps = {'Mw 9.0': 'mw9_map'}
+
+gauge_map1 = 'dampier_gauges_1.jpg'
+gauge_map2 = 'dampier_gauges_2.jpg'
+
+# Create sections and graphs for each designated production directory
+latex_output = []
+report_name = 'latexoutput_cipma'
+
+if is_parallel == True:
+
+    for i in range(nodes):
+        print 'Sending node %d of %d' %(i,nodes)
+        swwfiles = {}
+        reportname = report_name + '_%s' %(i)
+        for label_id in production_dirs.keys():
+            file_loc = project.output_dir + label_id + sep
+            sww_extra = '_P%s_%s' %(i,nodes)
+            swwfile = file_loc + project.scenario_name + sww_extra + '.sww'
+            swwfiles[swwfile] = label_id
+
+            texname, elev_output = sww2timeseries(swwfiles,
+                                                  project.gauge_filename,
+                                                  production_dirs,
+                                                  report = True,
+                                                  reportname = reportname,
+                                                  plot_quantity = ['stage', 'momentum'],
+                                                  surface = False,
+                                                  time_min = None,
+                                                  time_max = None,
+                                                  title_on = False,
+                                                  verbose = True)
+            
+            latex_output.append(texname)
+    
+else:
+    
+    for label_id in production_dirs.keys():
+
+        file_loc = project.output_dir + label_id + sep
+        swwfile = file_loc + project.scenario_name + '.sww'
+        swwfiles[swwfile] = label_id
+        
+    texname, elev_output = sww2timeseries(swwfiles,
+                                          project.gauge_filename,
+                                          production_dirs,
+                                          report = True,
+                                          reportname = report_name,
+                                          plot_quantity = ['stage', 'momentum'],
+                                          surface = False,
+                                          time_min = None,
+                                          time_max = None,
+                                          title_on = False,
+                                          verbose = True)
+
+# Start report generation
+# Future: generate_report(reportdir, scenario, report_title, 
+# project.gauge_filename, max_maps, damage_maps, production_dirs, latex_output)
+report_name = reportdir + scenario + '_report_CIPMA.tex'
+fid = open(report_name, 'w')
+
+s = """
+% This is based on an automatically generated file (by make_report.py).
+%
+% Manual parts are:
+% * an abstract must be written in abstract.tex
+% * an introduction must be written in introduction.tex; a basic outline and
+%   some of the core inputs are already in place
+% * outline of the modelling methodology provided in modelling_methodology.tex
+% * the tsunami-genic event should be discussed in tsunami_scenario.tex
+% * an computational_setup.tex file needs to be written for the particular scenario
+% * the interpretation of the results needs to be written to interpretation.tex
+% * maximum inundation maps need to be included in HAT_map.tex and LAT_map.tex etc.
+% * damage modelling maps need to be included in HAT_damage and LAT_damage etc.
+% * a summary must be written into summary.tex
+% * metadata for the scenario data to be included in metadata.tex
+
+\documentclass{article}
+
+\usepackage{ae} % or {zefonts}
+\usepackage[T1]{fontenc}
+\usepackage[ansinew]{inputenc}
+\usepackage{amsmath}
+\usepackage{amssymb}
+\usepackage{graphicx}
+\usepackage{color}
+\usepackage[colorlinks]{hyperref}
+\usepackage{lscape} %landcape pages support
+\usepackage{setspace}
+\usepackage{rotating}
+\usepackage{pdfpages}
+\include{appendix}
+\setstretch{1.25}
+\\topmargin 0pt
+\oddsidemargin 0pt
+\evensidemargin 0pt
+\marginparwidth 0.5pt
+\\textwidth \paperwidth
+\\advance\\textwidth -2in
+
+"""
+fid.write(s)
+
+s = """
+\date{\\today}
+%\\author{Geoscience Australia}
+
+\\begin{document}
+\\title{
+\\begin{figure}[hbt]
+  \centerline{ \includegraphics[scale=0.4]{../report_figures/GAlogo.jpg}}
+\end{figure}
+"""
+fid.write(s)
+s = '%s} '%report_title
+fid.write(s)
+s = """
+  \maketitle
+
+    \section{Executive Summary}
+    \label{sec:execsum}
+  \input{execsum}
+
+  \\tableofcontents
+  
+   \section{Modelling methodology}
+    \label{sec:methodology}
+    \input{modelling_methodology}
+
+  \section{Data sources}
+    \label{sec:data}
+    \input{data}
+    
+   %\section{Inundation model}
+   % \label{sec:anuga}
+   % \input{anuga}
+   % \input{computational_setup}
+        
+  \section{Inundation modelling results}
+     \label{sec:results}
+         
+"""
+fid.write(s)
+
+# Generate latex output for location points
+##s = '\\begin{table} \\begin{center} \n'
+##fid.write(s)
+##s = '\caption{Defined point locations for %s study area.}' %scenario_name
+##fid.write(s)
+##s = """
+##\label{table:locations}
+##\\begin{tabular}{|l|l|l|l|}\hline
+##\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline
+##"""
+##fid.write(s)
+##
+##gauges, locations, elevation = get_gauges_from_file(project.gauge_filename)
+##
+##for name, gauges, elev in zip(locations, gauges, elevation):
+##    east = gauges[0]
+##    north = gauges[1]
+##    s = '%s & %.2f & %.2f & %.2f \\\\ \hline \n' %(name.replace('_',' '), east, north, elev)
+##    fid.write(s)
+##
+##s = '\\end{tabular} \n  \end{center} \n \end{table} \n \n'
+##fid.write(s)
+##
+
+s = '\\begin{figure}[hbt] \centerline{\includegraphics[scale = 0.5]{../report_figures/%s}}' %gauge_map2
+fid.write(s)
+s  = """
+\caption{Point locations used for Dampier study; surrounding the facility.
+See Appendix \\ref{sec:timeseries} for stage and momentum results.
+}  
+\label{fig:points2}
+\end{figure}
+"""
+fid.write(s)
+s = '\\begin{figure}[hbt] \centerline{\includegraphics[scale=0.5]{../report_figures/%s}}' %gauge_map3
+fid.write(s)
+s  = """
+\caption{Point locations used for Dampier study; section
+of the pipeline towards the boundary of the region and other boundary gauges.
+See Appendix \\ref{sec:timeseries} for stage and momentum results.}  
+\label{fig:points3}
+\end{figure}
+"""
+fid.write(s)
+    
+#s = '\input{interpretation} \n'
+#fid.write(s)
+
+# Assign titles to each production section 
+# Must specify one name per section
+for i, name in enumerate(production_dirs.keys()):
+
+    s = '\input{%s} \n \clearpage \n \n' %max_maps[production_dirs[name]]
+    fid.write(s)
+
+# Closing
+
+
+s = """
+
+     \section{Summary}
+     \label{sec:summary}
+     \input{summary}
+
+     %\section{Acknowledgements}
+     %\input{acknowledgements}
+     
+    \input{references}
+
+    \\appendix
+
+   \section{ANUGA modelling parameters}
+\label{sec:anugasetup}
+\input{anuga_setup}
+
+   \section{Metadata}
+     \label{sec:metadata}
+     \input{metadata}
+     
+\section{ANUGA modelling parameters}
+\label{sec:anugasetup}
+\input{anuga_setup}
+
+\clearpage
+
+   \section{Time series}
+     \label{sec:timeseries}
+"""
+fid.write(s)
+
+for i in range(len(latex_output)):
+    if latex_output[i] <> '':
+        s = '\input{%s} \n \clearpage \n \n' %latex_output[i]   
+        fid.write(s)
+
+s="""
+\end{document}
+"""
+fid.write(s)
Index: /anuga_work/production/dampier_2006/make_report_cipma.py
===================================================================
--- /anuga_work/production/dampier_2006/make_report_cipma.py	(revision 4133)
+++ /anuga_work/production/dampier_2006/make_report_cipma.py	(revision 4134)
@@ -27,5 +27,4 @@
 * Tsunami scenario
 * Data sources
-* Inundation model
 * Inundation modelling results
 * Impact modelling
@@ -34,4 +33,5 @@
 * References
 * Appendix: Metadata
+* Appendix: Inundation model metadata
 * Appendix: Time series outputs
 
@@ -236,5 +236,5 @@
 ##\label{table:locations}
 ##\\begin{tabular}{|l|l|l|l|}\hline
-##\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation}\\\\ \hline
+##\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline
 ##"""
 ##fid.write(s)
@@ -298,5 +298,9 @@
 
     \\appendix
-    
+
+   \section{ANUGA modelling parameters}
+\label{sec:anugasetup}
+\input{anuga_setup}
+
    \section{Metadata}
      \label{sec:metadata}
Index: /anuga_work/production/dampier_2006/report/anuga_setup.tex
===================================================================
--- /anuga_work/production/dampier_2006/report/anuga_setup.tex	(revision 4133)
+++ /anuga_work/production/dampier_2006/report/anuga_setup.tex	(revision 4134)
@@ -78,11 +78,11 @@
 \begin{tabular}{|l|l|l|}\hline
 Mesh & & \hline
-& resolution in Region 1 & 500 m$^2$ \hline
-& resolution in Region 2 & 2000 m$^2$ \hline
-& resolution in Region 3 & 2000 m$^2$ \hline
-& remaining resolution & 100 000 m$^2$ \hline
-Model parameters & & \hline
-& friction & 0
-& minimum stored height & 0.1 m \hline
+& resolution in Region 1 & 500 m$^2$  \\ \hline
+& resolution in Region 2 & 2000 m$^2$ \\ \hline
+& resolution in Region 3 & 2000 m$^2$ \\ \hline
+& remaining resolution & 100 000 m$^2$ \\ \hline
+Model parameters & & \\ \hline
+& friction & 0 \\ \hline
+& minimum stored height & 0.1 m \\ \hline
 \end{tabular}
 \end{center}
Index: /anuga_work/production/onslow_2006/make_report.py
===================================================================
--- /anuga_work/production/onslow_2006/make_report.py	(revision 4133)
+++ /anuga_work/production/onslow_2006/make_report.py	(revision 4134)
@@ -25,5 +25,5 @@
 * Tsunami scenario
 * Data sources
-* Inundation model
+
 * Inundation modelling results
 * Impact modelling
@@ -32,4 +32,5 @@
 * References
 * Appendix: Metadata
+* Appendix: Inundation model metadata
 * Appendix: Time series outputs
 
@@ -53,5 +54,5 @@
 from os import getcwd, sep, altsep, mkdir, access, F_OK
 import project
-from anuga.pyvolution.util import sww2timeseries, get_gauges_from_file
+from anuga.abstract_2d_finite_volumes.util import sww2timeseries, get_gauges_from_file
 
 # Derive scenario name
@@ -103,4 +104,5 @@
 texname, elev_output = sww2timeseries(swwfiles,
                                       project.gauge_filename,
+                                      #project.gauge_filename_bindi,
                                       production_dirs,
                                       report = True,
@@ -202,9 +204,4 @@
     \label{sec:data}
     \input{data}
-    
-   \section{Inundation model}
-    \label{sec:anuga}
-    \input{anuga}
-    \input{computational_setup}
         
   \section{Inundation modelling results}
@@ -222,5 +219,5 @@
 \label{table:locations}
 \\begin{tabular}{|l|l|l|l|}\hline
-\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation}\\\\ \hline
+\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline
 """
 fid.write(s)
@@ -289,13 +286,20 @@
 
     \\appendix
-    
+
+    \section{ANUGA modelling parameters}
+    \label{sec:anugasetup}
+    \input{anuga_setup}
+
+\clearpage
+
    \section{Metadata}
      \label{sec:metadata}
      \input{metadata}
 
-\pagebreak
+\clearpage
 
    \section{Time series}
      \label{sec:timeseries}
+     \input{timeseriesdiscussion}
 """
 fid.write(s)
Index: /anuga_work/production/onslow_2006/project.py
===================================================================
--- /anuga_work/production/onslow_2006/project.py	(revision 4133)
+++ /anuga_work/production/onslow_2006/project.py	(revision 4134)
@@ -74,4 +74,5 @@
 #for MOST
 gauge_filename = gaugedir + 'gauge_location_onslow.csv'
+gauge_filename_bindi = gaugedir + 'gauge_location_bindi.csv'
 gauges50 = gaugedir + '50_gauges.xya'
 gauge_comparison = gaugedir + 'MOST_comparison_gauges.xya'
Index: uga_work/production/onslow_2006/report/anuga.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/anuga.tex	(revision 4133)
+++ 	(revision )
@@ -1,46 +1,0 @@
-The software tool, ANUGA \cite{ON:modsim}, has been used to estimate 
-the inundation extent
-and associated water level at various points in space and time.
-ANUGA has been developed by GA and the Australian National University
-(ANU) to solve the nonlinear shallow water
-wave equation using the finite volume technique\footnote{The finite volume
-technique belongs to the class of computational fluid dynamic (CFD)
-methods which is based on discretizing the study area in
-control ''volumes''. The method ensures that conservation
-of mass and horizontal momentum is satisfied for
-each control volume.
-An advantage of this technique is that the discretization
-can be changed
-according to areas of interest and that wetting and drying
-is treated robustly as part of the numerical scheme.}. 
-ANUGA is continually being refined and validated to ensure
-the modelling approximations are as accurate as possible.
-However, model sensitivity to errors in bathymetric data, 
-frictional resistance of the seafloor and the size of the tsunamigenic event are not well understood and the topic of ongoing research. 
-As such, the current results are preliminary.
-
-The following information is required to undertake the 
-inundation modelling;
-
-\begin{itemize}
-\item onshore and offshore elevation data (topographic and bathymetric data, 
-see Section \ref{sec:data}),
-\item initial conditions, such as initial water levels (e.g. determined by tides),
-\item boundary conditions (the tsunami source as described in 
-Section \ref{sec:methodology}), and
-\item computational requirements relating to the mesh construction.
-\end{itemize}
- 
-The initial conditions used for this scenario are MSL, HAT and LAT which were
-defined in Section \ref{sec:data}.
-The dynamics of
-tidal effects (that is, the changes in water height over time for
-the entire study area) are not currently modelled.
-Sea floor friction will generally provide resistance to the water flow
-and thus reduce the impact somewhat. However, limited
-research has been carried out to determine 
-the friction coefficients, and 
-thus it has not been incorporated
-in the scenario. The 
-results are therefore likely to be over estimates.
-
Index: /anuga_work/production/onslow_2006/report/anuga_setup.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/anuga_setup.tex	(revision 4134)
+++ /anuga_work/production/onslow_2006/report/anuga_setup.tex	(revision 4134)
@@ -0,0 +1,99 @@
+The following information is required to undertake the 
+inundation modelling;
+
+\begin{itemize}
+\item onshore and offshore elevation data (topographic and bathymetric data, 
+see Section \ref{sec:data}),
+\item initial conditions, such as initial water levels (e.g. determined by tides),
+\item boundary conditions (the tsunami source as described in 
+Section \ref{sec:methodology}), and
+\item computational requirements relating to the mesh construction.
+\end{itemize}
+ 
+The initial conditions used for this scenario are MSL, HAT and LAT which were
+defined in Section \ref{sec:data}.
+The dynamics of
+tidal effects (that is, the changes in water height over time for
+the entire study area) are not currently modelled.
+Sea floor friction will generally provide resistance to the water flow
+and thus reduce the impact somewhat. However, limited
+research has been carried out to determine 
+the friction coefficients, and 
+thus it has not been incorporated
+in the scenario. The 
+results are therefore likely to be over estimates.
+
+To set up a model for the tsunami scenario, a study area is first 
+determined. Preliminary investigations have indicated that the
+output from MOST or URS should be input to ANUGA 
+at the 100m water depth. Historical run-up heights are 
+of the order of 10 m and we would expect that a tsunami wave
+would penetrate no higher for this scenario, hence we have
+bounded our study region at 10m elevation.
+Current computation requirements define a coastline
+extent of around 100 km. Therefore, the study area of around 6300 km$^2$ 
+covers approximately 100 km of
+coastline and extends offshore to the 100m contour line and inshore to 
+approximately 10m elevation. 
+
+The finite volume technique relies on the construction of a triangular mesh which covers the study region. 
+This mesh can be altered to suit the needs of the scenario in question. The mesh can be refined in areas of 
+interest, particularly in the coastal region where complex behaviour is likely to occur. 
+In setting up the model, the user defines the area of the triangular cells in each region of interest\footnote{Note that the cell 
+area will be the maximum cell area within the defined region and that each 
+cell in the region does not necessarily have the same area.}.
+The cell areas should not be too small as this will cause unrealisticly long computational time, 
+and not too great as this may inadequately capture important behaviour. 
+%There are no gains in choosing the area to be less than the supporting data.
+Figure \ref{fig:onslow_area} shows the study area with regions of difference cell areas. The total number
+of cells is 401939. 
+Lateral accuracy refers to the distance at which we are confident in stating a region is inundated. 
+Figure \ref{fig:onslow_area} shows the maximum triangular cell area and lateral accuracy for each region. 
+Therefore we can only be confident in the calculated inundation extent in the Onslow town centre to within 30 m.
+
+
+\begin{figure}[hbt]
+
+  \centerline{ \includegraphics[scale=0.15]{../report_figures/onslow_resolution_zones.jpg}}
+
+  \caption{Study area for Onslow scenario highlighting four regions of increased refinement.
+Region 1: Surrounds Onslow town centre with a cell area of 500 m$^2$ (lateral accuracy 30 m).
+Region 2: Surrounds the coastal region with a cell area of 2500 m$^2$ (lateral accuracy 70 m).
+Region 3: Water depths to the 50m contour line (approximately) with a cell area of 20000 m$^2$ (lateral accuracy 200 m).
+Region 4: Water depths to the boundary (approximately 100m contour line) with a cell area of 100000 m$^2$ (lateral accuracy 445 m).
+}
+  \label{fig:onslow_area}
+\end{figure}
+
+
+The final item to be addressed to complete the model setup is the 
+definition of the boundary condition. As 
+discussed in Section \ref{sec:methodology}, a range of events with
+different return periods will modelled. It must be 
+remembered that tsunamis are typically made up of a series
+of waves with different amplitudes which is affected by the energy
+and style of the event as well as the bathymetry whilst it travels
+from its source to Onslow. The amplitude and velocity of each of these
+waves are then provided to ANUGA as boundary conditions and propagated
+inshore.
+
+Sea floor friction will generally provide resistance to the water flow
+and thus reduce the impact somewhat. Initial experimental results 
+indicate that the friction coefficient shoudl be 0.017.
+The following table summarises the modelling parameters;
+
+\begin{table}
+\begin{center}
+\caption{Parameters used in ANUGA for the Onslow scenario.}
+\begin{tabular}{|l|l|l|}\hline
+Mesh & & \hline
+& resolution in Region 1 & 500 m$^2$ \hline
+& resolution in Region 2 & 2500 m$^2$ \hline
+& resolution in Region 3 & 20000 m$^2$ \hline
+& remaining resolution & 100 000 m$^2$ \hline
+Model parameters & & \hline
+& friction & 0.017
+& minimum stored height & 0.1 m \hline
+\end{tabular}
+\end{center}
+\end{table}
Index: uga_work/production/onslow_2006/report/computational_setup.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/computational_setup.tex	(revision 4133)
+++ 	(revision )
@@ -1,98 +1,0 @@
-To set up a model for the tsunami scenario, a study area is first 
-determined. Preliminary investigations have indicated that the
-output from MOST should be input to ANUGA 
-at the 100m water depth\footnote{ 
-Preliminary investigations indicate that MOST and ANUGA compare
-well at a water depth of 100 m. In addition, the resolution for
-the MOST modelling indicate that it can theoretically model a
-tsunami wave with a wavelength of 20-30 km. The wavelength of
-the tsunami wave at the boundary in this scenario is approximately 20 km. A much
-higher model resolution will be used in developing the probabilistic
-models for further studies so that tsunami waves with shorter wavelengths
-can be captured.}. Historical run-up heights are 
-of the order of 10 m and we would expect that a tsunami wave
-would penetrate no higher for this scenario, hence we have
-bounded our study region at 10m elevation.
-Current computation requirements define a coastline
-extent of around 100 km. Therefore, the study area of around 6300 km$^2$ 
-covers approximately 100 km of
-coastline and extends offshore to the 100m contour line and inshore to 
-approximately 10m elevation. 
-
-The finite volume technique relies on the construction of a triangular mesh which covers the study region. 
-This mesh can be altered to suit the needs of the scenario in question. The mesh can be refined in areas of 
-interest, particularly in the coastal region where complex behaviour is likely to occur. 
-In setting up the model, the user defines the area of the triangular cells in each region of interest\footnote{Note that the cell 
-area will be the maximum cell area within the defined region and that each 
-cell in the region does not necessarily have the same area.}.
-The cell areas should not be too small as this will cause unrealisticly long computational time, 
-and not too great as this may inadequately capture important behaviour. 
-%There are no gains in choosing the area to be less than the supporting data.
-Figure \ref{fig:onslow_area} shows the study area with regions of difference cell areas. The total number
-of cells is 401939. 
-Lateral accuracy refers to the distance at which we are confident in stating a region is inundated. 
-Figure \ref{fig:onslow_area} shows the maximum triangular cell area and lateral accuracy for each region. 
-Therefore we can only be confident in the calculated inundation extent in the Onslow town centre to within 30 m.
-
-
-\begin{figure}[hbt]
-
-  \centerline{ \includegraphics[scale=0.15]{../report_figures/onslow_resolution_zones.jpg}}
-
-  \caption{Study area for Onslow scenario highlighting four regions of increased refinement.
-Region 1: Surrounds Onslow town centre with a cell area of 500 m$^2$ (lateral accuracy 30 m).
-Region 2: Surrounds the coastal region with a cell area of 2500 m$^2$ (lateral accuracy 70 m).
-Region 3: Water depths to the 50m contour line (approximately) with a cell area of 20000 m$^2$ (lateral accuracy 200 m).
-Region 4: Water depths to the boundary (approximately 100m contour line) with a cell area of 100000 m$^2$ (lateral accuracy 445 m).
-}
-  \label{fig:onslow_area}
-\end{figure}
-
-%\begin{figure}[hbt]
-%
-%  \centerline{ \includegraphics[width=100mm, height=75mm]
-%              {../report_figures/mesh.jpg}}
-
-%  \caption{Computational mesh for Onslow study area where the
-%cell areas increase in resolution; 500 m$^2$, 2500 m$^2$, 20000 
-%m$^2$ and 100000 m$^2$.}
-%  \label{fig:mesh_onslow}
-%\end{figure}
-
-The final item to be addressed to complete the model setup is the 
-definition of the boundary condition. As 
-discussed in Section \ref{sec:methodology}, a Mw 9 event provides
-the tsunami source. The resultant tsunami wave is made up of a series
-of waves with different amplitudes which is affected by the energy
-and style of the event as well as the bathymetry whilst it travels
-from its source to Onslow. The amplitude and velocity of each of these
-waves are then provided to ANUGA as boundary conditions and propagated
-inshore.
-%To complete the model setup, we illustrate the
-%tsunami wave from the earthquake source described
-%in Section \ref{sec:tsunamiscenario} which is used as the boundary condition,
-%as described in Section \ref{sec:methodology}. 
-%MOST was used to initiate the event and propagate the wave in deep water. 
-%ANUGA uses the MOST wave amplitude and velocity at
-%the boundary (the 100m contour line as shown in Figure \ref{fig:onslow_area})
-%and continues to propagate the wave in shallow water and onshore. 
-%To illustrate the form of the tsunami wave, we show the 
-%tsunami wave moving through the point locations shown in 
-%Figure \ref{fig:MOSTsolution} as a surface showing the wave's
-%amplitude as a function of its spatial location and time.
-%This figure shows how the wave has been affected by the bathymetry in 
-%arriving at these locations as the amplitude is variable. It is also
-%important to note that the tsunami is made up of a series of
-%waves with different amplitudes.
-
-%\begin{figure}[hbt] 
-% \centering 
-% \begin{tabular}{cc} 
-%\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/point_line_3d.png}& 
-%\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/solution_surfaceMOST.png}\\ 
-%(a) & (b) \\
-%\end{tabular} 
-% \caption{Point locations used to illustrate the form of the tsunami wave and the
-%corresponding surface function.} 
-% \label{fig:MOSTsolution} 
-% \end{figure} 
Index: /anuga_work/production/onslow_2006/report/damage_inputs.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/damage_inputs.tex	(revision 4133)
+++ /anuga_work/production/onslow_2006/report/damage_inputs.tex	(revision 4134)
@@ -58,5 +58,5 @@
 spinal column injuries, or crush syndrome. \\ \hline
 Fatal
-Severity 4	&Instantaneously killed or mortally injured, \\ \hline
+Severity 4	&Instantaneously killed or mortally injured. \\ \hline
 \end{tabular*}
 \end{center}
Index: /anuga_work/production/onslow_2006/report/execsum.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/execsum.tex	(revision 4133)
+++ /anuga_work/production/onslow_2006/report/execsum.tex	(revision 4134)
@@ -1,5 +1,5 @@
 This report is being provided to the Fire and Emergency Services Authority
 (FESA) as part of the Collaborative Research Agreement (CRA) 
-with Geoscience Australia (GA).
+with Geoscience Australia (GA), Tsunami Impact Modelling for WA.
 FESA has recognised the potential vulnerability of the Western Australia
 coastline to tsunami originating from earthquakes on
@@ -11,7 +11,8 @@
 threat and develop detailed response plans for a range of plausible events.
 
-This report describes the modelling methodology and initial results
-for a specific tsunami-genic event as it impacts the Onslow township
-and its surrounds. In particular, maximum inundation maps are shown
+This report describes the modelling methodology and results
+for a number of tsunami-genic events with varied return periods
+as they impact the Onslow region.
+In particular, maximum inundation maps are shown
 and discussed 
 for the event occurring at mean sea level as well as 
@@ -20,15 +21,9 @@
 the numbers of persons affected. The Onslow township has approximately
 350 residential structures and a population of around 800.
-For this specific event at high tide, approximately
-100 houses are inundated with two of those collapsing. Approximately
-15-20\% of the population will sustain injuries, including fatalities.
 
-Future studies
-will present a series of scenarios for a range of return periods to 
-assist FESA in developing appropriate plans for a range of event impacts.
-This will also allow an assessment of the relative tsunami risk 
+The results of this study will allow an assessment of the relative tsunami risk 
 to communities along the NW Shelf of WA. 
 This report and the decision support tool are the 
-June 2006 deliverables of the Collaborative Research Agreement,
-Tsunami Impact Modelling for WA, between FESA and GA.
+June 2007 deliverables of the Collaborative Research Agreement
+between FESA and GA.
 
Index: /anuga_work/production/onslow_2006/report/interpretation.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/interpretation.tex	(revision 4133)
+++ /anuga_work/production/onslow_2006/report/interpretation.tex	(revision 4134)
@@ -1,146 +1,8 @@
-The main features of the
-tsunami wave and resultant inundation ashore is described in this section. 
-We have
-chosen a number of locations to illustrate the features
-of the tsunami as it approaches and impacts Onslow.
-These locations have been chosen as we believe they would 
-either be critical
-in an emergency situation, (e.g. the hospital and power station) or
-effect recovery efforts, (e.g. the airport and docks). These locations
-are described in Table \ref{table:locations} and shown in 
-Figure \ref{fig:points}. The water's stage and speed 
-at each of these locations are shown
-as a function of time in the series of graphs shown in
-Appendix \ref{sec:timeseries}. It is assumed that the earthquake is
-generated at the beginning of the simulation, i.e. time = 0 minutes.
-Stage is defined as the absolute
-water level (in metres) relative to AHD 
-\footnote{For an offshore location such as Beadon Bay West,
-the initial water level will be that of the tidal scenario. In the 
-case of MSL, this water level will be 0. As the tsunami wave moves
-through this point, the water height may grow and thus the stage will 
-represent the amplitude of the wave. For an onshore location such as the 
-Light Tower, the actual water depth will be the difference between
-the stage and the elevation at that point. Therefore, at the beginning
-of the simulation, there will be no water onshore and therefore
-the stage and the elevation will be identical.}. Both stage and speed 
-(in metres/second) for
-each scenario (HAT, MSL and LAT) are shown
-on consistent scales to allow comparison between point locations.
-As a useful benchmark, Table \ref{table:speedexamples}
-describes typical examples for a range of speeds found in the 
-simulations. 
-
-\begin{table}[h]
-\label{table:speedexamples}
-\caption{Examples of a range of velocities.}
-\begin{center}
-\begin{tabular}{|l|l|}\hline
-{\bf Velocity (m/s)} & {\bf Example} \\ \hline
-1 & leisurely stroll pace\\ \hline
-1.5 & average walking pace \\ \hline
-%2 & 100m Olympic male freestyle \\ \hline
-%3 & mackeral \\ \hline
-4 & average person can maintain running for 1000m \\ \hline
-%5 & blue whale \\ \hline
-10 & 100m Olympic male sprinter \\ \hline
-16 & car travelling in urban zones (60 km/hr) \\ \hline
-\end{tabular}
-\end{center}
-\end{table}
-
-A tsunami wave typically has a small amplitude and typically travels at 100's of kilometres per hour. 
-The low amplitude complicates the ability to detect
-the wave. As the water depth decreases, 
-the speed of the wave 
-decreases and the amplitude grows. Another important feature of tsunamis 
-is drawdown. This means that the water is seen to retreat from the beaches
-before a tsunami wave 
-impacts that location. Other features 
-include reflections (where the wave is redirected due to the 
-influence 
-of the coast) and shoaling (where the wave's amplitude is amplified
-close to the coast due to wave interactions).
-These features are seen in these scenarios, and are consistent
-for HAT, MSL and LAT. 
-There is a small wave, followed
-by a large drawdown and then a large secondary wave. 
-
-These 
-features are illustrated in Figure \ref{fig:gaugeBeadonBayeast}
-where a small wave can be seen at around 200 mins. For the HAT
-case (shown in blue), the amplitude
-of the wave at this location is around 0.8 m\footnote{In this
-scenario, the initial water level is 1.5 m, which means that
-the actual amplitude is the difference between the stage value
-and the initial water level; 2.3 - 1.5}.
-The drawdown of around 4.3 m (i.e. 2.3 - -2) then occurs at around 230 mins
-(i.e. 3.8 hours after the event has been generated), before
-the second wave arrives 
-with an amplitude of around 3.6 m (i.e. 4.1 - 1.5). A further wave
-is then evident a short time later (around 255 mins) 
-which further increases the amplitude to around 5 m (i.e. 6.6 - 1.5).
-These features are replicated at each of the offshore points (those
-points with negative elevation as shown in Table \ref{table:locations}).
-
-The wave amplitude is typically greater
-for those locations which are in the shallowest water. For example,
-the maximum wave amplitude at the Beadon Bay East location 
-(Figure \ref{fig:gaugeBeadonBayeast}) is over
-4.5m where the water depth would normally be 3.56 m. In the 
-Beadon Bay West location (Figure \ref{fig:gaugeBeadonBaywest})
-where the water depth would normally be 4.62 m, 
-the maximum wave amplitude is much less (around 3 m). The wave amplitude
-at the West of Groyne location (Figure \ref{fig:gaugeWestofGroyne})
-is not greater than that seen
-at the Beadon Bay East location, even though the water depth is 
-much less, at 2.11m. This is probably due to its proximity
-to the groyne\footnote{A groyne is a man made structure to combat
-coastal erosion.}
-which has impeded the tsunami wave to some degree. However, the 
-maximum speed found amongst the locations is at the West of Groyne
-point which is in the shallowest water.
-
-The speed of the tsunami sharply increases as it moves onshore. There
-is minimal inundation found at the locations chosen, with the Bindi Bindi
-community receiving the greatest inundation for all tidal scenarios.
-At HAT, the community would receive over 1 m of inundation with 
-the water moving through the community at approximately 16 m/s. Referring
-to Table \ref{table:speedexamples}, a person in this location could
-not outrun this water movement. A small amount of water is found
-at the hospital (10 cm). Whilst this seems minimal, the water is moving 
-at around 6 m/s which could dislodge some items if the water was able to enter the hospital.
- 
-The geography of the Onslow area has played a role in offering
-some protection to the Onslow community. The tsunami wave is 
-travelling from the north west of the area. Most of
-the inundation along the coast is that which is open to this
-direction.  
-The sand dunes west of Onslow 
-appear to have halted this tsunami wave
-(see Figure \ref{fig:MSL_max_inundation}) with limited
-inundation found on the town's side of the dunes.
-The inundation within the community has occurred due to the
-wave reflecting from the beach area west of the creek and
-returning towards the Onslow town itself.  
-There are also sand dunes east of the creek which have also
-halted inundation beyond them. 
-Currently, we do not model changes 
-to the bathymetry or topography due to effects of the water flow. 
-Therefore, we do not know whether these sand dunes would withstand the
-transmitted energy of the tsunami wave. 
-
-Water features such as rivers, creeks and estuaries also play a role
-in the inundation extent.
-The tsunami wave penetrates the creek east of Onslow with a wave height 
-over 2 m at the mouth 
-(Figure \ref{fig:gaugeBeadonCreekmouth}) for the HAT scenario.
-Inundation exceeds 1 m at the Beadon Creek south of dock location (Figure
-\ref{fig:gaugeBeadonCreeksouthofdock}) suggesting that the wave's 
-energy dissipates as inundation overflows from the creek. A large
-tidal flat region surrounds the southern parts of the creek and
-it is evident that the inundation is essentially caught in this
-area.
-
+The inundation extent calculated at Onslow will be described in this section with
+impact assessments following in Section \ref{sec:impact}.
+% there will need to be something in here for when doing a range of events for each return period.
+Figures \ref{fig:HAT_max_inundation}, \ref{fig:MSL_max_inundation} and 
+\ref{fig:LAT_max_inundation} illustrate the maximum inundation extent
+for the Mw 9 event occurring at HAT, MSL and LAT respectively. 
 As expected, there is greater inundation at HAT with increased
 extent. The major road
@@ -163,2 +25,41 @@
 free of inundation for each tidal scenario. Section \ref{sec:impact}
 details the impact estimates to the residential infrastructure.
+
+The geography of the Onslow area has played a role in offering
+some protection to the Onslow community. The tsunami wave is 
+travelling from the north west of the area. Most of
+the inundation along the coast is that which is open to this
+direction.  
+The sand dunes west of Onslow 
+appear to have halted this tsunami wave
+(see Figure \ref{fig:MSL_max_inundation}) with limited
+inundation found on the town's side of the dunes.
+The inundation within the community has occurred due to the
+wave reflecting from the beach area west of the creek and
+returning towards the Onslow town itself.  
+There are also sand dunes east of the creek which have also
+halted inundation beyond them. 
+Currently, we do not model changes 
+to the bathymetry or topography due to effects of the water flow. 
+Therefore, we do not know whether these sand dunes would withstand the
+transmitted energy of the tsunami wave. 
+Water features such as rivers, creeks and estuaries also play a role
+in the inundation extent.
+The tsunami wave penetrates the creek east of Onslow, however it
+is evident that the inundation is essentailly maintained in the
+large tidal flat region surrounding the southern parts of the creek.
+
+In addition to describing the maximum inundation extent,
+we have
+chosen a number of locations to illustrate the features
+of the tsunami as it approaches and impacts Onslow.
+These locations have been chosen as we believe they would 
+either be critical
+in an emergency situation, (e.g. the hospital and power station) or
+effect recovery efforts, (e.g. the airport and docks). These locations
+are described in Table \ref{table:locations} and shown in 
+Figure \ref{fig:points}. The water's stage and speed 
+at each of these locations are shown
+as a function of time in the series of graphs shown in
+Appendix \ref{sec:timeseries}. Discussion of the main features of the
+tsunami wave is also described in Appendix \ref{sec:timeseries}.
Index: /anuga_work/production/onslow_2006/report/latexoutput.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/latexoutput.tex	(revision 4133)
+++ /anuga_work/production/onslow_2006/report/latexoutput.tex	(revision 4134)
@@ -1,22 +1,2 @@
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonPointLoadingBerthstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonPointLoadingBerthspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Beadon Point Loading Berth location (elevation -8.70m)} 
- \label{fig:gaugeBeadonPointLoadingBerth} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeHospitalstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeHospitalspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Hospital location (elevation 6.02m)} 
- \label{fig:gaugeHospital} 
- \end{figure} 
- 
 \begin{figure}[hbt] 
  \centering 
@@ -25,116 +5,6 @@
 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBindiBindiCommunityspeed.png}\\ 
 \end{tabular} 
- \caption{Time series for stage and speed at Bindi Bindi Community location (elevation 4.08m)} 
+ \caption{Time series for depth and speed at Bindi Bindi Community location (elevation 4.08m)} 
  \label{fig:gaugeBindiBindiCommunity} 
  \end{figure} 
  
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugePowerStationstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugePowerStationspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Power Station location (elevation 5.81m)} 
- \label{fig:gaugePowerStation} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeAirportRunwaystage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeAirportRunwayspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Airport Runway location (elevation 5.49m)} 
- \label{fig:gaugeAirportRunway} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekDocksstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekDocksspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Beadon Creek Docks location (elevation 1.65m)} 
- \label{fig:gaugeBeadonCreekDocks} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeWestofGroynestage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeWestofGroynespeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at West of Groyne location (elevation -2.11m)} 
- \label{fig:gaugeWestofGroyne} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekmouthstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekmouthspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Beadon Creek mouth location (elevation -2.90m)} 
- \label{fig:gaugeBeadonCreekmouth} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreeksouthofdockstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreeksouthofdockspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Beadon Creek south of dock location (elevation -1.81m)} 
- \label{fig:gaugeBeadonCreeksouthofdock} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeCentredamwallstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeCentredamwallspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Centre dam wall location (elevation 2.09m)} 
- \label{fig:gaugeCentredamwall} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeDamoverflowstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeDamoverflowspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Dam overflow location (elevation 0.88m)} 
- \label{fig:gaugeDamoverflow} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeLightTowerstage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeLightTowerspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Light Tower location (elevation 3.86m)} 
- \label{fig:gaugeLightTower} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayweststage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBaywestspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Beadon Bay west location (elevation -4.62m)} 
- \label{fig:gaugeBeadonBaywest} 
- \end{figure} 
- 
-\begin{figure}[hbt] 
- \centering 
- \begin{tabular}{cc} 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayeaststage.png}& 
-\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayeastspeed.png}\\ 
-\end{tabular} 
- \caption{Time series for stage and speed at Beadon Bay east location (elevation -3.56m)} 
- \label{fig:gaugeBeadonBayeast} 
- \end{figure} 
- 
Index: /anuga_work/production/onslow_2006/report/modelling_methodology.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/modelling_methodology.tex	(revision 4133)
+++ /anuga_work/production/onslow_2006/report/modelling_methodology.tex	(revision 4134)
@@ -29,13 +29,10 @@
 evenly spaced along the Sunda Arc subduction zone. 
 The assessment used the Method of Splitting Tsunamis (MOST)
-\cite{VT:MOST} model. 
-%The maximal wave height at a fixed contour line near the coastline
-%(e.g. 50m) is then reported as the hazard to communities ashore. 
-%Models such as Method of Splitting Tsunamis (MOST) \cite{VT:MOST} and the 
-%URS Corporation's
-%Probabilistic Tsunami Hazard Analysis  
-%\cite{somerville:urs} follow this paradigm.
+\cite{VT:MOST} model. A detailed probabilistic hazard map has now been completed
+for the WA coastline \cite{prob:fesa} which is based on the paradigm used by the URS corporation's 
+Probabilistic Tsunami Hazard Analysis \cite{somerville:urs}.
 
-While MOST is suitable for generating and propagating the tsunami wave from its source, it is not adequate to
+While MOST and URS are suitable for generating and propagating the tsunami wave from its source, 
+they are not adequate to
 model the wave's impact on communities ashore.  
 To capture the \emph{impact} of a tsunami to a coastal community, 
@@ -52,87 +49,67 @@
 details of the wave and its interactions, a much finer resolution is 
 required than that of the hazard model. As a result, ANUGA simulations concentrate 
-on specific coastal communities. MOST by contrast uses a 
-coarser resolution and covers often vast areas. To develop the impact 
+on specific coastal communities. MOST and URS by contrast use a 
+coarser resolution and cover often vast areas. To develop the impact 
 from an earthquake event from a distant source, we adopt a hybrid approach of 
-modelling the event itself with MOST and modelling the impact with ANUGA. 
-In this way, the output from MOST serves as an input to ANUGA. 
-In modelling terms, the MOST output is a boundary condition for ANUGA.
+modelling the event itself with MOST or URS and modelling the impact with ANUGA. 
+In this way, the output from MOST or URS serve as an input to ANUGA. 
+In modelling terms, the MOST or URS output is a boundary condition for ANUGA.
+Further details
+regarding the inundation modelling requirements for this study can be found in 
+Appendix \ref{anugasetup}.
 
-\bigskip %FIXME (Ole): Should this be a subsection even? 
-The risk of a given tsunami scenario cannot be determined until the 
-likelihood of the tsunami is known. GA is currently building a 
-complete probabilistic hazard map which is due for completion 
-in late 2006. We therefore report on the impact of a single 
-tsunami event only. When the hazard map is completed, the impact 
-will be assessed for a range of events which will ultimately 
-determine a tsunami risk assessment for the NW shelf.
-
-FESA is interested in the ``most frequent worst case scenario''. Whilst
-we currently cannot determine exactly what that event may be, the
-preliminary hazard assessment suggested that the maximum
-magnitude of earthquakes off Java was considered to be at 
-least 8.5 and could potentially be as high as 9. Therefore, 
-the Mw 9 event 
-provides a plausible worst case scenario and is used as the tsunami
-source in this report.
-Figure \ref{fig:mw9} shows the maximum wave height of a tsunami initiated 
-by a Mw 9 event off
-the coast of Java. This event provides the source and
-boundary condition to the
-inundation model presented in Section \ref{sec:anuga}. 
+The risk of a given tsunami scenario can only be determined when the likelihood
+of the event is known. The probabilistic hazard map for WA \cite{prob:fesa}
+calculates which events pose the most threat to an identified region. Figure
+\ref{fig:probonslow} describes the probability of each event generated along the Java trench
+impacting Onslow. For example, an event generated at point A end would have a
+smaller chance of impacting Onslow than an event generated at point B.
 
 \begin{figure}[h]
 
-  \centerline{ \includegraphics[width=140mm, height=100mm]
-{../report_figures/mw9.jpg}}
+%\centerline{ \includegraphics[width=140mm, height=100mm]{../report_figures/}}
 
-  \caption{Maximum wave height (in cms) for a Mw 9 event off the 
-coast of Java}
-  \label{fig:mw9}
+  \caption{}
+  \label{fig:probonslow}
 \end{figure}
 
-%To model the
-%details of tsunami inundation of a community one must therefore capture %what is 
-%known as non-linear effects and use a much higher resolution for the 
-%elevation data. 
-%Linear models typically use data resolutions of the order
-%of hundreds of metres, which is sufficient to model the tsunami waves
-%in deeper water where the wavelength is longer. 
-%Non-linear models however require much finer resolution in order to %capture
-%the complexity associated with the water flow from offshore 
-%to onshore. By contrast, the data
-%resolution required is typically of the order of tens of metres.
-%The model ANUGA \cite{ON:modsim} is suitable for this type of non-linear
-%modelling.
-%Using a non-linear model capable of resolving local bathymetric effects 
-%and runup using detailed elevation data will require more computational 
-%resources than the typical hazard model making it infeasible to use it 
-%for the entire, end-to-end, modelling. 
+To prepare a tsunami risk assessment, a number of events will be chosen
+for a range of probabilities (or return periods). As Figure \ref{fig:probonslow}
+shows, for a given probability, a number of events are possible. The resulting
+impact to Onslow would then vary depending on the source of the event.  
 
-%We have adopted a hybrid approach whereby the output from the  
-%hazard model MOST is used as input to ANUGA at the seaward boundary of its %study area. 
-%In other words, the output of MOST serves as boundary condition for the 
-%ANUGA model. In this way, we restrict the computationally intensive part %only to 
-%regions where we are interested in the detailed inundation process.  
+% used for the 2005 report when looking at one event
+%\bigskip %FIXME (Ole): Should this be a subsection even? 
+%The risk of a given tsunami scenario cannot be determined until the 
+%likelihood of the tsunami is known. GA is currently building a 
+%complete probabilistic hazard map which is due for completion 
+%in late 2006. We therefore report on the impact of a single 
+%tsunami event only. When the hazard map is completed, the impact 
+%will be assessed for a range of events which will ultimately 
+%determine a tsunami risk assessment for the NW shelf.
 
-%Furthermore, to avoid unnecessary computations ANUGA works with an 
-%unstructured triangular mesh rather than the rectangular grids 
-%used by e.g.\ MOST. The advantage of an unstructured mesh 
-%is that different regions can have different resolutions allowing 
-%computational resources to be directed where they are most needed. 
-%For example, one might use very high resolution near a community 
-%or in an estuary, whereas a coarser resolution may be sufficient 
-%in deeper water where the bathymetric effects are less pronounced. 
-%Figure \ref{fig:refinedmesh} shows a mesh of variable resolution.
+%FESA is interested in the ``most frequent worst case scenario''. Whilst
+%we currently cannot determine exactly what that event may be, the
+%preliminary hazard assessment suggested that the maximum
+%magnitude of earthquakes off Java was considered to be at 
+%least 8.5 and could potentially be as high as 9. Therefore, 
+%the Mw 9 event 
+%provides a plausible worst case scenario and is used as the tsunami
+%source in this report.
+%Figure \ref{fig:mw9} shows the maximum wave height of a tsunami initiated 
+%by a Mw 9 event off
+%the coast of Java. This event provides the source and
+%boundary condition to the
+%inundation model presented in Section \ref{sec:anuga}. 
 
-%\begin{figure}[hbt]
-%
-%  \centerline{ \includegraphics[width=100mm, height=75mm]
-%             {../report_figures/refined_mesh.jpg}}
-%
-%  \caption{Unstructured mesh with variable resolution.}
-%  \label{fig:refinedmesh}
+%\begin{figure}[h]
+
+%  \centerline{ \includegraphics[width=140mm, height=100mm]
+%{../report_figures/mw9.jpg}}
+
+%  \caption{Maximum wave height (in cms) for a Mw 9 event off the 
+%coast of Java}
+%  \label{fig:mw9}
 %\end{figure}
-   
 
 
Index: /anuga_work/production/onslow_2006/report/onslow_2006_report.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/onslow_2006_report.tex	(revision 4133)
+++ /anuga_work/production/onslow_2006/report/onslow_2006_report.tex	(revision 4134)
@@ -63,16 +63,11 @@
     \input{modelling_methodology}
     
-%  \section{Tsunami scenarios}
-%    \label{sec:tsunamiscenario}
-%    \input{tsunami_scenario}
+  \section{Tsunami scenarios}
+    \label{sec:tsunamiscenario}
+    \input{tsunami_scenario}
 
   \section{Data sources}
     \label{sec:data}
     \input{data}
-    
-   \section{Inundation model}
-    \label{sec:anuga}
-    \input{anuga}
-    \input{computational_setup}
         
   \section{Inundation modelling results}
@@ -83,5 +78,5 @@
 \label{table:locations}
 \begin{tabular}{|l|l|l|l|}\hline
-\bf{Point Name} & \bf{Easting} & \bf{Northing} & \bf{Elevation}\\ \hline
+\bf{Point Name} & \bf{Easting} & \bf{Northing} & \bf{Elevation (m)}\\ \hline
 Beadon Point Loading Berth & 302986.51 & 7607334.65 & -8.70 \\ \hline 
 Hospital & 304973.04 & 7605500.42 & 6.02 \\ \hline 
@@ -103,5 +98,5 @@
  
 \begin{figure}[hbt] 
- \centerline{ \includegraphics[width=150mm, height=100mm]{../report_figures/onslow_dli_gauge.jpg}}
+ \centerline{ \includegraphics[width=\paperwidth]{../report_figures/onslow_dli_gauge.jpg}}
 \caption{Point locations used for Onslow study.}  
 \label{fig:points}
@@ -122,7 +117,7 @@
      \input{damage}
 
-%   \section{Impact due to data accuracy}
-%     \input{discussion}
-%     \label{sec:issues}
+  % \section{Impact due to data accuracy}
+  %   \input{discussion}
+  %   \label{sec:issues}
 
      \section{Summary}
@@ -136,13 +131,20 @@
 
     \appendix
-    
+
+    \section{ANUGA modelling parameters}
+    \label{sec:anugasetup}
+    \input{anuga_setup}
+
+\clearpage
+
    \section{Metadata}
      \label{sec:metadata}
      \input{metadata}
 
-\pagebreak
+\clearpage
 
    \section{Time series}
      \label{sec:timeseries}
+     \input{timeseriesdiscussion}
 \input{latexoutput} 
  \clearpage 
Index: /anuga_work/production/onslow_2006/report/references.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/references.tex	(revision 4133)
+++ /anuga_work/production/onslow_2006/report/references.tex	(revision 4134)
@@ -12,7 +12,11 @@
 Tsunami (MOST) model, NOAA Technical Memorandum ERL PMEL-112.
 
-%\bibitem{somerville:urs} Somerville, P., Thio, H.K. and Ichinose, G. (2005)
-%Probabilistic Tsunami Hazard Analysis. Report delivered to Geoscience
-%Australia 2005.
+\bibitem{prob:fesa} Burbidge, D. and Cummins, P. (2006) Probabilistic
+Tsunami Hazard Assessment of Western Australia. Report to the 
+Fire and Emergency Services Authority of Western Australia.
+
+\bibitem{somerville:urs} Somerville, P., Thio, H.K. and Ichinose, G. (2005)
+Probabilistic Tsunami Hazard Analysis. Report delivered to Geoscience
+Australia 2005.
 
 \bibitem{matsuyama:1999}
@@ -46,3 +50,5 @@
 HAZUS-MH User Manual, Washington DC, USA.
 
+\bibitem{uq:friction} Duncan - do you have a reference for this?
+
 \end{thebibliography}
Index: /anuga_work/production/onslow_2006/report/timeseriesdiscussion.tex
===================================================================
--- /anuga_work/production/onslow_2006/report/timeseriesdiscussion.tex	(revision 4134)
+++ /anuga_work/production/onslow_2006/report/timeseriesdiscussion.tex	(revision 4134)
@@ -0,0 +1,99 @@
+It is assumed that the earthquake is
+generated at the beginning of the simulation, i.e. time = 0 minutes.
+Stage is defined as the absolute
+water level (in metres) relative to AHD 
+\footnote{For an offshore location such as Beadon Bay West,
+the initial water level will be that of the tidal scenario. In the 
+case of MSL, this water level will be 0. As the tsunami wave moves
+through this point, the water height may grow and thus the stage will 
+represent the amplitude of the wave. For an onshore location such as the 
+Light Tower, the actual water depth will be the difference between
+the stage and the elevation at that point. Therefore, at the beginning
+of the simulation, there will be no water onshore and therefore
+the stage and the elevation will be identical.}. Both stage and speed 
+(in metres/second) for
+each scenario (HAT, MSL and LAT) are shown
+on consistent scales to allow comparison between point locations.
+As a useful benchmark, Table \ref{table:speedexamples}
+describes typical examples for a range of speeds found in the 
+simulations. 
+
+\begin{table}[h]
+\label{table:speedexamples}
+\caption{Examples of a range of velocities.}
+\begin{center}
+\begin{tabular}{|l|l|}\hline
+{\bf Velocity (m/s)} & {\bf Example} \\ \hline
+1 & leisurely stroll pace\\ \hline
+1.5 & average walking pace \\ \hline
+%2 & 100m Olympic male freestyle \\ \hline
+%3 & mackeral \\ \hline
+4 & average person can maintain running for 1000m \\ \hline
+%5 & blue whale \\ \hline
+10 & 100m Olympic male sprinter \\ \hline
+16 & car travelling in urban zones (60 km/hr) \\ \hline
+\end{tabular}
+\end{center}
+\end{table}
+
+A tsunami wave typically has a small amplitude and typically travels at 100's of kilometres per hour. 
+The low amplitude complicates the ability to detect
+the wave. As the water depth decreases, 
+the speed of the wave 
+decreases and the amplitude grows. Another important feature of tsunamis 
+is drawdown. This means that the water is seen to retreat from the beaches
+before a tsunami wave 
+impacts that location. Other features 
+include reflections (where the wave is redirected due to the 
+influence 
+of the coast) and shoaling (where the wave's amplitude is amplified
+close to the coast due to wave interactions).
+These features are seen in these scenarios, and are consistent
+for HAT, MSL and LAT. 
+There is a small wave, followed
+by a large drawdown and then a large secondary wave. 
+
+These 
+features are illustrated in Figure \ref{fig:gaugeBeadonBayeast}
+where a small wave can be seen at around 200 mins. For the HAT
+case (shown in blue), the amplitude
+of the wave at this location is around 0.8 m\footnote{In this
+scenario, the initial water level is 1.5 m, which means that
+the actual amplitude is the difference between the stage value
+and the initial water level; 2.3 - 1.5}.
+The drawdown of around 4.3 m (i.e. 2.3 - -2) then occurs at around 230 mins
+(i.e. 3.8 hours after the event has been generated), before
+the second wave arrives 
+with an amplitude of around 3.6 m (i.e. 4.1 - 1.5). A further wave
+is then evident a short time later (around 255 mins) 
+which further increases the amplitude to around 5 m (i.e. 6.6 - 1.5).
+These features are replicated at each of the offshore points (those
+points with negative elevation as shown in Table \ref{table:locations}).
+
+The wave amplitude is typically greater
+for those locations which are in the shallowest water. For example,
+the maximum wave amplitude at the Beadon Bay East location 
+(Figure \ref{fig:gaugeBeadonBayeast}) is over
+4.5m where the water depth would normally be 3.56 m. In the 
+Beadon Bay West location (Figure \ref{fig:gaugeBeadonBaywest})
+where the water depth would normally be 4.62 m, 
+the maximum wave amplitude is much less (around 3 m). The wave amplitude
+at the West of Groyne location (Figure \ref{fig:gaugeWestofGroyne})
+is not greater than that seen
+at the Beadon Bay East location, even though the water depth is 
+much less, at 2.11m. This is probably due to its proximity
+to the groyne\footnote{A groyne is a man made structure to combat
+coastal erosion.}
+which has impeded the tsunami wave to some degree. However, the 
+maximum speed found amongst the locations is at the West of Groyne
+point which is in the shallowest water.
+
+The speed of the tsunami sharply increases as it moves onshore. There
+is minimal inundation found at the locations chosen, with the Bindi Bindi
+community receiving the greatest inundation for all tidal scenarios.
+At HAT, the community would receive over 1 m of inundation with 
+the water moving through the community at approximately 16 m/s. Referring
+to Table \ref{table:speedexamples}, a person in this location could
+not outrun this water movement. A small amount of water is found
+at the hospital (10 cm). Whilst this seems minimal, the water is moving 
+at around 6 m/s which could dislodge some items if the water was able to enter the hospital.
Index: /anuga_work/production/pt_hedland_2006/make_report.py
===================================================================
--- /anuga_work/production/pt_hedland_2006/make_report.py	(revision 4133)
+++ /anuga_work/production/pt_hedland_2006/make_report.py	(revision 4134)
@@ -211,5 +211,5 @@
 \label{table:locations}
 \\begin{tabular}{|l|l|l|l|}\hline
-\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation}\\\\ \hline
+\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline
 """
 fid.write(s)
@@ -271,10 +271,16 @@
 
     \\appendix
-    
+
+    \section{ANUGA modelling parameters}
+    \label{sec:anugasetup}
+    \input{anuga_setup}
+
+\clearpage
+
    \section{Metadata}
      \label{sec:metadata}
      \input{metadata}
 
-\pagebreak
+\clearpage
 
    \section{Time series}
Index: /anuga_work/production/sydney_2006/report/anuga_setup.tex
===================================================================
--- /anuga_work/production/sydney_2006/report/anuga_setup.tex	(revision 4133)
+++ /anuga_work/production/sydney_2006/report/anuga_setup.tex	(revision 4134)
@@ -48,6 +48,6 @@
 Region 1: Surrounding major population centre of Wollongong with a cell area of 500 m$^2$ (lateral accuracy 30 m).
 Region 2: Surrounding Lake Illawarra (south of Wollongong) with a cell area of 500 m$^2$ (lateral accuracy 30 m).
-Region 3: Surrounds the coastal regions with a cell area of 25000 m$^2$ (lateral accuracy  m).
-Region 4: The remaining area is given a cell area of 1000000 m$^2$ (lateral accuracy  m).
+Region 3: Surrounds the coastal regions with a cell area of 25000 m$^2$ (lateral accuracy 220 m).
+Region 4: The remaining area is given a cell area of 1000000 m$^2$ (lateral accuracy 1400 m).
 }
   \label{fig:regionA}
@@ -58,8 +58,8 @@
   \centerline{ \includegraphics[scale=0.5]{../report_figures/regionBmodel.jpg}}
 
-  \caption{Study area for region B highlighting four regions of increased refinement.
+  \caption{Study area for region B highlighting three regions of increased refinement.
 Region 1: Surrounding the entrance to Sydney harbour, Northern Beaches and Botany Bay with a cell area of 500 m$^2$ (lateral accuracy 30 m).
-Region 2: Surrounding the coastal regions with a cell area of 50000 m$^2$ (lateral accuracy m).
-Region 3: The remaining area is given a cell area of 250000 m$^2$ (lateral accuracy m).
+Region 2: Surrounding the coastal regions with a cell area of 50000 m$^2$ (lateral accuracy 315 m).
+Region 3: The remaining area is given a cell area of 250000 m$^2$ (lateral accuracy 700 m).
 }
   \label{fig:regionB}
@@ -70,8 +70,8 @@
   \centerline{ \includegraphics[scale=0.5]{../report_figures/regionCmodel.jpg}}
 
-  \caption{Study area for region C highlighting four regions of increased refinement.
-Region 1: Surrounding the major population centre of Newcastle with a cell area of 1000 m$^2$ (lateral accuracy m).
-Region 2: Surrounding the coastal regions with a cell area of 50000 m$^2$ (lateral accuracy  m).
-Region 3: The remaining area is given a cell area of 500000 m$^2$ (lateral accuracy  m).
+  \caption{Study area for region C highlighting three regions of increased refinement.
+Region 1: Surrounding the major population centre of Newcastle with a cell area of 1000 m$^2$ (lateral accuracy 45 m).
+Region 2: Surrounding the coastal regions with a cell area of 50000 m$^2$ (lateral accuracy  315 m).
+Region 3: The remaining area is given a cell area of 500000 m$^2$ (lateral accuracy  1000 m).
 }
   \label{fig:regionC}
@@ -110,2 +110,42 @@
 \end{center}
 \end{table}
+
+
+\begin{table}[h]
+\begin{center}
+
+\caption{Output directories for model simulations.}
+\begin{tabular}{|l|l|l|}\hline
+{\bf Region} & {\bf Slide}  & {\bf Location:}
+$\backslash {\rm inundation \backslash data \backslash new\_south\_wales
+\backslash }$ \\ \hline
+A & Bulli &  $ {\rm wollongong \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash 20061211\_060105}$ \\ \hline
+A & Shovel & ${\rm wollongong \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash
+20061212\_012715}$ \\ \hline
+A & Yacaaba & ${\rm wollongong \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash
+20061212\_064735}$ \\ \hline
+B & Bulli &  ${\rm sydney \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash
+20061211\_071516}$ \\ \hline
+B & Shovel & ${\rm sydney \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash
+20061212\_012705}$ \\ \hline
+B & Yacaaba & $\rm{ sydney \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash
+20061212\_064807}$ \\ \hline
+C & Bulli &  ${\rm newcastle \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash
+20061211\_073709}$ \\ \hline
+C & Shovel &  ${\rm newcastle \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash
+20061212\_012802}$ \\ \hline
+C & Yacaaba & ${\rm newcastle \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash
+20061212\_064757}$ \\ \hline
+\end{tabular}
+
+\end{center}
+\end{table}
Index: /anuga_work/production/sydney_2006/report/interpretation.tex
===================================================================
--- /anuga_work/production/sydney_2006/report/interpretation.tex	(revision 4133)
+++ /anuga_work/production/sydney_2006/report/interpretation.tex	(revision 4134)
@@ -5,7 +5,7 @@
 each model with respect to each other and highlights the location
 of the slides modelled in this study. See Appendix \ref{sec:anugasetup}
-showing the detail for each study region.
+for the detailed inundation modelling setup for each study region.
 
-\begin{figure}
+\begin{figure}[h]
 \centerline{\includegraphics[scale=0.4]{../report_figures/overallmodel.jpg}}
 \caption{NSW region showing the location of each study region and the location
@@ -26,5 +26,5 @@
 
 The time to impact after the initiation of the slide is approximately 25 mins 
-for region A, 25 mins for region B and mins for region C. See Appendix \ref{sec:timeseries}
+for region A and B and approximately 35 mins for region C. See Appendix \ref{sec:timeseries}
 for time series outputs at Fairy Meadow (for region A), Manly Beach (for region B) and 
 Stockton Beach (for region C).
@@ -68,6 +68,6 @@
 resultant amplitude are approximately equal.
 
-\begin{figure}
-\centerline{\includegraphics[scale=0.4]{../report_figures/depthvsampdataslide.png}}
+\begin{figure}[h]
+\centerline{\includegraphics[scale=0.25]{../report_figures/depthvsampdataslide.png}}
 \caption{Relationship between characteristic 3D amplitude and depth
 for each slide volume.}
Index: /anuga_work/production/sydney_2006/report/latexoutput20061211071516.tex
===================================================================
--- /anuga_work/production/sydney_2006/report/latexoutput20061211071516.tex	(revision 4133)
+++ /anuga_work/production/sydney_2006/report/latexoutput20061211071516.tex	(revision 4134)
@@ -9,2 +9,12 @@
  \end{figure} 
  
+\begin{figure}[hbt] 
+ \centering 
+ \begin{tabular}{cc} 
+\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeOffshoreManlystage20061211071516.png}& 
+\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeOffshoreManlymomentum20061211071516.png}\\ 
+\end{tabular} 
+ \caption{Time series for depth,  and momentum at Offshore Manly location (elevation 0.75m)} 
+ \label{fig:20061211071516gaugeOffshoreManly} 
+ \end{figure} 
+ 
Index: /anuga_work/production/sydney_2006/report/sydney_2006_report.tex
===================================================================
--- /anuga_work/production/sydney_2006/report/sydney_2006_report.tex	(revision 4133)
+++ /anuga_work/production/sydney_2006/report/sydney_2006_report.tex	(revision 4134)
@@ -72,9 +72,24 @@
          
 \input{interpretation} 
-\input{shovel_A_map} 
+\input{bulli_A_map} 
+ \clearpage 
+ \input{shovel_A_map} 
+ \clearpage 
+ \input{yacaaba_A_map} 
+ \clearpage 
+ \input{bulli_B_map} 
+ \clearpage 
+ \input{shovel_B_map} 
+  \clearpage 
+ 
+ \input{yacaaba_B_map} 
+ \clearpage 
+\input{bulli_C_map} 
+ \clearpage 
+\input{shovel_C_map} 
+ \clearpage 
+\input{yacaaba_C_map} 
  \clearpage 
  
-
-
      \section{Summary}
      \label{sec:summary}
@@ -95,4 +110,6 @@
 \label{sec:anugasetup}
 \input{anuga_setup}
+
+\pagebreak
 
 \section{Submarine Mass Failure Model}
