Changeset 4062


Ignore:
Timestamp:
Dec 6, 2006, 5:15:13 PM (17 years ago)
Author:
sexton
Message:

updates for usermanual - cairns demo

Location:
anuga_core/documentation/user_manual
Files:
2 added
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/documentation/user_manual/anuga_user_manual.tex

    r4040 r4062  
    979979supporting data is found in the ASCII grid, \code{cairns.asc}, which
    980980has been sourced from the publically available Australian Bathymetry
    981 and Topography Grid 2005, \cite{grid250}.
     981and Topography Grid 2005, \cite{grid250}. The required resolution
     982for inundation modelling will depend on the underlying topography and
     983bathymetry; as the terrain becomes more complex, the desired resolution
     984would decrease to the order of tens of metres.
    982985
    983986\begin{figure}[hbt]
     
    10961099The function returns the water displacement for all \code{x} and
    10971100\code{y} in the domain. The water displacement is a double Gaussian
    1098 function that depends on the characteristics of the slump (length,
    1099 thickness, slope, etc), its location (origin) and the depth at that
    1100 location.
    1101 
     1101function that depends on the characteristics of the slide (length,
     1102width, thickness, slope, etc), its location (origin) and the depth at that
     1103location. For this example, we choose to apply the slide function
     1104at a specified time into the simulation.
    11021105
    11031106\subsubsection{Friction}
     
    11611164
    11621165With the basics established, the running of the `evolve' step is
    1163 very similar to the corresponding step in \file{runup.py}. Here,
    1164 the simulation is run for five hours (10000 seconds) with
    1165 the output stored every two minutes (120 seconds).
     1166very similar to the corresponding step in \file{runup.py}. For the slide
     1167scenario,
     1168the simulation is run for 5000 seconds with the output stored every five seconds.
     1169For this example, we choose to apply the slide at 60 seconds into the simulation.
    11661170
    11671171{\small \begin{verbatim}
    11681172    import time t0 = time.time()
    11691173
    1170     for t in domain.evolve(yieldstep = 120, duration = 10000):
    1171         print domain.timestepping_statistics()
    1172         print domain.boundary_statistics(tags = 'bottom')
    1173 
    1174     print 'That took %.2f seconds' %(time.time()
     1174   
     1175    for t in domain.evolve(yieldstep = 5, finaltime = 60):
     1176            domain.write_time()
     1177            domain.write_boundary_statistics(tags = 'ocean_east')     
     1178           
     1179        # add slide
     1180        thisstagestep = domain.get_quantity('stage')
     1181        if allclose(t, 60):
     1182            slide = Quantity(domain)
     1183            slide.set_values(tsunami_source)
     1184            domain.set_quantity('stage', slide + thisstagestep)
     1185   
     1186        # save every two mins leading up to wave approaching land
     1187        for t in domain.evolve(yieldstep = 5, finaltime = 5000,
     1188                               skip_initial_step = True):
     1189            domain.write_time()
     1190        domain.write_boundary_statistics(tags = 'ocean_east')
     1191\end{verbatim}}
     1192
     1193For the fixed wave scenario, the simulation is run to 10000 seconds,
     1194with the first half of the simulation stored at two minute intervals,
     1195and the second half of the simulation stored at ten second intervals.
     1196This functionality is especially convenient as it allows the detailed
     1197parts of the simulation to be viewed at higher time resolution.
     1198
     1199
     1200{\small \begin{verbatim}
     1201
     1202# save every two mins leading up to wave approaching land
     1203    for t in domain.evolve(yieldstep = 120, finaltime = 5000):
     1204        domain.write_time()
     1205        domain.write_boundary_statistics(tags = 'ocean_east')     
     1206
     1207    # save every 30 secs as wave starts inundating ashore
     1208    for t in domain.evolve(yieldstep = 10, finaltime = 10000,
     1209                           skip_initial_step = True):
     1210        domain.write_time()
     1211        domain.write_boundary_statistics(tags = 'ocean_east')
     1212       
    11751213\end{verbatim}}
    11761214
     
    11801218As described earlier, the user may run animate to view a three-dimensional representation
    11811219of the simulation.
     1220
     1221The user may also be interested in a maximum inundation map. This simply shows the
     1222maximum water depth over the domain and is achieved with the function sww2dem (described in
     1223Section \label{sec:basicfileconversions}).
     1224\file{ExportResults.py} demonstrates how this function can be used:
     1225
     1226\verbatiminput{demos/cairns/ExportResults.py}
     1227
     1228The script generates an maximum water depth ASCII grid at a defined
     1229resolution (here 100 m$^2$) which can then be viewed in a GIS environment, for
     1230example. The parameters used in the function are defined in \file{project.py}.
     1231Figures \ref{fig:maxdepthcairnsslide} and \ref{fig:maxdepthcairnsfixedwave} show
     1232the maximum water depth within the defined region for the slide and fixed wave scenario
     1233respectively.
     1234The user could develop a maximum absolute momentum or other expressions which can be
     1235derived from the quantities.
     1236
     1237\begin{figure}[hbt]
     1238\centerline{\includegraphics[scale=0.5]{graphics/slidedepth.jpg}}
     1239\caption{Maximum inundation map for the Cairns side scenario.}
     1240\label{fig:maxdepthcairnsslide}
     1241\end{figure}
     1242
     1243\begin{figure}[hbt]
     1244\centerline{\includegraphics[scale=0.5]{graphics/fixedwavedepth.jpg}}
     1245\caption{Maximum inundation map for the Cairns fixed wave scenario.}
     1246\label{fig:maxdepthcairnsfixedwave}
     1247\end{figure}
    11821248
    11831249The user may also be interested in interrogating the solution at a particular spatial
     
    12061272verbatiminput{demos/cairns/GetTimeseries.py}
    12071273
    1208 The user may also be interested in a maximum inundation map. This simply shows the
    1209 maximum water depth over the domain and is achieved with the function sww2dem (described in
    1210 Section \label{sec:basicfileconversions}).
    1211 \file{ExportResults.py} demonstrates how this function can be used:
    1212 
    1213 \verbatiminput{demos/cairns/ExportResults.py}
    1214 
    1215 Here, the maximum water depth is shown for the region defined in \file{project.py} and
    1216 is shown in Figures \ref{fig:maxdepthcairnsslump} and \ref{fig:maxdepthcairnsfixedwave}.
    1217 The user could develop a maximum absolute momentum or other expressions which can be
    1218 derived from the quantities.
     1274Here, the time series for the quantities stage and speed will be generated for
     1275each gauge defined in the gauge file defined in file{\project.py}. As an example output,
     1276Figure \ref{fig:trinitybeach} shows the time series for the quantity stage for the
     1277Trinity Beach location.
    12191278
    12201279\begin{figure}[hbt]
    12211280%\centerline{\includegraphics[scale=0.5]{graphics/}}
    1222 \caption{Maximum inundation map for the Cairns slump scenario.}
    1223 \label{fig:maxdepthcairnsslump}
    1224 \end{figure}
    1225 
    1226 \begin{figure}[hbt]
    1227 %\centerline{\includegraphics[scale=0.5]{graphics/}}
    1228 \caption{Maximum inundation map for the Cairns fixed wave scenario.}
    1229 \label{fig:maxdepthcairnsfixedwave}
     1281\caption{Time series information of the quantity stage for the Trinity Beach location.}
     1282\label{fig:trinitybeach}
    12301283\end{figure}
    12311284
Note: See TracChangeset for help on using the changeset viewer.