In this tutorial model we will complete a simplistic flood impact assessment. We will start with the model developed in Module 4 and modify this to represent a "developed" scenario. This will involve changing the geometry in the 2D model domain to:
- Remove a portion of the area from the floodplain (fill for development); and
- Lower an area to increase conveyance.
Once the model has been simulated, the pre and post development cases will be compared to determine the changes in peak water level and flood behaviour caused by the development.
In this module, we will also use an alternative 2D output type which stores all the results in a single output file (as opposed to individual data files) and use a folder structure to make the results files easier to navigate.
In a later tutorial module we will look at a more detailed impact assessment in which a detailed design surface is included in the model.
To take advantage of the elevation output in this tutorial, you will need to be running TUFLOW version 2012-05-AE or later. If the ZH (elevation) output type is used (as is the case with earlier versions), this will be written out at each output time as opposed to once for the entire simulation.
We have provided the location of the proposed development in the module data under Module_Data\Module_05\. Two GIS layers have been provided, Fill_Area and Cut_Area. For this impact assessment the fill area is to be removed from the floodplain for the 100 year Average Recurrence Interval (ARI) event being modelled. The cut area will have 1 m removed from the existing ground surface to increase the conveyance in the area. To make the GIS changes, please select your GIS package from the list below:
There are two methods for running multiple scenarios in TUFLOW:
- Create separate control files for each simulation; or
- Run multiple simulations from the same control file (.tcf) using simple logic statements
For this tutorial we will use the second approach. We will describe the logic as we progress through this tutorial, however, there is more detail in Section 4.12.2 of the 2010 TUFLOW manual. For this tutorial we will run two versions of the simulation, called EXG and DEV. The "EXG" case represents the existing or base case scenario. The models have the same boundaries, extents and DEM; the only difference will be that the "DEV" case has an additional layer to represent the proposed changes to the topography.
While we are going to run the model using the same control file, we still want individual results outputs and check files for each scenario. When we run the model, ~s1~ included in .tcf filename will be replaced with the scenario name (EXG or DEV).
For MapInfo users add the following
If Scenario == DEV
- Read GIS Z Shape == mi\2d_zsh_M05_development_001.mif
- Pause == Scenario not recognised
For other GIS users add the following
If Scenario == DEV
- Read GIS Z Shape == gis\2d_zsh_M05_development_001_R.shp
- Pause == Scenario not recognised
The elevation section of the geometry control file will now look something like this:
Set Zpts == 75
Read GRID Zpts == grid\DEM_M01.asc
Read GIS Z Shape == mi\2d_zsh_M03_Rd_Crest_001.mif
Read GIS Z HX Line == mi\2d_bc_M04_hx_001.mif | mi\2d_bc_M04_hx_001_P.mif
If Scenario == DEV
- Read GIS Z Shape == gis\2d_zsh_M05_development_001.mif
- Pause == Scenario not recognised
If the model is simulated with the DEV scenario, the additional 2d_zsh GIS layer will be read into the model. Otherwise, it will be omitted for the EXG scenario. In the above set of commands, if a scenario other than "DEV" or "EXG" is specified, the Pause command will pause the simulation and output the message Scenario not recognised. Although an optional step, including these commands will ensure that any errors in defining scenarios is captured. This is highly recommended when using nested If Scenario command blocks.
As with all TUFLOW inputs the scenarios are case insensitive (or should that be cASe InSeNsItIvE)!
Next we will update a number of other options in the TUFLOW Control File .tcf.
- All results are stored in the same file
- Results are much faster to access because indexing information is included; and
- Use of a folder structure makes viewing results easier
Map Output Format == xmdf ASC
ASC Map Output Interval == 0
The second command applies to ASC map output file formats only as denoted by the 'ASC' prefix. Specifying an interval of 0 will output only the maximums.
Map Output Data Types == h V q d MB1 MB2 ZH
Run the models
We are now ready to run the model. Now that we have two scenarios, we can run them from a batch file, specifying the scenario switch, or from a text editor / windows explorer with the scenario included in the .tcf file. Both methods are described below.
From a batch file
Batch files are a powerful way of running multiple simulations. There is a page on the wiki dedicated to batch files. To run the simulations in series (one after the other), the batch file should look like the below:
Start "TUFLOW" /wait "C:\TUFLOW\Releases\2012-05\w64\TUFLOW_iSP_w64.exe" -b -s1 EXG M05_5m_~s1~_001.tcf Start "TUFLOW" /wait "C:\TUFLOW\Releases\2012-05\w64\TUFLOW_iSP_w64.exe" -b -s1 DEV M05_5m_~s1~_001.tcf
From left to right the inputs in the batch file do the following:
- Start "TUFLOW" opens a new window and gives it the title "TUFLOW".
- The /wait argument prevents the second line from running until the first simulation has finished. (To run the second line at the same time simply remove the /wait argument.)
- The next argument is the TUFLOW executable to be used. In our case this is C:\TUFLOW\Releases\2012-05\w64\TUFLOW_iSP_w64.exe. This filepath will need to be updated.
- The -b switch tells TUFLOW to run in "batch" mode. This suppresses the message boxes at the end of the simulation and allows subsequent runs to start without user intervention.
- We need to input the scenario that we wish to simulate. The -s1 switch indicates a scenario argument. -s1 DEV indicates that the "DEV" scenario should be simulated. -s1 EXG indicates that the "EXG" scenario should be simulated. Up to 9 scenarios can be used in a single simulation.
- The final argument M05_5m_~s1~_001.tcf is the TUFLOW control file to be used for the simulation.
From Text Editor / Windows
In the batch file method above we specified the scenario to run. When we are creating a model or have only a small number of scenarios to run it is useful to be be able to start a simulation without using a batch file. To do this we need to include the scenario in the .tcf. To do so add the following line to the M05_5m_~s1~_001.tcf:
Model Scenarios == EXG
Save the control file and start the simulation from your text editor. Update the scenario to:
Model Scenarios == DEV
Resave the file and start the simulation. There should now be two simulations running. This method provides a quick way to start the simulations from your text editor, but when you have multiple scenarios and events batch files can streamline the process!
Once the simulations have finished processing we can compare the differences between the results for the "developed" and "existing" cases.
There are a number of different methods for comparing results datasets. In this tutorial we will outline this process in three ways:
- Using SMS
- Using GIS to subtract the grids
- Using the TUFLOW utilities
Care needs to be taken as not all methods will compare results when values exist in only one of the datasets. For example, in GIS a typical grid (raster) calculation will only output a difference if both grids have values. This can exclude areas from the difference grid where there has been an increase or decrease in flood extent.
SMS allows the user to easily compare datasets for the temporal and maximums data. To step through the process please see the following page.
The TUFLOW to GIS utility can be used to export the results into GIS format as a raster dataset. Within your GIS package a grid calculation is performed to calculate the difference in peak water level. The first step in this process is to use the TUFLOW to GIS utility to create the gridded results files. Use of the utility was described in tutorial module 1. There is also a separate page discussing the TUFLOW to GIS utility.
- Create a new text file in the TUFLOW\results\M05\2d\ folder. Save this as convert_to_GIS.bat
- Open this file in a text editor and add the following text. Note that in the first line, you will need to update the filepath to TUFLOW_to_GIS.exe to match the location of the utility on your computer.
set t2g="C:\TUFLOW\utilities\TUFLOW_to_GIS.exe" %t2g% -asc -b -typeH -max -out M05_5m_exg_001_hMax M05_5m_exg_001.xmdf %t2g% -asc -b -typeH -max -out M05_5m_dev_001_hMax M05_5m_dev_001.xmdf
For this tutorial we will convert the peak water level results to raster grids using a batch file:
To perform the GIS calculation please select your GIS package below:
There are two methods for calculating the flood afflux using the TUFLOW utilities:
Both of these methods are outlined below.
This method compares the results at each of the model output points and requires an identical number of nodes in each results file. If the model extent or number of cells has changed between models, you will be unable to compare using this method.
- Create a new text file in the TUFLOW\Results\M05\2d\ folder and save as afflux_res_to_res.bat
- Open the file in a text editor and add the following lines:
"C:\TUFLOW\utilities\res_to_res.exe" -dif -typeH -t99999 -out dH_M05_5m_001.xmdf M05_5m_dev_001.xmdf M05_5m_exg_001.xmdf
- Double click on the batch file in Windows explorer to begin the utility.
The above line will calculate the flood afflux and output this into a new file named dH_M05_5m_001.xmdf. The -t99999 flag specifies that the afflux is to be calculated at the maximums. If this were omitted the afflux would be calculated at each output timestep.
This file can now be loaded into SMS or alternatively can be converted to GIS using the TUFLOW_to_GIS utility. The command for doing this is:
"C:\TUFLOW\utilities\TUFLOW_to_GIS_w64.exe" -asc -s1 -2dm M05_5m_Exg_001.2dm dH_M05_5m_001.xmdf
The results are stored in a folder "difference" in the .xmdf file with the name "diff_M05_5m_dev_001 - M05_5m_exg_001". When converting to .asc instead of using the -typeH argument -s1 is used instead to convert the 1st (and in this case only) scalar dataset. The -2dm M05_5m_Exg_001.2dm is used to specify the 2D mesh file to be associated with the .xmdf file. For the xmdf file output from TUFLOW the .2dm file has the same filename as the .xmdf file just a different extension. In this case a new .xmdf file has been created and so we need to specify the correct .2dm file.
The asc_to_asc utility performs a comparison of each of the raster cells in a pair of raster grids. For this reason, the output grid files must have the same extent and cell size. By default the utility will produce two output files, the first is the difference in locations where both models are wet. The second contains the areas where there has been a change in inundation extent (eg. was dry now wet). The first step in using the asc_to_asc utility is to convert the TUFLOW outputs into .asc format using TUFLOW_to_GIS.exe:
- Create a new text file in the TUFLOW\results\M05\2d\ folder. Save this as afflux_calc_asc_to_asc.bat
- Open this file in a text editor and add the following text. Note that in the first two lines you will need to update the filepaths to the TUFLOW_to_GIS and Asc_to_Asc utilities to match their locations on your computer.
set t2g="C:\TUFLOW\utilities\TUFLOW_to_GIS.exe" set asc2asc="C:\TUFLOW\utilities\asc_to_asc.exe" %t2g% -asc -b -typeH -max -out M05_5m_exg_001_hMax.asc M05_5m_exg_001.xmdf %t2g% -asc -b -typeH -max -out M05_5m_dev_001_hMax.asc M05_5m_dev_001.xmdf %asc2asc% -dif -out dhMax_M05_5m_001.asc M05_5m_dev_001_hMax.asc M05_5m_exg_001_hMax.asc
The first two lines set the locations of the utilities on your computer. Lines three and four convert the existing and developed case results to .asc format. The final line creates the difference grids using the asc_to_asc utility. The resulting grids can be opened / imported into your GIS package and contoured and a colouring style applied.
We have used the powerful topographic modification features in TUFLOW to run a 100 year ARI design flood for an existing and a proposed topography. We have compared the predicted peak water levels for the two simulations to create a flood afflux result grid. A variety of methods for calculating the flood afflux were provided.
If the TUFLOW simulation fails to start TUFLOW will output the error in a number of locations. Firstly, check the Console Window or the TUFLOW Log File (.tlf) (located in the TUFLOW\runs\log\ folder with the name M05_5m_001.tlf). This file can be opened in a text editor and the error is generally located at the end of the file. You can however search for "Error" if you can not see the error. In most cases there is also a spatial location for the error message (if the error reported in the log file is prefixed by XY:). To check the location of the geographic errors, open the M05_5m_001_messages.mif or M05_5m_001_messages_P.shp file in your GIS package.
This section contains links to some possible problems that may occur when progressing through the fifth tutorial module. If you experience an issue that is not detailed below, please add and describe the problem on the discussion page or email firstname.lastname@example.org.