# 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

dir_input = "data/input"
dir_output = "data/output"

[time]
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

[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.variables]
vegetation_canopy_water__depth = "canopystorage"

soil_water_sat-zone__depth = "satwaterdepth"
soil_layer_water_unsat-zone__depth = "ustorelayerdepth"

river_water__instantaneous_depth = "h_river"
river_water__instantaneous_volume_flow_rate = "q_river"

land_surface_water__instantaneous_volume_flow_rate = "q_land"
land_surface_water__instantaneous_depth = "h_land"

subsurface_water__instantaneous_hydraulic_head = "head"

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

# these are not directly part of the model
basin__local_drain_direction = "wflow_ldd"
river_location__mask = "wflow_river"
subbasin_location__count = "wflow_subcatch"

[input.forcing]
atmosphere_water__precipitation_volume_flux = "P"
land_surface_water__potential_evaporation_volume_flux = "PET"

[input.static]
soil_surface_water__vertical_saturated_hydraulic_conductivity = "kv"
soil_water__vertical_saturated_hydraulic_conductivity_scale_parameter.value = 0.001
soil_layer_water__brooks-corey_exponent.value = [10.0, 10.0, 10.0, 10.0]
soil_water__residual_volume_fraction = "thetaR"
soil_water__saturated_volume_fraction = "thetaS"
"soil~compacted__area_fraction" = "PathFrac"
soil__thickness = "soilthickness"

vegetation_canopy_water__mean_evaporation-to-mean_precipitation_ratio = "EoverR"
vegetation_root__depth = "rootingdepth"
vegetation_canopy__gap_fraction = "canopygapfraction"
vegetation_water__storage_capacity = "cmax"

river__length = "wflow_riverlength"
river_water_flow__manning_n_parameter = "N_river"
river__slope = "RiverSlope"
river__width = "wflow_riverwidth"

land_surface_water_flow__manning_n_parameter = "N"
land_surface__elevation = "wflow_dem"
land_surface__slope = "Slope"

subsurface_surface_water__horizontal_saturated_hydraulic_conductivity = "k"
"model_boundary_condition~constant_hydraulic_head" = "constant_head"
land_drain_location__mask = "drain"
land_drain__conductance = "cond_drain"
land_drain__elevation = "elev_drain"
subsurface_water__specific_yield = "specific_yield"
river_bottom__elevation = "river_bottom"
river_water__infiltration_conductance = "infiltration_conductance"
river_water__exfiltration_conductance = "infiltration_conductance"

[model]
constanthead__flag = true
drain__flag = true
kinematic_wave__adaptive_time_step_flag = true
cold_start__flag = true
cell_length_in_meter__flag = true
snow__flag = false
soil_layer__thickness = [100, 300, 800]
type = "sbm_gwf"

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

[output.netcdf_grid.variables]
vegetation_canopy_water__depth = "canopystorage"
soil_water_sat-zone__depth = "satwaterdepth"
soil_layer_water_unsat-zone__depth = "ustorelayerdepth"
land.soil.parameters.soilthickness = "soilthickness"
river_water__volume_flow_rate = "q_av"
subsurface_water__hydraulic_head = "head"
subsurface_water_sat-zone_top__net_recharge_volume_flow_rate = "rate"
"land_drain_water~to-subsurface__volume_flow_rate" = "drain_flux"
"river_water~to-subsurface__volume_flow_rate" = "flux"

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

[[output.csv.column]]
header = "Q_av"
index = 5
parameter = "river_water__volume_flow_rate"

[[output.csv.column]]
header = "head"
index = 5
parameter = "subsurface_water__hydraulic_head"