# 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"

[logging]
loglevel = "info"

[time]
calendar = "proleptic_gregorian"
endtime = 2000-02-01T00:00:00
starttime = 2000-01-01T00:00:00
time_units = "days since 1900-01-01 00:00:00"
timestepsecs = 86400

[state]
path_input = "instates-moselle.nc"
path_output = "outstates-moselle.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_surface__temperature = "tsoil"
soil_layer_water_unsat-zone__depth = "ustorelayerdepth"

"snowpack~dry__leq-depth" = "snow"
"snowpack~liquid__depth" = "snowwater"

river_water__instantaneous_depth = "h_river"
river_water__instantaneous_volume_flow_rate = "q_river"

reservoir_water_surface__instantaneous_elevation = "waterlevel_reservoir"

subsurface_water__volume_flow_rate  = "ssf"

land_surface_water__instantaneous_volume_flow_rate = "q_land"
land_surface_water__instantaneous_depth = "h_land"

[input]
path_forcing = "forcing-moselle.nc"
path_static = "staticmaps-moselle.nc"

# these are not directly part of the model
basin__local_drain_direction = "wflow_ldd"
river_location__mask = "wflow_river"
reservoir_area__count = "wflow_reservoirareas"
reservoir_location__count = "wflow_reservoirlocs"
subbasin_location__count = "wflow_subcatch"
river_gauge__count = "wflow_gauges_grdc"

[input.forcing]
atmosphere_water__precipitation_volume_flux = "precip"
land_surface_water__potential_evaporation_volume_flux = "pet"
atmosphere_air__temperature = "temp"

[input.static]
atmosphere_air__snowfall_temperature_threshold = "TT"
atmosphere_air__snowfall_temperature_interval = "TTI"

"land~water-covered__area_fraction" = "WaterFrac"

snowpack__melting_temperature_threshold = "TTM"
snowpack__degree-day_coefficient = "Cfmax"
snowpack__liquid_water_holding_capacity =  "WHC"

soil_layer_water__brooks-corey_exponent = "c"
soil_surface_water__infiltration_reduction_parameter = "cf_soil"
soil_water__vertical_saturated_hydraulic_conductivity_scale_parameter = "f"
"soil~compacted_surface_water__infiltration_capacity" = "InfiltCapPath"
soil_water__residual_volume_fraction = "thetaR"
soil_water__saturated_volume_fraction = "thetaS"
soil_water_sat-zone_bottom__max_leakage_volume_flux = "MaxLeakage"
"soil~compacted__area_fraction" = "PathFrac"
"soil_root~wet__sigmoid_function_shape_parameter" = "rootdistpar"
soil__thickness = "SoilThickness"

vegetation_canopy_water__mean_evaporation-to-mean_precipitation_ratio = "EoverR"
vegetation_canopy__light-extinction_coefficient = "Kext"
vegetation__specific-leaf_storage = "Sl"
vegetation_wood_water__storage_capacity = "Swood"
vegetation_root__depth = "RootingDepth"

river__length = "wflow_riverlength"
river_water_flow__manning_n_parameter = "N_River"
river__slope = "RiverSlope"
river__width = "wflow_riverwidth"
river_bank_water__depth = "RiverDepth"
river_bank_water__elevation = "RiverZ"

land_surface_water_flow__manning_n_parameter = "N"
land_surface__slope = "Slope"

reservoir_surface__area = "reservoir_area"
"reservoir_water_demand~required~downstream__volume_flow_rate" = "ResDemand"
reservoir_water_release-below-spillway__max_volume_flow_rate = "ResMaxRelease"
reservoir_water__max_volume = "ResMaxVolume"
"reservoir_water~full-target__volume_fraction" = "ResTargetFullFrac"
"reservoir_water~min-target__volume_fraction" = "ResTargetMinFrac"
reservoir_water_surface__initial_elevation = "waterlevel_reservoir"
reservoir_water__rating_curve_type_count = "outflowfunc"
reservoir_water__storage_curve_type_count = "storfunc"

subsurface_water__horizontal-to-vertical_saturated_hydraulic_conductivity_ratio = "KsatHorFrac"

[input.cyclic]
vegetation__leaf-area_index = "LAI"

[input.static.soil_surface_water__vertical_saturated_hydraulic_conductivity]
netcdf.variable.name = "KsatVer"
scale = 1.0
offset = 0.0

[model]
kinematic_wave__adaptive_time_step_flag = true
snow_gravitional_transport__flag = true
cold_start__flag = true
reservoir__flag = true
snow__flag = true
soil_layer__thickness = [100, 300, 800]
type = "sbm"
river_streamorder__min_count = 6
land_streamorder__min_count = 5

[output.netcdf_grid]
path = "output_moselle.nc"

[output.netcdf_grid.variables]
soil_water_sat-zone__depth = "satwaterdepth"
soil_surface__temperature = "tsoil"
soil_layer_water_unsat-zone__depth = "ustorelayerdepth"
"snowpack~dry__leq-depth" = "snow"
"snowpack~liquid__depth" = "snowwater"
river_water__depth = "h_av_river"
river_water__volume_flow_rate = "q_av_river"
reservoir_water__volume = "storage_reservoir"
subsurface_water__volume_flow_rate  = "ssf"
land_surface_water__volume_flow_rate = "q_av_land"
land_surface_water__depth = "h_av_land"
land.interception.variables.canopy_storage = "canopystorage"

[output.netcdf_scalar]
path = "output_scalar_moselle.nc"

[[output.netcdf_scalar.variable]]
name = "Q"
map = "river_gauge__count"
parameter = "routing.river_flow.variables.q"

[[output.netcdf_scalar.variable]]
coordinate.x = 6.255
coordinate.y = 50.012
name = "temp_coord"
location = "temp_bycoord"
parameter = "atmosphere_air__temperature"

[[output.netcdf_scalar.variable]]
location = "temp_byindex"
name = "temp_index"
index.x = 100
index.y = 264
parameter = "atmosphere_air__temperature"

[output.csv]
path = "output_moselle.csv"

[[output.csv.column]]
header = "Q"
parameter = "river_water__volume_flow_rate"
reducer = "maximum"

[[output.csv.column]]
header = "storage"
index = 1
parameter = "reservoir_water__volume"

[[output.csv.column]]
coordinate.x = 6.255
coordinate.y = 50.012
header = "temp_bycoord"
parameter = "atmosphere_air__temperature"

[[output.csv.column]]
coordinate.x = 6.255
coordinate.y = 50.012
header = "vwc_layer2_bycoord"
parameter = "land.soil.variables.vwc"
layer = 2

[[output.csv.column]]
header = "temp_byindex"
index.x = 100
index.y = 264
parameter = "atmosphere_air__temperature"

[[output.csv.column]]
header = "Q"
map = "river_gauge__count"
parameter = "river_water__volume_flow_rate"

[[output.csv.column]]
header = "recharge"
map = "subbasin_location__count"
parameter = "soil_water_sat-zone_top__net_recharge_volume_flux"
reducer = "mean"

[API]
variables = [
  "river_water__volume_flow_rate",
   "soil_water_unsat-zone__depth",
   "soil_water__transpiration_volume_flux",
   "soil_layer~1_water_unsat-zone__depth",
   "soil_layer~2_water_unsat-zone__depth",
   "soil_layer~3_water_unsat-zone__depth",
   "soil_layer~4_water_unsat-zone__depth",
]