Tutorial M09

From TUFLOW FV Wiki
Revision as of 17:27, 1 February 2023 by Michaelbarry (talk | contribs) (Refixes model download URL)
Jump to navigation Jump to search

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.

The tutorial is broken up into the following subsections:


Set up TUFLOW FV

This tutorial recommends the use of TUFLOW FV 2023.0 or later, available on the TUFLOW Downloads Page. Please download TUFLOW FV 2023.0 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\2023.0\). 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 2023.0 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. The folders and most of their contents are the same as those within Tutorial 05’s “Complete_Model” folder, with some additional folders. 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 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” (μg/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, 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” (μg/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, 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 the TUFLOW Viewer Plugin. To start this process:

1. Open QGIS
2. Make sure the TUFLOW Viewer Plugin is installed

T09 7.0.0.01.png

3. Within the TUFLOW Viewer Plugin, 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. This tutorial will initially use WQ_000_WQ.nc

T09 7.0.0.02.png

4. 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

5. 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:

1. Make sure the WQ_DISS_OXYGEN_MG_L field is selected as above
2. 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

3. Click again on the same down arrow to exit the drop down menu
4. Click once on the TUFLOW Viewer Vertical Profile button itself

T09 7.0.1.02.png

5. 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

6. Right click on the Plot window and select Customise

T09 7.0.1.05.png

7. 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

8. 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

9. Click on the Freeze Axis Limits button to keep the ranges selected throughout time

T09 7.0.1.08.png

10. 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

11. 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:

1. Make sure the WQ_DISS_OXYGEN_MG_L field is selected as above
2. 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

3. Click again on the same down arrow to exit the drop down menu
4. Click on the TUFLOW Viewer Curtain Plot button itself

T09 7.0.2.02.png

5. 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, as shown below. 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

6. Customise the plot as previously described, and click OK

T09 7.0.2.05.png

7. The colourbar scale, name and map can be adjusted by clicking on the Images, etc. tab

T09 7.0.2.06.png

8. 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

9. 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

10. 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:

1. 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
2. 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

3. Click on the QGIS Selection tool

T09 7.0.3.02.png

4. 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

5. 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

6. Click again on the same down arrow to exit the drop down menu
7. Select ‘Layer Selection’ from the dropdown menu under Map Output Plotting

T09 7.0.3.05.png

8. Click on the Curtain button

T09 7.0.3.07.png

9. A curtain plot should appear in the ‘Cross Section / Long Profile’ plot window as follows

T09 7.0.3.08.png

10. Make any changes to the plot axes or colour bar as required (see above)
11. 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

12. Overturning events will also occur during the simulation to restore dissolved oxygen concentrations at depth

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:

1. 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

2. 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

3. Make sure the shapefile "TUFLOWFV\model\gis\2d_zln_Curtain_000_L.shp" is loaded into the QGIS Map window as above
4. Click on the QGIS Selection tool

T09 7.1.03.png

5. 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

6. 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

7. Click again on the same down arrow to exit the drop down menu
8. Select ‘Layer Selection’ from the dropdown menu under Map Output Plotting

T09 7.1.06.png

9. Click on the Curtain button

T09 7.1.07.png

10. A curtain plot should appear in the ‘Cross Section / Long Profile’ plot window as follows

T09 7.1.08.png

11. Make any changes to the plot axes or colourbar as required (see above)
12. 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:

1. 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

2. 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

3. Make sure the shapefile "TUFLOWFV\model\gis\2d_zln_Curtain_000_L.shp" is loaded into the QGIS Map window as above
4. Click on the QGIS Selection tool

T09 7.2.03.png

5. 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

6. 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

7. Click again on the same down arrow to exit the drop down menu
8. Select ‘Layer Selection’ from the dropdown menu under Map Output Plotting

T09 7.2.06.png

9. Click on the Curtain button

T09 7.2.07.png

10. A curtain plot should appear in the ‘Cross Section / Long Profile’ plot window as follows

T09 7.2.08.png

11. Make any changes to the plot axes or colourbar as required (see above)
12. 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

13. 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 WQ Module simulation. To complete more tutorials or learn more tips and tricks (including how to take advantage of the WQ Module's diagnostic variables in your model calibration process), please return to the TUFLOW FV Wiki Mainpage.