Section 11 Outputs

Prior to carrying out a simulation the desired output needs to be configured and customised. For viewing outputs, see Chapter 15. There are also a variety of post-processing utilities available, see Chapter 17.

The chapter includes discussion on:

  • Setting the folder location of where outputs are written (Section 11.1).
  • Configuring map output options including 1D results in a 2D form and the powerful output zones feature (Section 11.2).
  • Options for setting up time-series output from 1D and 2D domains at specified locations (Section 11.3).
  • Specialised outputs such as time and depth of inundation along evacuation routes (Section 11.4).
  • Options for customising the check and log file outputs (Section 11.5).
  • List of commands available for customising all the different output options, including output to the simulation console window, check files and folder locations (Section 11.6).

11.1 Output Location

To set the location of where the outputs are written to, the Output Folder command is used. For example “Output Folder == <folder>” redirects all TUFLOW output data except the log and summary files to the specified folder.

Outputs can be written in a different location to where the model sits. It is typically recommended to write outputs to your local C: or D: drive instead of filling up the network drive, or to keep results separate to the input data. A URL can be used (e.g. \myserver), which is useful for running simulations on other computers, but with the output directed to your local drive (your drive will need to be shared) or other location.

The Output Drive command can be used to change the drive letter of any output files with a full path specified. For example “Output Drive == K” changes any output filepaths that have a drive letter to the K drive.

11.2 Map Outputs

TUFLOW has a range of commands that allows the user to control the formats, data types, frequency, sub-areas (output zones) or resolution (SGS high-resolution) of a simulation’s output. The options are highly flexible and can be customised to be different for different output formats and/or output zones.

For example, time-based map output may be produced using one format (e.g. XMDF), and the peak flood level surface for another output (e.g. TIF). Output over a township maybe at a higher output interval frequency than for the whole model so as to produce a higher quality animation of flooding through the township, or if using SGS, at a higher resolution (e.g. HRTIF).

The main commands to control map outputs are:

For example, application of the above commands may look like:

Output Folder == ..\results\ ! Location of the output files
Map Output Format == XMDF TIF ! Produce map output in XMDF and TIF formats
Map Output Data Types == h V d ! Output water level, velocities and depths
Start Map Output == 10 ! Start map output at 10 hours
Map Output Interval (s) == 300 ! Write map output every 5 minutes (300s)
TIF Map Output Interval (s) == 0 ! Output maximums only for TIF format

11.2.1 Output Time Controls

To control when and how frequent outputs are written, time control commands are used. They can be set as a global setting or customised for each output format. A variety of formats can be output for a single simulation, the supported formats are discussed in Section 11.2.2.

To set default settings the followings commands can be used:

To customise the time control commands to be different to the default or global settings, the map output format acronym is included at the start of the command. For example, “XMDF Map Output Interval == 360” will set the output frequency to 360s just for the XMDF format. The format acronym must be identical to that used in Map Output Format. The following commands can be applied in this manner.

To apply the same setting to more than one format, the command needs to be repeated for each format. If a format is not customised the default setting or the setting applied to the whole of the model is used. The order of commands is important, ensure that commands defining the map output settings for all format types are read in prior to any commands specific to a certain format type. For example:

  • The following commands set a Map Output Interval of 120 seconds for the XMDF format and 3600 seconds for the TIF format:
    Map Output Format == XMDF TIF
    Map Output Interval == 120
    TIF Map Output Interval == 3600

  • Swapping the order of the second and third command lines will set a Map Output Interval of 120 seconds for all output formats as the third command overwrites the TIF Map Output Interval command.
    Map Output Format == XMDF TIF
    TIF Map Output Interval == 3600
    Map Output Interval == 120

Of note is that each output is tracked every timestep for its maximum if Maximums and Minimums is set to “ON” or “ON MAXIMUMS ONLY” (the default) to ensure that the peak result (water level, depth, hazard etc.), independent of the time it occurred during a simulation, is recorded.

11.2.2 Map Output Formats

TUFLOW offers a wide range of map output formats, all non-proprietary, to cater for a range of GIS and GUI software. There are no constraints over how many output formats a single simulation can produce. Different formats can have different output settings, and by using Output Zones different regions of the model can output in different formats.

Map output is offered in the following forms:

  • Mesh Based: The output is based on a quadrilateral and/or triangular mesh of the 1D and 2D domains. This output is the closest reproduction of the hydraulic calculations, with minimal interpolation from the 1D and 2D computational points. These formats include: XMDF, DAT, TMO and WRB.
  • Grid Based: The output is over a regular north-south grid in a similar manner to a raster DEM. The hydraulic output at each output grid cell is interpolated from the computational mesh using the Mesh Based output above. These formats include: TIF, GPKG, FLT, ASC, NC, TGO, WRR.
  • High Resolution Grid Based: If using SGS, the retained sampled elevations can be used to produce high-resolution (HR) grid outputs. All supported Grid Based outputs can be high-resolution: HRTIF, HRGPKG, HRFLT, HRASC, HRNC, HRTGO, HRWRR.
  • Combination Mesh and Grid: This format includes: WRC.
  • GIS Based: GIS layers written as the simulation proceeds (similar to using the post-processing TUFLOW_to_GIS utility). The format is set by GIS Format.

The selection of formats is controlled by the TCF command Map Output Format. One or more formats may be specified for the whole model or for an Output Zone (refer to Section 11.2.5).

If no output format has been specified, the results are written by default using the XMDF format. The following sections describe the supported map output formats.

11.2.2.1 Mesh Based Formats

Mesh based map outputs use a mesh of quadrilaterals and/or triangles of the 2D cells and any 1D WLLs to manage and store the model results.

The recommended mesh format is the XMDF.

Table 11.1: Map Output Format Options - Mesh Based Output Formats
Format Description
XMDF

The default mesh format.

XMDF (.xmdf) was developed by Aquaveo as a faster and more space efficient replacement to the DAT format. The XMDF format complies with the HDF5 standard. The advantages of the XMDF format are:

  • all specified Map Output Data Types are in the one file;
  • clear labelling of output types and maximum outputs;
  • an internal folder structure is available for use by GUIs; and
  • supports compression (see XMDF Output Compression).

For more information on the XMDF structure for standard and specialised outputs see the TUFLOW Wiki.

The XMDF format is supported by most GUIs, the TUFLOW QGIS Viewer and the TUFLOW Utilities. For these reasons XMDF is the default.
DAT The DAT (.dat) format is a legacy format now superseded by XMDF and no longer recommended for use. For more info on the DAT format see the 2018 TUFLOW Manual.
TMO

The TMO (.tmo) output format is utilised by 12D Solutions for their TUFLOW GUI interface. If using 12D to display/view results this format will need to be specified using Map Output Format. This format contains 2D domain cell centred outputs from the model.

Note, the Output Zones feature (Section 11.2.5) is currently not yet available for the TMO format.
WRB

WaterRIDE by Worley Parsons is commercial software for visualising and post-processing hydraulic modelling results. TUFLOW supports the WRB (.wrb), WRC (.wrc) (Section 11.2.2.1.1) and WRR (.wrr) (Section 11.2.2.4) WaterRIDE formats. WaterRIDE triangulation format for visualising and post-processing hydraulic modelling results. Results are only output in the cell centred triangular mesh arrangement as WaterRIDE does not support quadrilateral elements in a mesh.

One .wrb file is produced for each simulation that contains the model’s ground/bathymetric elevations, water levels, velocities (scalar and vector), and optionally the Z0 (VxD product) and one hazard category. The .wrb format is restricted to these data types, other data types specified using Map Output Data Types are ignored for WaterRIDE output. Also, WaterRIDE can only display a single hazard category (the first one is used if more than one is specified). If Z0 and/or a hazard category are not specified for WRB output, WaterRIDE can optionally post-process these hazard values noting that the maximum of these hazard values, and the actual values, may differ from those directly output from TUFLOW due to post-processing interpolation effects. Other data types such as depth are also post-processed by WaterRIDE (and may differ from those directly output from TUFLOW).

If maximums are tracked (see Maximums and Minimums) these are also added to the .wrb file for the data types mentioned above. Note that if WaterRIDE is used to post-process maximums the values will be different to those provided by TUFLOW. TUFLOW tracks maximums every timestep, whilst WaterRIDE post-processes maximums using the values in the .wrb file, which only occur every Map Output Interval.
CC A cell centred map output format. The output writes to NetCDF file format and uses the 2D domain’s cell size and orientation. This feature is only compatible with single 2D domain models and does not include any 1D output via WLLs. The CC option allows for rotated TUFLOW model grids to be output at the 2D cell size resolution without interpolation to a north-south aligned raster. The file format of the output is described on the TUFLOW Wiki.

Additional time outputs are available using Time Output Cutoff Depths or Time Output Cutoff Hazards to output maps of the duration of inundation and the time of first inundation above specified depth(s) or VxD(s).

The options in the table below only apply to the XMDF format. Only one (or none) can be specified.

Note that the Map Output Format SMS HIGH RES option is a legacy feature no longer supported and is, therefore, not included in the table below. For details on the SMS HIGH RES feature search the 2018 TUFLOW Manual.

Table 11.2: Map Output Format Options - Options for XMDF Format
Format Description
SMS This is the default SMS option and outputs ground elevations and results at the cell corners only.
SMS TRIANGLES

Outputs 2D cells as four triangles rather than as a quadrilateral if DAT and/or XMDF are specified. The triangles are constructed so that the 2D cell centre is a common vertex to all four triangles. This means that the mesh is entirely constructed of triangles (four triangles per 2D cell and any 1D WLL triangles).

Note: For the formats that rely on this triangle only mesh (e.g. ASC, FLT, NC, T3, TGO, WRB, WRR), the SMS TRIANGLES option does not need to be invoked. Only specify SMS TRIANGLES if you require your XMDF or DAT output to be based on a triangle only mesh rather than the default mesh of quadrilaterals for 2D cells and triangles for 1D WLLs. See Section 11.2.2.1.1.
11.2.2.1.1 Mesh Configurations

Quadrilateral and Triangle Mesh Option

The default mesh used by the XMDF format uses quadrilaterals (squares) for the 2D cells and triangles for representing any 1D WLLs as discussed in Section 11.2.4. Each 2D domain and the 1D WLL triangles are treated as separate meshes, although they will appear as one when viewing the .2dm file.

The advantage of this format is that the amount of data output is optimised keeping file sizes small with little loss of accuracy in translating results from the computational mesh.

Only the XMDF format supports the combination of quadrilaterals and triangles.

Triangular Mesh Option

Except for the XMDF format, all the other mesh-based formats use a triangular mesh to represent the computational 1D and 2D domains. The triangular mesh is also used by the grid-based formats (TIF, FLT, ASC, NC, TGO, WRR) to interpolate from the triangular mesh the values at the grids cells’ centres (this is the cell centres of the output grid, not the 2D domain grid cells).

To utilise this format using the XMDF format, specify SMS TRIANGLES anywhere in the arguments for the Map Output Format command.

The triangular mesh option incorporates output at the cell centres, so the exact water level calculated by TUFLOW at the 2D domain cell centres is used when translating results, thereby providing a slightly more representative surface of the hydraulic calculations. Each 2D cell is represented as four triangles with a common vertex at the cell centre, giving a higher resolution spatial representation than just using the 2D cell corners.

However, the number of nodes in the mesh for the 2D domains increases by 20% and the number of elements by a factor of 4, therefore, the output file sizes will be larger and the viewing and processing times possibly slower. Any triangular elements from 1D WLLs are not affected and remain the same in either mesh arrangement.

2D Cell Corner Interpolation/Extrapolation

Irrespective of the mesh option being used, the Map Output Data Type values at the 2D cell corners need to be interpolated (or extrapolated if at the wet-dry interface). Two methods were developed that can be set using the command Map Output Corner Interpolation == METHOD B or METHOD C. Method C (the default) is recommended as it is simpler and extensive testing has indicated it resolves issues associated with prior methods. METHOD B does largely resolve earlier issues, and has the added advantage that the effect of thin breaklines is better handled when they are dry or flow is upstream controlled.

For model results where the water has risen upwards (e.g. river flooding), the Method B or C approaches should cause no significant changes in results (i.e. fractions of a mm). Where the flow is downwards over steep slopes, some changes in results at cell corners will occur, but usually only slightly. However, maximum hazard values on very steep slopes may experience a more significant change.

11.2.2.2 Grid Based Formats

Grid-based formats output the results over a regular grid in similar formats to DEMs and other raster-based data. The output grid resolution and origin is not necessarily the same as the 2D hydraulic modelling grid, and will also use any 1D WLL triangulations to include map based output of 1D domains. Grid based formats are highly suited for use in GIS, and especially for handling output from very large models (e.g. GPU models).

The available grid formats are listed in Table 11.3. The output grid format is set with the Map Output Format command, for example:

Map Output Format == TIF

Alternatively, instead of including a specific format (e.g. TIF or FLT) in the Map Output Format list, this can also be specified using the Map Output Format == Grid command. In this case the grids are written by default in TIF format. This can be changed to other grid formats by specifying the specific format in the .tcf command Grid Format. If Write Check Files is used, the DEM_Z and DEM_M grid outputs will be of the format specified by the Grid Format command (or TIF by default). Note that not all formats support compression, so those that don’t can produce considerably larger file sizes.

The output grids are placed in a sub-folder called “grids”. For TUFLOW Classic the default output grid resolution is half the smallest 2D cell size (considering that multiple 2D domains may exist). For TUFLOW HPC the default is half the 2D cell size. For HPC models using Quadtree, the output resolution is the smallest 2D cell size. The resolution of the output grid can be controlled by the Grid Output Cell Size command. Only one output grid resolution is possible. If multiple commands exist, the last occurrence of the command will prevail. The origin of the output grid is rounded to the nearest cell size. This can be set to the exact model origin using the Grid Output Origin == Model Origin.

The TUFLOW grid output formats (listed in Table 11.3) support all available Map Output Data Types (listed in Table 11.4 and Table 11.5). The minor exceptions are noted in Table 11.3.

Grids outputs store data on a grid by grid basis and in a north-south direction, therefore interpolation must occur to convert the data from the TUFLOW calculation points. The raster data is interpolated from the values at the cell corners and centres.

For example, the image below shows:

  • A _grd_check file (black lines)
  • A _zpt_check showing the TUFLOW calculation (ZC, ZU, ZV and ZH) points;
  • An output water level TIF grid (underlying grey grid); and
  • Orange dotted lines showing the triangulation between the cell corners and centres.

The raster data (grid centre indicated by red cross) obtains its value from the triangulation between A, B, and E.

To more closely align the grid output with the TUFLOW calculation points:

The <format> Map Output Interval can be set to zero as follows to trigger only outputting the maximums for those formats. For example:

TIF Map Output Interval == 0 ! Only output result maximums for TIF grids
FLT Map Output Interval == 0 ! Only output result maximums for FLT grids

All time outputs are supported by the grid map output formats. The grid file extensions are:

  • _TDur_<cutoff> for duration of inundation;
  • _TExc_<cutoff> for time <cutoff> is exceeded;
  • _TMax_h for time of peak water level; and
  • _TMax_V for time of peak velocity.
Table 11.3: Map Output Format Options - Grid Based Output Formats
Format Description
TIF

The GeoTIFF raster compressible format is supported from the 2023-03 Release and onwards and is the default format. A projection can be set for the output GeoTIFF rasters by using the TIF Projection command. Individual files are created for each output type, time, maximums, etc.

TUFLOW supports several compression methods for GeoTIFF, see TIF Compression, the default is the “deflate” method. A compression predictor “horizontal differencing” is used to improve the compression ratio, see TIF Compression Predictor. TUFLOW will default to using all available CPU cores when writing GeoTIFF files which can speed up processing when using compression. This can be changed by specifying the number of threads using the command line argument “-nt[thread count]”.
FLT

ESRI binary (float) version of the ASC format. The data is uncompressed. The file header containing the dimensions of the grid is output to a .hdr (text) file and contains the same header as for an .asc file. The remainder of the output, the 3D surface values, is written to a .flt file as a binary dump rather than as a text file.

Advantages of this format are that it is very simple and is much faster to write and open than the ASC format. Main limitations are that file sizes are large, and a separate file needs to be written for every Map Output Data Type specified and for every output time (that can be a lot of files!).
GPKG

The GeoPackage grid format is supported from the 2023-03 Release onwards. The format uses a tiled structure to make rendering and loading faster by enabling the ability to only process the required tiles. The tiled structure also makes pyramids (sometimes referred to as ‘overviews’) inherently available for the format. Individual files are created for each output type, time, maximums, etc.

The GPKG raster outputs will be grouped if the Spatial Database command is set to “grouped” (the default).

The GPKG raster format supports LZW compression of the data, see GPKG Compression. A compression predictor “horizontal differencing” is used to improve the compression ratio, see GPKG Compression Predictor.
ASC

ESRI ASCII (.asc) uncompressed grid format, a long-established industry standard format often used for transferring 3D surfaces between GIS software.

Advantages of this format are that it is very simple and can be viewed in a text editor. Limitations are that for large grids the file maybe slow to write, slow to open and work with, file sizes can be very large, and a separate file needs to be written for every Map Output Data Type specified and for every output time (that can be a lot of files!).
GRID Outputs grid files in the default grid output (.tif) unless Grid Format is specified in the .tcf.
NC

The NETCDF (Network Common Data Format) is a commonly used compressible format for storing modelling and scientific data. A single file is created that contains all output types, output times and static outputs (e.g. maximums).

A number of NetCDF specific commands are supported as listed below. The TUFLOW Wiki page TUFLOW NetCDF Raster Format provides additional information.

If using FEWS to view NC output, please ensure that NetCDF Output Format == FEWS has been set.
TGO The TGO format is utilised by 12D Solutions for their TUFLOW interface. The output is a north-south aligned raster and includes outputs from multiple domains and water level lines. Unlike the mesh .tmo format, all Map Output Data Types are supported, and the resolution of the output grids can be set using the command Grid Output Cell Size. Individual files are created for each output type. As of the 2020-10 release the TGO velocity angle is in radians, previously this was degrees.
WRR A WaterRIDE uncompressed format that contains the time varying grid output using a north-south aligned raster and includes outputs from multiple domains and 1D WLLs. All Map Output Data Types are supported, and the resolution of the output grids can be set using the command Grid Output Cell Size. A single file is created for all output types, times, etc. Due to no compression, file sizes can be very large.

11.2.2.3 High-Resolution Grids

When using the Method C sub-grid sampling (SGS) approach the sampled elevations are retained including topography modifiers such as breaklines. At the end of the geometry processing this allows a high-resolution grid to be written and used for high resolution map outputs and also DEM check grids (DEM_Z_HR). All grid formats in Table 11.3 are supported. Currently, high resolution raster outputs are only available for depth (d) and water level (h) output data types.

To produce high-resolution output, add “HR” to the grid format. For example “HRTIF” or “HRFLT” to produce a high resolution TIF or FLT respectively. Note, the Map Output Interval will also have to be set, for example, “HRTIF Map Output Interval == 0” will output just the maximums using the HRTIF format.

The high resolution grid output uses the following approach:

  • The water level at each HR raster output cell is interpolated from the computed 2D water levels.
  • The depth is the difference between the interpolated water level and the terrain elevation interpolated from the surrounding SGS sampled elevations.

The regular grid output interpolates depth from cell centres/corners with a default output resolution of half the cell size. The difference between the regular and the HR grid output is illustrated in Figure 11.1.

**Regular Grid Depth Output (Left) and HR Grid Depth Output (Right)**

Figure 11.1: Regular Grid Depth Output (Left) and HR Grid Depth Output (Right)

When modelling breaklines in TUFLOW, “thin” breaklines modify the cell face elevations but do not modify the cell storages. When outputting the high-resolution outputs, the user can set whether the cell face elevations are included using the .tcf command (the default is ON, to use face elevations):

HR Grid Output Use Face Elevations == {ON} | OFF

The water level and depth interpolation approaches for the high resolution outputs can be altered depending on the type of model (e.g. hydrology inflows or direct rainfall) and specifically how TUFLOW outputs the high resolution results around thin breaklines. These output options are discussed on the HR Output TUFLOW Wiki page.

HR grid outputs are also compatible with the output zone functionality, see Section 11.2.5.

11.2.2.4 Combined Format

Combination of a mesh and grid-based approach. The only available combined format is the WRC format. The WRC format produces a master (.wrc) file and one or more WRR and WRB files. The approach adopted by TUFLOW is as follows:

  • Each 2D domain is output as a rotated WRR format grid with the cell size equal to the TUFLOW cell size. The WRR format is significantly faster than than the WRB format due to its grid-based formatting. The output values are the cell centre values for each 2D cell.
  • Any 1D WLL triangulations are output as a separate WRB file.
  • The WRC “master” file is output to the specified results folder, while the WRR and WRB files are written to a “waterRIDE” sub folder.

11.2.2.5 GIS Based Format

GIS based map outputs use the format specified by GIS Format to write out map output as a series of GIS layers. Gridded output format as GIS layers can be written directly from TUFLOW during the simulation by including “GIS” in the Map Output Format command. This offers a similar functionality to that using the TUFLOW_to_GIS utility via, for example, the –shp option. For scalar outputs, these are output as a point GIS layer with a separate file for each output time. For vector outputs, this can be either as a point or region GIS layer. Specific commands to this output format are:

The commands above can be applied to all vector outputs or can be specific to the data type by prefixing with a “v”, “q” or “W”, for velocity, unit flow or wind respectively. The example below sets the scale factor to 1 for all outputs except unit flow, which has a smaller factor of 0.1.

GIS Grid Vector SF == 1.0 ! Set scale factor for all GIS vector output to 1.0
q GIS Grid Vector SF == 0.1 ! Set scale factor for unit flow GIS output to 0.1

11.2.3 Map Output Data Types

TUFLOW can output a wide range of output types in map format. Table 11.4 describes all the non-hazard map output types, while Table 11.5 contains all flood hazard category (Z) map output types. The map output types produced by a simulation are controlled using the .tcf command Map Output Data Types.

The map output types’ flags are listed in the first column of the tables and are used to denote the type(s) to be output. They can occur in any combination or order and are not case-sensitive. For example, to output water level, velocity and unit flow, enter the following line in the .tcf file:

Map Output Data Types == h V q

Although optional, it is strongly recommended that spaces are used between each data type for clarity.

The output types are available in a wide range of Map Output Formats and can be varied for different formats (Section 11.2.2) and between output zones (Section 11.2.5). Not all Map Output Data Types are available for all Map Output Formats due to limitations or constraints of the type/format. The supported formats for each type are documented in Table 11.4 and Table 11.5.

It is possible to get different output types for different output formats, as discussed in Section 11.2.2. For example:

TIF Map Output Data Types == h V Z0
XMDF Map Output Data Types == h V q d ci z0

Table 11.4: Map Output Types (Excluding Hazard (Z) Types)
Flag Map Output
Data Type
Description Supported Formats
AP Atmospheric Pressure Atmospheric pressure in hPa. Atmospheric Pressure is only available if using the Read GIS Cyclone or Read GIS Hurricane commands. Maximum and minimum output is not available. All formats excluding TMO, WRB and HR.
BSS Bed Shear Stress

Bed Shear Stress in N/m2 (or lbf/ft2 in English Units) is given by the equation below where ρ is density, g gravity, V velocity, n Manning’s n, and y depth:

Metric Units: \[\tau_{bed} = \frac{\rho g V^2 n^2}{y^{\frac{1}{3}}} \hspace{3em} (N/m^2)\]
English Units: \[\tau_{bed} = \frac{\rho g V^2 n^2}{2.208 y^{\frac{1}{3}}} \hspace{3em} (lbf/ft^2)\]

The Bed Shear Stress map output can be misleading at very shallow depths as the BSS formula divides by the depth. The BSS and SP outputs are linearly reduced to zero once the depth is below a threshold (by default, 0.1m). This threshold can be changed using the .tcf command BSS Cutoff Depth.

Note, prior to the 2017 release BSS output in English Units were in Poundals per square foot (pdl/ft2).
All formats excluding TMO, WRB and HR.
CI Cumulative Infiltration

The cumulative infiltration over the entire simulation in mm or inches when a soil infiltration method has been used (see Section 7.3.7). Also see the IR (infiltration rate) map output type below. Maximum and minimum output is not applicable/available as it is a cumulative output.

Note, the CI output is restricted to models for which the value is cumulative, specificially this means models that have:

  • A single vertical soil layer.
  • No horizontal movement of ground water.
  • No drying of the soil layer through negative rainfall.
When any of the three options (multiple vertical layers, horizontal movement or negative rainfall) are used, the water content in the layer is no longer a “cumulative value”.
All formats excluding TMO, WRB and HR.
Cr Courant Number

Classic Only.

Courant number (2D domains only). Maximum and minimum output is not available.
All formats excluding TMO, WRB and HR.
CWF Cell Width Factor No longer supported - previously available via the legacy SMS HIGH RES option. N/A
d Depth

Water depths in m or ft. For the cell cornered results formats (see Section 11.2.2.1.1) the depths are calculated as the interpolated water level at the nodes (see _h below) less the ZH value. The interpolated water level may occasionally lie below the ZH value, in which case a negative depth may result, which is set to zero by default (see Zero Negative Depths). Both maximum and minimum output are available.

For maximum depth output this is calculated at the end of the simulation based on the maximum water level and the ground elevation. For models that utilise varying ground elevations (using the Read GIS Variable Z Shape or variable geometry (VG) boundaries), care should be taken when interpreting maximum depth outputs. Hazard outputs (based on velocity and depth) are tracked at each timestep, and the maximum for these is the maximum at any timestep during the model.
All formats excluding WRB.
dGW Depth to Groundwater

Depth to groundwater (from the ground surface) over time in metres or feet when a groundwater depth or level has been defined (see Section 7.3.7).

If using multiple sub-surface layers in TUFLOW HPC (see Section 7.4.5 it is the distance from the ground surface to the groundwater level of the layer in question.

Maximum and minimum output is not available.
All formats excluding TMO, WRB and HR.
dt Minimum Timestep

HPC Only.

The dt map output is a grid map output that displays the calculated minimum timestep at each grid cell. It does not necessarily align with the timestep adopted in the model (the timestep in the model could be less if the previous timesteps were smaller). This helps identify which cells in the model are controlling the model timestep. For more information see the TUFLOW Wiki.
All formats excluding TMO, WRB and HR.
E Energy

Scalar data file containing the energy levels at the element nodes (cell corners). The energy levels are based on the interpolation of water levels and dynamic head (V2/2g) at the output location. As both the 1D and both 2D solutions use a staggered computational scheme, the interpolation of the water levels and/or velocities to the output location may occasionally cause an “increase” in energy to be apparent in the output, therefore, energy output should be treated with caution.

For 1D areas, this output should be treated with additional caution as it is derived from interpolation of water levels and approximations of the channel velocities across the WLLs, which can be problematic in 1D channels with high velocities. The energy output for 1D nodes is available as part of the plotting output (Section 15.3.1).

Maximum energy levels is for when the maximum water level occurs (Note: This may cause undulations in the maximum energy due to variations in the time of the maximum water level).
All formats excluding TMO, WRB and HR.
F Froude Number Froude number output. No maximum and minimum output is available at this stage. All formats excluding TMO, WRB and HR.
FLC Form Loss Coefficient No longer supported - previously available via the legacy SMS HIGH RES option. N/A
GWd Groundwater Depth

Depth of water within each sub-surface layer(s) when layered interflow has been defined (see Section 7.4.5). Determined by dividing the cumulative infiltration by porosity.

Reported in metres or feet.
All formats excluding TMO, WRB and HR.
GWh Groundwater Level

Elevation of the groundwater surface (water table) for each sub-surface layer(s) when layered interflow has been defined (see Section 7.4.5).

Reported in metres or feet.
All formats excluding TMO, WRB and HR.
GWm Groundwater Moisture Groundwater moisture output. Dimensionless number in the range of zero to one representing a “fraction full” for each sub-surface layer(s) when layered interflow has been defined (see Section 7.4.5). All formats excluding TMO, WRB and HR.
GWq Groundwater Unit Flow

Groundwater unit flow (m2/s, flow per unit width) at the nodes (cell corners) of the sub-surface layer(s) when layered interflow has been defined (see Section 7.4.5).

Reported in m2/hr or ft2/hr.
All formats excluding TMO, WRB and HR.
GWv Groundwater Velocity

Groundwater flow velocity when layered interflow has been defined (see Section 7.4.5). The resulting groundwater velocity vector is calculated from the surrounding u and v-points.

Reported in m/hr or ft/hr.
All formats excluding TMO, WRB and HR.
h Water Level Water level output. For the cell cornered results formats (see Section 11.2.2.1.1) the water levels are interpolated from the water levels calculated at the cell centres. Both maximum and minimum outputs are available. All formats.
IR Infiltration Rate

The infiltration rate in mm/hr or inches/hr over time when a soil infiltration method has been used (see Section 7.3.7). See also the CI (cumulative infiltration) map output type above. Maximum and minimum output is not available.

If using multiple sub-surface layers (see Section 7.4.5) the infiltration rate is reported for the first (top) layer only.
All formats excluding TMO, WRB and HR.
MB1 Mass Balance

TUFLOW Classic only (for HPC, see the ‘dt’ output).

Measure of the convergence level of the solution. The measure is a cumulative value since the last output time, therefore is an effective way of identifying problem areas in a model that repeatedly have poor convergence and most likely mass error. Very useful for identifying problem areas within a model.

This output does not include 1D output from WLLs.
All formats excluding TMO, WRB and HR.
MB2 Mass Balance

TUFLOW Classic only (for HPC, see the ‘dt’ output).

Same as MB1 above but is accumulated over the entire simulation.

This output does not include 1D output from WLLs.
All formats excluding TMO, WRB and HR.
n Manning’s n Manning’s n values. The n values only vary over time for materials using the Manning’s n varying with depth feature. The n values at the cell corners in the _n.xmdf file are interpolated from the surrounding four cell mid-sides. Maximum and minimum output is not available. All formats excluding TMO, WRB and HR.
Nc Shallow Wave Celerity Number

HPC Only.

The Shallow Wave Celerity Number. One of the three controls TUFLOW HPC uses to determine the maximum timestep to maintain stability. See the HPC Adaptive Timestepping Wiki Page for discussion on adaptive timestepping and the Nc number.
All formats excluding TMO, WRB and HR.
Nd Diffusion Number

HPC Only.

The Diffusion Number. One of the three controls TUFLOW HPC uses to determine the maximum timestep to maintain stability. See the HPC Adaptive Timestepping Wiki Page for discussion on adaptive timestepping and the Nd number.
All formats excluding TMO, WRB and HR.
Nu Courant Number

HPC Only.

The Courant Number. One of the three controls TUFLOW HPC uses to determine the maximum timestep to maintain stability. See the HPC Adaptive Timestepping Wiki Page for discussion on adaptive timestepping and the Nu number.
All formats excluding TMO, WRB and HR.
q Vector Unit Flow

Unit flow (m2/s, flow per unit width) at the nodes (cell corners). The resulting flow vector is calculated from the surrounding u and v-point velocities and the depth.

Unit flow may also be used as a measure of flood hazard (i.e. velocity by depth or VxD).

Note: The maximum unit flow is not tracked for the q output, the Z0 scalar hazard value option can be used, as this output is tracked at each timestep.
All formats excluding TMO, WRB and HR.
R Flow Regime

TUFLOW Classic only.

Flow regime. The output value is 0 (zero) for normal (sub-critical flow with momentum); greater than 1 for upstream controlled friction flow (e.g. supercritical flow); ‑1.5 for broad-crested weir flow; and ‑1 for flow through a flow constriction when the deck is submerged. No maximum and minimum output is available at this stage. The flow regime at the cell corners and centre is a weighted average (using unit flow) of the flow regimes at the four surrounding cell mid-sides, therefore, where there’s different flow regimes at any of the four cell faces the R value can be misleading due to the averaging.
All formats excluding TMO, WRB and HR.
RC Route Category

The route category output over time for evacuation routes. The definition and number of categories is based on the values specified within the Cut_Off_Values attribute of the 2d_zshr GIS layer (see Section 11.4.2). The RC values are output as an integer representing the closure category specified by the user.

The maximum RC category value is tracked every timestep and output (if tracking maximums is switched on, which is the default).
All formats excluding TMO, WRB and HR.
RFC Cumulative Rainfall

The cumulative rainfall in mm or inches over time when direct rainfall has been applied to the model (refer to Section 8.5.3). See also the RFR (rainfall rate) map output type below. Both the RFC and RFR outputs (see next item) are inclusive of any boundary adjustments (e.g. in the boundary database) and rainfall losses applied in the materials file. Soil infiltration is applied once the rainfall has been applied to the cells, so this is not accounted for in the rainfall outputs, see also CI (cumulative infiltration) and IR (infiltration rate) output types.

Maximum and minimum output is not applicable/available, as it is a cumulative output.
All formats excluding TMO, WRB and HR.
RFML Material Based Rainfall Loss The output contains the total rainfall losses applied due the initial and continuing rainfall losses specified in the Read Materials File (.tmf or .csv) file. The RFML option can be used to track the rainfall based material losses that have been applied spatially. The RFC and RFR map output data types can be used to output the cumulative rainfall and rainfall rate. All formats excluding TMO, WRB and HR.
RFR Rainfall Rate The rainfall rate in mm/hr or inches/hr over time when direct rainfall has been applied to the model (refer to Section 8.5.3). See also the RFC (cumulative rainfall) map output type above. Maximum and minimum output is not available. All formats excluding TMO, WRB and HR.
SP Stream Power

Stream Power as given by the equation below where τbed is bed shear stress (see BSS above) and V is velocity.

\[Stream Power = |V|\tau_{bed} \hspace{3em} (W/m^2)\]

The Stream Power map output can be misleading at very shallow depths as the BSS formula divides by the depth. The BSS and SP outputs are linearly reduced to zero once the depth is below a threshold (by default, 0.1m). This threshold can be changed using the .tcf command BSS Cutoff Depth.

Prior to the 2017 release SP output in English Units were in Poundals per square foot (pdl/ft2). From the 2017 release onwards the units are in Pounds Force per square foot (lbf/ft2), therefore, the SP values are 32.174 times smaller than for releases prior to 2017.
All formats excluding TMO, WRB and HR.
SS Sink / Source Flow

The net source/sink inflows. Note the flow rate for a cell is shown at the ZH point (top right of the cell).

Maximum and minimum output is not available.
All formats excluding TMO, WRB and HR.
t Viscosity Coeff Eddy viscosity coefficient. This is useful for checking the Smagorinsky coefficient values. No maximum and minimum output is available at this stage. All formats excluding TMO, WRB and HR.
tau Shear stress

This output contains the shear stress values applied via the external stress file (.tesf).

The output values are in Newtons per square metre (N/m2) for SI units and pound-force per square foot (lbf/ft2) for US customary (English) units.
All formats excluding TMO, WRB and HR.
V Vector Velocity

Flow velocity. The resulting velocity vector is calculated from the surrounding u and v-points.

Note: The maximum and minimum velocities are tracked over time, however, be careful interpreting maximum velocities displayed as vectors as the flow patterns may appear to be opposing each other - this is due to the maximum velocities not all occuring at the same instant in time. By default the maximum velocities are tracked over 0.1m depth, below this depth the velocity at maximum water level is used. See the Maximum Velocity Cutoff Depth command for more information.
All formats excluding HR.
WI10 Wind Vector Wind Vector output in m/s at a height of 10m. Available when a cyclone/hurricane boundary is used, see Section 8.7. All formats excluding TMO, WRB and HR.
ZH Bathymetry

Elevations at the cell corners (ZH points). This information is already contained in the .2dm file, however, this option is useful if the model’s bathymetry varies over time because of variable geometry (2d_vzsh or VG boundaries) or for morphological modelling. This output is very useful if you are comparing two or more runs that have different topography (e.g. before and after scenarios), and you wish to easily view or compare the topography for each scenario.

If the topography in the model does not change over time (i.e. no variable Z shapes or morphological changes), for the default .xmdf output format the ZH Zpt values are output once, rather than every timestep, thereby not consuming disk space unnecessarily. The ZH map output will appear under a XMDF folder “Fixed”. This feature is only available if using the XMDF format, for other output formats, the bathymetry will be output at each output interval.

Currently, no maximum and minimum output is available.
All formats excluding TMO, WRB and HR.

11.2.3.1 Hazard Data Types

Table 11.5 presents the hazard category outputs. Of note is that each hazard is tracked every timestep for its maximum if Maximums and Minimums is set to “ON” or “ON MAXIMUMS ONLY” (the default) to ensure that the peak hazard category is recorded during the simulation. Up to ten (10) different hazard categories per simulation can be specified for map output. Grid map output hazard categories are output as integer grids (i.e. values are rounded to the nearest integer when a grid output cell centre is located at a change in category) when using Map Output Data Types except for output Z0, ZUK0 and ZUK2, which are output as real numbers.

Users have the option to customise hazard outputs based on a .csv file input using the “ZUD1” hazard type. This feature allows users to define a number of depth, velocity and depth-velocity product thresholds to create custom hazard categories. Like other hazards these are assessed at each computational timestep and the maximum is tracked every timestep. To specify a user defined hazard, add output type “ZUD1” in the “Map Output Data Types” command, e.g.:

Map Output Data Types == h d V ZUD1

The thresholds are defined through reading a .csv file using the following command:

Read Hazard File == <file_name>.csv

The csv file should contain three (3) columns defining the thresholds for depth, velocity, and depth-velocity product respectively. The example below produces the same output as for “ZNZ1” described in Table 11.5.

**Example Hazard File for User Defined Hazard**

Figure 11.2: Example Hazard File for User Defined Hazard

  • Category 1: D < 0.1. Note: the velocity, and depth-velocity product thresholds are effectively not applied by using a large threshold value such as 99999.
  • Category 2: D < 0.5 and V < 1.0. Note: the depth-velocity product threshold is effectively not applied by using the large threshold value of 99999.
  • Category 3: D < 1.0 and V < 2.0 and D*V < 1.0
  • Category 4: Otherwise
  • Dry points are assigned Category 0.

A maximum of 10 categories can be specified by the .csv file, but just one User Defined Hazard is currently allowed per simulation. This functionality is compatible with both TUFLOW Classic and HPC solvers.

Please email if you would like to code a unique hazard output that is currently not included as an output type option or possible using the user defined hazard type (ZUD1).

Table 11.5: Map Output Hazard (Z) Types
Flag Supported Formats Description
Z0 All formats Velocity x Depth product
Z1 All formats

Flood hazard category based on the Australian NSW Floodplain Management Manual (NSW Goverment, 2005b). The output is an integer number from 1 to 3 as follows and as illustrated in the figure below.

  • 1: Low Hazard
  • 2: Intermediate Hazard (dependent on site conditions)
  • 3: High Hazard
Z2 All formats

V > 2.0: Category 5 (Extreme Hazard)
D*V > 1.0: Category 4 (High Hazard Floodway)
D > 1.0: Category 3 (High Hazard Depth)
V + 3.3*D > 2.64: Category 2 (High Hazard)
Otherwise Category 1 (Low Hazard)
Dry points are assigned Category 0

Z3 All formats

V > 2.4: Category 7 (Extreme Hazard)
D*V > 2.0: Category 6 (Extreme Hazard Floodway)
V > 2.0 OR V*D > 1.0: Category 5 (High Hazard Floodway)
D > 2.0: Category 4 (Extreme Hazard Depth)
D > 1.0: Category 3 (High Hazard Depth)
V + 3.3*D > 2.64: Category 2 (High Hazard)
Otherwise Category 1 (Low Hazard)
Dry points are assigned Category 0

Z4 All formats

Flood hazard mapping approach based on the Australian Guidelines (CSIRO, 2000) using the following logic:

  • D ≤ 0.3 and V ≤ 0.38: Category 1 (Low Hazard)
  • D ≤ 0.6 and V ≤ 0.8 and D + 0.64*V ≤ 0.82: Category 2 (Medium Hazard)
  • D ≤ 1.2 and V ≤ 1.5 and D + 0.69*V ≤ 1.38: Category 3 (High Hazard)
  • Otherwise Category 4 (Extreme Hazard)
  • Dry points are assigned Category 0
Z7 All formats

Based on Figure L1 of the NSW Floodplain Risk Development Manual (NSW Goverment, 2005b), using the following logic:

  • Category 1 (low hazard): V ≤ 2 and V + 6.667 x D ≤ 3
  • Category 2 (vehicles unstable) : V ≤ 2 and V + 3.333 x D ≤ 2.667
  • Category 3 (wading unsafe) : V ≤ 2 and D ≤ 2 and V x D ≤ 1
  • Category 4 (extreme) : Otherwise
Z9 All formats

Flood hazard mapping approach based on the draft storm tide hydraulic hazard categorisation developed for Moreton Bay Regional Council (GHD, 2011). Z9 output values are 0 (zero) for no hazard and 1 to 5 for H1 to H5 respectively.

ZAEM1 All formats

Flood hazard category as outlined by Australian Emergency Management Institute in 2014 (Australian Emergency Management Institute, 2014). ZAEM1 output values are 0 (zero) for no hazard and 1 to 6 for H1 to H6 respectively.

ZGC1 All formats

Flood Intensity output used by Gold Coast City Council (Australia). The logic applied is as follows:

  • D < 0.3, V < 0.4 and D*V < 0.25: Category 1 (Low Intensity)
  • D < 0.6, V < 0.8 and D*V < 0.4: Category 2 (Medium Intensity)
  • D < 1.2, V < 1.5 and D*V < 0.6: Category 3 (High Intensity)
  • Otherwise Category 4 (Extreme Intensity)
  • Dry points are assigned Category 0 (No Intensity).
ZMBRC All formats Flood hazard output used by Moreton Bay Regional Council (Australia). Where:

V > 2.5 or D > 2.5 or V*D > 2.5: Category 5 (H5)
V > 2.0 or D > 2.0 or V*D > 1.0: Category 4 (H4)
V > 3.2 - 4D: Category 3 (H3)
V > 0.5 or D > 0.3: Category 2 (H2)
Otherwise Category 1 (H1)
Dry points are assigned Category 0

Where:
H1: Hydraulically suitable for parked or moving cars.
H2: Hydraulically suitable for parked or moving heavy vehicles and wading by able-bodied adults.
H3: Hydraulically suitable for light construction (e.g. Timber frame and brick veneer).
H4: Hydraulically suitable for heavy construction (e.g. steel frame and reinforced concrete).
H5: Generally unsuitable
ZMW1 All formats No longer recommended by Melbourne Water. For details on this hazard category, see the 2018 TUFLOW Manual.
ZMW2 All formats No longer recommended by Melbourne Water. For details on this hazard category, see the 2018 TUFLOW Manual.
ZMW3 All formats

Based on Melbourne Water (Australia) FM&M Technical Specifications (Melbourne Water, 2016) to quantify Safety Risk in Roads. The logic used is as follows:

D*V < 0.4 and D < 0.4: Low Risk
D*V < 0.8 and D < 0.8: Medium Risk
Otherwise: High Risk
Dry points are assigned Category 0
ZNZ1 All formats

The flood hazard category based on the Hamilton City Council (New Zealand) Flood Hazard Report [HamiltonFloodMapping].

The categories are defined as:

  • D < 0.1: Category 1 (Insignificant)
  • D < 0.5 and V < 1.0: Category 2 (Low)
  • D < 1.0 and V < 2.0 and D*V < 1.0: Category 3 (Medium)
  • Otherwise Category 4 (High)
  • Dry points are assigned Category 0.
ZNZ2 All formats

Hazard category used by Auckland Council, New Zealand.

  • Category 1 (Low Safety Risk): D*V < 0.4, V < 3.0 m/s and D < 0.5 m
  • Category 2 (Moderate Safety Risk): D*V < 0.8, V < 3.0 m/s and D < 1.2 m
  • Category 3 (High Safety Risk): D*V ≥ 0.8, V ≥ 3.0 m/s and D ≥ 0.1m, and D ≥ 1.2 m
  • Dry points are assigned Category 0.
ZPA All formats People Hazard category “Hazard to Adults” based on the Australian Rainfall and Runoff (ARR) Project 10 Stage One Report (Cox et al., 2010).

The values within the ZPA output are:

  • 0 = Safe (no hazard)
  • 1 = Low Hazard
  • 2 = Moderate Hazard
  • 3 = Significant Hazard
  • 4 = Extreme Hazard

It is possible to specify cut-off depth/s representing when the Safe category applies by using the .tcf command ZP Hazard Cutoff Depth. Up to three values can be defined, which are the cut-off depths for ZPA, ZPC and ZPI respectively.

The relevant tables and figures are reproduced below:


ZPC All formats People Hazard category “Hazard to Children” based on the ARR Project 10 Stage One Report (Cox et al., 2010). Refer to the description of output values for ZPA.
ZPI All formats People Hazard category “Hazard to Infants and frail/elderly People” based on the ARR Project 10 Stage One Report (Cox et al., 2010). Refer to the description of output values for ZPA.
ZPMH All formats

The flood hazard category based on the Port Macquarie-Hastings Council Flood Policy (Port Macquarie-Hastings Council, 2018).

The logic used is as follows:

  • D < 0.4 and V < 0.5 and D*V < 0.2: Category 1 (Low Risk)
  • D < 0.8 and V < 2.0 and D*V < 0.5: Category 2 (Medium Risk)
  • D < 1.8 and V < 2.0 and D*V < 1.5: Category 3 (High Risk)
  • V < 4.0 and D*V < 2.5: Category 4 (Very High Risk)
  • Otherwise Category 5 (Extreme Risk)
  • Dry points are assigned Category 0.
ZQRA All formats Hazard categories for the Queensland Reconstruction Authority (Queensland Reconstruction Authority, 2012).
ZTMR All formats

Hazard category for the Queensland Department of Transport and Main Roads. Areas of a model with the Material ID set to a value of 100 (roads) will be tested for their submergence /closure status:

  • Material ID ≠ 100: No hazard reporting (0)
  • D <= 0: Road is immune (1)
  • D+V2/2g < 0.3: Road is submerged (2)
  • Other: Road Closed (3)
ZUD1 All formats Customisable hazard output based on a .csv file input read in using the Read Hazard File command. This feature allows users to define a number of depth, velocity and depth-velocity product thresholds to create custom hazard categories. The csv file should contain three (3) columns defining the thresholds for depth, velocity, and depth-velocity product respectively.

The example below produces the same output as for “ZNZ1”:

This hazard type is further detailed in Section 11.2.3.1.
ZUK0 All formats The (real) value of the UK Hazard formula based on Environment Agency (2006) (see UK Hazard Formula, UK Hazard Land Use and UK Hazard Debris Factor).
ZUK1 All formats

The UK Hazard category (integer) based on Environment Agency (2006) (see UK Hazard Formula, UK Hazard Land Use and UK Hazard Debris Factor).

ZUK1 outputs a .xmdf file containing an integer, as per:

  • 0 = No Hazard (H \(\le\) 0)
  • 1 = Low Hazard (H \(\le\) 0.75)
  • 2 = Moderate Hazard (H \(\le\) 1.25)
  • 3 = Significant Hazard (H \(\le\) 2.5)
  • 4 = Extreme Hazard (H \(\ge\) 2.5)
ZUK2 All formats The (real) value of the UK Hazard formula based on UK FD2320 Technical Report (Environment Agency, 2008).

\[Flood Hazard Rating = ((v+0.5)\:*\:D) + DF\]

Where:
v = velocity (m/s)
D = depth (m)
DF = debris factor
ZUK3 All formats

The UK Hazard category (integer) based on UK FD2320 Technical Report (Environment Agency, 2008).

Where:

  • 1 = Low
  • 2 = Moderate
  • 3 = Significant
  • 4 = Extreme

11.2.4 1D Map Outputs

1D domain results can be output in combination with 2D domain(s) by using the 1d_wll GIS layer and the Read GIS WLL 1D command. 1d_wll GIS layer(s) are used to define and control the 1D map output. The layer(s) contain lines (called Water Level Lines or WLLs) that cross or snap to 1D channels and/or nodes. A WLL is essentially a line of horizontal water level, and should be digitised on this basis (i.e. perpendicular to the flow direction).

The direction of WLLs is important. They must be digitised from left to right looking in the positive direction of the digitised channel.

When viewing the results, if the 1D WLLs and 2D domains overlap, the 1D results are displayed on top of the 2D results. However, depending on the viewing platform, when observing the scalar and vector magnitudes as the cursor is moved around, the 2D values maybe given precedence over the 1D where the overlap occurs.

Two WLL options are available. The preferred method can be specified by using the command WLL Approach. Of the two options, Method A and Method B, Method B allows more advanced and accurate mapping of 1D results in map formats and is the default and recommended method, and that described below. For documentation on Method A, refer to the TUFLOW 2010 Manual.

It should be noted that water level lines do NOT change the 1D hydraulic computations, they are purely used in order to display the 1D results in plan (2D) formats.

11.2.4.1 Water Level Lines

Ground elevations and optionally material (Manning’s n) values can be assigned to points along a WLL. A more accurate representation of 1D domain velocity and flood hazard can be mapped using this approach. The velocity at a point on the WLL is estimated by carrying out a parallel channel analysis along the WLL using the flow in the channel the WLL is associated with as explained further below.

The parallel channel analysis assumes the water surface slope is that of the associated channel, and the water level at the WLL is linearly interpolated between the upstream and downstream nodes based on the digitised channel length.

WLLs can have any number of vertices. The association of WLLs to channels and nodes is as follows:

  • At the channel ends, to lock in the water level at a node, use a 3 vertex line with the middle vertex snapped to the 1D node. If you use a 3 vertex line across a channel, the channel “thalweg” is taken at the middle vertex, otherwise, for 2 vertex lines the mid-point of the WLL line is used.
  • If a WLL crosses two or more channels, the channel closest to the middle vertex (3 point line) or half-way point (2 point line) is used.
  • If a WLL middle vertex snaps to a node with two or more channels on the upstream side, the channel that is closest in angle to the WLL’s perpendicular (based on the WLL’s two end points) is used.
  • For 4 or more vertices, one of the middle vertices (i.e. not an end vertex) must snap to a vertex on the channel line.

There is one 1d_wll attribute required as described in 11.6. The attribute, Dist_for_Add_Points, is the minimum distance in metres along which to generate points for that WLL. If Dist_for_Add_Points is zero, only the vertices along the WLL are generated.

Estimates of the elevation, depth, velocity and other hydraulic outputs from the parallel channel analysis are produced at each WLL point including any automatically added ones. The user can subsequently assign elevations and materials to each point to improve the topographic representation if so desired (see Section 11.2.4.2).

An automatically generated triangulation (TIN) of the WLL points is used to display the outputs in map based formats.

Table 11.6: 1D Water Level Line (1d_wll) Attribute Descriptions
No Default GIS Attribute Name Description Type
1 Dist_for_Add_Points The minimum distance interval along the WLL to generate elevation and material sampling points (WLL Points). These points form the corners of the triangulation.
If set to zero, no additional points are generated (i.e. only the existing vertices along the WLL are used).
Float

Use the 1D command Read GIS WLL to specify the 1d_wll layer and automatically create 1D map output for TUFLOW 1D domains. For Flood Modeller or XPSWMM linked TUFLOW models, use Read GIS ISIS WLL or Read GIS XP WLL in the .tcf file to read the 1d_wll layer(s). The Flood Modeller units or XPSWMM links will also need to be built into a GIS layer and read using Read GIS ISIS Network or Read GIS XP Network from the .tcf file.

Several 1d_wll layers can be specified covering different sections of the 1D domain(s) if required. The TIN of triangular elements is created from the WLL points and included in the .2dm mesh file that can be viewed using the QGIS TUFLOW Viewer or in SMS to check they have been created correctly. In particular, look for overlapping or strange triangulations that arise from WLLs not being digitised correctly.

The default approach uses the processed cross-section data (height–width data) from the channel for setting elevations at each WLL point. For Flood Modeller and XPSWMM, these data are automatically fed through to TUFLOW and is based on the cross-section information entered into the Flood Modeller/XPSWMM model.

If a WLL is snapped to a node, the processed data used for setting any bed elevations is from the higher channel unless it is a bridge in which case it uses the bridge processed data.

WLL Automatic == CULVERTS is particularly useful for automatically generating WLLs for pipe network systems. WLL No Weirs is useful for not assigning WLLs to weir channels where they are in parallel to another channel (this is particularly useful for BW, CW and RW channels).

11.2.4.2 Water Level Line Points

If Write Check Files in the .tcf file is specified, two GIS check layers are created from the WLLs. These are labelled 1d_WLLo and 1d_WLLp. 1d_WLLo (Water Level Line Objects) reproduces the WLLs with attributes containing the channel and node the WLL is allocated to for cross-checking purposes.

1d_WLLp layers (Water Level Line Points) contain all of the elevation points generated based on the Dist_for_Add_Points attribute. This layer can then be used to allocate elevations (first attribute) to each point from a DEM (in the same manner that 2D Zpts can be assigned elevations).

A second attribute, RR, contains the relative resistance of each point (which will have a value of 1 when first generated). The RR attribute can be replaced by the integer material value at each point by using GIS to assign values from material polygons. The material value must exist in the .tmf file or a materials.csv file (see Read Materials File).

The attributes of a 1d_WLLp layer created by Write Check Files and used in Read GIS WLL Points, Read GIS ISIS WLL Points or Read GIS XP WLL Points are listed in Table 11.7.

The parallel channel analysis carried out for each WLL in future simulations will then use the WLL point data to estimate velocities and other outputs across the WLL instead of using the default approach of using processed hydraulic data as previously described. Note, the WLLs do NOT change the 1D hydraulic computations, they are purely used in order to display the 1D results in plan (2D) formats.

For Flood Modeller and XP-SWMM, the layers are essentially the same, but are named using xWLLo and xWLLp.

Note: If using Read GIS WLL Points or Read GIS X1D WLL Points, this layer must be a copy of the 1d_WLLp layer produced by Write Check Files. Points from this layer can be deleted, but not added. At deleted points, the default of estimating an elevation from the channel’s processed data is used. If the 1d_WLL layer is modified or any of the Dist_for_Add_Points attribute values changed, the 1d_WLLp layer needs to be regenerated and repopulated again.

For culvert channels (R and C channel type), only the end and mid vertices are used along the WLL, and the elevations are set to the culvert invert irrespective of the number of points along the WLL or the Dist_for_Add_Points value.

Table 11.7: 1D Water Level Line Point (1d_wllp) Attribute Descriptions
No Default GIS Attribute Name Description Type
1 Z Ground elevation of the point. Automatically generated from the channel cross-section processed data or updated based on a point inspection of a DEM or other approach. Float
2 RR or Material

In the 1d_WLLp check file, the relative resistance of the point. A value of 1 is assigned if the elevation was estimated from the channel’s processed data. If the elevation was provided through a point using Read GIS WLL Points, RR is the material Manning’s n value divided by the channel’s n value.

In a 1d_WLLp layer being used in Read GIS WLL Points or Read GIS X1D WLL Points, this column should either be set to an integer material value (normally sourced from a GIS layer of material polygons) – the material value must exist in the .tmf file (see Read Materials File).
Float

Elevation values along WLLs for bridge channels are always based on the processed data (i.e. any WLLp Z values are overridden) to ensure that the bridge deck underside is correctly represented. This has benefits when using the post-processing utility TUFLOW_to_GIS.exe (see Section 17.3) when extracting obverts of structures for longitudinal profiles.

A useful tip at a junction of 1D channels is to use a connector for any side channels (Type = “X” – see Section 5.9.3). Separate WLLs can then be allocated to the side channel and main channel removing the confusion that sometimes occurs in generating the triangulation between WLLs at junctions.

11.2.4.3 Adding Triangles into the 1d_WLL Layer

When using Method B, triangular regions can also be included in the 1d_WLL layer (as a separate region or polygon layer if using the .gpkg or .shp formats). The triangles are inserted into the WLL triangulation. This is particularly useful at junctions or between parallel channels to enhance the WLL TIN by facilitating a smooth transition of the water level surface at, for example, junctions or across meanders to infill areas and avoid steps in the output surfaces.

The triangles must snap to the ends of WLL line objects. In the example below, the WLLs are shown in red and the triangles are shown as yellow. These are connected (snapped) to the ends of the WLLs. The map output surfaces will interpolate over the triangle based on the hydraulic output values at the snapped WLLs.

If the region object is not correctly snapped an ERROR 1311 message occurs, pointing to the vertex on the triangle that is not snapped.

**Adding Triangles into 1d_WLL Layer to Infill Areas**

Figure 11.3: Adding Triangles into 1d_WLL Layer to Infill Areas

11.2.5 Output Zones

Output Zones are a powerful feature that allows users to generate map and check file outputs for sub-regions of a model. The region is defined using a polygon feature using the 2d_oz_empty template file. No attributes are associated with the polygon. One polygon per output zone and GIS layer is recommended. Note, multiple polygons within a single layer are not supported.

Different Output Zones can have different output formats, start/end times, frequency of output, and output data types. Any number of Output Zones can be defined for a model, and all, a selection of, or none of the Output Zones can be activated when the model is run.

Output zones are very useful, particularly for large models, where higher frequency map output is required for a portion of a model, for example to create an animationfor an urban region using a shorter map output interval. Another useful purpose is to create an Output Zone for a problematic section of model that requires closer examination, by generating output on a fine time interval during the period in question. Another benefit is simply reducing the size of the output files by only outputting where information is required, and disabling the output for the whole model using Map Output Entire Model == OFF. Output zones currently only apply to TUFLOW’s map and GIS check file outputs (i.e. does not apply to time series and other plot outputs).

Note: the Model Output Zones command is used to select which Output Zones are to be used (activated). If this command is omitted, no output from the Output Zones is written. Separate multiple Output Zones using a “|” (vertical bar or pipe character). For example, to output from zones ZoneA and ZoneC specify:

Model Output Zones == ZoneA | ZoneC

If there are numerous Output Zones, it is recommended that the definitions are placed in one or more separate files, and use the Read File command to reference these file(s). This will keep the size of the .tcf file to a minimum. This, however, is not a requirement.

Each Output Zone is defined using a definition block as follows:

Define Output Zone == <ZoneA>

…..

End Define

The following commands can be used within an Output Zone definition. With the exception of Read GIS Output Zone, all commands are optional. The optional commands can be used to change the default setting or the setting applied outside the Output Zone definitions.

Read GIS Output Zone
Map Output Format
Start Map Output
End Map Output
Map Output Interval
Map Output Data Types
Output Folder
Maximums and Minimums
Maximums and Minimums Only for Grid
Write Check Files

The Read GIS Output Zone reads a GIS layer containing one or more polygons that define the regions to be output. The attributes of the layer are not used.

Note: If one of the commands above is not specified within the Output Zone’s definition, the latest setting of that command prior to the Output Zone’s definition, or the global default if the command has not been used prior to that location in the .tcf file, will be used. If, for example, all output is to be in the .xmdf format, only specify Map Output Format == XMDF once and prior to any Output Zone definitions.

Map Output Entire Model can be used outside an Output Zone definition block to turn on or off map output for the entire model (the default is ON). If set to OFF only map output for any Output Zones is written. Map output commands that occur outside Output Zone definitions apply to the entire model output.

Example: Defining an Output Zone

The example below defines Output Zone “ZoneA”. In the example, a DEM_Z check grid will be the only check layer written for the Output Zone, and Grid, WRB and XMDF outputs will be produced. The Grid output will consist of water level (h) and VxD (Z0) results, starting at time 0.5 hours and ending at 4 hours at an interval of 0.5 hours. The WRB and XMDF output files will contain water levels (h), velocities (V) and the ZMBRC hazard categories starting at the simulation beginning and finishing at 6 hours at an interval of 6 minutes.

! Entire Model settings
Map Output Format == TIF XMDF
Map Output Interval == 3600 ! Default output interval is hourly
TIF Map Output Interval == 0 ! Maximums only for .tif output
Map Output Data Types == h v q d z0 ! Default output types

Model Output Zones == ZoneA | ZoneB ! Output both of the zones below

Define Output Zone == ZoneA

  Read GIS Output Zone == ..\2d_oz_ZoneA.shp
  Write Check Files Include == DEM_Z
  Map Output Formats == TIF XMDF WRB
  End Map Output == 6.0
  Map Output Interval == 360
  TIF Start Map Output == 0.5
  TIF End Map Output == 4.0
  TIF Map Output Interval == 1800
  XMDF Map Output Data Types == h V ZMBRC
  TIF Map Output Data Types == h Z0

End Define

Define Output Zone == ZoneB

  Read GIS Output Zone == ..\2d_oz_ZoneB.shp
  Map Output Formats == XMDF
  Map Output Interval == 360
  XMDF Map Output Data Types == h V ZMBRC

End Define

11.2.6 Gauge Level Output

The 2d_glo GIS layer uses the .tcf command Read GIS GLO and writes map-based output data when the water level at the gauge reaches user defined levels.

Note: Map output based on reaching gauge levels replaces the conventional approach of using a Start Map Output time and a Map Output Interval.

The 2d_glo GIS layer is used to define the location of the gauge within the modelled extent. The gauge is digitised as a point object within a 2d_glo layer with the attributes as described in Table 11.8 and referenced within the .tcf using the command Read GIS GLO.

For TUFLOW Classic only, the water level at the gauge location is tracked by TUFLOW throughout the simulation and output to the Console DOS Window (see Section 14.2) preceded by the letters “GL”.

When the water level reaches a user specified value, the map output results at that simulation time will be written. The results are therefore useful for mapping the predicted area of inundation for specified gauge heights. The range of water levels at the gauge for which the results are written may be defined using one of two methods:

  • By referencing a .csv file within the 2d_glo GIS layer using the first attribute. The .csv file contains a single column of levels, and comment lines are allowed using the “!” or “#” character. If a .csv file is specified in the first attribute, the remaining attributes are not used, but dummy or empty values for these attributes must exist.
  • By populating the remaining attributes of the 2d_glo GIS layer as described below.

When using the attributes of the 2d_glo GIS layer to define the gauge heights, the map output results are first written when the water level at the gauge reaches the specified “Bottom_elevation”. Subsequent results are written as the water level at the gauge rises based on the value of the “increment” attribute. The map output ceases once the water level at the gauge reaches the “Top_elevation”.

For example, if the 2d_glo GIS layer has been defined as shown below, map output results will be written when the water level at the gauge reaches 41m, 41.5m, 42m, 42.5m and 43m.

Only one gauge location may be specified per model simulation. If more than one object exists within the 2d_glo layer, the gauge that is monitored by TUFLOW will be the last digitised point object. Similarly, if the Read GIS GLO command is used more than once, only the last occurrence of the command will be used.

Table 11.8: 2D Gauge Level Output (2d_glo) Attribute Descriptions
No Default GIS Attribute Name Description Type
1 Datafile

Filename (and relative or full path if needed) of the file containing the gauge levels to trigger map output. Must be a comma or space delimited text file such as a .csv file. Only the first column is used, and this column must contain the gauge levels. Comment lines using a leading “!” or “#” can be used at any point within the file.

If this attribute is blank, the following three attributes are used to define the gauge levels to trigger map output. If this attribute is not blank, the following three attributes are not used, but they must exist and can simply be populated with empty or default values.
Char (254)
2 Bottom_Elevation The water level in m above datum at the gauge at which the writing of map output results will commence. Float
3 Top_Elevation The water level in m above datum at the gauge at which writing of the map output results will cease. Float
4 Increment The water level increments in metres between the Bottom_Elevation and Top_Elevation at which the map output results will be written. Float

11.3 Time-Series Outputs

Time-series results (sometimes referred to as Plot Output) produces output for graphing in charts and profiles. The data output are specified using the Map Output Data Types command, at an interval specified by the Time Series Output Interval command.

Time-Series data can be output for 1D domains (Section 11.3.1) and 2D domains (Section 11.3.2). If the total flow across the floodplain is required for a 1D/2D model where the river is in 1D and the floodplain in 2D, flows can be combined using Reporting Locations (Section 11.3.3). The Structure Reporting feature outputs time-series and summary data for structures (Section 11.3.4).

There are two formats available for time-series outputs: comma-separated values (.csv) (the default) and NetCDF (.nc). One of the advantages of NetCDF is all the timeseries output is in a single compressed .nc file, rather than multiple uncompressed .csv files. To change the format use Time Series Output Format command.

It is possible to write 1D and 2D time-series outputs as the simulation progresses, using the Write PO Online command.

11.3.1 1D Time-Series Output

Time-series data output from 1D domains is available for a range of hydraulic parameters. The Output Data Types ECF command or the 1D Output Data Types TCF command controls the types to output. The options are:

  • A: flow area at channels (m2);
  • E: energy at nodes (m);
  • H: water level at nodes (m);
  • Q: flow at channels (m3/s);
  • QI: integral flow at channels (m3);
  • S: structure and grouped structure output (see Section 11.3.4);
  • V: velocity at channels (m/s); and
  • Vol: volume at nodes (m3).

To review the 1D time-series data see 15.3.

11.3.2 2D Time-Series Output

Time-series data output from 2D domains is available for a range of hydraulic parameters (as listed in Table 11.9). Output takes the form of time-series hydrographs (referred to as PO – Plot Output) or longitudinal profiles (LP) over time.

This is carried out by creating one or more GIS layers containing points, lines and regions that define the locations of PO and LP output. Figure 11.4 illustrates how 2d_po objects are interpreted.

The start time for PO and LP output and the output interval is set using Start Time Series Output and Time Series Output Interval. If no start time is specified the simulation start time is used. If no output interval is specified the simulation will stop with ERROR 0046 to prevent excessive amounts of memory and disk space from being used.

The output is written to a .csv file and also to the _TS layer (refer to Section 15.3.4). 2D domain time-series (PO) output is synchronised with 1D domain output by default. This allows both 1D and 2D time-series to be placed in the _TS layer. Set Output Times Same as 2D to OFF in the .ecf file if 1D and 2D time-series data is not to be synchronised. In this case, no 2D PO is written to the _TS layer.

Maximums and minimums are output to four additional rows near the top of the _PO.csv file, and columns in the _LP.csv files, containing the Maximum, Time of Maximum, Minimum, and Time of Minimum values. The _TS GIS layer also contains the tracked values. For TUFLOW Classic this information is tracked every computational timestep. For TUFLOW HPC the maximum/minimum values are post-processed at the end of the simulation based on the Time Series Output Interval, not every computational timestep. Tracking of maximums and minimums can be disabled by setting the Maximums and Minimums Time Series command to OFF in the .tcf file.

**Interpretation of PO Objects and Map Output**

Figure 11.4: Interpretation of PO Objects and Map Output

11.3.2.1 Plot Output

The data types available for PO are as listed in Table 11.9, noting that only supported object types (i.e. Point, Line and Region) for the data type are documented in the table. PO data is read into a simulation using the Read GIS PO .tcf command. This is carried out by creating one or more GIS layers containing points, lines and regions that define the locations of PO and LP output. Figure 11.4 illustrates how 2d_po objects are interpreted.

Table 11.9: Time-Series (PO) Data Types - Point, Line or Region
Flag Description Comments
CI Cumulative Infiltration Point: Represents the current water content of each layer (mm or in of pure water). This value may increase or decrease depending on the flows into and out of the cell (it is not the cumulative value).
D_ Depth

Point: Depth of the nearest cell.

Line: The average depth of all wet cells along the line. Note: In TUFLOW Classic if a line with more than 2 vertices (i.e. a polyline) is used, the average water level along each line segment is output, therefore, use of polylines is not recommended for this output type at present. TUFLOW HPC, will take the average depth along the entire polyline.
G_ Gauge Level Point: Water level at the cell center of the nearest cell. If the cell is dry, the ground level (ZC) is output. Used for Read GIS Objects to record gauge levels when a receptor is first inundated. Refer to Section 11.3.3.
GWd Groundwater Depth Point: Depth of water within sub-surface layer(s) (cumulative infiltration divided by porosity).
GWh Groundwater Level Point: Elevation of groundwater surface (water table) within sub-surface layer(s), reported in m or ft.
GWm Groundwater Moisture Point: A dimensionless number in range 0 - 1 representing “fraction full” of the sub-surface layer(s).
GWq Groundwater Unit Flow

Point: Unit flow magnitude within the sub-surface layer(s).

Line: Total flow passing through the given line within each sub-surface layer(s). The positive flow convention is left to right looking downstream (same as the Q_ type).
Gwqa Groundwater Unit Flow Angle Point: Angle of unit flow vector within each layer. Reported in degrees clockwise from north (i.e. a compass bearing).
Gwqu Groundwater Unit Flow
U-Component
Point: Unit flow u-component within each sub-surface layer(s).
GWqv Groundwater Unit Flow
V-Component
Point: Unit flow v-component within each sub-surface layer(s).
GWVol Groundwater Volume Region: Total volume of water within the polygon for each layer (m3 or ft3).
H_ Water Level (Head) Point: Water level of the nearest cell. If the cell is dry, the ground level (ZC) is output.
HAvg Average Water Level Region: The average water level within the region (wet cells only).
HD HU Downstream and Upstream Structure Water Levels

Point:
Used to define the upstream (HU) and downstream (HD) water levels of a QS structure (see QS below and Section 11.3.4). The 2D water level is used to populate the upstream and downstream water level data in the new structure output feature. For a point the water level at the nearest 2D cell centre is used.

To associate the HU and HD objects with the QS line, all three (QS, HU and HD) must have the same ID for the 2d_po Label attribute (see Table 11.10).

Note that the water levels over time are output to a 2D PO .csv file and the summary information at the flood peak to the new _SHmx.csv file. If the point or line is completely dry, -99999 is output to the .csv files. The _SHmx.csv file is currently not produced for 2D only models.
HMax Maximum Water Level Region: The maximum water level within the region.
Q_ Flow or discharge

Line: The flow crossing the line.

The flow is determined by summing the flow across cell sides whose perpendiculars intersect the line (see Figure 11.4).

The flow is positive if the water is flowing away from you when looking in a direction with the start of the PO line on your left and the end of the line on your right.

If digitising a flow line across a 1D channel that is carved through the 2D domain, ensure that the line is digitised so that it crosses the 1D channel where there is a change in colour in linked 2D HX cells as shown in the 1d_to_2d_check or _TSMB1d2d layers (thematically map using the Primary Node if not using MapInfo). The 1D flow can then be added manually in a spreadsheet. However, note that Read GIS Reporting Location is now the preferred method for cumulating flows across a line for 1D/2D models. See Section 11.3.3 for discussion on Reporting Locations.
QA Flow Area Line: The flow area is calculated using the same cell sides as for Q_. An adjustment for oblique lines is made.
QI Integral Flow Line: Integrates the flow (as determined for Q_ above) over time (i.e. the area under a Q_ time-series curve) (i.e. cumulative volume). If Write PO Online is set to ON, the integral flow is not calculated until the simulation is complete.
Qin Flow In Region: The flow into a region.
Qout Flow Out Region: The flow out of a region.
QS Structure Flow Line: Same as Q_ above, but also used to set up a 2D structure output (see Section 11.3.4) that will include in addition to the 2D flow any flows from intersected 1D structures and the split between below and above deck flows. Note the flow output to the 2D PO .csv files is only the 2D flow, while that to the new _SQ.csv file is the combined 1D/2D structure flow. The _SQ.csv file is currently not produced for 2D only models.
QX Flow in X-direction. Line: The X component of Q_ (i.e. the sum of the flows at the u-points).
QY Flow in Y-direction. Line: The Y component of Q_ (i.e. the sum of the flows at the v-points).
SS Sink Source Region: Sink / source flows applied within the region (rainfall, infiltration, source area inflow and SX flows).
V_ Velocity

Point: The magnitude of the resolved vector based on the two u-points and two v-points of the cell in which the point falls. Exactly which cell is selected may be uncertain if the point falls exactly on a cell’s side.

Line: Velocity as Q_/QA (i.e. the depth and width averaged velocity along the line). Prior to release 2020-10-AA it used the cell in which the line starts.
VA Velocity Angle Point: The angle of V_ (degrees relative to east where east is zero, north is 90, etc.).
Vol Volume Region: Total volume within the region.
Vu or Uu u-point velocity Point: The magnitude of the u-point velocity (i.e. across the right hand side of the cell).
Vv v-point velocity Point: The magnitude of the v-point velocity (i.e. across the top side of the cell).
VX Velocity in X-direction Point: The magnitude of the average of the u-point velocities (i.e. across the left and right hand sides of the cell).
VY Velocity in Y-direction Point: The magnitude of the average of the v-point velocities (i.e. across the bottom and top sides of the cell).

Table 11.10 describes the GIS attributes of the 2d_po layer.

Table 11.10: 2D Plot Output (2d_po) Attribute Descriptions
No Default GIS Attribute Name Description Type
1 Type Any combination of the two letter flags listed in Column 1 of Table 11.9 (limit of 10 flags per entry). For example, to output water level and flow time-series for the same line, enter “H_Q_”. Char (20)
2 Label

Label up to 30 characters long defining the name of the time-series. The label appears at the top of the columns of data in the _PO.csv file. Spaces are permitted, but do not use commas.

Read GIS Reporting Location is the preferred method for cumulating flows in 2D1D models. See Section 11.3.3 for discussion on Reporting Locations.
Char (30)
3 Comment Optional field for entering comments. Not used. Char (250)

11.3.2.2 Long Profile Output

H_ (water level) and V_ (velocity) are the only data types available for long profile (LP) outputs. LP locations and data types are initiated for a simulation using the Read GIS LP .tcf command.

Table 11.11 describes the GIS attributes of the 2d_lp layer. The 2d_lp layer(s) contain lines defining where the profile data are to be generated. Each line is given a label to uniquely define the profile in the output. The starting vertex of the line will set the distance origin for the profile.

The advantage of having TUFLOW generate the profiles directly rather than post-processing them is the outputs will be slightly more accurate due to no post-processing interpolation rounding, plus if the profile(s) are repeatedly being plotted, the plotting process can be streamlined via python scripts that use the LP outputs.

Table 11.11: 2D Longitudinal Profiles (2d_lp) Attribute Descriptions
No Default GIS Attribute Name Description Type
1 Type Specify “H_” to output water level. Specify “V_” to output velocity. Char (20)
2 Label Label up to 30 characters defining the name of the longitudinal profile. The label appears at the top of the columns of data in the _LP.csv file. Spaces are permitted. Commas are not permitted. Char (30)
3 Comment Optional field for entering comments. Not used. Char (250)

11.3.3 Reporting Locations

Read GIS Reporting Location allows for plotting of time-series results that automatically combines 1D and 2D outputs. For example, it is possible to digitise a reporting location line that extends across 1D and 2D domains, including multiple 2D domains, and TUFLOW will sum the flow across any 1D channels intersected by the line and all the 2D cells. This will save the need to post-process time-series output from 1D and 2D domains to accumulate the flow.

Reporting location lines are digitised into a 0d_RL GIS layer containing only a single attribute, the name of the reporting location, as outlined in Table 11.12. Points, lines and regions can be used. Points will be treated as water level output, lines as flow output and regions as volume output. For a point snapped to a 1D node, the 1D water level is used, if no 1D node is snapped a 2D water level is output.

The flow line can cross 1D and 2D sections of the model, for the 1D channels it does not have to snap to any vertices on the channels, it just needs to intersect them.

Table 11.12: 0d_RL Reporting Location Attributes
No Default GIS Attribute Name Description Type
1 Name The name of the reporting location. Lines and points can share the same name. Char (32)

The RL outputs are written to the “plot\csv\” folder. The following files are produced:

  • _RLL_Q.csv - flow time-series;
  • _RLL_Qmx.csv - maximum flow information;
  • _RLP_H.csv – water level time series;
  • _RLP_Hmx.csv – maximum water level information;
  • _RLR_Vol.csv - volume time-series; and
  • _RLR_Volmx.csv - maximum volume information.

As well as the maximum water level and flow information, the time that these occur, the water level at maximum flow and vice versa, and the maximum change between timesteps are also output to the mx.csv files.

The RLs are also output to the plot\gis PLOT GIS layers and can be viewed and their time-series data displayed using the TUFLOW Viewer using the QGIS TUFLOW Plugin as illustrated in Figure 11.5.

An example model using reporting locations is available in the Example Models on the TUFLOW Wiki. In addition, see Section 15.4.3 for information on plotting of Reporting Location results.

**Example of the QGIS TUFLOW Plugin for a Reporting Location**

Figure 11.5: Example of the QGIS TUFLOW Plugin for a Reporting Location

11.3.4 Structure Output

The Structure Reporting feature outputs time-series and summary data for single and grouped structures. 1D and 2D structures are all output together to give a complete set of results. The summary output at the flood peak also produces the flow split between below and above deck, along with other information such as the head drop. The structure output is particularly useful in the reporting of hydraulic structure flows and afflux.

Structures are classified according to the following logic:

1D Structures:

  • 1D structures that are in parallel (i.e. two structures that link to the same upstream and downstream nodes) are automatically grouped together and treated as a single structure for this output. The ID assigned to the group is the structure with the lowest bed elevation. Note that the directions of the digitised channels is important, that is, to form a group all channels must be digitised in the same direction.
  • For 1D structures that have no parallel 1D structures, these are also included in the output so as to provide a complete set of results for all 1D structures.
  • The flow split between below and above deck is based on the structure geometry, except weirs contributing to the below deck flow and/or above deck flow, depending on the configuration.

2D and 1D/2D Structures:

  • To create a structure output that includes 2D flow, and optionally any 1D structures, the “QS” PO line is digitised in a 2d_po layer (see Section 11.3.2.1). All 2D flow across this line and any 1D structures that intersect this line are grouped together. The 1D structure’s 1d_nwk line does not have to snap with the QS line; they only have to cross over each other. The ID assigned to the structure output group is the 2d_po QS Label (see Table 11.10).
  • If the 2d_po QS line selects cell sides that are modified by the Layered Flow Constriction or 2D Bridge feature, the summary output will split the flow into a below and above “deck” component based on Layers 1 to 2 being below “deck” and Layer 3 and 4 above.
  • 2d_po HU and HD lines or points can be used to define the upstream and downstream water levels of the structure. HU objects should be located upstream of the structure and HD downstream. The average 2D water level along a line object will be used to populate the upstream and downstream water level data in the output. Lines can have more than two vertices (i.e. polylines are accepted). To associate the HU and HD objects with the QS line, all three (QS, HU and HD) must have the same ID for the 2d_po Label attribute Label (see Table 11.10). If a QS line has no HU and/or HD objects associated with it, output that cannot be produced, such as the water level drop across the structure, is given a -99999 value in the _SHmx.csv output file described below. HU and HD inputs are necessary for 2d_lfcsh and 2d_bg bridges.

An example model of the structure output for a bridge represented by a 2d_lfcsh is available in the Example Models on the TUFLOW Wiki.

The structure output is located in the plot/csv folder and includes:

  • _SQ.csv file that contains time-series data of the flow through the structure. This file is similar to other time-series .csv output, but it only contains 1D and/or 2D structures as described above.
  • _SHmx.csv file that contains a summary of each structure when the upstream water level reached its maximum. To generate this output the flow and other information is tracked every timestep for grouped structures. The output columns include: flow, area and average velocity for below and above deck; total flow, area and average velocity for the whole structure; upstream and downstream water levels; the head drop across the structure (i.e. upstream minus downstream water level); and the time these data were recorded (e.g. the time the upstream water level peaked).
  • Time-series output for the upstream and downstream water levels are available through the 1d_H.csv, 2d_HD.csv and 2d_HU.csv files. Note that for the 2d_HD.csv and 2d_HU.csv files, if all 2D cells are dry a -99999 is output.

See Section 15.4.3 for further information on plotting of grouped structure results.

Two structure group check files are output if a model contains any structure groups (either automatically created, or via a “QS” type line in a 2d_po layer). The check files are both .csv files as follows:

  • <simulation name>_Str_Grp_All.csv, contains information for all structure groups, including single 1D structures.
  • <simulation name>_Str_Grp_Multi.csv, contains information for structure groups that are comprised of more than one 1D channel or are generated from a 2d_po “QS” line.

For more information on the structure group check files, please see the Check File Page on the TUFLOW Wiki.

11.4 Specialised Outputs

Specialised outputs include:

  • Receptors: records the flood level and time at one or more gauge(s) when the receptor is inundated above the specified trigger level, see Section 11.4.1.
  • Evacuation Routes: provides additional information (e.g. first point of closure) along routes, see Section 11.4.2.

11.4.1 Receptors

Read GIS Objects RECORD GAUGE DATA records the flood level and simulation time at one or more gauge(s) when receptors are first inundated above their trigger inundation levels (e.g. floor levels). Read GIS Receptors can be used as an alias to Read GIS Objects. Level output associates a flood level at one or more gauges with the time of first inundation at properties, buildings, or other areas of interest within the modelled extent.

When the water level at the property reaches a user-defined trigger inundation level, the gauge height and simulation time are recorded and tagged to the receptor. This is particularly useful for flood warning and forecasting studies where property specific information on the likelihood of a property being inundated for a given gauge height can be generated. An example of translating gauge data information to catchment receptors can be found in the TUFLOW Wiki Example Models.

Gauges are defined as a point within a 2d_po GIS layer with type “G_” (see Section 11.3.2 and Table 11.9). The levels from all gauges are recorded at each receptor once inundated.

Receptors must be GIS point and/or polygon objects located in one or more GIS layers nominally prefixed by 2d_obj or 2d_rec. Each object within the layer represents a receptor, property or other object of interest. For information on the attributes of the GIS layer see the Read GIS Objects .tgc command. The command also includes options to set the Zpt elevations to the receptor level or the first attribute in the layer (for example, to set the Zpts to the floor level of the buildings), or to alternatively use the existing ZC elevations.

Once simulated, a GIS layer is written to the location as specified using the .tcf command Output Folder, and has a _GDO extension standing for Gauge Data Output. The layer contains point objects (for regions the centroid is used). The attribute data for the layer are described in Table 11.13. Figure 11.6 shows an example of how the output GIS layers may be used to illustrate the flooding of properties in relation to the water level at a gauge.

As well as outputting the water levels at gauge objects (as described above), the water levels at all point reporting locations and flows for all line reporting locations are also recorded. The reporting locations can be 1D and/or 2D locations. See Section 11.3.3 for a full description of the reporting location functionality. For example, flows past a gauge can be recorded, including combined flows from 1D channels and the 2D domain(s).

Table 11.13: 2d_GDO_ Gauge Data Output Attributes
GIS Attribute Description
<gauge_1> The water level in metres or feet above datum at <gauge_1> (as defined in a 2d_po layer) when the water level first reaches the Trigger_Level at the property/building/object.
<gauge_n>… If a second, third, fourth… gauge exists, the water level in metres above datum at each gauge.
Time The simulation time in hours when the water level at the receptor first reaches the Trigger_Level.
Trigger_Level The trigger inundation level at the receptor. The level will be either the first attribute in the GIS layer (typically either the ground level of the property or the floor level of a building), or the lowest ZC value within the property or building polygon if the “USE ZC” option is specified (i.e. Read GIS Objects USE ZC).
Use_ZC “Y” if the ZC 2D cell elevation was used for the Trigger_Level; “N” if not.
**Example Use of Gauge Data Output Layer**

Figure 11.6: Example Use of Gauge Data Output Layer

11.4.2 Evacuation Routes

Evacuation routes can be specified to output information relating to:

  • Evacuation route suitability;
  • Warning times;
  • Risks;
  • Degree of route inundation; and
  • Duration of inundation.

The .tgc Read GIS Z Shape Route is used to define the routes and, by default, also adjust the Zpt elevations along the route using the standard Z Shape options. The 2d_zshr layer is the same as a 2d_zsh layer, but with three additional attributes, as shown in Table 11.14. Note, if using the TUFLOW HPC SGS functionality (Section 7.4.3, the cells/faces selected by this command are assumed flat (e.g. the traditional sampling approach using a single cell centre/face elevation).

The commands Set Route Cut Off Values and Set Route Cut Off Type can be used in the .tcf and/or .tgc files. If used in the .tcf file, this sets the default values for the 2d_zshr Cut_Off_Type and Cut_Off_Values attributes if these attributes are left blank. The default values can be changed between different Read GIS Z Shape Route commands in the .tgc file by repeat usage of the same commands.

The _RCP output layer is a layer of points showing where a route’s cut off value(s) were first exceeded (e.g. first point of closure). The layer contains the attributes as shown in Table 11.15.

The _RC.xmdf output file contains the Route Categories over time and is written when adding RC to the .tcf command Map Output Data Types. This file can be used to view and animate the route category values.

A range of evacuation route examples are available in the Example Models on the TUFLOW Wiki.

Table 11.14: 2D Z-Shape Route (2d_zshr) Attribute Descriptions
No Default GIS Attribute Name Description Type
1 Z Unless the BRIDGE Shape_Options is specified (see below), the Z Shape lines adjust the Zpts as described in Table 7.5 for the same attribute. Float
2 dZ Unless the BRIDGE Shape_Options is specified (see below), the Z Shape lines adjust the Zpts as described in Table 7.5 for the same attribute. Float
3 Shape_Width_or_dMax Refer to the same attribute in Table 7.5. Float
4 Shape_Options BRIDGE: Does not use the Z attribute to adjust the Zpts, instead uses it to assign evacuation route categories. This allows for a Z Shape Route Line to cross a river without physically blocking it. Because a route can be comprised of more than one line (as long as each line is given the same Route_Name) they will all be regarded as being part of the same route. Therefore, to represent a bridge, split the route line either side of the bridge and give all three lines the same Route_Name, with the line representing the bridge or river crossing assigned a Shape_Options of “BRIDGE”. Char(20)
5 Route_Name Used to label the evacuation route. A route can be split into several lines if required, provided all the lines have the same Route_Name attribute. This can be useful where the route is more easily sourced or digitised as several lines, or if using the BRIDGE option described above. Char(40)
6 Cut_Off_Type

Options are as follows:

  • Blank: The value from the latest Set Route Cut Off Type command is applied at that point in the .tgc file.
  • Depth: Depth is used to define the cut-off (default setting).
  • VxD, Z0 or Hazard: The velocity x depth hazard product is used.
  • V or Velocity: The velocity is used.
  • Energy: Energy depth is used. The energy-depth is calculated as:

\[Energy Depth = d + \frac{V^2}{2g}\]

More than one 2d_zshr layer may be used if different cut-off types are required.
Char(40)
7 Cut_Off_Values A comma delimited list of one or more values (e.g. depths) used to categorise the severity of the overtopping along the routes. For example, if “0.1, 0.3, 0.7” is specified for a depth Cut_Off_Type, then where the water depth exceeds 0.1 m, these sections of the route are assigned a Category 1; above 0.3 m deep Category 2; and above 0.7 m Category 3. The values must be in ascending order. Elsewhere the route is assigned Category 0 (i.e. no or less than Category 1 overtopping). The Route Categories are output over time in the _RC.dat file and summarised in the _RCP layer. Multiple cut-off values allow for the assessment of different risk levels (e.g. shallow depths would be acceptable for most vehicles and people to safely negotiate, while deeper depths would only be acceptable for higher set vehicles). Char(80)
Table 11.15: RCP Output (2d_zshr) Attribute Description
GIS Attribute Description
Route_Name The name of the route.
Cut_Off_Value The cut off value applied.
First_Cut_Off_Time The simulation time in hours the Cut_Off_Value was first exceeded.
Last_Cut_Off_Time The simulation time in hours when the Cut_Off_Value was last exceeded.
Duration_Cut_Off The duration in hours that the Cut_Off_Value was exceeded – not necessarily the difference between the first and last cut off times if the route reopened during this time.

11.5 Check and Log Files

The options to control check and log file outputs are discussed in Chapter 14, with common commands to control these outputs listed in Table 11.16.

11.6 Output Control Commands

A wide range and variety of commands allow the user to configure TUFLOW output to be different to the default settings. Table 11.16 lists these commands in different categories followed by a brief description.

Table 11.16: Commands used to Control TUFLOW Outputs
Requirement Command Description Category
Mandatory GPKG Projection, SHP Projection, or MI Projection Define the projection for the output of GIS layers. Not all projection commands are mandatory, i.e. only specify projection settings for the formats required. Commands are: GPKG Projection, MI Projection, SHP Projection, GIS Projection Check, SHP Projection Check Method, TIF Projection. GIS
Mandatory Map Output Format Specify one or more map output formats to be written. Map
Mandatory Map Output Interval Mandatory command that controls the frequency of map output. Alternatively, set this value to zero to only output maximums. Map
Mandatory Time Series Output Interval Mandatory command that controls the frequency of time-series output. Time-Series
Recommended Write Check Files Customise which check files/layers to output. Check
Recommended Log Folder Customise the location for the .tlf and other files. File Management
Recommended Output Folder Customise the parent folder for the 2D based output files. If the command is place in the .ecf or 1D Domain controls the folder for the 1D based output files. File Management
Recommended GIS Format Controls the format for output GIS layers. GIS
Recommended Map Output Data Types Controls the map data types to be output. Default is just water level and velocity, but there are many more! Map
Optional Pit Channel Offset Change the default distance that a pit and its results are displayed relative to the 1D node the pit flows into. Check
Optional Write X1D Check Files Set to ON to output additional check files for linked external 1D schemes. Check
Optional Mass Balance Output Switch mass balance tracking and outputting on or off. Console
Optional Mass Balance Output Interval Set the time interval for outputting mass balance information. Console
Optional Screen/Log Display Interval Controls how often information is displayed and written to the .tlf file whilst a simulation is underway. Console
Optional Output Drive Override the drive letter for all output files. File Management
Optional Output Files Control which output files are to be written. File Management
Optional Simulations Log Folder Customise the location for logging simulations. File Management
Optional Write Empty GIS Files Output empty GIS layers. Very useful command when setting up a new model. File Management
Optional Write Restart File at Time Controls how and when restart files are output. File Management
Optional Write Restart File Interval Sets the interval at which to write a restart file. File Management
Optional Write Restart File Version Controls the version of the restart file feature to output. File Management
Optional Write Restart Filename Controls the approach on whether to keep overwriting or give a unique name for each restart file. File Management
Optional XF Files Global setting on whether or not to use and output XF files. File Management
Optional XF Files Include in Filename Set a unique text to be included in XF filenames. File Management
Optional Read GIS Objects Add gauge levels and other information to GIS receptor layer(s). GIS
Optional Write GIS Domain Output a GIS layer of the 2D domains processed thus far in the .tgc file. GIS
Optional Write GIS Zpts Output a GIS layer of the current 2D Zpts elevations at that point in the .tgc file. GIS
Optional Spatial Database Output Sets how the outputs are written when using GeoPackages. The grouped (default) option will group outputs by folder location. For example all check file outputs will be contained within one database. GIS
Optional Grid Format Controls the format for output GIS grid layers. Grid
Optional Grid Output Cell Size Sets the cell size of output grids, including check grids. This can be different to the 2D cell size(s). Grid
Optional Grid Output Origin Change how the grid output origin is determined. Grid
Optional Write GIS Grid Output a GIS layer of the current active 2D cells and their attributes at that point in the .tgc file. Grid
Optional UK Hazard Debris Factor
UK Hazard Formula
UK Hazard Land Use
Change the default settings on how the UK flood hazard output is calculated. Hazard
Optional ZP Hazard Cutoff Depth Change the cutoff depth for outputting People Hazard categories (ZPA, ZPC and ZPI). Hazard
Optional BSS Cutoff Depth Change the default cutoff depth setting for bed shear stress output. Map
Optional End Map Output Change the end time for map output to be different to the simulation end time. Useful for output zones where a shorter output period is required. Map
Optional Map Cutoff Depth Change the default cutoff depth for setting whether a 2D cell is wet or dry in the map output. Useful for direct rainfall modelling. Only used for map output, ie. it does not affect the hydraulic computations like Cell Wet/Dry Depth does. Map
Optional Map Output Entire Model Switch on or off the map output for the entire model. Useful if using output zones and map output for the whole model is not needed. Map
Optional Maximum Velocity Cutoff Depth Change the default depth that the approach to tracking maximum velocities is changed. Map
Optional Read GIS GLO Change from regular time-based map output to map output produced at different water levels at a gauge. GLO stands for Gauge Level Output. Map
Optional Read GIS WLL Points
Read GIS X1D WLL Points
Change the ground elevations and materials at WLL points along WLLs. Map
Optional

Read GIS WLL

Read GIS X1D WLL
Include the results from 1D elements in the map output using WLLs. Map
Optional Start Map Output Change the start of map output to be different to the simulation start time. Map
Optional Time Output Cutoff Add time of exceedance and time of duration for one or more depth or hazard values to the map output data set. Map
Optional Meshparts Adds flags to the .2dm file for use by some map output viewers such as SMS to control the display of different parts of a model. Other
Optional Reveal 1D Nodes Add the hidden 1D nodes used in 2D/2D links to the 1D output. Other
Optional Set Route Cut Off Type Controls the parameter to use when determining whether a route is closed. For example, depth or hazard. Other
Optional Set Route Cut Off Values Define one or more values for route closure categories. Other
Optional CSV Time Change the units of the time column. Time-Series
Optional Excel Start Date Start the time column at a specific date for use in Excel. Time-Series
Optional Order Output Controls whether 1D nodes and channels are output in alphanumeric order or in the order they are processed during input. Time-Series
Optional Read GIS LP Set up 2D longitudinal profile output. Time-Series
Optional Read GIS PO Set up 2D time-series output locations. Time-Series
Optional Read GIS Reporting Location Set up combined 1D and 2D time-series output locations. Time-Series
Optional Start Time Series Output Change the start of time-series output to be different to the simulation start time. Time-Series
Optional Write PO Online Write time-series output at each time map output is written (rather than just at the end of the simulation). Time-Series
Optional Maximums and Minimums Change the default setting on whether to track maximums and/or minimums. Tracking Max/Min
Optional Maximums and Minimums Only for Grids Controls whether grid based map output is only written out as maximums, ie. switch on or off the grid temporal map output. Tracking Max/Min
Optional Maximums and Minimums Time Series Controls whether to track maximums for time-series outputs. Tracking Max/Min
Advanced Display Water Level Display a water level on the console window whilst running. Console
Advanced Force File IO Display Show file opening information to the console window. Console
Advanced Verbose Controls how much information is displayed and written to the .tlf file. Console
Advanced Read File Insert a file into a control file. File Management
Advanced Set Variable Define user variables for use within the control files. Useful for setting output folder names based on a scenario. File Management
Advanced TSF Update Interval Change when and how often the .tsf file is to be updated. File Management
Advanced Calibration Points GIS File Add the maximum water levels to a GIS layer of calibration or gauge points. GIS
Advanced Map Output Corner Interpolation Apply different approaches to interpolating results to 2D cell corners. Primarily provided for backward compatibility for legacy models. Map
Advanced WLLp Interpolate Bed Change how the bed level along a WLL is assigned. Map
Advanced Define Output Zone

End Define
Set up a definition block of commands for an output zone. Output Zones
Advanced Model Output Zones Controls which output zones are output. Output Zones
Advanced Read GIS Output Zone Defines the region for the output zone. Output Zones

References

Australian Emergency Management Institute. (2014). Australian Emergency Management Handbook 7: Managing Floodplain Best Practice in Flood Risk Management in Australia. https://knowledge.aidr.org.au/media/3521/adr-handbook-7.pdf
Cox, R., Shard, T., & Blacka, M. (2010). Australian Rainfall and Runoff Revision Project 10: Appropriate Safety Criteria for People. Austalian Institute of Engineers. https://arr.ga.gov.au/__data/assets/pdf_file/0006/40578/ARR_Project_10_Stage1_report_Final.pdf
CSIRO. (2000). Floodplain Management in Australia Best Practice Principles and Guidelines. CSIRO Publishing.
Environment Agency. (2006). UK FD2321 Technical Report - Flood Risks to People: Phase 2 Project Record. UK Government; https://assets.publishing.service.gov.uk/media/602bbb768fa8f50386a7f8aa/Flood_risks_to_people_-_Phase_2_Project_Record.pdf.
Environment Agency. (2008). Explanatory Note for FD2320 and FD2321 Project Record. https://assets.publishing.service.gov.uk/media/602bbdcfe90e070561b31432/Explanatory_note_for_FD2320_and_FD2321_project_record.pdf.
GHD. (2011). MBRC Storm Tide Management Study: Stage 1 – Scoping Report.
Melbourne Water. (2016). Flood Mapping Projects Guidelines and Technical Specifications (Version 1).
NSW Goverment. (2005b). Floodplain Management Manual. NSW Department of Land; Water Conservation (NSW Government).
Port Macquarie-Hastings Council. (2018). PORT MACQUARIE-HASTINGS COUNCIL FLOOD POLICY. https://www.pmhc.nsw.gov.au/files/assets/public/v/1/document-files/your-council/publications/policies/port-macquarie-hastings-council-flood-adopted-2018-12-12.pdf.
Queensland Reconstruction Authority. (2012). Planning for Stronger, More Resilient Fl Oodplains: Part 2. https://www.qra.qld.gov.au/sites/default/files/2018-10/resilient-floodplains-part2-full_0.pdf.