# This is a TOML configuration file for Wflow.
# Relative file paths are interpreted as being relative to this TOML file.
# Wflow documentation https://deltares.github.io/Wflow.jl/dev/
# TOML documentation: https://github.com/toml-lang/toml

calendar = "proleptic_gregorian"
endtime = 2000-06-30T00:00:00
starttime = 2000-05-31T00:00:00
time_units = "days since 1900-01-01 00:00:00"
timestepsecs = 86400
dir_input = "data/input"
dir_output = "data/output"

[state]
path_input = "instates-example-sbm-gwf.nc"
path_output = "outstates-example-sbm-gwf.nc"

# if listed, the variable must be present in the NetCDF or error
# if not listed, the variable can get a default value if it has one

[state.vertical.interception.variables]
canopy_storage = "canopystorage"

[state.vertical.soil.variables]
satwaterdepth = "satwaterdepth"
ustorelayerdepth = "ustorelayerdepth"

[state.lateral.river]
h = "h_river"
h_av = "h_av_river"
q = "q_river"

[state.lateral.land]
h = "h_land"
h_av = "h_av_land"
q = "q_land"

[state.lateral.subsurface.flow.aquifer]
head = "head"

[input]
path_forcing = "forcing-sbm-groundwater-part*.nc"
path_static = "staticmaps-sbm-groundwater.nc"

# these are not directly part of the model
ldd = "wflow_ldd"
river_location = "wflow_river"
subcatchment = "wflow_subcatch"
altitude = "wflow_dem"

# specify the internal IDs of the parameters which vary over time
# the external name mapping needs to be below together with the other mappings
forcing = [
  "vertical.atmospheric_forcing.precipitation",
  "vertical.atmospheric_forcing.potential_evaporation",
]

[input.vertical.vegetation_parameter_set]
rootingdepth = "rootingdepth"
canopygapfraction = "canopygapfraction"
cmax = "cmax"

[input.vertical.interception.parameters]
e_r = "EoverR"

[input.vertical.atmospheric_forcing]
potential_evaporation = "PET"
precipitation = "P"

[input.vertical.soil.parameters]
infiltcapsoil = "InfiltCapSoil"
kv_0 = "kv"
pathfrac = "PathFrac"
soilthickness = "soilthickness"
theta_r = "thetaR"
theta_s = "thetaS"

[input.lateral.river]
length = "wflow_riverlength"
n = "N_river"
slope = "RiverSlope"
width = "wflow_riverwidth"

[input.lateral.land]
n = "N"
slope = "Slope"

[input.lateral.subsurface]
conductivity = "k"
constant_head = "constant_head"
drain = "drain"
drain_conductance = "cond_drain"
drain_elevation = "elev_drain"
exfiltration_conductance = "exfiltration_conductance"
infiltration_conductance = "infiltration_conductance"
river_bottom = "river_bottom"
specific_yield = "specific_yield"

[model]
constanthead = true
drains = true
kin_wave_iteration = true
reinit = true
sizeinmetres = true
snow = false
thicknesslayers = [100, 300, 800]
type = "sbm_gwf"

[output]
path = "output_example-sbm-gwf.nc"

[output.vertical.interception.variables]
canopy_storage = "canopystorage"

[output.vertical.soil.variables]
satwaterdepth = "satwaterdepth"
ustorelayerdepth = "ustorelayerdepth"

[output.vertical.soil.parameters]
soilthickness = "soilthickness"

[output.lateral.river]
q = "q"

[output.lateral.subsurface.flow.aquifer]
head = "head"

[output.lateral.subsurface.recharge]
rate = "rate"

[output.lateral.subsurface.drain]
flux = "drain_flux"

[output.lateral.subsurface.river]
flux = "flux"

[csv]
path = "output_example-sbm-gwf.csv"

[[csv.column]]
header = "Q_av"
index = 5
parameter = "lateral.river.q_av"

[[csv.column]]
header = "head"
index = 5
parameter = "lateral.subsurface.flow.aquifer.head"