Is your feature request related to a problem? Please describe.
Request from @jdelsman, he has some prototype code for this.
It happens that large-scale simulations are run, but only states were stored. Later there is a request for fluxes, but recomputing everything, just to store fluxes, is too costly.
@Huite tried this before with iMOD-WQ but ran into the issue that that kernel saves output as float32, which resulted in grave errors when computing fluxes. This isn't the problem for MODFLOW6, as it stores everything float64.
Rough outline:
- Add
compute_conductances method to NodePropertyFlow or function.
- Add function that takes conductances, heads and, optionally, concentrations and buyoancy package, to compute fluxes.