Index: /anuga_work/publications/boxing_day_validation_2008/patong_validation.tex
===================================================================
--- /anuga_work/publications/boxing_day_validation_2008/patong_validation.tex	(revision 7209)
+++ /anuga_work/publications/boxing_day_validation_2008/patong_validation.tex	(revision 7210)
@@ -1,3 +1,8 @@
 \documentclass[a4paper]{article}
+% to be added when submitted to ocean dynamics
+%\documentclass[smallcondensed,draft]{svjour3}
+%\usepackage{mathptmx}
+%\journalname{Ocean Dynamics}
+
 \usepackage{graphicx}
 \usepackage{hyperref}
@@ -6,18 +11,19 @@
 \newcommand{\doi}[1]{\url{http://dx.doi.org/#1}}
 
+
 %----------title-------------%
-%\title{Inundation Modelling of the December 2004 Indian Ocean Tsunami}
 \title{Benchmarking Tsunami Models using the December 2004 Indian Ocean Tsunami and its Impact at Patong Beach}
 
 %-------authors-----------
-\author{J.~D. Jakeman\thanks{The Australian National University, Canberra, \textsc{Australia}.
-\protect\url{mailto:john.jakeman@anu.edu.au}}
-\and O.Nielsen\thanks{Geoscience Australia, Canberra, \textsc{Australia}}
-\and R. Mleczko\footnotemark[2]
-\and D. Burbidge\footnotemark[2]
-\and K. VanPutten\footnotemark[2]
-\and N. Horspool\footnotemark[2]
-}
-
+\author{J.~D. Jakeman \and O. Nielsen \and R. Mleczko \and D. Burbidge \and K. VanPutten \and N. Horspool}
+
+% to be added when submitted to ocean dynamics
+%\institute{J.~D. Jakeman \at 
+%	The Australian National University, Canberra, \textsc{Australia}\
+%	\email{john.jakeman@anu.edu.au}
+%	\and
+%	O. Nielsen \and R. Mleczko \and D. Burbidge \and K. VanPutten \and N. Horspool \at 
+%	Geoscience Australia, Canberra, \textsc{Australia}
+%}
 %================Start of Document================
 \begin{document}
@@ -25,5 +31,7 @@
 %------Abstract--------------
 \begin{abstract}
-In this paper a new benchmark for tsunami model validation is proposed. The benchmark is based upon the 2004 Indian Ocean tsunami, which provides a uniquely large amount of observational data for model comparison. Unlike the small number of existing benchmarks, the proposed test validates all three stages of tsunami evolution - generation, propagation and inundation. Specifically we use geodetic measurements of the Sumatra--Andaman earthquake to validate the tsunami source, altimetry data from the \textsc{jason} satellite to test open ocean propagation, eye-witness accounts to assess near shore propagation and a detailed inundation survey of Patong Bay, Thailand to compare model and observed inundation. Furthermore we utilise this benchmark to further validate the hydrodynamic modelling tool \textsc{anuga}  which is used to simulate the tsunami inundation and run rain-induced floods.
+In this paper a new benchmark for tsunami model validation is proposed. The benchmark is based upon the 2004 Indian Ocean tsunami, which provides a uniquely large amount of observational data for model comparison. Unlike the small number of existing benchmarks, the proposed test validates all three stages of tsunami evolution - generation, propagation and inundation. Specifically we use geodetic measurements of the Sumatra--Andaman earthquake to validate the tsunami source, altimetry data from the \textsc{jason} satellite to test open ocean propagation, eye-witness accounts to assess near shore propagation and a detailed inundation survey of Patong Bay, Thailand to compare model and observed inundation. Furthermore we utilise this benchmark to further validate the hydrodynamic modelling tool \textsc{anuga}  which is used to simulate the tsunami inundation and run rain-induced floods. Important buildings and other structures were incorporated into the underlying computational mesh and shown to have a large influence of inundation extent. Sensitivty analysis also showed that the model predictions are comparatively insensitive to large changes in friction and small perturbations in wave weight at the 100m depth contour.
+% to be added when submitted to ocean dynamics
+%\keywords{Tsunami \and modelling \and validation and verification \and benchmark}
 \end{abstract}
 
@@ -32,11 +40,13 @@
 
 \section{Introduction}
-Tsunami are a potential hazard to coastal communities all over the world. A number of recent large events have increased community and scientific awareness of the need for effective tsunami hazard mitigation. Tsunami modelling is major component of hazard mitigation, which involves detection, forecasting, and emergency preparedness~\cite{synolakis05}. Accurate models can be used to provide information that increases the effectiveness of action undertaken before the event to minimise damage (early warning systems, breakwalls etc.) and protocols put in place to be followed when the flood waters subside.
-
-Several approaches are currently used to model tsunami propagation and inundation. These methods differ in both the formulation used to describe the evolution of the tsunami and the numerical methods used to solve the governing equations. The shallow water wave equations, linearised shallow water wave equations, and Boussinesq-type equations are frequently used to simulate tsunami propagation. The nonlinear nature of these equations, the highly variable nature of the phenomena that they describe and the complex reality of the geometry they operate in, necessitate the use of numerical models. These models are typically used to predict quantities such as arrival times, wave speeds and heights and inundation extents which are used to develop efficient hazard mitigation plans. Inaccuracies in model prediction can result in inappropriate evacuation plans and town zoning which may result in loss of life and large financial losses. Consequently tsunami models must undergo sufficient testing to increase scientific and community confidence in the model predictions.
-
-Complete confidence in a model of a physical system frequently in general cannot be established.  One can only hope to state under what conditions the model hypothesis holds true. Specifically the utility of a model can be assessed through a process of verification and validation. Verification assesses the accuracy of the numerical method used to solve the governing equations and validation is used to investigate whether the model adequately represents the physical system~\cite{bates01}. Together these processes can be used to establish the likelihood that that a model is a legitimate hypothesis.
-
-The sources of data used to validate and verify a model can be separated into three main categories; analytical solutions, scale experiments and field measurements. Analytical solutions of the governing equations of a model, if available, provide the best means of verifying any numerical model. However analytical solutions to the governing equations used here are frequently limited to a small set of idealised examples that do not completely capture the more complex behaviour of `real' events. For tsunami inundation, scale experiments, typically in the form of wave-tank experiments, provide a much more realistic source of data that better captures the complex dynamics of natural tsunami, whilst allowing control of the event and much easier and accurate measurement of the tsunami properties. However comparison of numerical predictions with field data provides the most stringent test of model veracity. The use of field data increases the generality and significance of conclusions made regarding model utility. However the use of field data also significantly increases the uncertainty of the validation experiment that may constrain the ability to make unequivocal statements~\cite{bates01}.
+Tsunami are a potential hazard to coastal communities all over the world. A number of recent large events have increased community and scientific awareness of the need for effective detection, forecasting, and emergency preparedness. Probabilistic, geological, hydrodynamic, and economic models are required to predict the location and likelihood of an event, the initial sea floor deformation and subsequent propagation and inundation of the tsunami, the effectiveness of hazard mitigation procedures and the economic impact of such measures and the event itself. Here we focus on modelling of the physical processes. For discussion on economic and decision based models refer to~\cite{} and the references therein.
+
+Various approaches are currently used to assess the potential impact of tsunami. These methods differ in both the formulation used to describe the evolution of the tsunami and the numerical methods used to solve the governing equations. However any legitimate model must address each of the three distinct stages of tsunami evolution--- generation, inundation and propagation. Geological models must be used to provide estimates of initial sea floor and ocean surface deformation. The complexity of these models range from empirical to non-linear three-dimensional mechanical models. The shallow water wave equations, linearised shallow water wave equations, and Boussinesq-type equations are frequently used to simulate tsunami propagation. These models are typically used to predict quantities such as arrival times, wave speeds and heights, and inundation extents which are used to develop efficient hazard mitigation plans. 
+
+Inaccuracies in model prediction can result in inappropriate evacuation plans and town zoning which may result in loss of life and large financial losses. Consequently tsunami models must undergo sufficient end-to-end testing to increase scientific and community confidence in the model predictions.
+
+Complete confidence in a model of a physical system cannot be established.  One can only hope to state under what conditions the model hypothesis holds true. Specifically the utility of a model can be assessed through a process of verification and validation. Verification assesses the accuracy of the numerical method used to solve the governing equations and validation is used to investigate whether the model adequately represents the physical system~\cite{bates01}. Together these processes can be used to establish the likelihood that that a model is a legitimate hypothesis.
+
+The sources of data used to validate and verify a model can be separated into three main categories; analytical solutions, scale experiments and field measurements. Analytical solutions of the governing equations of a model, if available, provide the best means of verifying any numerical model. However analytical solutions are frequently limited to a small set of idealised examples that do not completely capture the more complex behaviour of `real' events. Scale experiments, typically in the form of wave-tank experiments, provide a much more realistic source of data that better captures the complex dynamics of natural tsunami, whilst allowing control of the event and much easier and accurate measurement of the tsunami properties. Comparison of numerical predictions with field data provides the most stringent test. The use of field data increases the generality and significance of conclusions made regarding model utility. However it must be noted that the use of field data also significantly increases the uncertainty of the validation experiment that may constrain the ability to make unequivocal statements~\cite{bates01}.
 
 Currently the extent of tsunami related field data is limited. The cost of tsunami monitoring programs, bathymetry and topography surveys prohibits the collection of data in many of the regions in which tsunamis pose greatest threat. The resulting lack of data has limited the number of field data sets available to validate tsunami models. Synolakis et. al~\cite{synolakis07} have developed a set of standards, criteria and procedures for evaluating numerical models of tsunami. They propose three analytical solutions to help identify the validity of a model and  five scale comparisons (wave-tank benchmarks) and two field events to assess model veracity. 
@@ -44,7 +54,9 @@
 The first field data benchmark introduced by Synolakis compares model results against observed data from the Hokkaido-Nansei-Oki tsunami that occurred around Okushiri Island, Japan on the 12th of July 1993. This tsunami provides an example of extreme runup generated from reflections and constructive interference resulting from local topography and bathymetry. The benchmark consists of two tide gauge records and numerous spatially distributed point sites at which modelled maximum runup elevations can be compared. The second benchmark is based upon the Rat Islands Tsunami that occurred off the coast of Alaska on the 17th of November 2003. The Rat island tsunami provides a good test for real-time forecasting models since tsunami was recorded at three tsunameters. The test requires matching the propagation model data with the DART recording to constrain the tsunami source model and using a propagation model to to reproduce the tide gauge record at Hilo.
 
-The evolution of earthquake-generated tsunamis has three distinctive stages: generation, propagation and run-up~\cite{titov97a}. To accurately model the evolution of a tsunami all three stages must be dealt with. In this section we present the data necessary to implement the proposed benchmark corresponding to each of the three stages of the tsunami's evolution.  In this paper we develop a field data benchmark to be used in conjunction with the other tests proposed by Synolakis et al.~\cite{synolakis07} to validate and verify tsunami models. Unlike the aforementioned tests, the proposed benchmark allows one to estimate the error in a model prediction for all three distinctive stages of the evolution of a tsunami: generation, propagation and run-up. The benchmark comprises of geodetic measurements of the Sumatra--Andaman earthquake to validate the tsunami source, altimetry data from the JASON satellite to test open ocean propagation, eye-witness accounts to assess near shore propagation and a detailed inundation survey of Patong Bay, Thailand to compare model and observed inundation. A description of this data is give in Section~\ref{sec:data}.
-
-An associated aim of this paper is to illustrate the use of this new benchmark to validate an operational tsunami model called \textsc{anuga}. A description of \textsc{anuga} is given in Secion~\ref{sec:models} and the validation results are given in Secion~\ref{sec:results}. 
+In this paper we develop a field data benchmark to be used in conjunction with the other tests proposed by Synolakis et al.~\cite{synolakis07} to validate and verify tsunami models. Unlike the aforementioned tests, the proposed benchmark allows evaluation of model structure during all three distinctive stages of the evolution of a tsunami. The benchmark consists of geodetic measurements of the Sumatra--Andaman earthquake which are used to validate the description of the tsunami source, altimetry data from the JASON satellite to test open ocean propagation, eye-witness accounts to assess near shore propagation and a detailed inundation survey of Patong Bay, Thailand to compare model and observed inundation. A description of the data required to construct the benchmark is give in Section~\ref{sec:data}.
+
+An associated aim of this paper is to illustrate the use of this new benchmark to validate an operational tsunami model called \textsc{anuga} used by Geoscience Australia. A description of \textsc{anuga} is given in Secion~\ref{sec:models} and the validation results are given in Secion~\ref{sec:results}. 
+
+The numerical models used to model tsunami are extremely computationally intensive. Full resolution models of the entire evolution process will often take a number of days to run. Consequently the uncertainty in model predictions is difficult to quantify. However model uncertainty should not be ignored. Section ~\ref{sec:sensitivity} provides a simple sensitivity analysis that can be used to investigate the sensitivity of model predictions to model parameters. 
 
 %================Section===========================
@@ -60,7 +72,7 @@
 All tsunami are generated from an initial disturbance of the ocean which develops into a low frequency wave that propagates outwards from the source. The initial deformation of the water surface is most commonly caused by coseismic displacement of the sea floor, but submarine mass failures, landslides, volcanoes or asteroids can also cause tsunami. In this section we detail the information we used in this study to validate models of the sea floor deformation generated by the 2004 Sumatra--Andaman earthquake.
 
-The 2004 Sumatra--Andaman tsunami was generated by severe coseismic displacement of the sea floor as a result of one of the largest earthquakes on record. The mega-thrust earthquake started on the 26 December 2004 at 0h58'53'' UTC (or just befor 8 am local time) approximately 70 km offshore North Sumatra (\url{http://earthquake.usgs.gov/eqcenter/eqinthenews/2004/usslav}). The rupture propagated 1000-1300 km along the Sumatra-Andaman trench time to the north at a rate of 2.5-3 km.s$^{-1}$ and lasted approximately 8-10 minutes~\cite{ammon05}. Estimates of the moment magnitude of this event range from about 9.1 to 9.3~\cite{chlieh07, stein07}.
-
-The unusually large surface deformation caused by this earthquakes means that there were a range of different geodetic measurements of the surface deformation available. These include field measurements of uplifted or subsided coral heads, continuous or campaign \textsc{GPS} measurments and remote sensing measurements of uplift or subsidence (see ~\cite{chlieh07} and references therein). Here we use the the near field estimates of vertical deformation in northwestern Sumatra and the Nicobar-Andaman islands collated by~\cite{chlieh07} to validate that our crustal deformation model of the 2004 Sumatra--Andaman earthquake is producing reasonable results. Note that the geodetic data used here is a combination of the vertical deformation that happened in the ~10 minutes of the earthquake plus the deformation that followed in the days following the earthquake before each particular measurement was actually made (typically of order days). Therefore some of the observations may not contain the purely co-seismic deformation but could include some post-seismic deformation as well~\cite{chlieh07}.
+The 2004 Sumatra--Andaman tsunami was generated by severe coseismic displacement of the sea floor as a result of one of the largest earthquakes on record. The mega-thrust earthquake started on the 26 December 2004 at 0h58'53'' UTC (or just before 8 am local time) approximately 70 km offshore North Sumatra (\url{http://earthquake.usgs.gov/eqcenter/eqinthenews/2004/usslav}). The rupture propagated 1000-1300 km along the Sumatra-Andaman trench time to the north at a rate of 2.5-3 km.s$^{-1}$ and lasted approximately 8-10 minutes~\cite{ammon05}. Estimates of the moment magnitude of this event range from about 9.1 to 9.3~\cite{chlieh07,stein07}.
+
+The unusually large surface deformation caused by this earthquakes means that there were a range of different geodetic measurements of the surface deformation available. These include field measurements of uplifted or subsided coral heads, continuous or campaign \textsc{GPS} measurements and remote sensing measurements of uplift or subsidence (see ~\cite{chlieh07} and references therein). Here we use the the near field estimates of vertical deformation in northwestern Sumatra and the Nicobar-Andaman islands collated by~\cite{chlieh07} to validate that our crustal deformation model of the 2004 Sumatra--Andaman earthquake is producing reasonable results. Note that the geodetic data used here is a combination of the vertical deformation that happened in the ~10 minutes of the earthquake plus the deformation that followed in the days following the earthquake before each particular measurement was actually made (typically of order days). Therefore some of the observations may not contain the purely co-seismic deformation but could include some post-seismic deformation as well~\cite{chlieh07}.
 
 %DAVID: I commented out the figure since we can combine it with the model result without obscuring it. That will keep the number of figures down.
@@ -93,5 +105,5 @@
 \begin{figure}[ht]
 \begin{center}
-\includegraphics[width=8.0cm,keepaspectratio=true]{nested_grids}
+\includegraphics[width=0.75\textwidth,keepaspectratio=true]{nested_grids}
 \caption{Nested grids of elevation data.}
 \label{fig:nested_grids}
@@ -136,12 +148,30 @@
 
 \subsubsection{Buildings and Other Structures}
-Human made build and structures can significantly effect tsunamni inundation. The location and size and number of floors of the buildings in Patong Bay were extracted from a GIS data set~\cite{}. The heights of these buildings was estimated assuming that each floor was 3m high.
+Human made build and structures can significantly effect tsunami inundation. The location and size and number of floors of the buildings in Patong Bay were extracted from a GIS data set~\cite{}. The heights of these buildings was estimated assuming that each floor was 3m high.
+
+\subsubsection{Inundation Survey}
+Tsunami run-up is often the cause of the largest financial and human losses yet run-up data that can be used to validate model runup predictions is scarce. Of the two field benchmarks proposed by Synolakis only the Okishiri benchmark facilitates comparison between modelled and observed run-up. One of the major strengths of the benchmark proposed here is that modelled runup can be compared to an inundation survey which maps the maximum run-up along an entire coast line rather than at a series of discrete sites. The survey map is shown in Figure~\ref{fig:patongescapemap} and plots the maximum run-up of the 2004 tsunami in Patong bay. Refer to Szczucinski et al~\cite{szczucinski06} for further details.
 
 \subsubsection{Eyewitness Accounts}
 Eyewitness accounts detailed in~\cite{papadopoulos06} report that most people at Patong Beach observed an initial retreat of the shoreline of more than 100m followed a few minutes later by a strong wave (crest). Another less powerful wave arrived another five or ten minutes later. Eyewitness statements place the arrival time of the strong wave between 2 hours and 55 minutes to 3 hours and 5 minutes after the source rupture (09:55am to 10:05am local time). 
 
-
-\subsubsection{Inundation Survey}
-Tsunami run-up is often the cause of the largest financial and human losses yet run-up data that can be used to validate model runup predictions is scarce. Of the two field benchmarks proposed by Synolakis only the Okishiri benchmark facilitates comparison between modelled and observed run-up. One of the major strengths of the benchmark proposed here is that modelled runup can be compared to an inundation survey which maps the maximum run-up along an entire coast line rather than at a series of discrete sites. The survey map is shown in Figure~\ref{fig:patongescapemap} and plots the maximum run-up of the 2004 tsunami in Patong bay. Refer to Szczucinski et al~\cite{szczucinski06} for further details.
+Two videos were sourced from the internet which include footage of the tsunami in Patong Bay on the day of the Indian Ocean Tsunami. Both videos show an
+already inundated town, they then show what is to be assumed as the second and third waves approaching and further flooding the town. The first video
+is in the very north filmed from what is believed to be the roof of the Novotel Hotel. The second video is in the very south
+filmed from a building next door to the Comort Resort. Figure~\ref{fig:video_flow} shows a video taken near the corner of Ruam Chai St. from which we made our estimates.
+
+\begin{figure}[ht]
+\begin{center}
+\includegraphics[width=6.0cm,keepaspectratio=true]{flow_rate_south_0_00sec.jpg}
+\includegraphics[width=6.0cm,keepaspectratio=true]{flow_rate_south_5_04sec.jpg}
+\includegraphics[width=6.0cm,keepaspectratio=true]{flow_rate_south_7_12sec.jpg}
+\includegraphics[width=6.0cm,keepaspectratio=true]{flow_rate_south_7_60sec.jpg}
+\caption{four frames from a video where flow rate could be estimated, circle indicates tracked debris, clockwise from top left: 0.0 sec, 5.0 s, 7.1 s, 7.6 s.}
+\label{fig:video_flow}
+\end{center}
+\end{figure}
+
+Crude flow rates can be estimated with landmarks found in satellite imagery and the use of a GIS and were found to be in
+the range of 5 to 7 metres per second (+/- 2 m/s) in the north and 0.5 to 2 metres per second (+/- 1 m/s) in the south. 
 
 \begin{figure}[ht]
@@ -159,7 +189,7 @@
  \item Reproduce the vertical deformation observed in north-western Sumatra and along the Nicobar--Andaman islands, see Section~\ref{sec:gen_data}.
  \item Reproduce the \textsc{jason} satellite altimetry sea surface anomalies, see Section~\ref{sec:data_jason}.
+ \item Reproduce the inundation survey map in Patong bay (Figure~\ref{fig:patongescapemap}).
  \item Simulate a leading depression followed by two distinct crests of decreasing magnitude.
- \item Predict the arrival time of the first crest should arrive at Patong beach between 2 hours and 55 minutes to 3 hours and 5 minutes after the initial rupture of the source. The subsequent crest arrive five to ten minutes later.
- \item Reproduce the inundation survey map in Patong bay (Figure~\ref{fig:patongescapemap}).
+ \item Predict the water depths and flow speeds, at the locations of the eye-witness videos, that fall within the bounds obtained from the videos.
 \end{itemize}
 
@@ -170,7 +200,7 @@
 \subsection{Generation}\label{sec:modelGeneration}
 
-There are various approaches to modelling the expected crustal deformation from an earthquake at depth. Most approaches model the earthquake as a dislocation in linear, elastic medium. Here we use the method of~\cite{wang03}. One of the main advantages of their method is that it allows the dislocation to be located in a stratified linear elastic half-space with an arbitary number of layers. Other methods (such as those based on Okada's equations) can only model the dislocation in a homogeneous elastic half space, or can only include a limited number of layers, and thus cannot model the effect of the depth dependence of the elasticity of the Earth~\cite{wang03}. The original versions of the codes described here are available from \url{http://www.iamg.org/CGEditor/index.htm}. The first program, \textsc{EDGRN}, calculates elastic Green's function for a set of point sources at a regular set of depths out to a specified distance. The equations controlling the deformation are solved by using a combination of Hankel's transform and Wang et al's implementation of the Thomson-Haskell propagator algorithm~\cite{wang03}. Once the Green's functions are calculated we use a slightly modified version of \textsc{EDCMP} to calculate the sea floor deformation for a specific subfault. This second code discretises the subfault into a set of unit sources and sums the elastic Green's functions calculated from \textsc{EDGRN} for all the unit sources on the fault plane in order to calculate the final static deformation caused by a two dimensional dislocation along the subfault. This step is possible because of the linearity of the governing equations. For this study, we have made minor modifications to \textsc{EDCMP} in order for it to output in a file format compatible with the propagation code in the following section but it is otherwise the similar to the original code.
-
-In order to calculate the crustal deformation using these codes we thus need to have a model describing the variation in elastic properties with depth and a slip model of the earthquake to describe the dislocation. The elastic parameters used for this study are the same as those in Table 2 of \cite{burbidge08}. For the slip model, there are many possible models for the 2004 Andaman--Sumatran earthquake to choose from ~\cite{chlieh07,asavanant08,arcas06,grilli07,ioualalen07}. Some are determined from various geological surveys of the site, others solve an inverse problem which calibrates the source based upon the tsunami wave signal, the seismic signal and/or the runup. The source parameters used here to simulate the 2004 Indian Ocean Tsunami were taken from model G-M9.15 from~\cite{chlieh07}. This model was created by inversion of wide range of geodetic and seismic data. The slip model consists of 686 20km x 20km subsegments each with a different slip, strike and dip angle. The dip subfaults go from $17.5^0$ in the north and $12^0$ in the south. Refer to~\cite{chlieh07} for a detailed discussion of this model and its derivation. Note that the geodetic data used in the validation was also included by~\cite{chlieh07} in the inversion used to find G-M9.15, thus the validation is not completely independent. However, a successful validation would still show that the crustal deformation and elastic properties model used here is at least as valid as the one used in~\cite{chlieh07} and can reproduce the observations just as accurately.
+There are various approaches to modelling the expected crustal deformation from an earthquake at depth. Most approaches model the earthquake as a dislocation in linear, elastic medium. Here we use the method of Wang et. al.~\cite{wang03}. One of the main advantages of their method is that it allows the dislocation to be located in a stratified linear elastic half-space with an arbitrary number of layers. Other methods (such as those based on Okada's equations) can only model the dislocation in a homogeneous elastic half space, or can only include a limited number of layers, and thus cannot model the effect of the depth dependence of the elasticity of the Earth~\cite{wang03}. The original versions of the codes described here are available from \url{http://www.iamg.org/CGEditor/index.htm}. The first program, \textsc{edgrn}, calculates elastic Green's function for a set of point sources at a regular set of depths out to a specified distance. The equations controlling the deformation are solved by using a combination of Hankel's transform and Wang et al's implementation of the Thomson-Haskell propagator algorithm~\cite{wang03}. Once the Green's functions are calculated we use a slightly modified version of \textsc{edcmp} to calculate the sea floor deformation for a specific subfault. This second code discretises the subfault into a set of unit sources and sums the elastic Green's functions calculated from \textsc{edgrn} for all the unit sources on the fault plane in order to calculate the final static deformation caused by a two dimensional dislocation along the subfault. This step is possible because of the linearity of the governing equations. For this study, we have made minor modifications to \textsc{edcmp} in order for it to output in a file format compatible with the propagation code in the following section but it is otherwise the similar to the original code.
+
+In order to calculate the crustal deformation using these codes we thus need to have a model describing the variation in elastic properties with depth and a slip model of the earthquake to describe the dislocation. The elastic parameters used for this study are the same as those in Table 2 of Burbidge~\cite{burbidge08}. For the slip model, there are many possible models for the 2004 Andaman--Sumatran earthquake to choose from ~\cite{chlieh07,asavanant08,arcas06,grilli07,ioualalen07}. Some are determined from various geological surveys of the site, others solve an inverse problem which calibrates the source based upon the tsunami wave signal, the seismic signal and/or the runup. The source parameters used here to simulate the 2004 Indian Ocean tsunami were taken from model G-M9.15 from Chlieh et. al.~\cite{chlieh07}. This model was created by inversion of wide range of geodetic and seismic data. The slip model consists of 686 20km x 20km subsegments each with a different slip, strike and dip angle. The dip subfaults go from $17.5^0$ in the north and $12^0$ in the south. Refer to Chlieh et. al.~\cite{chlieh07} for a detailed discussion of this model and its derivation. Note that the geodetic data used in the validation was also included by~\cite{chlieh07} in the inversion used to find G-M9.15, thus the validation is not completely independent. However, a successful validation would still show that the crustal deformation and elastic properties model used here is at least as valid as the one used by Chlieh et. al.~\cite{chlieh07} and can reproduce the observations just as accurately.
 
 \subsection{Propagation}\label{sec:modelPropagation}
@@ -191,5 +221,5 @@
 
 \subsection{Generation}\label{modelGeneration}
-The location and magnitude of the sea floor displacement associated with the 2004 Sumatra--Andaman tsunami calculated from the G-M9.15 model  of~\cite{chlieh07} is shown in Figure~\ref{fig:surface_deformation}. The magnitude of the sea floor displacement ranges from about $-3.0$ to $5.0$ metres. The region near the fault is predicted to uplift, while that further away from the fault subsides. Also shown in Figure~\ref{fig:surface_deformation} are the areas that were observed to uplift (arrows pointing up) or subside (arrows point down) during and immediately after the earthquake. Most of this data comes uplifted or subsided coral heads. The length of vector increases with the magniutude of the displacement, the length corresponding to 1m of observed motion is shown in the top right corner of the figure. As can be seen, the source model detailed in Section~\ref{sec:modelGeneration} produces a crustal deformation that matches the vertical displacements in the Nicobar-Andaman islands and Sumatra very well. Uplifted regions are close to the fault and subsided regions are further away. The crosses on Figure~\ref{fig:surface_deformation} show estimates of the pivot line from the remote sensing data~\cite{chlieh07} and they follow the predicted pivgot line quite accurately. The average difference between the observed motion and the predicted motion (including the pivot line points) is only 0.06m, well below the typical error of the observations of between 0.25 and 1.0m. However, the occasional point has quite a large error (over 1m), for example a couple uplifted/subsided points appear to be on a wrong side of the predicted pivot line~\ref{fig:surface_deformation}. These points are rare, most fit the predicted vertical motion very well. The excellence of the fit is not surprising, since the original slip model was chosen by~\cite{chlieh07} to fit this (and the seismic data) well. However, this does demonstrate that \textsc{edgrn} and our modified version of \textsc{edstat} can reproduce the correct pattern of vertical deformation very well when the slip distribution is well constrained and when reasonable values for the elastic properties are used.
+The location and magnitude of the sea floor displacement associated with the 2004 Sumatra--Andaman tsunami calculated from the G-M9.15 model  of~\cite{chlieh07} is shown in Figure~\ref{fig:surface_deformation}. The magnitude of the sea floor displacement ranges from about $-3.0$ to $5.0$ metres. The region near the fault is predicted to uplift, while that further away from the fault subsides. Also shown in Figure~\ref{fig:surface_deformation} are the areas that were observed to uplift (arrows pointing up) or subside (arrows point down) during and immediately after the earthquake. Most of this data comes uplifted or subsided coral heads. The length of vector increases with the magnitude of the displacement, the length corresponding to 1m of observed motion is shown in the top right corner of the figure. As can be seen, the source model detailed in Section~\ref{sec:modelGeneration} produces a crustal deformation that matches the vertical displacements in the Nicobar-Andaman islands and Sumatra very well. Uplifted regions are close to the fault and subsided regions are further away. The crosses on Figure~\ref{fig:surface_deformation} show estimates of the pivot line from the remote sensing data~\cite{chlieh07} and they follow the predicted pivot line quite accurately. The average difference between the observed motion and the predicted motion (including the pivot line points) is only 0.06m, well below the typical error of the observations of between 0.25 and 1.0m. However, the occasional point has quite a large error (over 1m), for example a couple uplifted/subsided points appear to be on a wrong side of the predicted pivot line~\ref{fig:surface_deformation}. These points are rare, most fit the predicted vertical motion very well. The excellence of the fit is not surprising, since the original slip model was chosen by~\cite{chlieh07} to fit this (and the seismic data) well. However, this does demonstrate that \textsc{edgrn} and our modified version of \textsc{edstat} can reproduce the correct pattern of vertical deformation very well when the slip distribution is well constrained and when reasonable values for the elastic properties are used.
 
 \begin{figure}[ht]
@@ -203,5 +233,5 @@
 
 \subsection{Propagation}
-The deformation results described in Section~\ref{sec:modelGeneration} was then used to provide an profile of the initial ocean surface displacement. This wave was used as an initial condition for \textsc{ursga} and was propagated the tsunami throughtout the Bay of Bengal. The rectangular computational domain of the largest grid extended from 90$^0$ to 100$^0$East and 0 to 15$^0$North and contained 1335$\times$1996 finite difference points. Inside this grid, a nested sequence of grids was used. The grid resolution of the nested grids went from 27 arc seconds in the coarsest grid, down to 9 arc seconds in the second grid, 3 arc seconds in the thrid grid and finally 1 arc second in the finest grid near Patong. The computational domain is shown in Figure\ref{fig:computational_domain}.
+The deformation results described in Section~\ref{sec:modelGeneration} were used to provide an profile of the initial ocean surface displacement. This wave was used as an initial condition for \textsc{ursga} and was propagated the tsunami throughtout the Bay of Bengal. The rectangular computational domain of the largest grid extended from 90$^0$ to 100$^0$East and 0 to 15$^0$North and contained 1335$\times$1996 finite difference points. Inside this grid, a nested sequence of grids was used. The grid resolution of the nested grids went from 27 arc seconds in the coarsest grid, down to 9 arc seconds in the second grid, 3 arc seconds in the third grid and finally 1 arc second in the finest grid near Patong. The computational domain is shown in Figure\ref{fig:computational_domain}.
 
 Figure \ref{fig:jasonComparison} provides a comparison of the \textsc{ursga} predicted surface elevation with the JASON satellite altimetry data. The \textsc{ursga} model replicates the amplitude and timing of the first peak and trough well. However the model does not resolve the double peak of the first wave. Also note that the \textsc{ursga} model prediction of the ocean surface elevation becomes out of phase with the JASON data at 3 to 7 degrees latitude. Chlieh et al~\cite{chlieh07} also observe these misfits and suggest it is caused by a reflected wave from the Aceh Peninsula that is not resolved in the model due to insufficient resolution of the computational mesh and bathymetry data. This is also a limitation of the model presented here, but probably could be improved by nesting grids near Aceh. 
@@ -235,7 +265,4 @@
 During the \textsc{anuga} simulation the tide was kept constant at $0.80$m. This value was chosen to correspond to the tidal height specified by the Thai Navy tide charts (\url{http://www.navy.mi.th/hydro/}) at the time the tsunami arrived at Patong Bay. Although the tsunami propagated for approximately 3 hours before it reach Patong Bay, the period of time during which the wave propagated through the \textsc{anuga} domain is much smaller. Consequently the assumption of constant tide height is reasonable
 
-FIXME(John): Need a commentary on the dynamics of what is being observed and whether it aligns with eye witness observations.
-Both the URS model and the \textsc{anuga} inundation model shows that the event comprises a train of waves some with preceding drawdown effects Add details of waveform with a graph from URL and a gauge from \textsc{anuga} and discuss. This will come from the work of Kristy (plots of water depth and speed at two locations on beach and one of shore) and Richard (approximating water depth and wave speeds from videos). 
-
 Maximum onshore inundation elevation was simulated throughout the entire Patong Bay region. Figure~\ref{fig:inundationcomparison1cm} shows very good agreement between the measured and simulated inundation. The \textsc{anuga} simulation determines a region to be inundated if at some point in time it was covered by at least 1cm of water. This precision in field measurements is impossible to obtain. The inundation boundary is determined by observing water marks and other signs left by the receding waters. The precision of the observed inundation map is, most likely, at least an order of magnitude worse than the \textsc{anuga} simulation. The simulated inundation based upon a 10cm threshold is shown in Figure~\ref{fig:inundationcomparison1cm}. An inundation threshold of 10cm was selected for all future simulations to reflect the likely accuracy of the survey and subsequently facilitate a more appropriate comparison between the modelled and observed inundation area. 
 
@@ -243,7 +270,4 @@
 
 %\url{https://datamining.anu.edu.au/anuga/attachment/wiki/AnugaPublications/patong_2004_indian_ocean_tsunami_ANUGA_animation.mov}.
-
-
-
 
 \begin{figure}[ht]
@@ -268,7 +292,5 @@
 Additional causes of the discrepancies between the survey data and the modelled inundated include: unknown distribution of surface roughness, inappropriate parameterisation of the source model, effect of humans structures on flow, as well as uncertainties in the elevation data, effects of erosion and deposition by the tsunami event, measurement errors, and missing data in the field survey data itself. The impact of some of these sources of uncertainties are is investigated in Section~\ref{sec:sensitivity}
 
-
-\section{Dynamic behaviour of the tsunami}
-
+\subsection{Eye-witness accounts}
 Figure \ref{fig:gauge_locations} shows four locations where time series have been extracted from the model. The two offshore timeseries are shown in Figure \ref{fig:offshore_timeseries} and
 the two onshore timeseries are shown in Figure \ref{fig:onshore_timeseries}. The latter coincide with locations where video footage from the event is available.
@@ -301,27 +323,9 @@
 \end{figure}
 
-\subsubsection{Video Evidence}
-Two videos were sourced from the internet which include footage of the tsunami in Patong Bay on the day of the Indian Ocean Tsunami. Both videos show an
-already inundated town, they then show what is to be assumed as the second and third waves approaching and further flooding the town. The first video
-is in the very north filmed from what is believed to be the roof of the Novotel Hotel. The second video is in the very south
-filmed from a building next door to the Comort Resort. Figure 11 shows a video taken near the corner of Ruam Chai St. from which we made our estimates.
-
-
-\begin{figure}[ht]
-\begin{center}
-\includegraphics[width=6.0cm,keepaspectratio=true]{flow_rate_south_0_00sec.jpg}
-\includegraphics[width=6.0cm,keepaspectratio=true]{flow_rate_south_5_04sec.jpg}
-\includegraphics[width=6.0cm,keepaspectratio=true]{flow_rate_south_7_12sec.jpg}
-\includegraphics[width=6.0cm,keepaspectratio=true]{flow_rate_south_7_60sec.jpg}
-\caption{four frames from a video where flow rate could be estimated, circle indicates tracked debris, clockwise from top left: 0.0 sec, 5.0 s, 7.1 s, 7.6 s.}
-\label{fig:video_flow}
-\end{center}
-\end{figure}
-
-
 Crude flow rates can be estimated with landmarks found in satellite imagery and the use of a GIS and were found to be in
 the range of 5 to 7 metres per second (+/- 2 m/s) in the north and 0.5 to 2 metres per second (+/- 1 m/s) in the south. This is in agreement
 with results from our simulations. Our modelled flow rates show maximum values in the order of 0.2 to 2.6 m/s in the south and 0.1 to
-3.3 m/s for the north as shown in the figures. Water depths could also be estimated from the videos by the level at which water rose up the sides of buildings such as shops. Our estimates are in the order of 1.5 to 2.0 metres (+/- 0.5 m). This is in the same range as our modelled maximum depths of 1.4 m in the north and 1.5 m in the south as seen in the figure. Fritz ~\cite{fritz06} did a detailed analysis of video frames taken around Banda Aceh and arrived at flow speeds in the range of 2 to 5 m/s.
+3.3 m/s for the north as shown in the figures. Water depths could also be estimated from the videos by the level at which water rose up the sides of buildings such as shops. Our estimates are in the order of 1.5 to 2.0 metres (+/- 0.5 m). This is in the same range as our modelled maximum depths of 1.4 m in the north and 1.5 m in the south as seen in the figure. Fritz ~\cite{fritz06} performed a detailed analysis of video frames taken around Banda Aceh and arrived at flow speeds in the range of 2 to 5 m/s.
+
 
 
@@ -330,5 +334,5 @@
 \section{Sensitivity Analysis}
 \label{sec:sensitivity}
-This section shows how model maximum inundation varies with: different values of Manning's friction coefficient; changing waveheight at the ANUGA boundary (where it was coupled with the URSGA model); and finally the presence and absence of buildings in the elevation dataset.
+This section investigates the effect of different values of Manning's friction coefficient, changing waveheight at the 100m depth contour, and the presence and absence of buildings in the elevation dataset on model maximum inundation.
 
 %========================Friction==========================%
@@ -336,5 +340,5 @@
 \label{sec:friction sensitivity}
 The first study investigated the impact of surface roughness on the predicted run-up. According to Schoettle~\cite{schoettle2007} appropriate values of Manning's coefficient range from 0.007 to 0.030 for tsunami propagation over a sandy sea floor.  Consequently we simulated the maximum onshore inundation using the a Manning's coefficient of 0.0003 and 0.03. The resulting run-up is shown in Figures
-\ref{fig:sensitivity_friction} and  the maximum flow speeds\ref{fig:sensitivity_friction_speed}. These figurers show that the on-shoer inundation extent decreases with increasing friction and that small perturbations in the friction cause bounded changes in the output. This is consistent with the conclusions of Synolakis~\cite{synolakis05} who states that the long wavelength of tsunami tends to mean that the friction is less important in comparison to the motion of the wave.
+\ref{fig:sensitivity_friction} and  the maximum flow speeds\ref{fig:sensitivity_friction_speed}. These figures show that the on-shore inundation extent decreases with increasing friction and that small perturbations in the friction cause bounded changes in the output. This is consistent with the conclusions of Synolakis~\cite{synolakis05} who states that the long wavelength of tsunami tends to mean that the friction is less important in comparison to the motion of the wave.
 
 %========================Wave-Height==========================%
@@ -370,7 +374,9 @@
 
 \section{Conclusion}
-This paper proposes an additional field data benchmark for the verification of tsunami inundation models. Currently, there is a scarcity of appropriate validation datasets due to a lack of well documented historical tsunami impacts. FIXME(John): Complete...
-
-From the sensitivity studies it appears that the presence or absence of buildings is the most important parameter followed by the right choice of friction whereas a small perturbation in the waveheight at the ANUGA boundary has comparatively little effect on the model results.
+This paper proposes an additional field data benchmark for the verification of tsunami inundation models. Currently, there is a scarcity of appropriate validation datasets due to a lack of well documented historical tsunami impacts. The benchmark proposed here utilises the uniquely large amount of observational data for model comparison obtained during, and immediately following, the Sumatra--Andaman tsunami of 26th December 2004. Unlike the small number of existing benchmarks, the proposed test validates all three stages of tsunami evolution - generation, propagation and inundation. In an attempt to provide higher visability and easier accessibility for tsunami benchmark problems the data used to construct the proposed benchmark is documented and freely available at \url{http://tinyurl.com/patong2004-data}. 
+
+An associated aim of this paper was to further validate the hydrodynamic modelling tool \textsc{anuga}  which is used to simulate the tsunami inundation and run rain-induced floods. Model predictions matched well geodetic measurements of the Sumatra--Andaman earthquake, altimetry data from the \textsc{jason}, eye-witness accounts of wave front arrival times and flow speeds and a detailed inundation survey of Patong Bay, Thailand. 
+
+A simple sensitivity analysis was performed to assess the influence of small changes in friction, wave-height at the 100m depth contour and the presence of buildings and other structures on the model predictions. The presence of buildings has the greatest influence on the simulated inundation extent. The value of friction and small perturbations in the waveheight at the ANUGA boundary have comparatively little effect on the model results.
 
 %================Acknowledgement===================
@@ -384,5 +390,5 @@
 \includegraphics[width=3.5cm,keepaspectratio=true]{sensitivity_minus10}
 \includegraphics[width=3.5cm,keepaspectratio=true]{sensitivity_plus10}
-\caption{Model results with wav height at ANUGA boundary artificially modified 
+\caption{Model results with wave height at ANUGA boundary artificially modified 
 to asses sensitivities. The first image is the reference inundation extent as reported in Section \protect \ref{sec:results} while the second and third show the inundation results if the wave at the ANUGA boundary is reduced or increased by 10cm respectively. The inundation severity varies in proportion to the boundary waveheight, but the model results are only slightly sensitive to this parameter for the range of values tested.}
 \label{fig:sensitivity_boundary}
@@ -449,31 +455,5 @@
 
 %====================Bibliography==================
-\bibliographystyle{plain}
+\bibliographystyle{spmpsci}
 \bibliography{tsunami07}
 \end{document}
-
-
-===================
-NOTES TO BE REMOVED
-
-Main source of uncertainty arises from inaccuracies in initial condition (source), inaccurate bathymetry data, to a lesser extent friction
-
-scale comparisons (laboratory benchmarking):
-Scale differences are not believed to be important. scale experiments generally do not have same bootom friction characteristics as real scenario but has not proven to be a problem. The long wavelength of tsunami tends to mean that the friction is less important in comparison to the motion of the wave
-Single wave on a simple beach
-Solitary wave on composite beach
-Conical island
-Monai Valley
-Landslide
-
-Field benchmarking:
-Most important verification process
-Hydrodynamic inversion to predict the source is an ill posed problem
-12 July 1993 Hokkaido-Nansei-Oki tsunami around Okushiri Island Japan extreme runup height of 31.7m was found at the tip of a narrow gully with the small cove at Monai
-17 November 2003 Rat Islands Tsunami
-
-Construction of more than one model can reveal biases in a single model. Two types of comparisons 1 between those that are conceptually similar and those that re different. In former case interested in how choice of numerical solver and discretisation effects results and the later can help determine the level of physical process representation necessary to represent an observed data set.
-
-Moving to field data increases the generality and significance of scientific evidence obtained. However we also significantly increase the uncertainty of the validation experiment that may constrain the ability to make unequivacol statements. E.g. in bathymetry source condition friction.
-
-The two field data benchmarks are very useful but only capture a small subset of possible tsunami behaviours and do not assess all three stages of tsunami evolution (generation,propagation and inundation) together. The type and size of a tsunami source, propagation extent, and local bathymetry and topography all affect the energy, waveform and subsequent inundation of a tsunami. Consequently additional field data benchmarks, such as the one proposed here, that further capture the variability and sensitivity of the real world system would be useful to allow model developers verify their models and subsequently use their models with greater confidence.
Index: /anuga_work/publications/boxing_day_validation_2008/svjour3.cls
===================================================================
--- /anuga_work/publications/boxing_day_validation_2008/svjour3.cls	(revision 7210)
+++ /anuga_work/publications/boxing_day_validation_2008/svjour3.cls	(revision 7210)
@@ -0,0 +1,1431 @@
+% SVJour3 DOCUMENT CLASS -- version 3.2 for LaTeX2e
+%
+% LaTeX document class for Springer journals
+%
+%%
+%%
+%% \CharacterTable
+%%  {Upper-case    \A\B\C\D\E\F\G\H\I\J\K\L\M\N\O\P\Q\R\S\T\U\V\W\X\Y\Z
+%%   Lower-case    \a\b\c\d\e\f\g\h\i\j\k\l\m\n\o\p\q\r\s\t\u\v\w\x\y\z
+%%   Digits        \0\1\2\3\4\5\6\7\8\9
+%%   Exclamation   \!     Double quote  \"     Hash (number) \#
+%%   Dollar        \$     Percent       \%     Ampersand     \&
+%%   Acute accent  \'     Left paren    \(     Right paren   \)
+%%   Asterisk      \*     Plus          \+     Comma         \,
+%%   Minus         \-     Point         \.     Solidus       \/
+%%   Colon         \:     Semicolon     \;     Less than     \<
+%%   Equals        \=     Greater than  \>     Question mark \?
+%%   Commercial at \@     Left bracket  \[     Backslash     \\
+%%   Right bracket \]     Circumflex    \^     Underscore    \_
+%%   Grave accent  \`     Left brace    \{     Vertical bar  \|
+%%   Right brace   \}     Tilde         \~}
+\NeedsTeXFormat{LaTeX2e}[1995/12/01]
+\ProvidesClass{svjour3}[2007/05/08 v3.2
+^^JLaTeX document class for Springer journals]
+\newcommand\@ptsize{}
+\newif\if@restonecol
+\newif\if@titlepage
+\@titlepagefalse
+\DeclareOption{a4paper}
+   {\setlength\paperheight {297mm}%
+    \setlength\paperwidth  {210mm}}
+\DeclareOption{10pt}{\renewcommand\@ptsize{0}}
+\DeclareOption{twoside}{\@twosidetrue  \@mparswitchtrue}
+\DeclareOption{draft}{\setlength\overfullrule{5pt}}
+\DeclareOption{final}{\setlength\overfullrule{0pt}}
+\DeclareOption{fleqn}{\input{fleqn.clo}\AtBeginDocument{\mathindent\z@}%
+\AtBeginDocument{\@ifpackageloaded{amsmath}{\@mathmargin\z@}{}}%
+\PassOptionsToPackage{fleqn}{amsmath}}
+%%%
+\DeclareOption{onecolumn}{}
+\DeclareOption{smallcondensed}{}
+\DeclareOption{twocolumn}{\@twocolumntrue\ExecuteOptions{fleqn}}
+\newif\if@smallext\@smallextfalse
+\DeclareOption{smallextended}{\@smallexttrue}
+\let\if@mathematic\iftrue
+\let\if@numbook\iffalse
+\DeclareOption{numbook}{\let\if@envcntsect\iftrue
+  \AtEndOfPackage{%
+   \renewcommand\thefigure{\thesection.\@arabic\c@figure}%
+   \renewcommand\thetable{\thesection.\@arabic\c@table}%
+   \renewcommand\theequation{\thesection.\@arabic\c@equation}%
+   \@addtoreset{figure}{section}%
+   \@addtoreset{table}{section}%
+   \@addtoreset{equation}{section}%
+  }%
+}
+\DeclareOption{openbib}{%
+  \AtEndOfPackage{%
+   \renewcommand\@openbib@code{%
+      \advance\leftmargin\bibindent
+      \itemindent -\bibindent
+      \listparindent \itemindent
+      \parsep \z@
+      }%
+   \renewcommand\newblock{\par}}%
+}
+\DeclareOption{natbib}{%
+\AtEndOfClass{\RequirePackage{natbib}%
+% Changing some parameters of NATBIB
+\setlength{\bibhang}{\parindent}%
+%\setlength{\bibsep}{0mm}%
+\let\bibfont=\small
+\def\@biblabel#1{#1.}%
+\newcommand{\etal}{et al.}%
+\bibpunct{(}{)}{;}{a}{}{,}}}
+%
+\let\if@runhead\iffalse
+\DeclareOption{runningheads}{\let\if@runhead\iftrue}
+\let\if@smartrunh\iffalse
+\DeclareOption{smartrunhead}{\let\if@smartrunh\iftrue}
+\DeclareOption{nosmartrunhead}{\let\if@smartrunh\iffalse}
+\let\if@envcntreset\iffalse
+\DeclareOption{envcountreset}{\let\if@envcntreset\iftrue}
+\let\if@envcntsame\iffalse
+\DeclareOption{envcountsame}{\let\if@envcntsame\iftrue}
+\let\if@envcntsect\iffalse
+\DeclareOption{envcountsect}{\let\if@envcntsect\iftrue}
+\let\if@referee\iffalse
+\DeclareOption{referee}{\let\if@referee\iftrue}
+\def\makereferee{\def\baselinestretch{2}}
+\let\if@instindent\iffalse
+\DeclareOption{instindent}{\let\if@instindent\iftrue}
+\let\if@smartand\iffalse
+\DeclareOption{smartand}{\let\if@smartand\iftrue}
+\let\if@spthms\iftrue
+\DeclareOption{nospthms}{\let\if@spthms\iffalse}
+%
+% language and babel dependencies
+\DeclareOption{deutsch}{\def\switcht@@therlang{\switcht@deutsch}%
+\gdef\svlanginfo{\typeout{Man spricht deutsch.}\global\let\svlanginfo\relax}}
+\DeclareOption{francais}{\def\switcht@@therlang{\switcht@francais}%
+\gdef\svlanginfo{\typeout{On parle francais.}\global\let\svlanginfo\relax}}
+\let\switcht@@therlang\relax
+\let\svlanginfo\relax
+%
+\AtBeginDocument{\@ifpackageloaded{babel}{%
+\@ifundefined{extrasenglish}{}{\addto\extrasenglish{\switcht@albion}}%
+\@ifundefined{extrasUKenglish}{}{\addto\extrasUKenglish{\switcht@albion}}%
+\@ifundefined{extrasfrenchb}{}{\addto\extrasfrenchb{\switcht@francais}}%
+\@ifundefined{extrasgerman}{}{\addto\extrasgerman{\switcht@deutsch}}%
+\@ifundefined{extrasngerman}{}{\addto\extrasngerman{\switcht@deutsch}}%
+}{\switcht@@therlang}%
+}
+%
+\def\ClassInfoNoLine#1#2{%
+   \ClassInfo{#1}{#2\@gobble}%
+}
+\let\journalopt\@empty
+\DeclareOption*{%
+\InputIfFileExists{sv\CurrentOption.clo}{%
+\global\let\journalopt\CurrentOption}{%
+\ClassWarning{Springer-SVJour3}{Specified option or subpackage
+"\CurrentOption" not found -}\OptionNotUsed}}
+\ExecuteOptions{a4paper,twoside,10pt,instindent}
+\ProcessOptions
+%
+\ifx\journalopt\@empty\relax
+\ClassInfoNoLine{Springer-SVJour3}{extra/valid Springer sub-package (-> *.clo)
+\MessageBreak not found in option list of \string\documentclass
+\MessageBreak  - autoactivating "global" style}{}
+\input{svglov3.clo}
+\else
+\@ifundefined{validfor}{%
+\ClassError{Springer-SVJour3}{Possible option clash for sub-package
+\MessageBreak "sv\journalopt.clo" - option file not valid
+\MessageBreak for this class}{Perhaps you used an option of the old
+Springer class SVJour!}
+}{}
+\fi
+%
+\if@smartrunh\AtEndDocument{\islastpageeven\getlastpagenumber}\fi
+%
+\newcommand{\twocoltest}[2]{\if@twocolumn\def\@gtempa{#2}\else\def\@gtempa{#1}\fi
+\@gtempa\makeatother}
+\newcommand{\columncase}{\makeatletter\twocoltest}
+%
+\DeclareMathSymbol{\Gamma}{\mathalpha}{letters}{"00}
+\DeclareMathSymbol{\Delta}{\mathalpha}{letters}{"01}
+\DeclareMathSymbol{\Theta}{\mathalpha}{letters}{"02}
+\DeclareMathSymbol{\Lambda}{\mathalpha}{letters}{"03}
+\DeclareMathSymbol{\Xi}{\mathalpha}{letters}{"04}
+\DeclareMathSymbol{\Pi}{\mathalpha}{letters}{"05}
+\DeclareMathSymbol{\Sigma}{\mathalpha}{letters}{"06}
+\DeclareMathSymbol{\Upsilon}{\mathalpha}{letters}{"07}
+\DeclareMathSymbol{\Phi}{\mathalpha}{letters}{"08}
+\DeclareMathSymbol{\Psi}{\mathalpha}{letters}{"09}
+\DeclareMathSymbol{\Omega}{\mathalpha}{letters}{"0A}
+%
+\setlength\parindent{15\p@}
+\setlength\smallskipamount{3\p@ \@plus 1\p@ \@minus 1\p@}
+\setlength\medskipamount{6\p@ \@plus 2\p@ \@minus 2\p@}
+\setlength\bigskipamount{12\p@ \@plus 4\p@ \@minus 4\p@}
+\setlength\headheight{12\p@}
+\setlength\headsep   {14.50dd}
+\setlength\topskip   {10\p@}
+\setlength\footskip{30\p@}
+\setlength\maxdepth{.5\topskip}
+%
+\@settopoint\textwidth
+\setlength\marginparsep {10\p@}
+\setlength\marginparpush{5\p@}
+\setlength\topmargin{-10pt}
+\if@twocolumn
+   \setlength\oddsidemargin {-30\p@}
+   \setlength\evensidemargin{-30\p@}
+\else
+   \setlength\oddsidemargin {\z@}
+   \setlength\evensidemargin{\z@}
+\fi
+\setlength\marginparwidth  {48\p@}
+\setlength\footnotesep{8\p@}
+\setlength{\skip\footins}{9\p@ \@plus 4\p@ \@minus 2\p@}
+\setlength\floatsep    {12\p@ \@plus 2\p@ \@minus 2\p@}
+\setlength\textfloatsep{20\p@ \@plus 2\p@ \@minus 4\p@}
+\setlength\intextsep   {20\p@ \@plus 2\p@ \@minus 2\p@}
+\setlength\dblfloatsep    {12\p@ \@plus 2\p@ \@minus 2\p@}
+\setlength\dbltextfloatsep{20\p@ \@plus 2\p@ \@minus 4\p@}
+\setlength\@fptop{0\p@}
+\setlength\@fpsep{12\p@ \@plus 2\p@ \@minus 2\p@}
+\setlength\@fpbot{0\p@ \@plus 1fil}
+\setlength\@dblfptop{0\p@}
+\setlength\@dblfpsep{12\p@ \@plus 2\p@ \@minus 2\p@}
+\setlength\@dblfpbot{0\p@ \@plus 1fil}
+\setlength\partopsep{2\p@ \@plus 1\p@ \@minus 1\p@}
+\def\@listi{\leftmargin\leftmargini
+            \parsep \z@
+            \topsep 6\p@ \@plus2\p@ \@minus4\p@
+            \itemsep\parsep}
+\let\@listI\@listi
+\@listi
+\def\@listii {\leftmargin\leftmarginii
+              \labelwidth\leftmarginii
+              \advance\labelwidth-\labelsep
+              \topsep    \z@
+              \parsep    \topsep
+              \itemsep   \parsep}
+\def\@listiii{\leftmargin\leftmarginiii
+              \labelwidth\leftmarginiii
+              \advance\labelwidth-\labelsep
+              \topsep    \z@
+              \parsep    \topsep
+              \itemsep   \parsep}
+\def\@listiv {\leftmargin\leftmarginiv
+              \labelwidth\leftmarginiv
+              \advance\labelwidth-\labelsep}
+\def\@listv  {\leftmargin\leftmarginv
+              \labelwidth\leftmarginv
+              \advance\labelwidth-\labelsep}
+\def\@listvi {\leftmargin\leftmarginvi
+              \labelwidth\leftmarginvi
+              \advance\labelwidth-\labelsep}
+%
+\setlength\lineskip{1\p@}
+\setlength\normallineskip{1\p@}
+\renewcommand\baselinestretch{}
+\setlength\parskip{0\p@ \@plus \p@}
+\@lowpenalty   51
+\@medpenalty  151
+\@highpenalty 301
+\setcounter{topnumber}{4}
+\renewcommand\topfraction{.9}
+\setcounter{bottomnumber}{2}
+\renewcommand\bottomfraction{.7}
+\setcounter{totalnumber}{6}
+\renewcommand\textfraction{.1}
+\renewcommand\floatpagefraction{.85}
+\setcounter{dbltopnumber}{3}
+\renewcommand\dbltopfraction{.85}
+\renewcommand\dblfloatpagefraction{.85}
+\def\ps@headings{%
+    \let\@oddfoot\@empty\let\@evenfoot\@empty
+    \def\@evenhead{\small\csname runheadhook\endcsname
+    \rlap{\thepage}\hfil\leftmark\unskip}%
+    \def\@oddhead{\small\csname runheadhook\endcsname
+    \ignorespaces\rightmark\hfil\llap{\thepage}}%
+    \let\@mkboth\@gobbletwo
+    \let\sectionmark\@gobble
+    \let\subsectionmark\@gobble
+    }
+% make indentations changeable
+\def\setitemindent#1{\settowidth{\labelwidth}{#1}%
+        \leftmargini\labelwidth
+        \advance\leftmargini\labelsep
+   \def\@listi{\leftmargin\leftmargini
+        \labelwidth\leftmargini\advance\labelwidth by -\labelsep
+        \parsep=\parskip
+        \topsep=\medskipamount
+        \itemsep=\parskip \advance\itemsep by -\parsep}}
+\def\setitemitemindent#1{\settowidth{\labelwidth}{#1}%
+        \leftmarginii\labelwidth
+        \advance\leftmarginii\labelsep
+\def\@listii{\leftmargin\leftmarginii
+        \labelwidth\leftmarginii\advance\labelwidth by -\labelsep
+        \parsep=\parskip
+        \topsep=\z@
+        \itemsep=\parskip \advance\itemsep by -\parsep}}
+% labels of description
+\def\descriptionlabel#1{\hspace\labelsep #1\hfil}
+% adjusted environment "description"
+% if an optional parameter (at the first two levels of lists)
+% is present, its width is considered to be the widest mark
+% throughout the current list.
+\def\description{\@ifnextchar[{\@describe}{\list{}{\labelwidth\z@
+          \itemindent-\leftmargin \let\makelabel\descriptionlabel}}}
+\let\enddescription\endlist
+%
+\def\describelabel#1{#1\hfil}
+\def\@describe[#1]{\relax\ifnum\@listdepth=0
+\setitemindent{#1}\else\ifnum\@listdepth=1
+\setitemitemindent{#1}\fi\fi
+\list{--}{\let\makelabel\describelabel}}
+%
+\newdimen\logodepth
+\logodepth=1.2cm
+\newdimen\headerboxheight
+\headerboxheight=163pt % 18 10.5dd-lines - 2\baselineskip
+\if@twocolumn\else\advance\headerboxheight by-14.5mm\fi
+\newdimen\betweenumberspace          % dimension for space between
+\betweenumberspace=3.33pt            % number and text of titles.
+\newdimen\aftertext                  % dimension for space after
+\aftertext=5pt                       % text of title.
+\newdimen\headlineindent             % dimension for space between
+\headlineindent=1.166cm              % number and text of headings.
+\if@mathematic
+   \def\runinend{} % \enspace}
+   \def\floatcounterend{\enspace}
+   \def\sectcounterend{}
+\else
+   \def\runinend{.}
+   \def\floatcounterend{.\ }
+   \def\sectcounterend{.}
+\fi
+\def\email#1{\emailname: #1}
+\def\keywords#1{\par\addvspace\medskipamount{\rightskip=0pt plus1cm
+\def\and{\ifhmode\unskip\nobreak\fi\ $\cdot$
+}\noindent\keywordname\enspace\ignorespaces#1\par}}
+%
+\def\subclassname{{\bfseries Mathematics Subject Classification
+(2000)}\enspace}
+\def\subclass#1{\par\addvspace\medskipamount{\rightskip=0pt plus1cm
+\def\and{\ifhmode\unskip\nobreak\fi\ $\cdot$
+}\noindent\subclassname\ignorespaces#1\par}}
+%
+\def\PACSname{\textbf{PACS}\enspace}
+\def\PACS#1{\par\addvspace\medskipamount{\rightskip=0pt plus1cm
+\def\and{\ifhmode\unskip\nobreak\fi\ $\cdot$
+}\noindent\PACSname\ignorespaces#1\par}}
+%
+\def\CRclassname{{\bfseries CR Subject Classification}\enspace}
+\def\CRclass#1{\par\addvspace\medskipamount{\rightskip=0pt plus1cm
+\def\and{\ifhmode\unskip\nobreak\fi\ $\cdot$
+}\noindent\CRclassname\ignorespaces#1\par}}
+%
+\def\ESMname{\textbf{Electronic Supplementary Material}\enspace}
+\def\ESM#1{\par\addvspace\medskipamount
+\noindent\ESMname\ignorespaces#1\par}
+%
+\newcounter{inst}
+\newcounter{auth}
+\def\authdepth{2}
+\newdimen\instindent
+\newbox\authrun
+\newtoks\authorrunning
+\newbox\titrun
+\newtoks\titlerunning
+\def\authorfont{\bfseries}
+
+\def\combirunning#1{\gdef\@combi{#1}}
+\def\@combi{}
+\newbox\combirun
+%
+\def\ps@last{\def\@evenhead{\small\rlap{\thepage}\hfil
+\lastevenhead}}
+\newcounter{lastpage}
+\def\islastpageeven{\@ifundefined{lastpagenumber}
+{\setcounter{lastpage}{0}}{\setcounter{lastpage}{\lastpagenumber}}
+\ifnum\value{lastpage}>0
+   \ifodd\value{lastpage}%
+   \else
+      \if@smartrunh
+         \thispagestyle{last}%
+      \fi
+   \fi
+\fi}
+\def\getlastpagenumber{\clearpage
+\addtocounter{page}{-1}%
+   \immediate\write\@auxout{\string\gdef\string\lastpagenumber{\thepage}}%
+   \immediate\write\@auxout{\string\newlabel{LastPage}{{}{\thepage}}}%
+   \addtocounter{page}{1}}
+
+\def\journalname#1{\gdef\@journalname{#1}}
+
+\def\dedication#1{\gdef\@dedic{#1}}
+\def\@dedic{}
+
+\let\@date\undefined
+\def\notused{~}
+
+\def\institute#1{\gdef\@institute{#1}}
+
+\def\offprints#1{\begingroup
+\def\protect{\noexpand\protect\noexpand}\xdef\@thanks{\@thanks
+\protect\footnotetext[0]{\unskip\hskip-15pt{\itshape Send offprint requests
+to\/}: \ignorespaces#1}}\endgroup\ignorespaces}
+
+%\def\mail#1{\gdef\@mail{#1}}
+%\def\@mail{}
+
+\def\@thanks{}
+
+\def\@fnsymbol#1{\ifcase#1\or\star\or{\star\star}\or{\star\star\star}%
+   \or \dagger\or \ddagger\or
+   \mathchar "278\or \mathchar "27B\or \|\or **\or \dagger\dagger
+   \or \ddagger\ddagger \else\@ctrerr\fi\relax}
+%
+%\def\invthanks#1{\footnotetext[0]{\kern-\bibindent#1}}
+%
+\def\nothanksmarks{\def\thanks##1{\protected@xdef\@thanks{\@thanks
+        \protect\footnotetext[0]{\kern-\bibindent##1}}}}
+%
+\def\subtitle#1{\gdef\@subtitle{#1}}
+\def\@subtitle{}
+
+\def\headnote#1{\gdef\@headnote{#1}}
+\def\@headnote{}
+
+\def\papertype#1{\gdef\paper@type{\MakeUppercase{#1}}}
+\def\paper@type{}
+
+\def\ch@ckobl#1#2{\@ifundefined{@#1}
+ {\typeout{SVJour3 warning: Missing
+\expandafter\string\csname#1\endcsname}%
+  \csname #1\endcsname{#2}}
+ {}}
+%
+\def\ProcessRunnHead{%
+    \def\\{\unskip\ \ignorespaces}%
+    \def\thanks##1{\unskip{}}%
+    \instindent=\textwidth
+    \advance\instindent by-\headlineindent
+    \if!\the\titlerunning!\else
+      \edef\@title{\the\titlerunning}%
+    \fi
+    \global\setbox\titrun=\hbox{\small\rmfamily\unboldmath\ignorespaces\@title
+                                \unskip}%
+    \ifdim\wd\titrun>\instindent
+       \typeout{^^JSVJour3 Warning: Title too long for running head.}%
+       \typeout{Please supply a shorter form with \string\titlerunning
+                \space prior to \string\maketitle}%
+       \global\setbox\titrun=\hbox{\small\rmfamily
+       Title Suppressed Due to Excessive Length}%
+    \fi
+    \xdef\@title{\copy\titrun}%
+%
+    \if!\the\authorrunning!
+    \else
+      \setcounter{auth}{1}%
+      \edef\@author{\the\authorrunning}%
+    \fi
+    \ifnum\value{inst}>\authdepth
+       \def\stripauthor##1\and##2\endauthor{%
+       \protected@xdef\@author{##1\unskip\unskip\if!##2!\else\ et al.\fi}}%
+       \expandafter\stripauthor\@author\and\endauthor
+    \else
+       \gdef\and{\unskip, \ignorespaces}%
+       {\def\and{\noexpand\protect\noexpand\and}%
+       \protected@xdef\@author{\@author}}
+    \fi
+    \global\setbox\authrun=\hbox{\small\rmfamily\unboldmath\ignorespaces
+                                 \@author\unskip}%
+    \ifdim\wd\authrun>\instindent
+    \typeout{^^JSVJour3 Warning: Author name(s) too long for running head.
+             ^^JPlease supply a shorter form with \string\authorrunning
+             \space prior to \string\maketitle}%
+    \global\setbox\authrun=\hbox{\small\rmfamily Please give a shorter version
+          with: {\tt\string\authorrunning\space and
+             \string\titlerunning\space prior to \string\maketitle}}%
+    \fi
+    \xdef\@author{\copy\authrun}%
+    \markboth{\@author}{\@title}%
+}
+%
+\let\orithanks=\thanks
+\def\thanks#1{\ClassWarning{SVJour3}{\string\thanks\space may only be
+used inside of \string\title, \string\author,\MessageBreak
+and \string\date\space prior to \string\maketitle}}
+%
+\def\maketitle{\par\let\thanks=\orithanks
+\ch@ckobl{journalname}{Noname}
+\ch@ckobl{date}{the date of receipt and acceptance should be inserted
+later}
+\ch@ckobl{title}{A title should be given}
+\ch@ckobl{author}{Name(s) and initial(s) of author(s) should be given}
+\ch@ckobl{institute}{Address(es) of author(s) should be given}
+\begingroup
+%
+    \renewcommand\thefootnote{\@fnsymbol\c@footnote}%
+    \def\@makefnmark{$^{\@thefnmark}$}%
+    \renewcommand\@makefntext[1]{%
+    \noindent
+    \hb@xt@\bibindent{\hss\@makefnmark\enspace}##1\vrule height0pt
+    width0pt depth8pt}
+%
+ \def\lastand{\ifnum\value{inst}=2\relax
+                 \unskip{} \andname\
+              \else
+                 \unskip, \andname\
+              \fi}%
+ \def\and{\stepcounter{auth}\relax
+          \if@smartand
+             \ifnum\value{auth}=\value{inst}%
+                \lastand
+             \else
+                \unskip,
+             \fi
+          \else
+             \unskip,
+          \fi}%
+ \thispagestyle{empty}
+ \ifnum \col@number=\@ne
+   \@maketitle
+ \else
+   \twocolumn[\@maketitle]%
+ \fi
+%
+ \global\@topnum\z@
+ \if!\@thanks!\else
+    \@thanks
+\insert\footins{\vskip-3pt\hrule\@width\if@twocolumn\columnwidth
+\else 38mm\fi\vskip3pt}%
+ \fi
+ {\def\thanks##1{\unskip{}}%
+ \def\iand{\\[5pt]\let\and=\nand}%
+ \def\nand{\ifhmode\unskip\nobreak\fi\ $\cdot$ }%
+ \let\and=\nand
+ \def\at{\\\let\and=\iand}%
+ \footnotetext[0]{\kern-\bibindent
+ \ignorespaces\@institute}\vspace{5dd}}%
+%\if!\@mail!\else
+%   \footnotetext[0]{\kern-\bibindent\mailname\
+%   \ignorespaces\@mail}%
+%\fi
+%
+ \if@runhead
+    \ProcessRunnHead
+ \fi
+%
+ \endgroup
+ \setcounter{footnote}{0}
+ \global\let\thanks\relax
+ \global\let\maketitle\relax
+ \global\let\@maketitle\relax
+ \global\let\@thanks\@empty
+ \global\let\@author\@empty
+ \global\let\@date\@empty
+ \global\let\@title\@empty
+ \global\let\@subtitle\@empty
+ \global\let\title\relax
+ \global\let\author\relax
+ \global\let\date\relax
+ \global\let\and\relax}
+
+\def\makeheadbox{{%
+\hbox to0pt{\vbox{\baselineskip=10dd\hrule\hbox
+to\hsize{\vrule\kern3pt\vbox{\kern3pt
+\hbox{\bfseries\@journalname\ manuscript No.}
+\hbox{(will be inserted by the editor)}
+\kern3pt}\hfil\kern3pt\vrule}\hrule}%
+\hss}}}
+%
+\def\rubric{\setbox0=\hbox{\small\strut}\@tempdima=\ht0\advance
+\@tempdima\dp0\advance\@tempdima2\fboxsep\vrule\@height\@tempdima
+\@width\z@}
+\newdimen\rubricwidth
+%
+\def\@maketitle{\newpage
+\normalfont
+\vbox to0pt{\if@twocolumn\vskip-39pt\else\vskip-49pt\fi
+\nointerlineskip
+\makeheadbox\vss}\nointerlineskip
+\vbox to 0pt{\offinterlineskip\rubricwidth=\columnwidth
+%%%%\vskip-12.5pt          % -12.5pt
+\if@twocolumn\else % one column journal
+   \divide\rubricwidth by144\multiply\rubricwidth by89 % perform golden section
+   \vskip-\topskip
+\fi
+\hrule\@height0.35mm\noindent
+\advance\fboxsep by.25mm
+\global\advance\rubricwidth by0pt
+\rubric
+\vss}\vskip19.5pt    % war 9pt
+%
+\if@twocolumn\else
+ \gdef\footnoterule{%
+  \kern-3\p@
+  \hrule\@width38mm     % \columnwidth  \rubricwidth
+  \kern2.6\p@}
+\fi
+%
+ \setbox\authrun=\vbox\bgroup
+     \if@twocolumn
+        \hrule\@height10.5mm\@width0\p@
+     \else
+        \hrule\@height 2mm\@width0\p@
+     \fi
+     \pretolerance=10000
+     \rightskip=0pt plus 4cm
+    \nothanksmarks
+%   \if!\@headnote!\else
+%     \noindent
+%     {\LARGE\normalfont\itshape\ignorespaces\@headnote\par}\vskip 3.5mm
+%   \fi
+    {\LARGE\bfseries
+     \noindent\ignorespaces
+     \@title \par}\vskip 17pt\relax
+    \if!\@subtitle!\else
+      {\large\bfseries
+      \pretolerance=10000
+      \rightskip=0pt plus 3cm
+      \vskip-12pt
+%     \noindent\ignorespaces\@subtitle \par}\vskip 11.24pt\relax
+      \noindent\ignorespaces\@subtitle \par}\vskip 17pt\relax
+    \fi
+    {\authorfont
+    \setbox0=\vbox{\setcounter{auth}{1}\def\and{\stepcounter{auth} }%
+                   \hfuzz=2\textwidth\def\thanks##1{}\@author}%
+    \setcounter{footnote}{0}%
+    \global\value{inst}=\value{auth}%
+    \setcounter{auth}{1}%
+    \if@twocolumn
+       \rightskip43mm plus 4cm minus 3mm
+    \else % one column journal
+       \rightskip=\linewidth
+       \advance\rightskip by-\rubricwidth
+       \advance\rightskip by0pt plus 4cm minus 3mm
+    \fi
+%
+\def\and{\unskip\nobreak\enskip{\boldmath$\cdot$}\enskip\ignorespaces}%
+    \noindent\ignorespaces\@author\vskip7.23pt}
+%
+    \small
+    \if!\@dedic!\else
+       \par
+       \normalsize\it
+       \addvspace\baselineskip
+       \noindent\@dedic
+    \fi
+ \egroup % end of header box
+ \@tempdima=\headerboxheight
+ \advance\@tempdima by-\ht\authrun
+ \unvbox\authrun
+ \ifdim\@tempdima>0pt
+    \vrule width0pt height\@tempdima\par
+ \fi
+ \noindent{\small\@date\if@twocolumn\vskip 7.2mm\else\vskip 5.2mm\fi}
+ \global\@minipagetrue
+ \global\everypar{\global\@minipagefalse\global\everypar{}}%
+%\vskip22.47pt
+}
+%
+\if@mathematic
+   \def\vec#1{\ensuremath{\mathchoice
+                     {\mbox{\boldmath$\displaystyle\mathbf{#1}$}}
+                     {\mbox{\boldmath$\textstyle\mathbf{#1}$}}
+                     {\mbox{\boldmath$\scriptstyle\mathbf{#1}$}}
+                     {\mbox{\boldmath$\scriptscriptstyle\mathbf{#1}$}}}}
+\else
+   \def\vec#1{\ensuremath{\mathchoice
+                     {\mbox{\boldmath$\displaystyle#1$}}
+                     {\mbox{\boldmath$\textstyle#1$}}
+                     {\mbox{\boldmath$\scriptstyle#1$}}
+                     {\mbox{\boldmath$\scriptscriptstyle#1$}}}}
+\fi
+%
+\def\tens#1{\ensuremath{\mathsf{#1}}}
+%
+\setcounter{secnumdepth}{3}
+\newcounter {section}
+\newcounter {subsection}[section]
+\newcounter {subsubsection}[subsection]
+\newcounter {paragraph}[subsubsection]
+\newcounter {subparagraph}[paragraph]
+\renewcommand\thesection      {\@arabic\c@section}
+\renewcommand\thesubsection   {\thesection.\@arabic\c@subsection}
+\renewcommand\thesubsubsection{\thesubsection.\@arabic\c@subsubsection}
+\renewcommand\theparagraph    {\thesubsubsection.\@arabic\c@paragraph}
+\renewcommand\thesubparagraph {\theparagraph.\@arabic\c@subparagraph}
+%
+\def\@hangfrom#1{\setbox\@tempboxa\hbox{#1}%
+      \hangindent \z@\noindent\box\@tempboxa}
+%
+\def\@seccntformat#1{\csname the#1\endcsname\sectcounterend
+\hskip\betweenumberspace}
+%
+% \newif\if@sectrule
+% \if@twocolumn\else\let\@sectruletrue=\relax\fi
+% \if@avier\let\@sectruletrue=\relax\fi
+% \def\makesectrule{\if@sectrule\global\@sectrulefalse\null\vglue-\topskip
+% \hrule\nobreak\parskip=5pt\relax\fi}
+% %
+% \let\makesectruleori=\makesectrule
+% \def\restoresectrule{\global\let\makesectrule=\makesectruleori\global\@sectrulefalse}
+% \def\nosectrule{\let\makesectrule=\restoresectrule}
+%
+\def\@startsection#1#2#3#4#5#6{%
+  \if@noskipsec \leavevmode \fi
+  \par
+  \@tempskipa #4\relax
+  \@afterindenttrue
+  \ifdim \@tempskipa <\z@
+    \@tempskipa -\@tempskipa \@afterindentfalse
+  \fi
+  \if@nobreak
+    \everypar{}%
+  \else
+    \addpenalty\@secpenalty\addvspace\@tempskipa
+  \fi
+% \ifnum#2=1\relax\@sectruletrue\fi
+  \@ifstar
+    {\@ssect{#3}{#4}{#5}{#6}}%
+    {\@dblarg{\@sect{#1}{#2}{#3}{#4}{#5}{#6}}}}
+%
+\def\@sect#1#2#3#4#5#6[#7]#8{%
+  \ifnum #2>\c@secnumdepth
+    \let\@svsec\@empty
+  \else
+    \refstepcounter{#1}%
+    \protected@edef\@svsec{\@seccntformat{#1}\relax}%
+  \fi
+  \@tempskipa #5\relax
+  \ifdim \@tempskipa>\z@
+    \begingroup
+      #6{% \makesectrule
+        \@hangfrom{\hskip #3\relax\@svsec}%
+          \raggedright
+          \hyphenpenalty \@M%
+          \interlinepenalty \@M #8\@@par}%
+    \endgroup
+    \csname #1mark\endcsname{#7}%
+    \addcontentsline{toc}{#1}{%
+      \ifnum #2>\c@secnumdepth \else
+        \protect\numberline{\csname the#1\endcsname\sectcounterend}%
+      \fi
+      #7}%
+  \else
+    \def\@svsechd{%
+      #6{\hskip #3\relax
+      \@svsec #8\/\hskip\aftertext}%
+      \csname #1mark\endcsname{#7}%
+      \addcontentsline{toc}{#1}{%
+        \ifnum #2>\c@secnumdepth \else
+          \protect\numberline{\csname the#1\endcsname}%
+        \fi
+        #7}}%
+  \fi
+  \@xsect{#5}}
+%
+\def\@ssect#1#2#3#4#5{%
+  \@tempskipa #3\relax
+  \ifdim \@tempskipa>\z@
+    \begingroup
+      #4{% \makesectrule
+        \@hangfrom{\hskip #1}%
+          \interlinepenalty \@M #5\@@par}%
+    \endgroup
+  \else
+    \def\@svsechd{#4{\hskip #1\relax #5}}%
+  \fi
+  \@xsect{#3}}
+
+%
+% measures and setting of sections
+%
+\def\section{\@startsection{section}{1}{\z@}%
+    {-21dd plus-8pt minus-4pt}{10.5dd}
+     {\normalsize\bfseries\boldmath}}
+\def\subsection{\@startsection{subsection}{2}{\z@}%
+    {-21dd plus-8pt minus-4pt}{10.5dd}
+     {\normalsize\upshape}}
+\def\subsubsection{\@startsection{subsubsection}{3}{\z@}%
+    {-13dd plus-8pt minus-4pt}{10.5dd}
+     {\normalsize\itshape}}
+\def\paragraph{\@startsection{paragraph}{4}{\z@}%
+    {-13pt plus-8pt minus-4pt}{\z@}{\normalsize\itshape}}
+
+\setlength\leftmargini  {\parindent}
+\leftmargin  \leftmargini
+\setlength\leftmarginii {\parindent}
+\setlength\leftmarginiii {1.87em}
+\setlength\leftmarginiv  {1.7em}
+\setlength\leftmarginv  {.5em}
+\setlength\leftmarginvi {.5em}
+\setlength  \labelsep  {.5em}
+\setlength  \labelwidth{\leftmargini}
+\addtolength\labelwidth{-\labelsep}
+\@beginparpenalty -\@lowpenalty
+\@endparpenalty   -\@lowpenalty
+\@itempenalty     -\@lowpenalty
+\renewcommand\theenumi{\@arabic\c@enumi}
+\renewcommand\theenumii{\@alph\c@enumii}
+\renewcommand\theenumiii{\@roman\c@enumiii}
+\renewcommand\theenumiv{\@Alph\c@enumiv}
+\newcommand\labelenumi{\theenumi.}
+\newcommand\labelenumii{(\theenumii)}
+\newcommand\labelenumiii{\theenumiii.}
+\newcommand\labelenumiv{\theenumiv.}
+\renewcommand\p@enumii{\theenumi}
+\renewcommand\p@enumiii{\theenumi(\theenumii)}
+\renewcommand\p@enumiv{\p@enumiii\theenumiii}
+\newcommand\labelitemi{\normalfont\bfseries --}
+\newcommand\labelitemii{\normalfont\bfseries --}
+\newcommand\labelitemiii{$\m@th\bullet$}
+\newcommand\labelitemiv{$\m@th\cdot$}
+
+\if@spthms
+% definition of the "\spnewtheorem" command.
+%
+% Usage:
+%
+%     \spnewtheorem{env_nam}{caption}[within]{cap_font}{body_font}
+% or  \spnewtheorem{env_nam}[numbered_like]{caption}{cap_font}{body_font}
+% or  \spnewtheorem*{env_nam}{caption}{cap_font}{body_font}
+%
+% New is "cap_font" and "body_font". It stands for
+% fontdefinition of the caption and the text itself.
+%
+% "\spnewtheorem*" gives a theorem without number.
+%
+% A defined spnewthoerem environment is used as described
+% by Lamport.
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+\def\@thmcountersep{}
+\def\@thmcounterend{}
+\newcommand\nocaption{\noexpand\@gobble}
+\newdimen\spthmsep \spthmsep=5pt
+
+\def\spnewtheorem{\@ifstar{\@sthm}{\@Sthm}}
+
+% definition of \spnewtheorem with number
+
+\def\@spnthm#1#2{%
+  \@ifnextchar[{\@spxnthm{#1}{#2}}{\@spynthm{#1}{#2}}}
+\def\@Sthm#1{\@ifnextchar[{\@spothm{#1}}{\@spnthm{#1}}}
+
+\def\@spxnthm#1#2[#3]#4#5{\expandafter\@ifdefinable\csname #1\endcsname
+   {\@definecounter{#1}\@addtoreset{#1}{#3}%
+   \expandafter\xdef\csname the#1\endcsname{\expandafter\noexpand
+     \csname the#3\endcsname \noexpand\@thmcountersep \@thmcounter{#1}}%
+   \expandafter\xdef\csname #1name\endcsname{#2}%
+   \global\@namedef{#1}{\@spthm{#1}{\csname #1name\endcsname}{#4}{#5}}%
+                              \global\@namedef{end#1}{\@endtheorem}}}
+
+\def\@spynthm#1#2#3#4{\expandafter\@ifdefinable\csname #1\endcsname
+   {\@definecounter{#1}%
+   \expandafter\xdef\csname the#1\endcsname{\@thmcounter{#1}}%
+   \expandafter\xdef\csname #1name\endcsname{#2}%
+   \global\@namedef{#1}{\@spthm{#1}{\csname #1name\endcsname}{#3}{#4}}%
+                               \global\@namedef{end#1}{\@endtheorem}}}
+
+\def\@spothm#1[#2]#3#4#5{%
+  \@ifundefined{c@#2}{\@latexerr{No theorem environment `#2' defined}\@eha}%
+  {\expandafter\@ifdefinable\csname #1\endcsname
+  {\global\@namedef{the#1}{\@nameuse{the#2}}%
+  \expandafter\xdef\csname #1name\endcsname{#3}%
+  \global\@namedef{#1}{\@spthm{#2}{\csname #1name\endcsname}{#4}{#5}}%
+  \global\@namedef{end#1}{\@endtheorem}}}}
+
+\def\@spthm#1#2#3#4{\topsep 7\p@ \@plus2\p@ \@minus4\p@
+\labelsep=\spthmsep\refstepcounter{#1}%
+\@ifnextchar[{\@spythm{#1}{#2}{#3}{#4}}{\@spxthm{#1}{#2}{#3}{#4}}}
+
+\def\@spxthm#1#2#3#4{\@spbegintheorem{#2}{\csname the#1\endcsname}{#3}{#4}%
+                    \ignorespaces}
+
+\def\@spythm#1#2#3#4[#5]{\@spopargbegintheorem{#2}{\csname
+       the#1\endcsname}{#5}{#3}{#4}\ignorespaces}
+
+\def\normalthmheadings{\def\@spbegintheorem##1##2##3##4{\trivlist\normalfont
+                 \item[\hskip\labelsep{##3##1\ ##2\@thmcounterend}]##4}
+\def\@spopargbegintheorem##1##2##3##4##5{\trivlist
+      \item[\hskip\labelsep{##4##1\ ##2}]{##4(##3)\@thmcounterend\ }##5}}
+\normalthmheadings
+
+\def\reversethmheadings{\def\@spbegintheorem##1##2##3##4{\trivlist\normalfont
+                 \item[\hskip\labelsep{##3##2\ ##1\@thmcounterend}]##4}
+\def\@spopargbegintheorem##1##2##3##4##5{\trivlist
+      \item[\hskip\labelsep{##4##2\ ##1}]{##4(##3)\@thmcounterend\ }##5}}
+
+% definition of \spnewtheorem* without number
+
+\def\@sthm#1#2{\@Ynthm{#1}{#2}}
+
+\def\@Ynthm#1#2#3#4{\expandafter\@ifdefinable\csname #1\endcsname
+   {\global\@namedef{#1}{\@Thm{\csname #1name\endcsname}{#3}{#4}}%
+    \expandafter\xdef\csname #1name\endcsname{#2}%
+    \global\@namedef{end#1}{\@endtheorem}}}
+
+\def\@Thm#1#2#3{\topsep 7\p@ \@plus2\p@ \@minus4\p@
+\@ifnextchar[{\@Ythm{#1}{#2}{#3}}{\@Xthm{#1}{#2}{#3}}}
+
+\def\@Xthm#1#2#3{\@Begintheorem{#1}{#2}{#3}\ignorespaces}
+
+\def\@Ythm#1#2#3[#4]{\@Opargbegintheorem{#1}
+       {#4}{#2}{#3}\ignorespaces}
+
+\def\@Begintheorem#1#2#3{#3\trivlist
+                           \item[\hskip\labelsep{#2#1\@thmcounterend}]}
+
+\def\@Opargbegintheorem#1#2#3#4{#4\trivlist
+      \item[\hskip\labelsep{#3#1}]{#3(#2)\@thmcounterend\ }}
+
+% initialize theorem environment
+
+\if@envcntsect
+   \def\@thmcountersep{.}
+   \spnewtheorem{theorem}{Theorem}[section]{\bfseries}{\itshape}
+\else
+   \spnewtheorem{theorem}{Theorem}{\bfseries}{\itshape}
+   \if@envcntreset
+      \@addtoreset{theorem}{section}
+   \else
+      \@addtoreset{theorem}{chapter}
+   \fi
+\fi
+
+%definition of divers theorem environments
+\spnewtheorem*{claim}{Claim}{\itshape}{\rmfamily}
+\spnewtheorem*{proof}{Proof}{\itshape}{\rmfamily}
+\if@envcntsame % all environments like "Theorem" - using its counter
+   \def\spn@wtheorem#1#2#3#4{\@spothm{#1}[theorem]{#2}{#3}{#4}}
+\else % all environments with their own counter
+   \if@envcntsect % show section counter
+      \def\spn@wtheorem#1#2#3#4{\@spxnthm{#1}{#2}[section]{#3}{#4}}
+   \else % not numbered with section
+      \if@envcntreset
+         \def\spn@wtheorem#1#2#3#4{\@spynthm{#1}{#2}{#3}{#4}
+                                   \@addtoreset{#1}{section}}
+      \else
+         \let\spn@wtheorem=\@spynthm
+      \fi
+   \fi
+\fi
+%
+\let\spdefaulttheorem=\spn@wtheorem
+%
+\spn@wtheorem{case}{Case}{\itshape}{\rmfamily}
+\spn@wtheorem{conjecture}{Conjecture}{\itshape}{\rmfamily}
+\spn@wtheorem{corollary}{Corollary}{\bfseries}{\itshape}
+\spn@wtheorem{definition}{Definition}{\bfseries}{\rmfamily}
+\spn@wtheorem{example}{Example}{\itshape}{\rmfamily}
+\spn@wtheorem{exercise}{Exercise}{\bfseries}{\rmfamily}
+\spn@wtheorem{lemma}{Lemma}{\bfseries}{\itshape}
+\spn@wtheorem{note}{Note}{\itshape}{\rmfamily}
+\spn@wtheorem{problem}{Problem}{\bfseries}{\rmfamily}
+\spn@wtheorem{property}{Property}{\itshape}{\rmfamily}
+\spn@wtheorem{proposition}{Proposition}{\bfseries}{\itshape}
+\spn@wtheorem{question}{Question}{\itshape}{\rmfamily}
+\spn@wtheorem{solution}{Solution}{\bfseries}{\rmfamily}
+\spn@wtheorem{remark}{Remark}{\itshape}{\rmfamily}
+%
+\newenvironment{theopargself}
+    {\def\@spopargbegintheorem##1##2##3##4##5{\trivlist
+         \item[\hskip\labelsep{##4##1\ ##2}]{##4##3\@thmcounterend\ }##5}
+     \def\@Opargbegintheorem##1##2##3##4{##4\trivlist
+         \item[\hskip\labelsep{##3##1}]{##3##2\@thmcounterend\ }}}{}
+\newenvironment{theopargself*}
+    {\def\@spopargbegintheorem##1##2##3##4##5{\trivlist
+         \item[\hskip\labelsep{##4##1\ ##2}]{\hspace*{-\labelsep}##4##3\@thmcounterend}##5}
+     \def\@Opargbegintheorem##1##2##3##4{##4\trivlist
+         \item[\hskip\labelsep{##3##1}]{\hspace*{-\labelsep}##3##2\@thmcounterend}}}{}
+%
+\fi
+
+\def\@takefromreset#1#2{%
+    \def\@tempa{#1}%
+    \let\@tempd\@elt
+    \def\@elt##1{%
+        \def\@tempb{##1}%
+        \ifx\@tempa\@tempb\else
+            \@addtoreset{##1}{#2}%
+        \fi}%
+    \expandafter\expandafter\let\expandafter\@tempc\csname cl@#2\endcsname
+    \expandafter\def\csname cl@#2\endcsname{}%
+    \@tempc
+    \let\@elt\@tempd}
+
+\def\squareforqed{\hbox{\rlap{$\sqcap$}$\sqcup$}}
+\def\qed{\ifmmode\else\unskip\quad\fi\squareforqed}
+\def\smartqed{\def\qed{\ifmmode\squareforqed\else{\unskip\nobreak\hfil
+\penalty50\hskip1em\null\nobreak\hfil\squareforqed
+\parfillskip=0pt\finalhyphendemerits=0\endgraf}\fi}}
+
+% Define `abstract' environment
+\def\abstract{\topsep=0pt\partopsep=0pt\parsep=0pt\itemsep=0pt\relax
+\trivlist\item[\hskip\labelsep
+{\bfseries\abstractname}]\if!\abstractname!\hskip-\labelsep\fi}
+\if@twocolumn
+% \if@avier
+%   \def\endabstract{\endtrivlist\addvspace{5mm}\strich}
+%   \def\strich{\hrule\vskip1ptplus12pt}
+% \else
+    \def\endabstract{\endtrivlist\addvspace{3mm}}
+% \fi
+\else
+\fi
+%
+\newenvironment{verse}
+               {\let\\\@centercr
+                \list{}{\itemsep      \z@
+                        \itemindent   -1.5em%
+                        \listparindent\itemindent
+                        \rightmargin  \leftmargin
+                        \advance\leftmargin 1.5em}%
+                \item\relax}
+               {\endlist}
+\newenvironment{quotation}
+               {\list{}{\listparindent 1.5em%
+                        \itemindent    \listparindent
+                        \rightmargin   \leftmargin
+                        \parsep        \z@ \@plus\p@}%
+                \item\relax}
+               {\endlist}
+\newenvironment{quote}
+               {\list{}{\rightmargin\leftmargin}%
+                \item\relax}
+               {\endlist}
+\newcommand\appendix{\par\small
+  \setcounter{section}{0}%
+  \setcounter{subsection}{0}%
+  \renewcommand\thesection{\@Alph\c@section}}
+\setlength\arraycolsep{1.5\p@}
+\setlength\tabcolsep{6\p@}
+\setlength\arrayrulewidth{.4\p@}
+\setlength\doublerulesep{2\p@}
+\setlength\tabbingsep{\labelsep}
+\skip\@mpfootins = \skip\footins
+\setlength\fboxsep{3\p@}
+\setlength\fboxrule{.4\p@}
+\renewcommand\theequation{\@arabic\c@equation}
+\newcounter{figure}
+\renewcommand\thefigure{\@arabic\c@figure}
+\def\fps@figure{tbp}
+\def\ftype@figure{1}
+\def\ext@figure{lof}
+\def\fnum@figure{\figurename~\thefigure}
+\newenvironment{figure}
+               {\@float{figure}}
+               {\end@float}
+\newenvironment{figure*}
+               {\@dblfloat{figure}}
+               {\end@dblfloat}
+\newcounter{table}
+\renewcommand\thetable{\@arabic\c@table}
+\def\fps@table{tbp}
+\def\ftype@table{2}
+\def\ext@table{lot}
+\def\fnum@table{\tablename~\thetable}
+\newenvironment{table}
+               {\@float{table}}
+               {\end@float}
+\newenvironment{table*}
+               {\@dblfloat{table}}
+               {\end@dblfloat}
+%
+\def \@floatboxreset {%
+        \reset@font
+        \small
+        \@setnobreak
+        \@setminipage
+}
+%
+\newcommand{\tableheadseprule}{\noalign{\hrule height.375mm}}
+%
+\newlength\abovecaptionskip
+\newlength\belowcaptionskip
+\setlength\abovecaptionskip{10\p@}
+\setlength\belowcaptionskip{0\p@}
+\newcommand\leftlegendglue{}
+
+\def\fig@type{figure}
+
+\newdimen\figcapgap\figcapgap=3pt
+\newdimen\tabcapgap\tabcapgap=5.5pt
+
+\@ifundefined{floatlegendstyle}{\def\floatlegendstyle{\bfseries}}{}
+
+\long\def\@caption#1[#2]#3{\par\addcontentsline{\csname
+  ext@#1\endcsname}{#1}{\protect\numberline{\csname
+  the#1\endcsname}{\ignorespaces #2}}\begingroup
+    \@parboxrestore
+    \@makecaption{\csname fnum@#1\endcsname}{\ignorespaces #3}\par
+  \endgroup}
+
+\def\capstrut{\vrule\@width\z@\@height\topskip}
+
+\@ifundefined{captionstyle}{\def\captionstyle{\normalfont\small}}{}
+
+\long\def\@makecaption#1#2{%
+ \captionstyle
+ \ifx\@captype\fig@type
+   \vskip\figcapgap
+ \fi
+ \setbox\@tempboxa\hbox{{\floatlegendstyle #1\floatcounterend}%
+ \capstrut #2}%
+ \ifdim \wd\@tempboxa >\hsize
+   {\floatlegendstyle #1\floatcounterend}\capstrut #2\par
+ \else
+   \hbox to\hsize{\leftlegendglue\unhbox\@tempboxa\hfil}%
+ \fi
+ \ifx\@captype\fig@type\else
+   \vskip\tabcapgap
+ \fi}
+
+\newdimen\figgap\figgap=1cc
+\long\def\@makesidecaption#1#2{%
+   \parbox[b]{\@tempdimb}{\captionstyle{\floatlegendstyle
+                                         #1\floatcounterend}#2}}
+\def\sidecaption#1\caption{%
+\setbox\@tempboxa=\hbox{#1\unskip}%
+\if@twocolumn
+ \ifdim\hsize<\textwidth\else
+   \ifdim\wd\@tempboxa<\columnwidth
+      \typeout{Double column float fits into single column -
+            ^^Jyou'd better switch the environment. }%
+   \fi
+ \fi
+\fi
+\@tempdimb=\hsize
+\advance\@tempdimb by-\figgap
+\advance\@tempdimb by-\wd\@tempboxa
+\ifdim\@tempdimb<3cm
+    \typeout{\string\sidecaption: No sufficient room for the legend;
+             using normal \string\caption. }%
+   \unhbox\@tempboxa
+   \let\@capcommand=\@caption
+\else
+   \let\@capcommand=\@sidecaption
+   \leavevmode
+   \unhbox\@tempboxa
+   \hfill
+\fi
+\refstepcounter\@captype
+\@dblarg{\@capcommand\@captype}}
+
+\long\def\@sidecaption#1[#2]#3{\addcontentsline{\csname
+  ext@#1\endcsname}{#1}{\protect\numberline{\csname
+  the#1\endcsname}{\ignorespaces #2}}\begingroup
+    \@parboxrestore
+    \@makesidecaption{\csname fnum@#1\endcsname}{\ignorespaces #3}\par
+  \endgroup}
+
+% Define `acknowledgement' environment
+\def\acknowledgement{\par\addvspace{17pt}\small\rmfamily
+\trivlist\if!\ackname!\item[]\else
+\item[\hskip\labelsep
+{\bfseries\ackname}]\fi}
+\def\endacknowledgement{\endtrivlist\addvspace{6pt}}
+\newenvironment{acknowledgements}{\begin{acknowledgement}}
+{\end{acknowledgement}}
+% Define `noteadd' environment
+\def\noteadd{\par\addvspace{17pt}\small\rmfamily
+\trivlist\item[\hskip\labelsep
+{\itshape\noteaddname}]}
+\def\endnoteadd{\endtrivlist\addvspace{6pt}}
+
+\DeclareOldFontCommand{\rm}{\normalfont\rmfamily}{\mathrm}
+\DeclareOldFontCommand{\sf}{\normalfont\sffamily}{\mathsf}
+\DeclareOldFontCommand{\tt}{\normalfont\ttfamily}{\mathtt}
+\DeclareOldFontCommand{\bf}{\normalfont\bfseries}{\mathbf}
+\DeclareOldFontCommand{\it}{\normalfont\itshape}{\mathit}
+\DeclareOldFontCommand{\sl}{\normalfont\slshape}{\@nomath\sl}
+\DeclareOldFontCommand{\sc}{\normalfont\scshape}{\@nomath\sc}
+\DeclareRobustCommand*\cal{\@fontswitch\relax\mathcal}
+\DeclareRobustCommand*\mit{\@fontswitch\relax\mathnormal}
+\newcommand\@pnumwidth{1.55em}
+\newcommand\@tocrmarg{2.55em}
+\newcommand\@dotsep{4.5}
+\setcounter{tocdepth}{1}
+\newcommand\tableofcontents{%
+    \section*{\contentsname}%
+    \@starttoc{toc}%
+    \addtocontents{toc}{\begingroup\protect\small}%
+    \AtEndDocument{\addtocontents{toc}{\endgroup}}%
+    }
+\newcommand*\l@part[2]{%
+  \ifnum \c@tocdepth >-2\relax
+    \addpenalty\@secpenalty
+    \addvspace{2.25em \@plus\p@}%
+    \begingroup
+      \setlength\@tempdima{3em}%
+      \parindent \z@ \rightskip \@pnumwidth
+      \parfillskip -\@pnumwidth
+      {\leavevmode
+       \large \bfseries #1\hfil \hb@xt@\@pnumwidth{\hss #2}}\par
+       \nobreak
+       \if@compatibility
+         \global\@nobreaktrue
+         \everypar{\global\@nobreakfalse\everypar{}}%
+      \fi
+    \endgroup
+  \fi}
+\newcommand*\l@section{\@dottedtocline{1}{0pt}{1.5em}}
+\newcommand*\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}}
+\newcommand*\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}}
+\newcommand*\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}}
+\newcommand*\l@subparagraph{\@dottedtocline{5}{10em}{5em}}
+\newcommand\listoffigures{%
+    \section*{\listfigurename
+      \@mkboth{\listfigurename}%
+              {\listfigurename}}%
+    \@starttoc{lof}%
+    }
+\newcommand*\l@figure{\@dottedtocline{1}{1.5em}{2.3em}}
+\newcommand\listoftables{%
+    \section*{\listtablename
+      \@mkboth{\listtablename}{\listtablename}}%
+    \@starttoc{lot}%
+    }
+\let\l@table\l@figure
+\newdimen\bibindent
+\setlength\bibindent{\parindent}
+\def\@biblabel#1{#1.}
+\def\@lbibitem[#1]#2{\item[{[#1]}\hfill]\if@filesw
+      {\let\protect\noexpand
+       \immediate
+       \write\@auxout{\string\bibcite{#2}{#1}}}\fi\ignorespaces}
+\newenvironment{thebibliography}[1]
+     {\section*{\refname
+        \@mkboth{\refname}{\refname}}\small
+      \list{\@biblabel{\@arabic\c@enumiv}}%
+           {\settowidth\labelwidth{\@biblabel{#1}}%
+            \leftmargin\labelwidth
+            \advance\leftmargin\labelsep
+            \@openbib@code
+            \usecounter{enumiv}%
+            \let\p@enumiv\@empty
+            \renewcommand\theenumiv{\@arabic\c@enumiv}}%
+      \sloppy\clubpenalty4000\widowpenalty4000%
+      \sfcode`\.\@m}
+     {\def\@noitemerr
+       {\@latex@warning{Empty `thebibliography' environment}}%
+      \endlist}
+%
+\newcount\@tempcntc
+\def\@citex[#1]#2{\if@filesw\immediate\write\@auxout{\string\citation{#2}}\fi
+  \@tempcnta\z@\@tempcntb\m@ne\def\@citea{}\@cite{\@for\@citeb:=#2\do
+    {\@ifundefined
+       {b@\@citeb}{\@citeo\@tempcntb\m@ne\@citea\def\@citea{,}{\bfseries
+        ?}\@warning
+       {Citation `\@citeb' on page \thepage \space undefined}}%
+    {\setbox\z@\hbox{\global\@tempcntc0\csname b@\@citeb\endcsname\relax}%
+     \ifnum\@tempcntc=\z@ \@citeo\@tempcntb\m@ne
+       \@citea\def\@citea{,\hskip0.1em\ignorespaces}\hbox{\csname b@\@citeb\endcsname}%
+     \else
+      \advance\@tempcntb\@ne
+      \ifnum\@tempcntb=\@tempcntc
+      \else\advance\@tempcntb\m@ne\@citeo
+      \@tempcnta\@tempcntc\@tempcntb\@tempcntc\fi\fi}}\@citeo}{#1}}
+\def\@citeo{\ifnum\@tempcnta>\@tempcntb\else
+            \@citea\def\@citea{,\hskip0.1em\ignorespaces}%
+  \ifnum\@tempcnta=\@tempcntb\the\@tempcnta\else
+   {\advance\@tempcnta\@ne\ifnum\@tempcnta=\@tempcntb \else \def\@citea{--}\fi
+    \advance\@tempcnta\m@ne\the\@tempcnta\@citea\the\@tempcntb}\fi\fi}
+%
+\newcommand\newblock{\hskip .11em\@plus.33em\@minus.07em}
+\let\@openbib@code\@empty
+\newenvironment{theindex}
+               {\if@twocolumn
+                  \@restonecolfalse
+                \else
+                  \@restonecoltrue
+                \fi
+                \columnseprule \z@
+                \columnsep 35\p@
+                \twocolumn[\section*{\indexname}]%
+                \@mkboth{\indexname}{\indexname}%
+                \thispagestyle{plain}\parindent\z@
+                \parskip\z@ \@plus .3\p@\relax
+                \let\item\@idxitem}
+               {\if@restonecol\onecolumn\else\clearpage\fi}
+\newcommand\@idxitem{\par\hangindent 40\p@}
+\newcommand\subitem{\@idxitem \hspace*{20\p@}}
+\newcommand\subsubitem{\@idxitem \hspace*{30\p@}}
+\newcommand\indexspace{\par \vskip 10\p@ \@plus5\p@ \@minus3\p@\relax}
+
+\if@twocolumn
+ \renewcommand\footnoterule{%
+  \kern-3\p@
+  \hrule\@width\columnwidth
+  \kern2.6\p@}
+\else
+ \renewcommand\footnoterule{%
+  \kern-3\p@
+  \hrule\@width.382\columnwidth
+  \kern2.6\p@}
+\fi
+\newcommand\@makefntext[1]{%
+    \noindent
+    \hb@xt@\bibindent{\hss\@makefnmark\enspace}#1}
+%
+\def\trans@english{\switcht@albion}
+\def\trans@french{\switcht@francais}
+\def\trans@german{\switcht@deutsch}
+\newenvironment{translation}[1]{\if!#1!\else
+\@ifundefined{selectlanguage}{\csname trans@#1\endcsname}{\selectlanguage{#1}}%
+\fi}{}
+% languages
+% English section
+\def\switcht@albion{%\typeout{English spoken.}%
+ \def\abstractname{Abstract}%
+ \def\ackname{Acknowledgements}%
+ \def\andname{and}%
+ \def\lastandname{, and}%
+ \def\appendixname{Appendix}%
+ \def\chaptername{Chapter}%
+ \def\claimname{Claim}%
+ \def\conjecturename{Conjecture}%
+ \def\contentsname{Contents}%
+ \def\corollaryname{Corollary}%
+ \def\definitionname{Definition}%
+ \def\emailname{E-mail}%
+ \def\examplename{Example}%
+ \def\exercisename{Exercise}%
+ \def\figurename{Fig.}%
+ \def\keywordname{{\bfseries Keywords}}%
+ \def\indexname{Index}%
+ \def\lemmaname{Lemma}%
+ \def\contriblistname{List of Contributors}%
+ \def\listfigurename{List of Figures}%
+ \def\listtablename{List of Tables}%
+ \def\mailname{{\itshape Correspondence to\/}:}%
+ \def\noteaddname{Note added in proof}%
+ \def\notename{Note}%
+ \def\partname{Part}%
+ \def\problemname{Problem}%
+ \def\proofname{Proof}%
+ \def\propertyname{Property}%
+ \def\questionname{Question}%
+ \def\refname{References}%
+ \def\remarkname{Remark}%
+ \def\seename{see}%
+ \def\solutionname{Solution}%
+ \def\tablename{Table}%
+ \def\theoremname{Theorem}%
+}\switcht@albion % make English default
+%
+% French section
+\def\switcht@francais{\svlanginfo
+%\typeout{On parle francais.}%
+ \def\abstractname{R\'esum\'e\runinend}%
+ \def\ackname{Remerciements\runinend}%
+ \def\andname{et}%
+ \def\lastandname{ et}%
+ \def\appendixname{Appendice}%
+ \def\chaptername{Chapitre}%
+ \def\claimname{Pr\'etention}%
+ \def\conjecturename{Hypoth\`ese}%
+ \def\contentsname{Table des mati\`eres}%
+ \def\corollaryname{Corollaire}%
+ \def\definitionname{D\'efinition}%
+ \def\emailname{E-mail}%
+ \def\examplename{Exemple}%
+ \def\exercisename{Exercice}%
+ \def\figurename{Fig.}%
+ \def\keywordname{{\bfseries Mots-cl\'e\runinend}}%
+ \def\indexname{Index}%
+ \def\lemmaname{Lemme}%
+ \def\contriblistname{Liste des contributeurs}%
+ \def\listfigurename{Liste des figures}%
+ \def\listtablename{Liste des tables}%
+ \def\mailname{{\itshape Correspondence to\/}:}%
+ \def\noteaddname{Note ajout\'ee \`a l'\'epreuve}%
+ \def\notename{Remarque}%
+ \def\partname{Partie}%
+ \def\problemname{Probl\`eme}%
+ \def\proofname{Preuve}%
+ \def\propertyname{Caract\'eristique}%
+%\def\propositionname{Proposition}%
+ \def\questionname{Question}%
+ \def\refname{Bibliographie}%
+ \def\remarkname{Remarque}%
+ \def\seename{voyez}%
+ \def\solutionname{Solution}%
+%\def\subclassname{{\it Subject Classifications\/}:}%
+ \def\tablename{Tableau}%
+ \def\theoremname{Th\'eor\`eme}%
+}
+%
+% German section
+\def\switcht@deutsch{\svlanginfo
+%\typeout{Man spricht deutsch.}%
+ \def\abstractname{Zusammenfassung\runinend}%
+ \def\ackname{Danksagung\runinend}%
+ \def\andname{und}%
+ \def\lastandname{ und}%
+ \def\appendixname{Anhang}%
+ \def\chaptername{Kapitel}%
+ \def\claimname{Behauptung}%
+ \def\conjecturename{Hypothese}%
+ \def\contentsname{Inhaltsverzeichnis}%
+ \def\corollaryname{Korollar}%
+%\def\definitionname{Definition}%
+ \def\emailname{E-Mail}%
+ \def\examplename{Beispiel}%
+ \def\exercisename{\"Ubung}%
+ \def\figurename{Abb.}%
+ \def\keywordname{{\bfseries Schl\"usselw\"orter\runinend}}%
+ \def\indexname{Index}%
+%\def\lemmaname{Lemma}%
+ \def\contriblistname{Mitarbeiter}%
+ \def\listfigurename{Abbildungsverzeichnis}%
+ \def\listtablename{Tabellenverzeichnis}%
+ \def\mailname{{\itshape Correspondence to\/}:}%
+ \def\noteaddname{Nachtrag}%
+ \def\notename{Anmerkung}%
+ \def\partname{Teil}%
+%\def\problemname{Problem}%
+ \def\proofname{Beweis}%
+ \def\propertyname{Eigenschaft}%
+%\def\propositionname{Proposition}%
+ \def\questionname{Frage}%
+ \def\refname{Literatur}%
+ \def\remarkname{Anmerkung}%
+ \def\seename{siehe}%
+ \def\solutionname{L\"osung}%
+%\def\subclassname{{\it Subject Classifications\/}:}%
+ \def\tablename{Tabelle}%
+%\def\theoremname{Theorem}%
+}
+\newcommand\today{}
+\edef\today{\ifcase\month\or
+  January\or February\or March\or April\or May\or June\or
+  July\or August\or September\or October\or November\or December\fi
+  \space\number\day, \number\year}
+\setlength\columnsep{1.5cc}
+\setlength\columnseprule{0\p@}
+%
+\frenchspacing
+\clubpenalty=10000
+\widowpenalty=10000
+\def\thisbottomragged{\def\@textbottom{\vskip\z@ plus.0001fil
+\global\let\@textbottom\relax}}
+\pagestyle{headings}
+\pagenumbering{arabic}
+\if@twocolumn
+   \twocolumn
+\fi
+%\if@avier
+%   \onecolumn
+%   \setlength{\textwidth}{156mm}
+%   \setlength{\textheight}{226mm}
+%\fi
+\if@referee
+   \makereferee
+\fi
+\flushbottom
+\endinput
+%%
+%% End of file `svjour3.cls'.
Index: /anuga_work/publications/boxing_day_validation_2008/tsunami07.bib
===================================================================
--- /anuga_work/publications/boxing_day_validation_2008/tsunami07.bib	(revision 7209)
+++ /anuga_work/publications/boxing_day_validation_2008/tsunami07.bib	(revision 7210)
@@ -68,5 +68,5 @@
 
 @article{zoppou03,
-  AUTHOR = {C., Zoppou and S.G Roberts},
+  AUTHOR = {Zoppou, C. and Roberts, S.G },
   TITLE = {Explicit schemes for dam-break simulations},
   YEAR = {2003},
@@ -80,5 +80,5 @@
 
 @article{zoppou99,
-  AUTHOR = {C., Zoppou and S.G Roberts},
+  AUTHOR = {Zoppou, C. and Roberts, S.G },
   TITLE = {Catastrophic collapse of water supply reserviours in urban areas},
   YEAR = {1999},
@@ -91,5 +91,5 @@
 
 @article{zoppou00,
-  AUTHOR = {C., Zoppou and S.G Roberts},
+  AUTHOR = {Zoppou, C. and Roberts, S.G },
   TITLE = {Numerical solution of the two-dimensional unsteady dam break},
   YEAR = {2000},
@@ -100,5 +100,5 @@
 
 @ARTICLE{kurganov01,
-  AUTHOR = {A., Kurganov and S. Noelle and G. Petrova},
+  AUTHOR = {Kurganov, A. and Noelle, S.  and Petrova, G. },
   TITLE = {Semidiscrete central-upwind schemes for hyperbolic conservation
 laws and {H}amilton-{J}acobi equations},
@@ -141,5 +141,5 @@
 
 @ARTICLE{synolakis05,
-AUTHOR = {C., Synolakis and E. Okal and E. Bernard},
+AUTHOR = {Synolakis, C. and Okal, E.  and  Bernard, E.},
 TITLE = {The megatsunami of {D}ecember 26 2004},
 YEAR = {2005},
@@ -634,5 +634,5 @@
 
 @ARTICLE{titov97a,
-AUTHOR = {V.V., Titov and F.I. Gonzalez},
+AUTHOR = {Titov, V.V. and  Gonzalez, F.I.},
 TITLE = {Implementation and testing of the method of splitting tsunami ({MOST}) model},
 YEAR = {1997},
@@ -641,5 +641,5 @@
 
 @ARTICLE{titov98,
-AUTHOR = {V.V. Titov and C. Synolakis},
+AUTHOR = {Titov, V.V.  and  Synolakis, C.},
 TITLE = {Numerical Modelling of tidal wave runup},
 YEAR = {1998},
@@ -715,5 +715,5 @@
 
 @Unpublished{roberts06,
-author = {S.G., Roberts and O.M. Nielsen and J.D. Jakeman},
+author = {Roberts, S.G.  and Nielsen, O.M.  and Jakeman, J.D. },
 title = {Simulation of Tsunami and Flash Flood},
 note = {Accepted for publication in the refereed proceedings of the International Conference on High Performance Scientific Computing: Modeling, Simulation and Optimization of Complex Processes, March 6-10, 2006, Hanoi Vietnam},
@@ -732,5 +732,5 @@
 
 @PHDTHESIS{ham06,
-AUTHOR = {D., Ham},
+AUTHOR = {Ham, D.},
 TITLE = {On techniques for modelling coastal and ocean flow with unstructed meshes},
 SCHOOL = {de Technische Universiteit Delft},
@@ -739,5 +739,5 @@
 
 @ARTICLE{anastasiou97,
-AUTHOR = {K., Anastasiou and C.T. Chan},
+AUTHOR = {Anastasiou, K. and Chan, C.T. },
 TITLE = {Solution of the 2D shallow water equations using the finite volume method on unstructured triangular meshes},
 JOURNAL = {International Journal for Numerical Methods in Fluids},
@@ -749,5 +749,5 @@
 
 @TECHREPORT{greensdale07,
-AUTHOR = {D., Greensdale and M . Simanjuntak and D. Burbidge and J. Chittleborough},
+AUTHOR = {Greensdale, D. and Simanjuntak, M .  and Burbidge, D.  and Chittleborough, J. },
 TITLE = {A first-generation real-time tsunami forecasting system for the {A}ustralian region},
 INSTITUTION = {Bureau of Meteorology Australia},
@@ -758,5 +758,5 @@
 
 @ARTICLE{titov97b,
-AUTHOR = {V.V., Titov and C.E. Synolakis},
+AUTHOR = {Titov, V.V. and Synolakis, C.E.},
 TITLE = {Extreme inundation flows during the Hokkaido-
 Nansei-Oki tsunami},
@@ -779,5 +779,5 @@
 
 @ARTICLE{bourgeois99,
-AUTHOR = {J., Bourgeois and C. Petroff and H. Yeh and V. Titov and C. Synolakis and B. Benson and J. Kuroiwa and J. Lander and E. Norabuena},
+AUTHOR = {Bourgeois, J.,  and Petroff, C.  and Yeh, H.  and Titov, V.  and Synolakis, C.  and Benson, B.  and Kuroiwa, J.  and Lander, J.  and Norabuena, E. },
 TITLE = {Geologic setting, field survey and modeling of the Chimbote,
 northern Peru, tsunami of 21 February 1996},
@@ -790,5 +790,5 @@
 
 @ARTICLE{gusiakov72,
-AUTHOR = {V.K., Gusiakov},
+AUTHOR = {Gusiakov, V.K. },
 TITLE = {Static displacement on the surface of an elastic space. Ill-posed problems of
 mathematical physics and interpretation of geophysical data},
@@ -800,5 +800,5 @@
 
 @ARTICLE{grilli06,
-AUTHOR = {S.T., Grilli and M. Ioualalen and J. Asavanant and F. Shi and J.T Kirby and P. Watts},
+AUTHOR = {Grilli, S.T. and Ioualalen, M.  and Asavanant, J.  and Shi, F.  and Kirby, J.T  and Watts, P. },
 TITLE = {Source Constraints and Model Simulation of the December 26, 2004 Indian
 Ocean Tsunami},
@@ -809,5 +809,5 @@
 
 @ARTICLE{titov05,
-AUTHOR = {V., Titov and F. Gonz\'{a}lez and E. Bernard and M. Eble and H. Mofjeld and J. Newman and A. Venturato},
+AUTHOR = {Titov, V.V. and Gonz\'{a}lez, F.  and Bernard, E.  and Eble, M. and  Mofjeld, H. and Newman, J.  and Venturato, A. },
 TITLE = {Real-Time Tsunami Forecasting: Challenges and Solutions},
 JOURNAL = {Natural Hazards},
@@ -819,5 +819,5 @@
 
 @ARTICLE{titov95,
-AUTHOR = {V., Titov and C. Synolakis},
+AUTHOR = {Titov, V.  and Synolakis, C. },
 TITLE = {Modeling of breaking and nonbreaking long wave
 evolution and runup using {VTCS}-2},
@@ -830,5 +830,5 @@
 
 @Article{vigny05,
-author = {C., Vigny and W.J.F. Simons and S. Abu and R. Bamphenyu and C. Satirapod, N. and Choosakul, C. Subarya and A. Socquet and K. Omar and H.Z. Abidin and B.A.C. Ambrosius},
+author = {Vigny, C.,  and Simons, W.J.F.  and Abu, S.  and Bamphenyu, R.  and Satirapod, C.  and Choosakul, N. and Subarya, C.  and  Socquet, A. and Omar, K.  and  Abidin, H.Z. and  Ambrosius, B.A.C.},
 title = {Insight into the 2004 {S}umatra-{A}ndaman earthquake from {GPS} measurements in southeast {A}sia},
 journal = {Nature},
@@ -839,5 +839,5 @@
 
 @Article{ammon05,
-author = {C.J., Ammon and C. Ji and H. Thio and D. Robinson and Sidao Ni and V. Hjorleifsdottir and H. and T. Lay and S. Das and D. Helmberger and G. Ichinose and J. Polet and D. Wald},
+author = {Ammon, C.J. and Ji, C. and Thio, H.  and Robinson, D.  and  Ni, S. and Hjorleifsdottir, V.  and Lay, H. and Lay. T. and Das, S.  and Helmberger, D.  and Ichinose, G.  and Polet, J.  and  Wald, D.},
 title = {Rupture Process of the 2004 Sumatra-Andaman Earthquake},
 journal = {Science},
@@ -848,5 +848,5 @@
 
 @TechReport{kawata05,
-author = {T. et XIV alia, Kawata},
+author = {Kawata, T. et XIV alia},
 title = {Comprehensive analysis of the damage and its impact on coastal zones by the 2004 {I}ndian {O}cean tsunami disaster},
 institution = {Disaster Prevention Research Institute},
@@ -856,5 +856,5 @@
 
 @Article{liu05,
-author = {P. L.-F., Liu and P. Lynett and H. Fernando and B.E. Jaffe and H. Fritz and B. Higman and R. Morton and J. Goff and C. Synolakis},
+author = {Liu, P. L.-F.,  and Lynett, P.  and Fernando, H.  and Jaffe, B.E.  and Fritz, H.  and Higman, B.  and Morton, R.  and Goff, J.  and Synolakis, C. },
 title = {Observations by the International Tsunami Survey Team in {S}ri {L}anka},
 journal = {Science},
@@ -865,5 +865,5 @@
 
 @Article{merrifield05,
-author = {M.A. et XXIII alia, Merrifield},
+author = {Merrifield, M.A. et XXIII alia},
 title = {Tide gauge observations of the {I}ndian {O}cean tsunami, {D}ecember 26, 2004},
 journal = {Geophysical Research Letters},
@@ -874,5 +874,5 @@
 
 @Article{roberts00,
-author = {S.G, Roberts and C. Zoppou},
+author = {Roberts, S.G and Zoppou, C. },
 title = {Robust and efficent solution of the 2D shallow water wave equation with domains containg dry beds},
 journal = {The ANZIAM Journal},
@@ -883,5 +883,5 @@
 }
 
-@Misc{ngdc,
+@Misc{ngdc, 
 author = {National Geophysical Data Center (NGDC)},
 title = {Indian Ocean December 26, 2004: Runups},
@@ -890,10 +890,10 @@
 
 @InProceedings{nielsen05,
-author = {O.M, Nielsen and  S.G Roberts and D. Gray and A. McPherson and A. Hitchman},
+author = {Nielsen, O.M  and  Roberts, S.G  and Gray, D.  and McPherson, A.  and Hitchman, A. },
 title = {Hydrodynamic modelling of coastal inundation},
 booktitle = {MODSIM 2005 International Congress on Modelling and Simulation},
 pages = {518-523},
 year = {2005},
-editor = {A. Zerger and R.M. Argent},
+editor = {Zerger, A.  and Argent, R.M. },
 Tmonth = {December},
 organization = {Modelling and Simulation Society of Australia and New Zealand},
@@ -902,5 +902,5 @@
 
 @Article{wei95,
-author = {G. Wei and J.T. Kirby and S.T. Grilli and R. Subramanya},
+author = {Wei, G.  and Kirby, J.T.  and Grilli, S.T.  and Subramanya, R. },
 title = {A fully nonlinear {B}oussinesq model for free surface waves. {P}art 1: Highly nonlinear unsteady waves},
 journal = {Journal of Fluid Mechanics},
@@ -932,5 +932,5 @@
 @article{satake95,
 title={Linear and nonlinear computations of the 1992 Nicaragua earthquake tsunami},
-author={Sataka K.},
+author={Satake, K.},
 journal={Pure and Applied Geophysics},
 year={1995},
@@ -941,5 +941,5 @@
 
 @Misc{schoettle2007,
-  author = 	 {Schoettle E. and Sakimoto S.},
+  author = 	 {Schoettle, E. and Sakimoto, S.},
   title = 	 {Modeling the Effects of Coral Reef Health on Tsunami Run-up 
                   with the Finite-element Model ADCIRC},
@@ -960,5 +960,5 @@
 
 @article{grilli07,
-author = {St\'{e}phan T. Grilli and Mansour Ioualalen and Jack Asavanant and Fengyan Shi and James T. Kirby and Philip Watts},
+author = {Grilli, S.T. and Ioualalen, M.  and Asavanant, J.  and Shi, F.  and  Kirby, J.T. and Watts, P.},
 title = {Source Constraints and Model Simulation of the December 26, 2004, Indian Ocean Tsunami},
 publisher = {ASCE},
@@ -1065,5 +1065,5 @@
 title = {New, improved version of Generic Mapping Tools released},
 journal = {EOS trans. AGU},
-year = {1998},
+year = {199  8},
 OPTvolume = {79},
 OPTpages = {579},
@@ -1091,9 +1091,9 @@
 
 @Article{fritz06,
-author = {Fritz, Hermann M. and Jose C. Borrero and Costas E. Synolakis and Jeseon Yoo.},
+author = {Fritz, H.~M. and Borrero, J.~C.  and Synolakis, C.~E. and  Yoo., J.},
 title = {2004 Indian Ocean tsunami flow velocity measurements from survivor videos},
 journal = {Geophys. Res. Lett.},
 year = {2006},
-number = {L24605},
+volume = {L24605},
 pages = {1--5},
 OPTdoi = {10.1029/2006GL026784},
