Tutorial M09: Difference between revisions

From TUFLOW FV Wiki
Jump to navigation Jump to search
Line 201: Line 201:


|-
|-
|   silicate min max == 0.0, 100.0
|     silicate min max == 0.0, 100.0
| Set min and max values for silicate
| Set min and max values for silicate


|-
|-
|   oxygen benthic == 4.2, 1.01
|     oxygen benthic == 4.2, 1.01
| Set silicate sediment release properties
| Set silicate sediment release properties


|-
|-
|   oxygen == on
|     oxygen == on
| Include oxygen in the simulation of silicate dynamics
| Include oxygen in the simulation of silicate dynamics


Line 229: Line 229:


|-
|-
|   ammonium min max == 0.0, 50.0
|     ammonium min max == 0.0, 50.0
| rowspan="2" | Set min and max values for ammonium and nitrate
| rowspan="2" | Set min and max values for ammonium and nitrate
|-
|-
|   ammonium min max == 0.0, 50.0
|     ammonium min max == 0.0, 50.0


|-
|-
|   ammonium benthic == 4.05, 1.06
|     ammonium benthic == 4.05, 1.06
| rowspan="2" | Set ammonium and nitrate sediment release properties
| rowspan="2" | Set ammonium and nitrate sediment release properties
|-
|-
|   nitrate benthic == 4.25, 1.10
|     nitrate benthic == 4.25, 1.10


|-
|-
|   nitrification == 0.05, 4.15, 1.01
|     nitrification == 0.05, 4.15, 1.01
| Set nitrification parameters
| Set nitrification parameters


|-
|-
|   denitrification == michaelis menten, 0.05, 4.01, 1.03
|     denitrification == michaelis menten, 0.05, 4.01, 1.03
| Set denitrification model and properties
| Set denitrification model and properties


|-
|-
|     oxygen == on
|     oxygen == on
| Include oxygen in the simulation of inorganic nitrogen dynamics
| Include oxygen in the simulation of inorganic nitrogen dynamics


|-
|-
|   atmospheric deposition == 5.0, 0.0, 0.5
|     atmospheric deposition == 5.0, 0.0, 0.5
| Set atmospheric deposition of inorganic nitrogen parameters
| Set atmospheric deposition of inorganic nitrogen parameters


|-
|-
|   anaerobic oxidation of ammonium == 0.005, 2.5, 2.5
|     anaerobic oxidation of ammonium == 0.005, 2.5, 2.5
| Set anammox parameters
| Set anammox parameters


|-
|-
|     diss nitrate reduction to ammonium == 0.005, 5.0
|     diss nitrate reduction to ammonium == 0.005, 5.0
| Set DRNA parameters
| Set DRNA parameters


Line 281: Line 281:


|-
|-
|   FRP min max == 0.0, 5.0
|     FRP min max == 0.0, 5.0
| rowspan="2" | Set min and max values for inorganic phosphorus and adsorbed inorganic phosphorus
| rowspan="2" | Set min and max values for inorganic phosphorus and adsorbed inorganic phosphorus
|-
|-
|   FRP ads min max == 0.0, 5.0
|     FRP ads min max == 0.0, 5.0


|-
|-
|   FRP benthic == 4.3, 1.095
|     FRP benthic == 4.3, 1.095
| Set inorganic phosphorus sediment release properties
| Set inorganic phosphorus sediment release properties


|-
|-
|   oxygen == on
|     oxygen == on
| Include oxygen in the simulation of inorganic phosphorus dynamics
| Include oxygen in the simulation of inorganic phosphorus dynamics


|-
|-
|   atmospheric deposition == 0.5, 0.0
|     atmospheric deposition == 0.5, 0.0
| Set atmospheric deposition of inorganic phosphorus parameters
| Set atmospheric deposition of inorganic phosphorus parameters


|-
|-
|   adsorption == linear, 0.005
|     adsorption == linear, 0.005
| Set adsorption model and parameter
| Set adsorption model and parameter


|-
|-
|   settling == 0.000
|     settling == 0.000
| Set settling rate of adsorbed inorganic phosphorus
| Set settling rate of adsorbed inorganic phosphorus


Line 323: Line 323:


|-
|-
|   min max == 0.05, 50.0
|     min max == 0.05, 50.0
| Set min and max values the phytoplankton group
| Set min and max values the phytoplankton group


|-
|-
|   temperature limitation == none
|     temperature limitation == none
| Set the temperature limitation model to none. Users may want to trial other models
| Set the temperature limitation model to none. Users may want to trial other models


|-
|-
|   salinity limitation == none
|     salinity limitation == none
| Set the salinity limitation model to be off
| Set the salinity limitation model to be off


|-
|-
|   light limitation == basic, 0.0, 20.0
|     light limitation == basic, 0.0, 20.0
| Set the light limitation model and parameters
| Set the light limitation model and parameters


|-
|-
|   nitrogen limitation == 0.01, 2.5
|     nitrogen limitation == 0.01, 2.5
| Set the nitrogen limitation model and parameters
| Set the nitrogen limitation model and parameters


|-
|-
|   phosphorus limitation == 0.005, 0.01
|     phosphorus limitation == 0.005, 0.01
| Set the phosphorus limitation model and parameters
| Set the phosphorus limitation model and parameters


|-
|-
|   uptake == 3.0, 0.6, 8.0
|     uptake == 3.0, 0.6, 8.0
| Set nutrient uptake parameters
| Set nutrient uptake parameters


|-
|-
|   primary productivity == 2.4, 1.2  
|     primary productivity == 2.4, 1.2  
| Set primary productivity parameters
| Set primary productivity parameters


|-
|-
|   respiration == 0.01, 1.0, 0.3, 0.4, 0.5  
|     respiration == 0.01, 1.0, 0.3, 0.4, 0.5  
| Set respiration parameters
| Set respiration parameters


|-
|-
|   carbon chla ratio == 27.8
|     carbon chla ratio == 27.8
| Set the carbon to chlorophyll a ratio
| Set the carbon to chlorophyll a ratio


|-
|-
|   nitrogen fixing == 0.001, 0.5
|     nitrogen fixing == 0.001, 0.5
| Turn on nitrogen fixing and set its parameters
| Turn on nitrogen fixing and set its parameters


|-
|-
|   settling == constant, -0.05
|     settling == constant, -0.05
| Set a constant settling rate model and its parameter
| Set a constant settling rate model and its parameter


Line 386: Line 386:
| rowspan="7" | Set default material properties for all constituents that are allowed sediment fluxes
| rowspan="7" | Set default material properties for all constituents that are allowed sediment fluxes
|-
|-
|   oxygen flux == -10.0  
|     oxygen flux == -10.0  
|-
|-
|   silicate flux == 10.0
|     silicate flux == 10.0
|-
|-
|   ammonium flux == 20.0
|     ammonium flux == 20.0
|-
|-
|   nitrate flux == 30.0  
|     nitrate flux == 30.0  
|-
|-
|   frp flux == 5.0  
|     frp flux == 5.0  
|-
|-
| end material
| end material
Line 402: Line 402:
| rowspan="3" | Set sediment oxygen properties for material 1 (unchanged from model 000)
| rowspan="3" | Set sediment oxygen properties for material 1 (unchanged from model 000)
|-
|-
|   oxygen flux == -550.0  
|     oxygen flux == -550.0  
|-
|-
| end material
| end material
Line 410: Line 410:
| rowspan="3" | Set sediment oxygen properties for material 2 (unchanged from model 000)
| rowspan="3" | Set sediment oxygen properties for material 2 (unchanged from model 000)
|-
|-
|   oxygen flux == 0.0  
|     oxygen flux == 0.0  
|-
|-
| end material
| end material
Line 418: Line 418:
| rowspan="3" | Set sediment oxygen properties for material 4 (unchanged from model 000)
| rowspan="3" | Set sediment oxygen properties for material 4 (unchanged from model 000)
|-
|-
|   oxygen flux == -3500.0  
|     oxygen flux == -3500.0  
|-
|-
| end material
| end material

Revision as of 15:08, 1 December 2022

Tutorial Description

Tutorial Module 09 demonstrates the development of a 3D water quality model using the TUFLOW FV Water Quality (WQ) Module. Starting with the final model of Tutorial 05, this tutorial will progressively increase in complexity, incorporating the simulation of dissolved oxygen, then inorganic nutrients, and finally capturing full organic nutrient cycling. Tutorial Module 09 provides a full set of working model files that can be re-run and experimented with, in a manner similar to the mass conservation model.

The tutorial is broken up into the following subsections:


Set up TUFLOW FV

This tutorial recommends the use of TUFLOW FV 2022.0.0 beta or later, available on the TUFLOW Downloads Page. Please download TUFLOW FV 2022.0.0 beta or later, and save to C:\TUFLOWFV\ or another preferred location. This location will be used in building scripts to run Tutorial 09's models.

This demonstration model dynamically couples TUFLOW FV with GOTM to model vertical turbulence exchanges. To do so, please download the TUFLOW FV External Turbulence API. Once downloaded, unzip the file .\GOTM\tuflowfv_external_turb.dll and replace the default tuflowfv_external_turb.dll that resides in the same directory as TUFLOWFV.exe (for example, C:\TUFLOWFV\2022.0.0_beta\). If this step has already been undertaken as part of the running of Tutorial 05, then it need not be repeated here.

TUFLOW FV WQ uses the dynamically linked library tuflowfv_wq.dll, which ships with TUFLOW FV 2022.0.0 beta or later, so no additional user set up is required to access the compiled binary necessary for water quality simulation.

Tutorial Data

The data for Tutorial 09 is available from the TUFLOW Website Tutorial Module Download. The folders and most of their contents are the same as those within Tutorial 05’s “Complete_Model” folder, with the addition of one new folder (named “wqm”) that houses TUFLOW FV WQ control files. The initial condition and boundary condition files in Tutorial 09 are slightly altered and expanded (to account for water quality specification) versions of those provided in Tutorial 05.

It is recommended that Tutorial 05 be completed before this water quality tutorial is attempted.

Complete Model Folder Contents

Three water quality models are provided in the Complete Model folder under the TUFLOWFV folder, and include:

  • .2dm model mesh file ‘hydraul_006.2dm’ (model\geo folder, unchanged from Tutorial 05)
  • 3D layer .csv file (model\csv folder, unchanged from Tutorial 05)
  • Supporting GIS files for nodestring and material definition (model\gis folder)
  • GIS empty files for adding features to the model (model\gis\empty folder, unchanged from Tutorial 05)
  • Boundary condition data, in .csv and .nc format (bc_dbase and bc_dbase\met\ folders, with the latter unchanged from Tutorial 05)
  • Meteorology FVC file (bc_dbase\met\ folder, unchanged from Tutorial 05)
  • Sediment file FVSED for modelling sediment transport (stm folder, unchanged from Tutorial 05)
  • Water quality control files (wqm folder)
  • Results are not provided as part of this tutorial data package. The simulations need to run to produce results.

Upstream And Downstream Boundary Condition Files

The upstream boundary condition file from Tutorial 05 is the same in this tutorial but has added columns (one for each water quality constituent simulated). The added columns have headers that are not keywords in themselves, but are used by TUFLOW FV to determine water quality boundary conditions for each constituent based on ordered listing of these same headers in “bc” blocks in each TUFLOW FV control file (see subsequent sections for a description of this ordering process). This upstream boundary’s location is set in the model mesh via association with a nodestring.

The downstream nodestring boundary condition of Tutorial 09’s model mesh has been removed so that it presents a barrier (such as a dam wall) to flow rather than an outflow. A new boundary condition has been created in its place that that is a polygon extraction that has outflows that approximately balance those entering the model from its upstream end. No salinity, temperature, sediment or water quality specifications (columns) are necessary for this extraction boundary condition because it represents outflows that have no influence on the model other than volumetric removal.

For the second and third simulations, a further inflow boundary condition is added (also as an inflow polygon) to represent a high nutrient wastewater treatment plant discharge.

Model Mesh File

This is the same model mesh as that used in Tutorial 05, in 3D (TUFLOW FV control file HYD_002.fvc).

Meteorological NetCDF Files

The meteorological data for the model is the same as that used in Tutorial 05, in 3D (TUFLOW FV control file HYD_002.fvc).

TUFLOW FV Hydrodynamic Model Setups

The TUFLOW FV control files for the models covered in this tutorial are located in the Complete_Model\TUFLOWFV\runs (TUFLOW FV control files) and Complete_Model\TUFLOWFV\wqm (TUFLOW FV WQ Module control files) directories. Descriptions are listed below.

WQ_000.fvc The same as HYD_002.fvc but with additional commands to allow for simulation of dissolved oxygen.
Companion water quality control file is WQ_000.fvwq (see following section)
WQ_001.fvc The same as HYD_002.fvc but with additional commands to allow for simulation of dissolved oxygen, inorganic nutrients and phytoplankton.
Companion water quality control file is WQ_001.fvwq (see following section)
WQ_002.fvc The same as HYD_002.fvc but with additional commands to allow for simulation of dissolved oxygen, inorganic nutrients, organic nutrients and phytoplankton.
Companion water quality control file is WQ_002.fvwq (see following section)


TUFLOW FV WQ Module Setups

The commands and syntax of the fvwq control files are provided below. Copy and paste commands into the search facility of the TUFLOW FV Water Quality Module manual for detailed descriptions of each.

Model 000: Dissolved oxygen

Dissolved oxygen simulation is often a good place to start when constructing water quality models. The model WQ_000.fvc provides this simulation capability, in 3D, and includes density effects, atmospheric heat exchange and meteorological forcing.

! SIMULATION CONTROLS
simulation class == do Set the overall simulation class to be "DO"
wq dt == 1800.0 Set the water quality module timestep. This is the frequency at which the WQ module is called from TUFLOW FV
wq units == mgl Set the units of simulation of the WQ module.
! CONSTITUENT MODEL SETTINGS
! Oxygen
oxygen model == O2 Set oxygen constituent model to be O2
    oxygen benthic == 0.5, 1.08 Set sediment oxygen consumption properties
end oxygen model End oxygen constituent model
! MATERIAL SPECIFICATIONS
material == default Set default material properties
    oxygen flux == -10.0
end material
material == 1 Set sediment oxygen properties for material 1
    oxygen flux == -550.0
end material
material == 2 ! Mudflats Set sediment oxygen properties for material 2
    oxygen flux == 0.0
end material
material == 4 ! Deep channel Set sediment oxygen properties for material 4
    oxygen flux == -3500.0
end material

The fvwq file has a set up that:

  • Deploys the “DO” Simulation Class, so therefore includes the water quality computed variables:
    • Dissolved oxygen (mg/L)
  • Uses the milligrams per litre units system, with TUFLOW FV calling the WQ Module every 30 minutes
  • Has sediment oxygen demand that varies spatially across various material types. These demands are greatest in the drowned river thalweg (material 4), and zero in the shallows

Model 001: Inorganics

Simulation of inorganics is often the next logical step after simulating dissolved oxygen. The model WQ_001.fvc provides this simulation capability, in 3D, and includes density effects, atmospheric heat exchange and meteorological forcing. The constituent model block for oxygen is unchanged from model 000. One phytoplankton group is simulated.

! SIMULATION CONTROLS
simulation class == inorganics Set the overall simulation class to be “inorganics”
wq dt == 1800.0 Set the water quality module timestep. This is the frequency at which the WQ module is called from TUFLOW FV
wq units == mgl Set the units of simulation of the WQ module.
! CONSTITUENT MODEL SETTINGS
! Silicate
silicate model == Si Set silicate constituent model to be Si
    silicate min max == 0.0, 100.0 Set min and max values for silicate
    oxygen benthic == 4.2, 1.01 Set silicate sediment release properties
    oxygen == on Include oxygen in the simulation of silicate dynamics
end silicate model End silicate constituent model
! Inorganic Nitrogen
Inorganic nitrogen model == ammoniumnitrate Set inorganic nitrogen constituent model to be ammoniumnitrate
    ammonium min max == 0.0, 50.0 Set min and max values for ammonium and nitrate
    ammonium min max == 0.0, 50.0
    ammonium benthic == 4.05, 1.06 Set ammonium and nitrate sediment release properties
    nitrate benthic == 4.25, 1.10
    nitrification == 0.05, 4.15, 1.01 Set nitrification parameters
    denitrification == michaelis menten, 0.05, 4.01, 1.03 Set denitrification model and properties
    oxygen == on Include oxygen in the simulation of inorganic nitrogen dynamics
    atmospheric deposition == 5.0, 0.0, 0.5 Set atmospheric deposition of inorganic nitrogen parameters
    anaerobic oxidation of ammonium == 0.005, 2.5, 2.5 Set anammox parameters
    diss nitrate reduction to ammonium == 0.005, 5.0 Set DRNA parameters
end inorganic nitrogen model End inorganic nitrogen constituent model
! Inorganic Phosphorus
Inorganic phosphorus model == frphsads Set inorganic phosphorus constituent model to be frphads, which includes simulation of free and adsorbed inorganic phosphorus
    FRP min max == 0.0, 5.0 Set min and max values for inorganic phosphorus and adsorbed inorganic phosphorus
    FRP ads min max == 0.0, 5.0
    FRP benthic == 4.3, 1.095 Set inorganic phosphorus sediment release properties
    oxygen == on Include oxygen in the simulation of inorganic phosphorus dynamics
    atmospheric deposition == 0.5, 0.0 Set atmospheric deposition of inorganic phosphorus parameters
    adsorption == linear, 0.005 Set adsorption model and parameter
    settling == 0.000 Set settling rate of adsorbed inorganic phosphorus
end inorganic phosphorus model End inorganic phosphorus constituent model
! Phytoplankton
phyto model == basic, green Set phytoplankton constituent model to be the basic model (i.e. assumed constant internal carbon to nutrient ratios), and named as “green”
    min max == 0.05, 50.0 Set min and max values the phytoplankton group
    temperature limitation == none Set the temperature limitation model to none. Users may want to trial other models
    salinity limitation == none Set the salinity limitation model to be off
    light limitation == basic, 0.0, 20.0 Set the light limitation model and parameters
    nitrogen limitation == 0.01, 2.5 Set the nitrogen limitation model and parameters
    phosphorus limitation == 0.005, 0.01 Set the phosphorus limitation model and parameters
    uptake == 3.0, 0.6, 8.0 Set nutrient uptake parameters
    primary productivity == 2.4, 1.2 Set primary productivity parameters
    respiration == 0.01, 1.0, 0.3, 0.4, 0.5 Set respiration parameters
    carbon chla ratio == 27.8 Set the carbon to chlorophyll a ratio
    nitrogen fixing == 0.001, 0.5 Turn on nitrogen fixing and set its parameters
    settling == constant, -0.05 Set a constant settling rate model and its parameter
end phyto model End phytoplankton constituent model
! MATERIAL SPECIFICATIONS
material == default Set default material properties for all constituents that are allowed sediment fluxes
    oxygen flux == -10.0
    silicate flux == 10.0
    ammonium flux == 20.0
    nitrate flux == 30.0
    frp flux == 5.0
end material
material == 1 Set sediment oxygen properties for material 1 (unchanged from model 000)
    oxygen flux == -550.0
end material
material == 2 ! Mudflats Set sediment oxygen properties for material 2 (unchanged from model 000)
    oxygen flux == 0.0
end material
material == 4 ! Deep channel Set sediment oxygen properties for material 4 (unchanged from model 000)
    oxygen flux == -3500.0
end material

The fvwq file has a set up that:

  • Deploys the “inorganics” Simulation Class, so therefore includes the water quality computed variables:
    • Dissolved oxygen (mg/L)
    • Silicate (mg/L)
    • Ammonium (mg/L)
    • Nitrate (mg/L)
    • Free reactive phosphorus (mg/L)
    • Adsorbed free reactive phosphorus (mg/L)
    • One phytoplankton group called “green”
  • Uses the milligrams per litre units system, with TUFLOW FV calling the WQ Module every 30 minutes
  • Has sediment oxygen demand that varies spatially across various material types, but constant behaviour for all other sediment releases.

Model 002: Organics

Simulation of organics is often the next logical step after simulating inorganics. The model WQ_002.fvc provides this simulation capability, in 3D, and includes density effects, atmospheric heat exchange and meteorological forcing. Constituent model blocks for oxygen, silicate, inorganic nitrogen, inorganic phosphorus and phytoplankton are unchanged from model 001.

! SIMULATION CONTROLS
simulation class == organics Set the overall simulation class to be “organics”
wq dt == 1800.0 Set the water quality module timestep. This is the frequency at which the WQ module is called from TUFLOW FV
wq units == mgl Set the units of simulation of the WQ module.
! CONSTITUENT MODEL SETTINGS
! Organics
organic matter model == refractory Set organics constituent model to be both labile and refractory
carbon min max == 0.0, 10.0, 0.0, 10.0 Set min and max values for labile particulate and labile dissolved organic carbon, nitrogen and phosphorus
nitrogen min max == 0.0, 10.0, 0.0, 10.0
phosphorus min max == 0.0, 10.0, 0.0, 10.0
organics benthic == 1.0, 1.08 Set labile organic matter hydrolysis parameters
mineralisation == 0.01, 1.2, 1.03, 0.5, 0.5 Set labile organic matter mineralisation parameters
self shading == 0.0, 0.0 Set labile organic matter self shading parameters
settling == stokes, 0.00000035, 2450.8 Set labile particulate settling model and parameters
ref carbon min max == 0.0, 10.0, 0.0, 10.0 Set min and max values for refractory particulate and refractory dissolved organic carbon, nitrogen and phosphorus
ref nitrogen min max == 0.0, 10.0, 0.0, 10.0
ref phosphorus min max == 0.0, 9.0, 0.0, 9.0
ref breakdown == 0.001, 0.1, 0.01 Set refractory organic matter breakdown parameters
ref activation == 0.01 Set refractory organic matter activation parameters
ref photolysis == 0.5 Set refractory organic matter photolysis parameter
ref self shading == 0.0, 0.0 Set refractory organic matter self shading parameters
ref settling == stokes, 0.00000045, 2450.8 Set refractory particulate settling model and parameters
end organic matter model End organics constituent model
! MATERIAL SPECIFICATIONS
material == default Set default material properties for all constituents that are allowed sediment fluxes
oxygen flux == -10.0
silicate flux == 10.0
ammonium flux == 20.0
nitrate flux == 30.0
frp flux == 5.0
doc flux == 20.0
don flux == 10.0
dop flux == 5.0
end material
material == 1 Set sediment oxygen properties for material 1 (unchanged from model 000)
oxygen flux == -550.0
end material
material == 2 ! Mudflats Set sediment oxygen properties for material 2 (unchanged from model 000)
oxygen flux == 0.0
end material
material == 4 ! Deep channel Set sediment oxygen properties for material 4 (unchanged from model 000)
oxygen flux == -3500.0
end material

The fvwq file has a set up that:

  • Deploys the “organics” Simulation Class, so therefore includes the water quality computed variables:
    • Dissolved oxygen (mg/L)
    • Silicate (mg/L)
    • Ammonium (mg/L)
    • Nitrate (mg/L)
    • Free reactive phosphorus (mg/L)
    • Adsorbed free reactive phosphorus (mg/L)
    • Particulate and dissolved labile organic carbon (mg/L)
    • Particulate and dissolved labile organic nitrogen (mg/L)
    • Particulate and dissolved labile organic phosphorus (mg/L)
    • Particulate and dissolved refractory organic carbon (mg/L)
    • Particulate and dissolved refractory organic nitrogen (mg/L)
    • Particulate and dissolved refractory organic phosphorus (mg/L)
    • One phytoplankton group called “green”
  • Uses the milligrams per litre units system, with TUFLOW FV calling the WQ Module every 30 minutes
  • Has sediment oxygen demand that varies spatially across various material types, but constant behaviour for all other sediment releases.


Running the Models

Results have not been provided with the data download package and you will need to run each model to obtain them. To do this, you can refer to Running TUFLOW FV for a detailed description of the various methods for running a TUFLOW FV simulation. A template batch file is provided in the Complete_Model\TUFLOWFV\runs folder. Within this batch file, you will need to repath to the version of TUFLOW FV you have downloaded and updated as defined in Set up TUFLOW FV.

Results are written to three netcdf files, where XXX is either 000 (“DO”), 001 (“inorganics”) or 002 (“organics”):

  • WQ_XXX_HD.nc (Hydrodynamics and sediment)
  • WQ_XXX_WQ.nc (All water quality variables)
  • WQ_XXX_WQ_diag.nc (All water quality diagnostic variables)


Inspecting Results

Model 000: Dissolved Oxygen

All results from this simulation can be viewed and interrogated in QGIS using TUFLOW Viewer. To start this process:

T09 7.0.0.01.png
  • Load the water quality netcdf results file, by navigating to the results folder of your simulation and selecting the netcdf file you wish to view. Initially, this tutorial will use WQ_000_WQ.nc
T09 7.0.0.02.png
  • Once loaded, click once on the WQ_DISS_OXYGEN_MG_L field in the Result Type window, and a plan view of the results will appear in the main map window. This is automatically set to be the surface layer of the three dimensional dataset
T09 7.0.0.03.png
T09 7.0.0.04.png
  • Click the TUFLOW Viewer mesh rendering button, if desired T09 7.0.0.05.png. This will show the mesh, but is not required to be displayed
T09 7.0.0.06.png

Given that the model is a stationary waterway, the dimension of primary interest is the vertical. It therefore makes sense to examine the dissolved oxygen predictions using vertical interrogation techniques. Different ways of doing this are presented below.

Vertical line profiles

The three dimensional netcdf dataset can be interrogated at a single x/y coordinate location, and a vertical profile over the water column at that location then plotted. This output is a line, and the time at which this line is plotted can be adjusted using the TUFLOW Viewer time controls. To do this:

  • Make sure the WQ_DISS_OXYGEN_MG_L field is selected as above
  • Click the down arrow next to the TUFLOW Viewer Vertical Profile button and select the WQ_DISS_OXYGEN_MG_L field checkbox
T09 7.0.1.01.png
  • Click again on the same down arrow to exit the drop down menu
  • Click once on the TUFLOW Viewer Vertical Profile button itself
T09 7.0.1.02.png
  • Click once anywhere within the model domain – somewhere near the deeper northern end of the model will be more interesting than the shallow areas. A pink dot will appear on the clicked location, and a profile of dissolved oxygen concentration at that location will appear in the TUFLOW Viewer Plot window, within the Vertical Profile tab. The initial profile will be uniform with depth because the water quality model was initialised with a single dissolved oxygen concentration everywhere within the domain
T09 7.0.1.03.png
T09 7.0.1.04.png
  • Right click on the Plot window and select Customise
T09 7.0.1.05.png
  • Select the desired depth (y-axis) and dissolved oxygen (x-axis) ranges, and rename axes labels if desired, and click OK. Make sure the ranges include all the data to be displayed, otherwise the plot will disappear off the axes. The ranges shown below are satisfactory in this example
T09 7.0.1.06.png
  • The Plot window should appear similar to the below. The vertical dimension (depth) is shown in the y-axis and dissolved oxygen concentration on the x-axis. The green line is the dissolved oxygen concentration at all model depths at the X/Y location previously chosen, at the initial timestep of the model simulation.
T09 7.0.1.07.png
  • Click on the Freeze axis limits button to keep the ranges selected throughout time
T09 7.0.1.08.png
  • Use the TUFLOW Viewer time navigation buttons to see how the (green line) vertical profile of dissolved oxygen at the chosen X/Y location changes in time. Check the Show as dates box to have the model time presented as dates
T09 7.0.1.09.png
  • It is expected that dissolved oxygen will decrease at depth (due to sediment consumption), and increase at the surface (due to atmospheric reaeration), to show profiles similar to the below. Overturning events where the water column vertically mixes are also expected
T09 7.0.1.10.png

Vertical curtain profiles – manual location

The three dimensional netcdf dataset can be interrogated concurrently at multiple x/y coordinate locations, and a continuous vertical curtain profile over the water column at those locations plotted. This is essentially a vertical slice through the three dimensional dataset, at a series of X/Y points that form a line over the model domain when viewed in plan. This output is a colour contoured curtain, and the time at which this curtain is plotted can be adjusted using the TUFLOW Viewer time controls. To do this:

  • Make sure the WQ_DISS_OXYGEN_MG_L field is selected as above
  • Click the down arrow next to the TUFLOW Viewer Curtain Plot button, and select the WQ_DISS_OXYGEN_MG_L field checkbox
T09 7.0.2.01.png
  • Click again on the same down arrow to exit the drop down menu
  • Click on the TUFLOW Viewer Curtain Plot button itself
T09 7.0.2.02.png
  • Click a path along the desired curtain X/Y location within the model domain – somewhere along the deeper drowned river thalweg from the northern end of the model to the south will be more interesting than the shallow areas. Double click to end the curtain path definition. A green line will appear along the clicked locations, and a curtain profile of dissolved oxygen concentration along that path will appear in the TUFLOW Viewer Plot window, within the Cross Section / Long Profile tab. The left hand end of the curtain colour contour will be the first location clicked in the green path
T09 7.0.2.03.png
T09 7.0.2.04.png
  • Customise the plot as previously described, and click OK
T09 7.0.2.05.png
  • The colourbar scale, name and map can be adjusted by clicking on the Images, etc. tab
T09 7.0.2.06.png
  • The curtain plot should then look similar to the following (it will have different data depending on the time selected)
T09 7.0.2.07.png
  • Use the TUFLOW Viewer time navigation buttons to see how the vertical curtain of dissolved oxygen along the chosen X/Y path changes in time. Check the Show as dates box to have the model time presented as dates
T09 7.0.2.08.png
  • It is expected that dissolved oxygen will decrease at depth (due to sediment consumption), and increase at the surface (due to atmospheric reaeration). Overturning events where the water column vertically mixes are also expected

Vertical curtain profiles – GIS location

Whilst the three dimensional netcdf dataset can be interrogated along a manually defined curtain as above, it is sometimes more convenient to use a predefined and repeatable curtain transect path. This can be achieved by using a shapefile to set the transect location rather than manual digitising. To do this:

  • Make sure the WQ_000_WQ.nc dataset is loaded into the TUFLOW Viewer as above, and the WQ_DISS_OXYGEN_MG_L field is selected
  • Drag the shapefile "TUFLOWFV\model\gis\2d_zln_Curtain_000_L.shp" onto the QGIS map window. The location of the transect will appear on top of the WQ_DISS_OXYGEN_MGL contours as a line
T09 7.0.3.01.png
  • Click on the QGIS selection tool
T09 7.0.3.02.png
  • Select the line from the loaded shapefile "TUFLOWFV\model\gis\2d_zln_Curtain_000_L.shp". The line will change colour
T09 7.0.3.03.png
  • Click the down arrow next to the TUFLOW Viewer Curtain Plot button, and select the WQ_DISS_OXYGEN_MG_L field checkbox
T09 7.0.3.04.png
  • Click again on the same down arrow to exit the drop down menu
  • Select ‘Layer Selection’ from the dropdown menu under Map Output Plotting
T09 7.0.3.05.png
  • Click on the Curtain button
T09 7.0.3.07.png
  • A curtain plot should appear in the ‘Cross Section / Long Profile’ plot window as follows
T09 7.0.3.08.png
  • Make any changes to the plot axes or colour bar as required (see above)
  • As above, use the TUFLOW Viewer time navigation buttons to see how the vertical curtain of dissolved oxygen along the chosen X/Y path changes in time. Check the Show as dates box (as above) to have the model time presented as dates. As the timeseries is clicked through, dissolved oxygen will decrease at depth, as expected
T09 7.0.3.09.png
  • Overturning events will also occur during the simulation to return the lower depth dissolved oxygen concentrations in part
T09 7.0.3.10.png

Model 001: Inorganics

The same curtain interrogation techniques can be used to examine the inorganics (001) simulation results. Using the same interrogation techniques described above, the inorganic phosphorus concentrations can be examined along the same curtain, especially in the location of the wastewater treatment plant. To do this:

  • Load the inorganics water quality netcdf results file, by navigating to the results folder of your simulation and selecting the netcdf file you wish to view. In this case, WQ_001_WQ.nc will be loaded
T09 7.1.01.png
  • Once loaded, click once on the WQ_FRP_MG_L field in the Result Type window, and a plan view of the results will appear in the main map window. This is automatically set to be the surface layer of the three dimensional dataset
T09 7.1.02.png
  • Make sure the shapefile "TUFLOWFV\model\gis\2d_zln_Curtain_000_L.shp" is loaded into the QGIS map window as above
  • Click on the QGIS selection tool
T09 7.1.03.png
  • Select the line from the loaded shapefile "TUFLOWFV\model\gis\2d_zln_Curtain_000_L.shp". The line will change colour
T09 7.1.04.png
  • Click the down arrow next to the TUFLOW Viewer Curtain Plot button, and select the WQ_FRP_MG_L field checkbox
T09 7.1.05.png
  • Click again on the same down arrow to exit the drop down menu
  • Select ‘Layer Selection’ from the dropdown menu under Map Output Plotting
T09 7.1.06.png
  • Click on the Curtain button
T09 7.1.07.png
  • A curtain plot should appear in the ‘Cross Section / Long Profile’ plot window as follows
T09 7.1.08.png
  • Make any changes to the plot axes or colourbar as required (see above)
  • As above, use the TUFLOW Viewer time navigation buttons to see how the vertical curtain of FRP along the chosen X/Y path changes in time. Check the Show as dates box (as above) to have the model time presented as dates. As the timeseries is clicked through, FRP will increase at the location of the wastewater treatment plant discharge, as expected
T09 7.1.09.png

Model 002: Organics

The same curtain interrogation techniques can be used to examine the organics (002) simulation results. Using the same interrogation techniques described above, the phytoplankton concentrations can be examined along the same curtain, especially in the location of the wastewater treatment plant. To do this:

  • Load the organics water quality netcdf results file, by navigating to the results folder of your simulation and selecting the netcdf file you wish to view. In this case, WQ_002_WQ.nc will be loaded
T09 7.2.01.png
  • Once loaded, click once on the WQ_PHYTO_GREEN_CONC_MICG_L field in the Result Type window, and a plan view of the results will appear in the main map window. This is automatically set to be the surface layer of the three dimensional dataset
T09 7.2.02.png
  • Make sure the shapefile "TUFLOWFV\model\gis\2d_zln_Curtain_000_L.shp" is loaded into the QGIS map window as above
  • Click on the QGIS selection tool
T09 7.2.03.png
  • Select the line from the loaded shapefile "TUFLOWFV\model\gis\2d_zln_Curtain_000_L.shp". The line will change colour
T09 7.2.04.png
  • Click the down arrow next to the TUFLOW Viewer Curtain Plot button, and select the WQ_PHYTO_GREEN_CONC_MICG_L field checkbox.
T09 7.2.05.png
  • Click again on the same down arrow to exit the drop down menu
  • Select ‘Layer Selection’ from the dropdown menu under Map Output Plotting
T09 7.2.06.png
  • Click on the Curtain button
T09 7.2.07.png
  • A curtain plot should appear in the ‘Cross Section / Long Profile’ plot window as follows
T09 7.2.08.png
  • Make any changes to the plot axes or colourbar as required (see above)
  • As above, use the TUFLOW Viewer time navigation buttons to see how the vertical curtain of phytoplankton concentration along the chosen X/Y path changes in time. Check the Show as dates box (as above) to have the model time presented as dates. As the timeseries is clicked through, phytoplankton grows and concentrations increase at the location of the wastewater treatment plant discharge (and elsewhere through advection), as expected
T09 7.0.3.06.png
  • Try plotting other TUFLOW FV WQ Module variables to see how they evolve in time


Conclusion

This tutorial has covered the basics of setting up, executing and interrogating a TUFLOW FV water quality module simulation. To complete more tutorials or learn more tips and tricks, please return to the TUFLOW FV Wiki Mainpage.