points_indices but for layers
Get the layer number for a given value of Z, similar to select.points_indices. Input DataArray with the Layer dimension and points with Z-coordinate. I think this would be a useful addition to the current selection capabilities in imod-python.
More elaborate, but also useful, is the capability to add well screens to a model based on geometry and transmissivity values, similar to iMOD MKWELLIPF. Returns a well screen fraction per well per model layer.