Deltares / Deltares/imod-python
[FEATURE] - Function to compute fluxes based on previously computed heads.
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 41
- Forks
- 12
- Avg merge
- 21h 8m
- Merged PRs (30d)
- 1
Description
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_conductancesmethod to NodePropertyFlow or function. - Add function that takes conductances, heads and, optionally, concentrations and buyoancy package, to compute fluxes.
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by reviewing NodePropertyFlow and the prototype code mentioned in the issue. Trace how conductances, heads, concentrations, and the buoyancy package are represented, then define the conductance and flux computation entry points. Done means fluxes can be computed from previously stored heads, with optional concentrations and buoyancy inputs, without rerunning the full simulation.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- backend
- Issue type
- Feature
- Difficulty
- 5/5
- Estimated time
- Over a week
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100