!> \mainpage Delft3D-WAQ documentation !> Use the first \"Modules\" item or tab to reach grouped tables of source code\n !> The following groups are distinguished: !> \li input processing !> \li functional modules for input processing !> \li advection diffusion solvers !> \li functional modules for simulation !> \li water quality processes !> !> Use the \"File List\" item or \"Files\" tab to get an alphabetic list of all source files !< !> \defgroup input input processing !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !>
source fileGroupFunction
dlwq00.fMain program of the DELWAQ input processing step
dlwq01.f1Reads the model identification and substances IDs
dlwq02.f2Reads integration method; monitoring areas/transects and timers
dlwq03.f3Reads grid layout; features and the computational volumes
dlwq04.f4Reads flow dimensions and pointers and all transport information
dlwq05.f5Reads all inputs associated with open boundaries
dlwq06.f6Reads all inputs associated with waste loads and withdrawals
dlwq07.f7Reads all model constants, variables and distributed variables/functions - old style
dlwq7a.f7Reads all model constants, variables and distributed variables/functions - new style
dlwq08.f8Reads initial conditions
dlwq09.f9Defines variables for output per available output file
dlwqs1.fsDefines process steering for statistical output processing
dlwqp1.fpDefines process steering for all water quality processing
space.fComputes and reports array space sizes and sets memory pointers at run time
!> \defgroup inp_funcs functional modules for input processing !> !> !> !> !> !>
source fileRole
opt0.f Top routine of system to read old style matrix input
dlwq5a.fNew style intuitive flexible input processing for boundaries and waste loads
dlwq7a.fNew style intuitive flexible input processing for constants and functions
!> \defgroup solvers advection diffusion solvers !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !>
methodsolverNumerical methodstability
delwaq2.fMain program of the simulations step
1 dlwqn1.fFirst order upwind in space and timeexplicit
!> !> !>
dlwq16.fFills derivative with advection diffusion processes
dlwq18.fSets an explicit time step from the derivative
2 dlwqn2.fFirst order upwind in space 2nd order Runge Kutta in timeexplicit
!> !> !> !> !>
dlwq16.fFills derivative with advection diffusion processes
dlwq20.fSets a preliminary step of dt to provide concentrations
dlwq16.fFills derivative with advection diffusion processes from these new concentrations
dlwq18.fSets a full final step of 2*dt
3 dlwqn3.fSecond order Lax Wendrof in space and timeexplicit
!> !> !>
dlwq30.fFills derivative with advection diffusion processes
dlwq18.fSets an explicit time step from the derivative
4 dlwqn4.fHorizontal ADI on regular grid, vertical implicit centralhalf implicit
!> !> !> !> !> !>
dlwq40.fFills derivative with advection diffusion processes in one horizontal explicit direction
dlwq42.fSets an explicit time step from the derivative
dlwq43.fFills derivative with advection diffusion processes in other horizontal direction, performs implicit double sweep
dlwq44.fUpdate the necessary arrays
dlwq46.fMass balances for this method
5 dlwqn5.fFlux corrected transport according to Boris and Book (FCT)explicit
!> !> !> !> !>
dlwq50.fFills derivative with advection processes only
dlwq18.fSets an explicit time step from the derivative
dlwq51.fPerform a flux correction step on this preliminary result
dlwq52.fUpdate the necessary arrays
6 dlwqn6.fStationary first order upwind in space, direct methodimplicit
!> !> !> !> !> !> !> !> !>
dlwq60.fScales derivative coming from processes routine
dlwq61.fFills the diagonal for 1 substance
dlwq62.fFills the sparce matrix with advection diffusion terms
dlwq67.fSets diagonal entries of approximately 0.0 at 1.0
delmat.fDirect matrix solver
dlwq63.fUpdate the necessary arrays
dlwq64.fMakes mass balances for this method
dlwq66.fUpdate the necessary arrays
7 dlwqn7.fStationary central in space, direct methodimplicit
!> !> !> !> !> !> !> !> !>
dlwq60.fScales derivative coming from processes routine
dlwq61.fFills the diagonal for 1 substance
dlwq70.fFills the sparce matrix with advection diffusion terms
dlwq67.fSets diagonal entries of approximately 0.0 at 1.0
delmat.fDirect matrix solver
dlwq63.fUpdate the necessary arrays
dlwq71.fMakes mass balances for this method
dlwq66.fUpdate the necessary arrays
8 dlwqn8.fStationary first order upwind in space, iterative methodobsolete

This method is obsolete

9 dlwqn9.fStationary central in space, iterative methodobsolete

This method is obsolete

10 dlwqnb.fFirst order upwind in space, implicit, direct methodimplicit
!> !> !> !> !> !> !> !> !> !> !>
dlwqb8.fRestores CONC array from MASS array
dlwqb3.fNew volumes from optional computational volumes
dlwqb1.fFills the diagonal for 1 substance
dlwqb6.fComputes right hand side
dlwqb7.fFills the sparce matrix with advection diffusion terms
dlwqb9.fFills the sparce matrix with advection diffusion terms
delmat.fDirect matrix solver
dlwqb2.fUpdate the necessary arrays
dlwqb5.fMakes mass balances for this method
dlwqb4.fUpdate the necessary arrays
11 dlwqnc.fFirst order upwind in space, implicit central verticalexplicit
!> !> !> !> !> !> !>
dlwq16.fFills derivative with advection diffusion processes
dlwq42.fSets an explicit time step from the derivative
dlwqd1.fSets an implicit time step for the central transport of the vertical
dlwqd2.fForester filter to avoid spurious extremes
dlwq44.fUpdate the necessary arrays
dlwq46.fMass balances for this method
12 dlwqnd.fFCT horizontal, central implicit verticalexplicit
!> !> !> !> !> !> !> !> !> !>
dlwq50.fFills derivative with advection processes only
dlwq18.fSets an explicit time step from the derivative
dlwq51.fPerform a flux correction step on this preliminary result
dlwq52.fUpdate the necessary arrays
dlwq42.fSets volumes on the diagonal
dlwqd1.fSets an implicit time step for the central transport of the vertical
dlwqd2.fForester filter to avoid spurious extremes
dlwq44.fUpdate the necessary arrays
dlwq46.fMass balances for this method
13 dlwqne.fFirst order upwind in space, upwind implicit verticalexplicit
!> !> !> !> !> !>
dlwq16.fFills derivative with advection diffusion processes
dlwq42.fSets volumes on the diagonal
dlwqe1.fSets an implicit time step for the upwind transport of the vertical
dlwq44.fUpdate the necessary arrays
dlwq46.fMass balances for this method
14 dlwqna.fFCT horizontal, upwind implicit verticalexplicit
!> !> !> !> !> !> !> !> !>
dlwq50.fFills derivative with advection processes only
dlwq18.fSets an explicit time step from the derivative
dlwq51.fPerform a flux correction step on this preliminary result
dlwq52.fUpdate the necessary arrays
dlwq42.fSets volumes on the diagonal
dlwqe1.fSets an implicit time step for the upwind transport of the vertical
dlwq44.fUpdate the necessary arrays
dlwq46.fMass balances for this method
15 dlwqnf.fFirst order upwind, iterative GMRES solverimplicit
!> !> !> !> !> !> !> !> !> !> !>
dlwqf1.fInitializes the necessary pointers for matrices and vectors for fast solver
dlwqb8.fRestores CONC array from MASS array
dlwqm7.fMakes mixing length array according to numerical settings
dlwqf2.f90Initializes the diagonal
dlwqf3.f90Fills the sparce matrix with advection and diffusion terms (one substance at a time)
dlwqf4.f90Computes Right Hand Side (one substance at a time)
sgmres.f90The Krilov subspace iterative sparce matrix solver
dlwqf6.f90Copies the solution for this substance in the total substances array
dlwqb5.fMakes mass balances for this method
dlwqb4.fUpdate the necessary arrays
16 dlwqng.fUpwind horizontally, central vertically, GMRES solverimplicit
!> !> !> !> !> !> !> !> !> !>
dlwqf1.fInitializes the necessary pointers for matrices and vectors for fast solver
dlwqb8.fRestores CONC array from MASS array
dlwqf2.f90Initializes the diagonal
dlwqg3.f90Fills the sparce matrix with advection and diffusion terms (one substance at a time)
dlwqf4.f90Computes Right Hand Side (one substance at a time)
sgmres.f90The Krilov subspace iterative sparce matrix solver
dlwqf6.f90Copies the solution for this substance in the total substances array
dlwqb5.fMakes mass balances for this method
dlwqb4.fUpdate the necessary arrays
17 dlwqnh.fFirst order upwind, iterative stationary GMRES solver implicit
!> !> !> !> !> !> !> !> !>
dlwqm7.fMakes mixing length array according to numerical settings
dlwqh1.fInitializes the diagonal
dlwqf3.f90Fills the sparce matrix with advection and diffusion terms (one substance at a time)
dlwqh3.fComputes Right Hand Side (one substance at a time)
sgmres.f90The Krilov subspace iterative sparce matrix solver
dlwqf6.f90Copies the solution for this substance in the total substances array
dlwq64.fMakes mass balances for this method
dlwq66.fUpdate the necessary arrays
18 dlwqni.fUpwind horizontally, central vertically, stationary GMRES solver implicit
!> !> !> !> !> !> !> !> !>
dlwqm7.fMakes mixing length array according to numerical settings
dlwqh1.fInitializes the diagonal
dlwqf3.f90Fills the sparce matrix with advection and diffusion terms (one substance at a time)
dlwqh3.fComputes Right Hand Side (one substance at a time)
sgmres.f90The Krilov subspace iterative sparce matrix solver
dlwqf6.f90Copies the solution for this substance in the total substances array
dlwq64.fMakes mass balances for this method
dlwq66.fUpdate the necessary arrays
19 dlwqnj.fADI solver of Delft3D-FLOW (difu) 3rd order, upwind vertically half implicit
!> !> !> !> !> !> !>
dlinit.fInitializes the Delft3D-Flow ADE solver arrays
dlconv.fConverts the hydrodynamics to Delft3D-Flow ADE full matrix arrays
dldifu.fAims to be identical to the Delft3D-Flow difu.f90 source file. Is difficult to maintain.
dlflux.fComputes fluxes for balances for this solver
dlmasb.fMakes the mass balances from these fluxes
dlback.fBack conversion of some arrays
20 dlwqnj.fADI solver of Delft3D-FLOW (difu) 3rd order, central vertically half implicit
!> !> !> !> !> !> !> !>
dlinit.fInitializes the Delft3D-Flow ADE solver arrays
dlconv.fConverts the hydrodynamics to Delft3D-Flow ADE full matrix arrays
dldifu.fAims to be identical to the Delft3D-Flow difu.f90 source file. Is difficult to maintain.
dlflux.fComputes fluxes for balances for this solver
dlmasb.fMakes the mass balances from these fluxes
dlback.fBack conversion of some arrays
dlwqd2.fForester filter to avoid spurious extremes
21 dlwqnm.fSelf adjusting theta method, Salezac limiter implicit
!> !> !> !> !> !> !> !> !> !> !> !> !>
dlwqf1.fInitializes the necessary pointers for matrices and vectors for fast solver
dlwqm7.fMakes mixing length array according to numerical settings
dlwqm0.fMakes specific flow and dispersion arrays for this substance
dlwqm1.fComputes space and time varyin theta coefficients
dlwq_output_theta.fSaves a theta field for output for the first substance only (often continuity)
dlwqm2.fConstructs the sparse advection diffusion matrix
dlwqm3.fConstructs the right hand side
sgmres.f90The Krilov subspace iterative sparce matrix solver
dlwqm4.fMakes the mass balance
dlwqm5.fFlux correction according to Salezac
dlwqb4.fUpdate the necessary arrays
proint.fIntegrate the fluxes at dump segments
22 dlwqnm.fSelf adjusting theta method, Boris and Book limiter implicit
!> !> !> !> !> !> !> !> !> !> !> !> !>
dlwqf1.fInitializes the necessary pointers for matrices and vectors for fast solver
dlwqm7.fMakes mixing length array according to numerical settings
dlwqm0.fMakes specific flow and dispersion arrays for this substance
dlwqm1.fComputes space and time varyin theta coefficients
dlwq_output_theta.fSaves a theta field for output for the first substance only (often continuity)
dlwqm2.fConstructs the sparse advection diffusion matrix
dlwqm3.fConstructs the right hand side
sgmres.f90The Krilov subspace iterative sparce matrix solver
dlwqm4.fMakes the mass balance
dlwqm8.fFlux correction according to Boris and Book
dlwqb4.fUpdate the necessary arrays
proint.fIntegrate the fluxes at dump segments
!> \defgroup run_funcs functional modules for simulation !> !> !> !> !> !> !> !> !> !> !>
source fileRole
dlwqi0.fInitialises the system
proces.fDeals with all water quality processes
dlwqo2.fDeals with all output of model results
dlwq17.fDeals with the resolving of open boundaries (dlwq_boundio should be included in here)
dlwq15.fDeals with the resolving of all waste loads
dlwqt0.fDeals with the resolving of all other external forcing
dlwqce.fMakes closure error corrections
dlwq13.fWrites the restart file
!> \defgroup waterquality water quality processes !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !> !>
groupsource filefunctionremark
age residense time
restim.fResidence time per volume, for advective transport only
weirch.fCalculates the residence time of a volume
watage.fAge of water through the tracer substances
statistics
staday.fPeriodic (day) average of a given substance
stadpt.fDepth-averaged, max and min value per timestep
stadsc.fMean, min, max, stdev of an output variable
stageo.fGeometric mean of a variable during a certian time span
staprc.fExceedence frequency, its complement and the mean
staqtl.fQuantiles for a given substance during a given period
bacteria
bacmrt.fMortality of bacteria depending on UV-light, salinity and temperature
groab.fAerobic growth of nitrifying autotrophic bact.
grohb.fGrowth of heterotrophic bact. on Cdom and Cfom
lysis.fLysis of bacteria
sedcar.fSedimentation routine used for IMx, OOC, algae, BOD pools, bacteria etc.
sumcol.fSum of all Coliform Bacteria
inorganic sediment
atmdep.fAtmosferic deposition and diffuse input of IMx, N, P, Org_us and Metals
burcar.fBurial of IMx, DetC, OOC, Green, Diat and AAP
calsnd.fTransport of non-cohesive sediment
calsed.fSedimentation velocity IMx, DetC OOC, BODC, all algea = f (Temp SS Sal)
digcar.fDigging of IMx, DetC, OOC, Green, Diat and AAP
dmvol.fVolume of dry matter in a segment
rescar.fResuspension of particulates by fraction * resdm, from S2 if S1 is exhausted
sedcar.fSedimentation routine used for IMx, OOC, algae, BOD pools, bacteria etc.
wkcomp.fComputes sum parameters from fractions (GEM)
BOD oxygen OC
bod.fCalculates the bod5 and bod-inf from available model substances
bodcod.fDecay of BOD, COD and NBOD and associated oxygen consumption
botmin.fMineralisation of organic substances and desorption of AAP in the bed (S1,S2) for C, N, P and Si.
burcar.fBurial of IMx, DetC, OOC, Green, Diat and AAP
calsed.fSedimentation velocity IMx, DetC OOC, BODC, all algea = f (Temp SS Sal)
d40blo.fBLOOM II algae module
d40swi.fSediment-water interaction, the SWITCH module
decbod.fOxydation of BOD-fractions with Monod kinetics for the TEWOR models
decfsn.fMineralisation of fast composing detritus and conversion of the C:N and C:P ratiosalmost identical to decfst
decfst.fMineralisation of fast composing detritus with equal rations for C, N and Pwhy 2 separate routines ?
decref.fMineralisation of refractory detritus C, N and Pwhy 2 separate routines ?
decren.fMineralisation of refractory detritus C, N and Pis identical to decref !!
decsln.fMineralisation of fast composing detritus and conversion of the C:N and C:P ratiosalmost identical to decslw
decslw.fMineralisation of fast composing detritus with equal rations for C, N and Pwhy 2 separate routines ?
densed.fDenitrification in sediment
denwat.fDenitrification in water column
digcar.fDigging of IMx, DetC, OOC, Green, Diat and AAP
gemmfb.fGEM microfytobenthos production
gemmnd.fGEM algae primary production module
gemmin.fGEM mineralisation
groab.fAerobic growth of nitrifying autotrophic bact.
makpoc.fDerive OOC from IM-fractions and percentage POM in IMx
methox.fOxidation of methane (new and generic)
minlim.fMineralisation organic C, N, P and Si
nitrif.fNitrification of ammonium + decay of CBOD
oxymin.fPotential daily mimimum dissolved oxygen concentration
posoxy.fReturns positive oxygen concentration or zero
pripro.fNett primary production and mortality DYNAMO algae
ptewor.fProduction fluxes for TEWOR+
rear.fReaeration of carbon dioxide and oxygen
rescar.fResuspension of particulates by fraction * resdm, from S2 if S1 is exhausted
sdppro.fTraditional algal growth module (DYNAMO)
satco2.fSaturation concentration carbon dioxide
satoxy.fSaturation concentration of oxygen
sedcar.fSedimentation routine used for IMx, OOC, algae, BOD pools, bacteria etc.
sednut.fSedimentation of nutrients in the organic carbon matrix (GEM)
sedox.fSediment oxygen demand
sedsod.fSedimentation of oxygen demand
strear.fAeration at weirs (Gameson and Nakasone) (input is array of structures)
swbur.fFluxes in the water bed due to burial and digging
varoxy.fVariation of oxygen due to variation in primary production within day
weirox.fRearation at weirs (input is a parameter indicating weir or not)
wkcomp.fComposition
nutrients
adspo4.fP-ad/desorption to particulate inorganic matter. 3 options for sorption formulation.
atmdep.fAtmosferic deposition and diffuse input of IMx, N, P, Org_us and Metals
botmin.fMineralisation of organic substances and desorption of AAP in the bed (S1,S2) for C, N, P and Si.
burcar.fBurial of IMx, DetC, OOC, Green, Diat and AAP
burnut.fBurial of nutrients in organicC from S2
consbl.fGrazing module
d40blo.fBLOOM II algae module
d40swi.fSediment-water interaction, the SWITCH module
decfsn.fMineralisation of fast composing detritus and conversion of the C:N and C:P ratiosalmost identical to decfst
decfst.fMineralisation of fast composing detritus with equal rations for C, N and Pwhy 2 separate routines ?
decref.fMineralisation of refractory detritus C, N and Pwhy 2 separate routines ?
decren.fMineralisation of refractory detritus C, N and Pis identical to decref !!
decsln.fMineralisation of fast composing detritus and conversion of the C:N and C:P ratiosalmost identical to decslw
decslw.fMineralisation of fast composing detritus with equal rations for C, N and Pwhy 2 separate routines ?
densed.fDenitrification in sediment
denwat.fDenitrification in water column
digcar.fDigging of IMx, DetC, OOC, Green, Diat and AAP
dignut.fDigging of nutrients in organicC
dissi.fDissolution of Si in opal (SWITCH defaults)
dmvol.fVolume of dry matter in a segment
explfl.fSimulation of explosive nutrient sediment release
ferdom.fFermentation of dissolved fermentable org. matter
gemcmp.fGEM computation of C, N and P from stoichiometry
gemmfb.fGEM microfytobenthos production
gemmin.fGEM mineralisation
gemmnd.fGEM algae primary production module
groab.fAerobic growth of nitrifying autotrophic bact.
grohb.fGrowth of heterotrophic bact. on Cdom and Cfom
hydpom.fHydrolysis of degradable particulate org. matter
lysis.fLysis of bacteria
mfbnut.fMicrofytobenthos nutrient concentration in the bottom GEM
minlim.fMineralisation organic C, N, P and Si
nh3fre.fCalculation conc. unionized ammonia
nitrif.fNitrification of ammonium + decay of CBOD
nralgs.fNutrient release of algae in S1 and S2
nutcnk.fNutrient cycle model by Nakata and Kuramoto
nutrel.fRelease (nutrients/detritus) by of mortality algae
nutupt.fUptake of nutrients by growth of algae
oyster.fOysterfarming as forcing function for primary consumption
phcomb.fCalculates total C, P, N, Si, Dm, Chlorophyll in algae from fractions in Bloom
ptewor.fProduction fluxes for TEWOR+
resant.fResuspension of nutrients in organic carbon matrix
rescar.fResuspension of particulates by fraction * resdm, from S2 if S1 is exhausted
resnut.fResuspension (from 2 layers) of nutrients in the organic carbon matrix
sdppro.fTraditional algal growth module (DYNAMO)
sedaap.fSedimentation flux and velocity for PAP and AAP (adsorbed PO4)
sedcar.fSedimentation routine used for IMx, OOC, algae, BOD pools, bacteria etc.
sednal.fSedimentation of nutrients in algae, diatoms
sednut.fSedimentation of nutrients in the organic carbon matrix (GEM)
swbura.fSedimentation of nutrients in organic carbon matrix (gx/m2/d => gC/m2/d)
swburn.fSedimentation of nutrients in organic carbon matrix (gC/m2/d => gx/m2/d)
swsedn.fSedimentation of nutrients in organic carbon matrix (SWITCH)
vertnk.fVert.transport TIM for nutr. cycle model (Nakata and Kuramoto)
vivia2.fDissolution/precipitation of P in vivianitean error(?) corrected in vivian
vivian.fDissolution/precipitation of P in vivianitewhy 2 separate routines ?
wkcomp.fComposition
zoodyn.fDynamic calculation of grazing and zooplankton biomass
algae growth
algmrt.fAlgea mortality of greens and diatoms in the bed (S1,S2) DYNAMO
bluetd.f
burcar.fBurial of IMx, DetC, OOC, Green, Diat and AAP
consbl.fGrazing module
d40blo.fBLOOM II algae module
depave.fAverage depth for a Bloom time step (typically a day)
digcar.fDigging of IMx, DetC, OOC, Green, Diat and AAP
dlalg.fDaylength function for algae DYNAMO
dmvol.fVolume of dry matter in a segment
extinc.fExtinction formulae ( inorganic , organic and algae )
gemmfb.fGEM microfytobenthos production
gemmnd.fGEM algae primary production module
gemnlm.fGEM algae nutrient limitation function
gemtmp.fGEM algae temperature limitation function
gmreed.fGrowth, mortality and decay of reed as helofyte filter
mfbllm.fMicrofytobenthos light limitation function
mndini.fInitialisation MND algae GEM
mndllm.fMND algae light limitation function GEM
nlalg.fNutrient limiation function for green algae
oyster.fOysterfarming as forcing function for primary consumption
phcomb.fCalculates total C, P, N, Si, Dm, Chlorophyll in algae from fractions in Bloom
phcomp.fComposition of phytoplankton by summing algae fractions - Dynamo - GEM
pocomp.fComposition of POC by summing algae fractions Dynamo, Bloom and GEM
pprlim.fLimitation (numerical) on primary production DYNAMO
pripro.fNett primary production and mortality DYNAMO algae
radalg.fLight efficiency function DYNAMO algae
rdbalg.fLight efficiency function DYNAMO algae
resalg.fResuspension of algae using resdm * fraction
ptewor.fProduction fluxes for TEWOR+
rescar.fResuspension of particulates by fraction * resdm, from S2 if S1 is exhausted
sdppro.fNett primary production and mortality diatoms
secchi.fExtinction of visible-light (370-680nm) by inorganic, organic and algae
sednal.fSedimentation of nutrients in algae, diatoms
sedphy.fCalculates the settling velocity of algae (BLOOM)
ssedph.fSum of sedimentation flux of algae Dynamo - Bloom - GEM
stvar.fDummy for consbl to save the state variables
swbur.fFluxes in the water bed due to burial and digging
swbura.fBurial nutrients in diatoms from sediment S1
tfalg.fTemperature functions for algae growth and mortality
ulfix.fFixation of BLOOM algae at the water bed (e.g. for Ulvae)
varoxy.fVariation of primary production within day
vtrans.fVertical distribution after a longer time span to correct 3D-BLOOM
wkcomp.fComposition
zoodyn.fGrazing module
Micro pollutants
burhm.fBurial of all heavy metals
buromv.fBurial of all organic micros from S1 and S2
cec.fCalculation Cation Exchange Capacity
degmp.fDegradation of organic micropolutants
digfix.fDigging flux of reversible sorbed HM/OMP
dighm.fDigging of all heavy metals
digomv.fDigging of all organic micros towards S1 and S2
ksorhm.fSlow adsorption and desorption of mive's. 3 compartment model instantaneous + 2 first order
ksorom.fSlow adsorption and desorption of omive's. 3 compartment model instantaneous + 2 first order
metres.fResuspension of metals for western scheldt estuary
norhm.fNormalisation of heavy metal content of the total solid
partmp.fPartitioning of micropollutants
resfix.fResuspension fluxes for HM/OMP are multiplied with fraction fixed HM/OMP
reshm.fResuspension of adsorbed heavy metal
resomv.fResuspension adsorbed organic micro pollutants
rfpart.fReprofunctions for HM partition coefficients
sedfix.fMultiplication sedimentatio fluxes HM/OMP with fraction fixed HM/OMP
sedhm.fSedimentation flux and velocity of adsorbed heavy metals
sedomv.fSedimentation flux and velocity of adsorbed organic micro pollutants
suloxi.fDecay of metal sulfides for western scheldt estuary
sulsed.fSedimentation of metal sulfides for Western Scheldt estuary
swoxy.fPartitioning switch in WC, S1 and S2 based on actual and critical oxygen concentration
trcoef.fGas and liquid exchange organic micro pollutants (Lyman and O'Connor)
vervlu.fAtmospheric exch. OMPs (volatilization/intake)
watmin.fGeneral mineralisation routine 0+1st order, temperature corrected
Physics
calchz.fCalculate chezy coefficient using roughness and depth
caltau.fCalculation of bottom friction
calvs.fCalculate individual sedimentation velocities
calwav.fWave characteristics
chlor.fCalculation of chloride from salinity
clcrad.fRadiation at segment upper and lower boundaries
dayl.fDaylength calculation
ddelt.fScaling of model DT to aux time step
ddepth.fDynamic calculation of the depth
dsptra.fDispersion/diffusion in the sediment
dsurf.fDynamic calculation of the horizontal surface area from volume and depth
emersi.femersion of segments in z-layers, set segment features accordingly
extina.fExtinction of light by algae and POC
extinc.fExtinction formulae ( inorganic , organic and algae )
grd.f90gradient in density or in the horizontal stream velocity
hdisp.f(1D) Horizontal dispersion as velocity dependent reprofunction
hdispv.f(2D) Horizontal dispersion as velocity dependent reprofunction
heatfl.fTotal heat flux for surface water absolute temperature model
intpol.fDepth where wave is created or wind fetch from wind direction
linpol.fLinear interpolation q-h relation
meteo.fProcess meteo from various meteo-stations
s2x.f90Conversion of variables per volume to per exchange by linear interpolation)
stox3d.fVertical dispersion (segment -> exchange)
temper.fExchange of excess temperature at the surface (Sweers)
totdep.fTotal depth water column
varsal.fsalinity in case of constant river discharge
vdisp.fVertical dispersion in a stratified water body using Richardson's formulae
veloc.f90Horizontal stream velocity in a segment based on flows at interfaces
velocv.f90Horizontal stream velocity in a segment based on velocities at interfaces
xtos3d.fConversion of units on interfaces to volumes for 3rd direction only
xveloc.fCalculation of velocity's on exchanges from flows
bed behaviour
advtra.fAdvective transport, velocities and fluxes, of solids in sediment
burial.fBurial total bottom mass (dry matter)
burcar.fBurial of IMx, DetC, OOC, Green, Diat and AAP
calsnd.fTransport of non-cohesive sediment
digcar.fDigging of IMx, DetC, OOC, Green, Diat and AAP
diggin.fDigging dry matter to sediment S1 and S2
rescar.fResuspension of particulates by fraction * resdm, from S2 if S1 is exhausted
sedcar.fSedimentation routine used for IMx, OOC, algae, BOD pools, bacteria etc.
sedcom.fComposition, thickness, total dry mass and density in sediment layers
somres.fSum of (in)organic substances resusp./buried/digged from S1/S2
somsed.fTotal of all sedimenting substances
trase2.fTotal of transport in sediment for 66 substances
chemistry
cselac.fConsumption of electron acceptors
d40cha.fCHARON chemical equilibrium module
ebuch4.fEbullition of methane
feadso.fComputes Fe adsorbens for P for WESTERN SCHELDT ESTUARY only!!
resdm.fResuspension total bottom material (dry mass)
salin.fConverts chloride concentration to salinity (Aquatic Chemistry p 567)
satch4.fMethane saturation concentration based on atmospheric methane pressure
simph.fSimple calculation of pH
sulfid.fSpeciation of dissolved sulphide (S= and HS-) in pore water
sulfox.fOxidation of dissolved sulphide (0 and 2nd order) (new, generic !)
sulfpr.fPrecipitation and dissolution of sulphide as first order process