subsurface_water__to_river_volume_flow_rate not recognized as an [output.netcdf_grid.variables]

[AquaPore]

Wflow version checks

• I have checked that this issue has not already been reported.

• I have checked that this bug exists on the latest version of Wflow.

Reproducible Example

Dear Wflow developer,

Based on Issue_889 which was partially fixed.

I thank you for fixing the error of not getting the output land_surface_water__to_river_volume_flow_rate, which after testing seems to be working.

Nevertheless the computation of the subsurface which flows in the river is not recognised
parameters_routing

TOML FILE

[output.netcdf_grid.variables] 
subsurface_water__to_river_volume_flow_rate = "subsurface_water__to_river_volume_flow_rate" #  [m3 / day] Part of subsurface flow that flows to the river"

Please note that I am running the latest Wflow model.

With kindness,
Joseph A.P. Pollacco

Current behaviour

Error message

 Error: Wflow simulation failed
│   exception =
│    Couldn't obtain a field from this model specified by 'subsurface_water__to_river_volume_flow_rate'.
│    Stacktrace:
│      [1] error(s::String)
│        @ Base .\error.jl:44
│      [2] get_field_in_model(model::@NamedTuple{land::Main.Wflow.LandHydrologySBM{Main.Wflow.NoDemandModel, Main.Wflow.NoAllocationLandModel}, routing::Main.Wflow.Routing{Main.Wflow.OverlandFlowModel{Main.Wflow.KinematicWave, Main.Wflow.LandFlowBC, Main.Wflow.ManningFlowParameters, Main.Wflow.OverLandFlowVariables}, Main.Wflow.RiverFlowModel{Main.Wflow.KinematicWave, Nothing, Main.Wflow.RiverFlowParameters, Main.Wflow.RiverFlowVariables, Main.Wflow.NoAllocationRiverModel}, Main.Wflow.GroundwaterFlowModel{Main.Wflow.SubsurfaceFlowBC{Main.Wflow.RechargeModel, Main.Wflow.GwfRiverModel, Nothing, Nothing}}}, mass_balance::Main.Wflow.NoMassBalance}, name::String; check_allow_dynamic_input::Bool)
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\standard_name\standard_name_utils.jl:188
│      [3] get_field_in_model
│        @ D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\standard_name\standard_name_utils.jl:172 [inlined]
│      [4] out_map(output_names_dict::Dict{String, String}, modelmap::@NamedTuple{land::Main.Wflow.LandHydrologySBM{Main.Wflow.NoDemandModel, Main.Wflow.NoAllocationLandModel}, routing::Main.Wflow.Routing{Main.Wflow.OverlandFlowModel{Main.Wflow.KinematicWave, Main.Wflow.LandFlowBC, Main.Wflow.ManningFlowParameters, Main.Wflow.OverLandFlowVariables}, Main.Wflow.RiverFlowModel{Main.Wflow.KinematicWave, Nothing, Main.Wflow.RiverFlowParameters, Main.Wflow.RiverFlowVariables, Main.Wflow.NoAllocationRiverModel}, Main.Wflow.GroundwaterFlowModel{Main.Wflow.SubsurfaceFlowBC{Main.Wflow.RechargeModel, Main.Wflow.GwfRiverModel, Nothing, Nothing}}}, mass_balance::Main.Wflow.NoMassBalance})
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\io.jl:673
│      [5] Main.Wflow.Writer(config::Main.Wflow.Config, modelmap::@NamedTuple{land::Main.Wflow.LandHydrologySBM{Main.Wflow.NoDemandModel, Main.Wflow.NoAllocationLandModel}, routing::Main.Wflow.Routing{Main.Wflow.OverlandFlowModel{Main.Wflow.KinematicWave, Main.Wflow.LandFlowBC, Main.Wflow.ManningFlowParameters, Main.Wflow.OverLandFlowVariables}, Main.Wflow.RiverFlowModel{Main.Wflow.KinematicWave, Nothing, Main.Wflow.RiverFlowParameters, Main.Wflow.RiverFlowVariables, Main.Wflow.NoAllocationRiverModel}, Main.Wflow.GroundwaterFlowModel{Main.Wflow.SubsurfaceFlowBC{Main.Wflow.RechargeModel, Main.Wflow.GwfRiverModel, Nothing, Nothing}}}, mass_balance::Main.Wflow.NoMassBalance}, domain::Main.Wflow.Domain, nc_static::NCDatasets.NCDataset{Nothing, Missing}; extra_dim::@NamedTuple{name::String, value::Array{Float64, 1}})
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\io.jl:715
│      [6] Main.Wflow.Model(config::Main.Wflow.Config, type::Main.Wflow.SbmGwfModel)
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\sbm_gwf_model.jl:37
│      [7] Main.Wflow.Model(config::Main.Wflow.Config)
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\Wflow.jl:186
│      [8] run(config::Main.Wflow.Config)
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\Wflow.jl:305
│      [9] (::Main.Wflow.var"#751#752"{IOStream, Main.Wflow.Config})()
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\Wflow.jl:286
│     [10] with_logstate(f::Main.Wflow.var"#751#752"{IOStream, Main.Wflow.Config}, logstate::Base.CoreLogging.LogState)
│        @ Base.CoreLogging .\logging\logging.jl:542
│     [11] with_logger
│        @ .\logging\logging.jl:653 [inlined]
│     [12] run(tomlpath::String; silent::Nothing)
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\Wflow.jl:282
│     [13] run(tomlpath::String)
│        @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\Wflow.jl:275
│     [14] top-level scope
│        @ REPL[9]:1
│     [15] eval(m::Module, e::Any)
│        @ Core .\boot.jl:489
│     [16] eval
│        @ .\Base_compiler.jl:146 [inlined]
│     [17] repleval(m::Module, code::Expr, ::String)
│        @ VSCodeServer c:\Users\jpollacco.local\.vscode\extensions\julialang.language-julia-1.219.2\scripts\packages\VSCodeServer\src\repl.jl:278
│     [18] #evalrepl##2
│        @ c:\Users\jpollacco.local\.vscode\extensions\julialang.language-julia-1.219.2\scripts\packages\VSCodeServer\src\repl.jl:241 [inlined]
│     [19] with_logstate(f::VSCodeServer.var"#evalrepl##2#evalrepl##3"{Module, Expr, REPL.LineEditREPL, REPL.LineEdit.Prompt}, logstate::Base.CoreLogging.LogState)
│        @ Base.CoreLogging .\logging\logging.jl:542
│     [20] with_logger
│        @ .\logging\logging.jl:653 [inlined]
│     [21] (::VSCodeServer.var"#evalrepl##0#evalrepl##1"{Module, Expr, REPL.LineEditREPL, REPL.LineEdit.Prompt})()
│        @ VSCodeServer c:\Users\jpollacco.local\.vscode\extensions\julialang.language-julia-1.219.2\scripts\packages\VSCodeServer\src\repl.jl:242
│     [22] (::VSCodeServer.var"#start_eval_backend##0#start_eval_backend##1")()
│        @ VSCodeServer c:\Users\jpollacco.local\.vscode\extensions\julialang.language-julia-1.219.2\scripts\packages\VSCodeServer\src\eval.jl:34
└ @ Main.Wflow D:\JOE\MAIN\MODELS\WFLOW\Wflow.jl\Wflow\src\Wflow.jl:293

Desired behaviour

Obtaining output subsurface_water__to_river_volume_flow_rate

Additional Context

I also tried:

routing.subsurface_flow.variables.to_river_average= "subsurface_water__to_river_volume_flow_rate"

which also failed

caused by: FieldError: type Main.Wflow.GroundwaterFlowVariables has no field `to_river_average`, available fields: `n`, `head`, `conductance`, `storage`, `q_net`, `q_net_cumulative`, `q_net_average`, `q_net_bnds`, `q_in_cumulative`, `q_in_average`, `q_cumulative`, `q_average`, `exfiltwater_cumulative`, `exfiltwater_average`

Many thanks for any help you may provide,
Joseph A.P. Pollacco

[Timh37]

I am also getting this error when trying to run the Moselle test case with a fresh install. I followed the exact steps on the Getting Started page. Help would be appreciated! Error log below:

julia> using Wflow

julia> toml_path = "sbm_config.toml"
"sbm_config.toml"

julia> Wflow.run(toml_path)
┌ Warning: 'wflow_version' is not recognized as a valid field of the [input] section in the TOML, this will be ignored.
└ @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\config_init.jl:81
[ Info: Wflow version v1.0.4
[ Info: Initialize model variables for model type sbm.
[ Info: Cyclic parameters are provided by data\input\staticmaps-moselle.nc.
[ Info: Forcing parameters are provided by data\input\forcing-moselle.nc.
[ Info: Set atmosphere_water__precipitation_volume_flux using netCDF variable precip as forcing parameter.
[ Info: Set atmosphere_air__temperature using netCDF variable temp as forcing parameter.
[ Info: Set land_surface_water__potential_evaporation_volume_flux using netCDF variable pet as forcing parameter.
[ Info: Set vegetation__leaf_area_index using netCDF variable LAI as cyclic parameter, with 12 timesteps.
┌ Info: General model settings
│   snow = true
│   gravitational_snow_transport = true
│   glacier = false
│   reservoirs = true
│   pits = false
└   water_demand = false
[ Info: Set subbasin_location__count using netCDF variable wflow_subcatch.
[ Info: Set basin__local_drain_direction using netCDF variable wflow_ldd.
[ Info: Set land_surface__slope using netCDF variable Slope.
[ Info: Set river_location__mask using netCDF variable wflow_river.
[ Info: Set land_water_covered__area_fraction using netCDF variable WaterFrac.
[ Info: Set reservoir_location__count using netCDF variable wflow_reservoirlocs.
[ Info: Set reservoir_area__count using netCDF variable wflow_reservoirareas.
[ Info: Set river__length using netCDF variable wflow_riverlength.
[ Info: Set river__width using netCDF variable wflow_riverwidth.
[ Info: Set river__slope using netCDF variable RiverSlope.
[ Info: Set vegetation_root__depth using netCDF variable RootingDepth.
[ Info: Set vegetation__crop_factor using default value 1.0.
[ Info: Set vegetation__specific_leaf_storage using netCDF variable Sl.
[ Info: Set vegetation_wood_water__storage_capacity using netCDF variable Swood.
[ Info: Set vegetation_canopy__light_extinction_coefficient using netCDF variable Kext.
[ Info: Set vegetation_canopy_water__mean_evaporation_to_mean_precipitation_ratio using netCDF variable EoverR.
[ Info: Set snowpack__degree_day_coefficient using netCDF variable Cfmax.
[ Info: Set atmosphere_air__snowfall_temperature_threshold using netCDF variable TT.
[ Info: Set atmosphere_air__snowfall_temperature_interval using netCDF variable TTI.
[ Info: Set snowpack__melting_temperature_threshold using netCDF variable TTM.
[ Info: Set snowpack__liquid_water_holding_capacity using netCDF variable WHC.
[ Info: Using 3 soil layers with the following thickness: [100, 300, 800]
[ Info: Set soil_surface_temperature__weight_coefficient using default value 0.1125.
[ Info: Set soil_surface_water__infiltration_reduction_parameter using netCDF variable cf_soil.
[ Info: Set soil_water__saturated_volume_fraction using netCDF variable thetaS.
[ Info: Set soil_water__residual_volume_fraction using netCDF variable thetaR.
[ Info: Set soil_surface_water__vertical_saturated_hydraulic_conductivity using netCDF variable KsatVer.
[ Info: Set soil_water__vertical_saturated_hydraulic_conductivity_scale_parameter using netCDF variable f.
[ Info: Set soil_water__air_entry_pressure_head using default value -10.0.
[ Info: Set vegetation_root__feddes_critical_pressure_head_h1 using default value 0.0.
[ Info: Set vegetation_root__feddes_critical_pressure_head_h2 using default value -100.0.
[ Info: Set vegetation_root__feddes_critical_pressure_head_h3_high using default value -400.0.
[ Info: Set vegetation_root__feddes_critical_pressure_head_h3_low using default value -1000.0.
[ Info: Set vegetation_root__feddes_critical_pressure_head_h4 using default value -16000.0.
[ Info: Set vegetation_root__feddes_critical_pressure_head_h1_reduction_coefficient using default value 1.0.
[ Info: Set soil__thickness using netCDF variable SoilThickness.
[ Info: Set compacted_soil_surface_water__infiltration_capacity using netCDF variable InfiltCapPath.
[ Info: Set soil_water_saturated_zone_bottom__max_leakage_volume_flux using netCDF variable MaxLeakage.
[ Info: Set soil_layer_water__brooks_corey_exponent using netCDF variable c.
[ Info: Set soil_layer_water__vertical_saturated_hydraulic_conductivity_factor using default value 1.0.
[ Info: Set compacted_soil__area_fraction using netCDF variable PathFrac.
[ Info: Set soil_wet_root__sigmoid_function_shape_parameter using netCDF variable rootdistpar.
[ Info: Set soil_water_saturated_zone_top__capillary_rise_max_water_table_depth using default value 2000.0.
[ Info: Set soil_water_saturated_zone_top__capillary_rise_averianov_exponent using default value 2.0.
┌ Info: Set subsurface_water__horizontal_to_vertical_saturated_hydraulic_conductivity_ratio using netCDF variable
└ KsatHorFrac.
[ Info: Set land_surface_water_flow__manning_n_parameter using netCDF variable N.
[ Info: Kinematic wave approach is used for land flow, adaptive timestepping = true.
[ Info: Set reservoir_surface__area using netCDF variable reservoir_area.
[ Info: Set reservoir_water_surface__initial_elevation using netCDF variable waterlevel_reservoir.
[ Info: Set reservoir_water__storage_curve_type_count using netCDF variable storfunc.
[ Info: Set reservoir_water__rating_curve_type_count using netCDF variable outflowfunc.
[ Info: Set reservoir_lower_location__count using default value 0.
[ Info: Set reservoir_location__count using netCDF variable wflow_reservoirlocs.
[ Info: Read 2 reservoir locations.
[ Info: Set reservoir_water_demand__required_downstream_volume_flow_rate using netCDF variable ResDemand.
[ Info: Set reservoir_water_release_below_spillway__max_volume_flow_rate using netCDF variable ResMaxRelease.
[ Info: Set reservoir_water__max_volume using netCDF variable ResMaxVolume.
[ Info: Set reservoir_water__target_full_volume_fraction using netCDF variable ResTargetFullFrac.
[ Info: Set reservoir_water__target_min_volume_fraction using netCDF variable ResTargetMinFrac.
[ Info: Set reservoir_water__outgoing_observed_volume_flow_rate using default value NaN.
[ Info: Set reservoir_water__external_inflow_volume_flow_rate using default value 0.0.
[ Info: Kinematic wave approach is used for river flow, adaptive timestepping = true.
[ Info: Set river_water_flow__manning_n_parameter using netCDF variable N_River.
[ Info: Set river_bank_water__depth using netCDF variable RiverDepth.
[ Info: Set river_water__external_inflow_volume_flow_rate using default value 0.0.
[ Info: Create an output netCDF file data\output\output_moselle.nc for grid data, using compression level 0.
┌ Error: State variable subsurface_water__volume_flow_rate not provided, please ensure that it is set correctly in
│ the model configuration
└ @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\states.jl:202
ERROR: ArgumentError: Not all required states are provided, these states are missing: `("subsurface_water__volume_flow_rate",)`
Stacktrace:
  [1] check_states(config::Wflow.Config)
    @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\states.jl:212
  [2] Wflow.Writer(config::Wflow.Config, modelmap::@NamedTuple{…}, domain::Wflow.Domain, nc_static::NCDatasets.NCDataset{…}; extra_dim::@NamedTuple{…})
    @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\io.jl:764
  [3] Wflow.Model(config::Wflow.Config, type::Wflow.SbmModel)
    @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\sbm_model.jl:32
  [4] Wflow.Model(config::Wflow.Config)
    @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\Wflow.jl:179
  [5] run(config::Wflow.Config)
    @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\Wflow.jl:280
  [6] (::Wflow.var"#669#670"{IOStream, Wflow.Config})()
    @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\Wflow.jl:261
  [7] with_logstate(f::Wflow.var"#669#670"{IOStream, Wflow.Config}, logstate::Base.CoreLogging.LogState)
    @ Base.CoreLogging .\logging\logging.jl:542
  [8] with_logger
    @ .\logging\logging.jl:653 [inlined]
  [9] run(tomlpath::String; silent::Nothing)
    @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\Wflow.jl:257
 [10] run(tomlpath::String)
    @ Wflow C:\Users\TST-WVL-hermanst\.julia\packages\Wflow\G0tmj\src\Wflow.jl:250
 [11] top-level scope
    @ REPL[38]:1
Some type information was truncated. Use `show(err)` to see complete types.

[visr]

@Timh37 it looks like you are running a TOML from the development branch (main/master) with the latest stable release v1.0.4. I think that is our bad, because looking at the stable docs https://deltares.github.io/Wflow.jl/stable/getting_started/download_example_models.html and searching for master it indeed gets them from there. If you replace that with v1.0.4, e.g. https://raw.githubusercontent.com/Deltares/Wflow.jl/v1.0.4/Wflow/test/sbm_config.toml you should get the right one. Does that help?

We should fix those links on the stable docs and rebuild.

[visr]

The stable docs are now fixed, and based on the v1.0.4 latest release + doc fixes on the release/v1.0 branch.

I believe that issue is different from the original by @AquaPore, so leaving this open. Sorry if we are not so quick to respond, it is summer holiday period here. @AquaPore can you clarify if by latest Wflow version, do you mean latest release v1.0.4, or the main branch? You link to the dev (development branch / main) docs, so I assume the latter? Note by the way that it may be good to update your local clone to rename the default branch, see #1021.