repro: LHM-vechtstromen

[evetion]

This model is the LHM-vechtstromen, send by @HendrikKok. It crashes with code unstable, and when run with saveat=0 it crashes on the waterbalance error (in the save step). I've isolated the basin in question, and uploaded it to repro-models/lhm-vecht which you can get with pixi run s3-download repro-models/lhm-vecht lhm-vecht (with the model ending up at models/lhm-vecht/.

In this model Basin #2189 falls dry, which triggers the balance error(?)

[evetion]

@SouthEndMusic Updated model at repro-models/vechtstromen_v2.

Details

➜ ribasim (7675bed0) ✗ julia --project=core -e 'using Ribasim; Ribasim.main(ARGS)' /Users/evetion/Downloads/vechtstromen_v2/ribasim.toml
┌ Warning: The Ribasim version in the TOML config file does not match the used Ribasim CLI version.
│   config.ribasim_version = "2025.2.0"
│   cli.ribasim_version = "2025.3.0"
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/main.jl:38
┌ Info: Starting a Ribasim simulation.
│   toml_path = "/Users/evetion/Downloads/vechtstromen_v2/ribasim.toml"
│   cli.ribasim_version = "2025.3.0"
│   starttime = 2020-01-01T00:00:00
└   endtime = 2021-01-01T00:00:00
Simulating 100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| Time: 0:03:12
[ Info: Computation time: 2 minutes, 46 seconds, 49 milliseconds
┌ Error: The model exited at model time 2020-01-03T17:44:00.212 with return code Unstable. See https://docs.sciml.ai/DiffEqDocs/stable/basics/solution/#retcodes
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/main.jl:70
➜ ribasim (7675bed0) ✗ julia --project=core -e 'using Ribasim; Ribasim.main(ARGS)' /Users/evetion/Downloads/vechtstromen_v2/ribasim.toml
┌ Warning: The Ribasim version in the TOML config file does not match the used Ribasim CLI version.
│   config.ribasim_version = "2025.2.0"
│   cli.ribasim_version = "2025.3.0"
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/main.jl:38
┌ Info: Starting a Ribasim simulation.
│   toml_path = "/Users/evetion/Downloads/vechtstromen_v2/ribasim.toml"
│   cli.ribasim_version = "2025.3.0"
│   starttime = 2020-01-01T00:00:00
└   endtime = 2021-01-01T00:00:00
┌ Error: Too large water balance error
│   id = Basin #1375
│   balance_error = 0.005510760481520104
│   relative_error = 0.05224829271657602
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1408
│   balance_error = -0.005676367391501675
│   relative_error = -0.6666666666666666
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1442
│   balance_error = -0.0024834107337819607
│   relative_error = -0.013252869482901315
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1445
│   balance_error = 0.0069755461950623335
│   relative_error = 0.0248832852450371
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1446
│   balance_error = -0.0023429950007473888
│   relative_error = -0.017051829868033774
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1469
│   balance_error = 0.0013660104001470552
│   relative_error = 0.028718139471442364
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1475
│   balance_error = -0.01702910217450504
│   relative_error = -0.13333333333333344
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1498
│   balance_error = -0.0058537538724861005
│   relative_error = -0.1377870563674321
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1518
│   balance_error = -0.0022783076151437383
│   relative_error = -0.012772503379585748
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1556
│   balance_error = -0.00141909184787542
│   relative_error = -0.03174603174603178
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1583
│   balance_error = 0.011574467884233882
│   relative_error = 0.31016042780748654
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1605
│   balance_error = 0.0030155701767352605
│   relative_error = 0.013704151551793612
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1615
│   balance_error = 0.004026018925366456
│   relative_error = 0.036103150359947024
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1621
│   balance_error = -0.003587015237797527
│   relative_error = -0.03964133215830546
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1634
│   balance_error = 0.005731800666809311
│   relative_error = 0.08776844070961723
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1635
│   balance_error = -0.017156598707712606
│   relative_error = -0.10926936044766913
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1652
│   balance_error = -0.00212863777181313
│   relative_error = -0.05128205128205132
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1656
│   balance_error = 0.0015438167173175976
│   relative_error = 0.033428357089272394
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1659
│   balance_error = 0.0029822070612430546
│   relative_error = 0.3024050468854004
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1672
│   balance_error = -0.0014700211695643085
│   relative_error = -0.05987652143235144
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1687
│   balance_error = 0.006247110510975557
│   relative_error = 0.4492531323907872
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1701
│   balance_error = -0.011627822411717448
│   relative_error = -0.0526607303965793
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1714
│   balance_error = -0.0021869430404084927
│   relative_error = -0.015293023898116438
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1755
│   balance_error = 0.006319393385070232
│   relative_error = 0.09970264124540855
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1760
│   balance_error = 0.0019527127294882117
│   relative_error = 0.03896778495121842
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1764
│   balance_error = -0.001510489573220225
│   relative_error = -0.02026196746027687
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1793
│   balance_error = -0.0017493589596210037
│   relative_error = -0.14306859340990627
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1807
│   balance_error = -0.005598553645458547
│   relative_error = -0.11612795202398984
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1824
│   balance_error = -0.01135273478300335
│   relative_error = -0.058823529411764705
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1873
│   balance_error = 0.0011512798365454226
│   relative_error = 0.016483339327124712
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1879
│   balance_error = 0.023751495327203193
│   relative_error = 0.04059512714006568
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1928
│   balance_error = -0.001374745227629308
│   relative_error = -0.052320675105485076
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #1955
│   balance_error = -0.0011777295870254523
│   relative_error = -0.34365711215207634
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #2075
│   balance_error = -0.0016188313604118612
│   relative_error = -0.05828971582953494
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #2127
│   balance_error = 0.0015964783288598457
│   relative_error = 0.024691358024691353
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #2147
│   balance_error = 0.001751691499721219
│   relative_error = 0.046375110067508045
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #2163
│   balance_error = -0.010643188859065456
│   relative_error = -0.015251652262328155
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #2166
│   balance_error = -0.0010490747351969698
│   relative_error = -0.04780700369446484
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #2180
│   balance_error = 0.004506756214920312
│   relative_error = 1.3166141757824001
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #2188
│   balance_error = -0.008514551087252513
│   relative_error = -0.03337969401947149
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Error: Too large water balance error
│   id = Basin #2256
│   balance_error = -0.0010953360218719978
│   relative_error = -0.03447956297514749
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/callback.jl:446
┌ Warning: Simulation crashed or interrupted.
└ @ Ribasim /Users/evetion/code/Ribasim/core/src/main.jl:51
ERROR: Too large water balance error(s) detected at t = 2020-01-03T17:44:00.212
Stacktrace:
  [1] error(s::String)
    @ Base ./error.jl:35
  [2] check_water_balance_error!(saved_flow::Ribasim.SavedFlow, integrator::OrdinaryDiffEqCore.ODEIntegrator{OrdinaryDiffEqBDF.QNDF{5, 0, ADTypes.AutoSparse{ADTypes.AutoForwardDiff{nothing, ForwardDiff.Tag{DiffEqBase.OrdinaryDiffEqTag, Float64}}, ADType

[visr]

repro-models/vechtstromen_v2 now runs in 30 seconds.