Deltares / Deltares/imod-python
Select horizontal cross-section
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 41
- Forks
- 12
- Avg merge
- 21h 8m
- Merged PRs (30d)
- 1
Description
In GitLab by @Huite on Jun 8, 2020, 23:25
Horizontal cross-sections are pretty easy to get in a voxel model, not so much in a layered model.
Xarray can be used with .sel(z=-10.0, method="nearest"), but these don't take cell boundaries into account.
The LayerRegridder can be used by using a thin fictitious thickness, but it's not exactly convenient.
A plane defined by z-coordinates, matching in y and x, seems like a decent generalized approach.
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 comparing the issue's xarray .sel(z=-10.0, method="nearest") approach with the LayerRegridder workaround. Clarify the intended plane defined by z-coordinates and how matching in y and x should work for layered models. Done means a documented, convenient generalized approach for selecting horizontal cross-sections that accounts for cell boundaries.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- data
- Issue type
- Feature
- Difficulty
- 5/5
- Estimated time
- Over a week
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100