File lake_storage.stat

[alicelg99]

This file is present is the folder init_conds in the Input folder but also in the Output folder. For the one in the folder init_conds we have positive values for the volume which are :

Lake volume status (for analysis or model re-start)
Subbasin	reservoir_size_class	volume[m^3]
15	1	       0.01
15	2	       0.08
15	3	       0.74
15	4	       2.91
15	5	      10.23
69	1	       0.01
69	2	       0.08
69	3	       0.74
69	4	       2.91
69	5	       0.00

But for each case the file with the same name in the Output folder has zero values for the volume column such as the file as follows :

Lake volume status (for analysis or model re-start)
Subbasin	reservoir_size_class	volume[m^3]
15	1	       0.00
15	2	       0.00
15	3	       0.00
15	4	       0.00
15	5	       0.00
69	1	       0.00
69	2	       0.00
69	3	       0.00
69	4	       0.00
69	5	       0.00

Also we note that the file lake_storage.stat_start in Output folder, organized in the same way, does not have zero values but it has the same values as the file lake_storage.out in the init_conds folder.

Note: The files lake_storage.stat_start and lake_storage.out are not present in the Output folder for the cases 5, 6 and 14. And they are empty for the case 15 (and 7 but all the files are empty cf the issue case7).

[TillF]

This is not the case when dosediment=.false.
I assume it is about silting up the small reservoirs. Did you find any way to check for the deposition in the small reservoir to see if this is the case?

[alicelg99]

What is the case corresponding to the small reservoir ?

[TillF]

I suggest you check https://tillf.github.io/WASA-SED/#4-3-output-of-the-reservoir-module and look for the output files for the small reservoirs ("lake_*.out") to see if any of them shows that the storage volume decreases. I just didn't check which file(s) is the best to check this, please try.

[TillF]

I confirm the problem is that the problem is that the small reservoirs completely silt up, thus, cannot retain any more water. This is now visible in lake_maxstorcap.out (needs updateing the executable, it was unavailable before).
However, I cannot tell yet if we have
a) an excessive unrealistiv erosion rate in our fake cases that just fills up the small reservoirs
or
b) the rates are ok, but there is a bug in the calculation of sedimentation.
You could try to check a), I will look into this after my holiday.