Huite

I think Sunny set it up for imod-python: https://deltares.github.io/imod-python/developing.html#how-to-release-a-new-version

For testing, do these give exactly the same interpolation? ```julia itp1 = SmoothedLinearInterpolation([-1.0, 0.0, 1.0], [0.0, 0.0, 10.0]; extrapolate = true) itp2 = SmoothedLinearInterpolation([-1.0, 0.0, 0.1], [0.0, 0.0, 0.1]; extrapolate...

Bad copy-pasting on my part yes, I've fixed it now.

A pragmatic solution for now: if you have the xugrid representation of the network, you can convert it with some associated data directly to a geodataframe in a single `to_geodataframe()`...

`y/n` license prompts are a horrible idea -- totally symbolic, nobody reads them, nobody presses `n` even if they were to read it. Where does this requirement come from?

This is a cute idea. The transformation can also be done ad hoc, wherever it's relevant, which has the benefit of things looking slightly less exotic. I think the main...

For anyone needing this on a recent Snakemake version: ```python def read_snakemake_rule(path, name: str) -> "snakemake.rules.Rule": """ Parameters ---------- path: str, pathlib.Path The path to the snakefile. name: str Name...

Sure; I think we'd want a TypeError along the lines of `TypeError: axis must be integer, not {type(axis).__name__}`?

Just had a thought: obviously if the initial stage is higher than the initial head estimate (for a steady-state computation), a linearized formulation is ill-posed since there's only fixed fluxes...

For anyone looking to do something similar in the mean time: chunking manually followed by `np.prange` works reasonably well. Instead of a list, use a dynamically growing array (like a...