WQ 000

From TUFLOW FV Wiki
Jump to: navigation, 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.