WQ 000
Jump to navigation
Jump to search
The syntax additions to HYD_002.fvc to include water quality modelling are detailed below .
! TUFLOW FV TUTORIAL | The first lines are a description of the model simulation. |
!Estuary Model – 3D Water Quality Model | |
GIS FORMAT == SHP | No syntax changes. |
SHP Projection == ..\model\gis\projection.prj | |
Tutorial Model == ON | |
! GENERAL PARAMETERS | Sediment has been included in the water quality model. To run the sediment model, a sediment control file is provided in the ST folder. |
Bottom Drag Model == ks | |
Include Salinity == 1,1 | |
Include Temperature == 1,1 | |
Spatial Order == 1,2 | |
Include Heat == 1 | |
Include sediment == 1,0 | |
! TIME COMMANDS | No syntax changes. |
time format == ISODATE | |
Start time == 01/05/2011 00:00:00 | |
End Time == 07/05/2011 00:00:00 | |
CFL == 0.6 | |
Timestep Limits == 0.0001,100 | |
Display dt == 900 | |
! MODEL PARAMETERS | An external water quality model is called (AED2). The files required to run this model are located in the WQ dir. The AED2 model is called every 900s by the command WQ Update DT == 900. |
stability limits == 10.0,10.0 | |
Momentum mixing model == Smagorinsky | |
Global horizontal eddy viscosity == 0.2 | |
vertical mixing model == External | |
External turbulence model directory == ..\GOTM\ | |
turbulence update dt == 1800 | |
global vertical eddy viscosity limits == 1.0e-4, 1.0 | |
global vertical scalar diffusivity limits == 0., 1.0 | |
Water quality model == EXTERNAL | |
External water quality model dir == ..\AED2\ | |
WQ Update DT == 900 | |
! GEOMETRY | No syntax changes. |
spherical == 1 | |
geometry 2d == ..\geo\hydraul_006.2dm | |
Read GIS Nodestring == ..\model\gis\2d_ns_OpenBoundaries_001_L.shp | |
echo geometry netcdf == 1 | |
!3D Geometry Commands | No syntax changes. |
Vertical Mesh Type == z | |
Layer Faces == ..\geo\3D_Z_Layers_003.csv | |
Sigma Layers == 4 | |
Min Bottom Layer Thickness == 0.5 | |
Cell 3D Depth == 1.0 | |
! COHESIVE SEDIMENT CONTROL FILE | A link to the fvsed sediment control file is given here. The key inputs to the fvsed file include material properties such as bottom roughness. |
sediment control file == .\WQ_000.fvsed | |
! MATERIAL PROPERTIES | No syntax changes. |
Set Mat == 1 | |
Read GIS Mat == ..\model\gis\2d_mat_Materials_002_R.shp | |
material == 1 | |
bottom roughness == 0.02 | |
end material | |
material == 2 | |
bottom roughness == 0.75 | |
end material | |
material == 3 | |
bottom roughness == 1.3 | |
end material | |
! INITIAL CONDITIONS | An initial sediment concentration value is specified. Additionally, initial values for the AED2 water quality model are provided here. |
Initial Temperature == 20. | |
Initial Water Level == 0.052 | |
Initial salinity == 10. | |
Initial sediment concentration == 25 | |
!DO,Si,NH3,NOx,FRP,Pads,DOC,POC,DON,PON,DOP,POP | |
initial WQ concentration == 294., 10, 2.6, 4.6, 1.3, 0, 100, 100, 16.5, 49., 0.5, 1.4, 6 | |
! BOUNDARY CONDITIONS | New bc files are used containing additional WQ information. |
bc == Q, 1, ..\bc\Upstream_Q_Temp_Sal_Sed_WQ_001.csv | |
bc header == time_hr,flow_m3s-1,sal_ppt,temp_degC,FineSed_mgL-1,Sand_mgL-1,DO,Si,NH3,NOx,FRP,P_ADS,DOC,POC,DON,PON,DOP,POP,PHY01 | |
end bc | |
bc == WL, 2, ..\bc\Downstream_H_Temp_Sal_Sed_WQ_001.csv | |
bc header == time_hr,tidelevel_m,sal_ppt,temp_degC,FineSed_mgL-1,Sand_mgL-1,DO,Si,NH3,NOx,FRP,P_ADS,DOC,POC,DON,PON,DOP,POP,PHY01 | |
end bc | |
!MET | |
include == ..\bc\met\MET_2011.fvc | |
! OUTPUT COMMANDS | Additional WQ outputs are included as a netcdf file. |
output dir == ..\Output\ | |
output == datv | |
output parameters == h,v,d, temp, sal, Air_temp,Rel_hum,LW_rad,SW_rad | |
output interval == 3600 | |
end output | |
output == netcdf | |
output parameters == h,v,d, temp, sal, Rhow | |
output interval == 3600 | |
end output | |
output == netcdf | |
suffix == WQ | |
output parameters == wq_all | |
output interval == 3600 | |
end output | |
output == flux | |
output interval == 900 | |
end output | |
output == mass | |
output interval == 900 | |
end output |
To return to Tutorial Module 5 please refer to Tutorial Module 5.