Skip to content

Conversation

@JoerivanEngelen
Copy link
Contributor

@JoerivanEngelen JoerivanEngelen commented Jan 10, 2025

Fixes #1376

Description

This aligns the way we find active cells, consistently with MODFLOW6:
idomain > 0 is now the active domain, instead of idomain == 1. I updated some model fixtures to have idomain = 2 to see if stuff like masking works now for idomain == 2. Furthermore, I tried my best to update all docstrings.

Checklist

  • Links to correct issue
  • Update changelog, if changes affect users
  • PR title starts with Issue #nr, e.g. Issue #737
  • Unit tests were added
  • If feature added: Added/extended example

imod/mf6/disv.py Outdated
bottom: array of floats (xu.UgridDataArray)
is the bottom elevation for each cell.
idomain: array of integers (xu.UgridDataArray)
Indicates the existence status of a cell. If the idomain value for a
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Don't we use active instead existence throughout our code?

Edit:
Ah nvm. In model.py you added clear definitions:

1 : active
=0 : inactive
<0 : vertical passthrough

@sonarqubecloud
Copy link

@JoerivanEngelen JoerivanEngelen added this pull request to the merge queue Jan 13, 2025
Merged via the queue into master with commit ae03531 Jan 13, 2025
5 checks passed
@JoerivanEngelen JoerivanEngelen deleted the issue_#1376_active_idomain_larger_than_zero branch January 13, 2025 11:34
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Find active cells with idomain > 0 instead of idomain == 1

3 participants