Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Issue #1395 Change default value for missing layers for angle1 #1396

Merged

Conversation

JoerivanEngelen
Copy link
Contributor

@JoerivanEngelen JoerivanEngelen commented Jan 29, 2025

Fixes #1395

Description

Use proper value for missing layers in projectfile.

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

@Huite
Copy link
Contributor

Huite commented Jan 29, 2025

See my comment in the issue: entering a value for ANGLE will likely slow MF6 down. Better to apply the trigonometry yourself.

imod/mf6/npf.py Outdated
@@ -510,7 +510,7 @@ def from_imod5_data(
angle1 = imod5_data["ani"]["angle"]
angle1 = 90.0 - angle1
angle1 = xr.where(angle1 < 0, 360.0 + angle1, angle1)
angle1 = fill_missing_layers(angle1, target_grid, 0)
angle1 = fill_missing_layers(angle1, target_grid, 90)
Copy link
Collaborator

Choose a reason for hiding this comment

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

Should be 90.0 to be consistent with the way other angles are written

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Nice catch, fixed!

@JoerivanEngelen
Copy link
Contributor Author

See my comment in the issue: entering a value for ANGLE will likely slow MF6 down. Better to apply the trigonometry yourself.

Created a separate issue for this #1397

@JoerivanEngelen JoerivanEngelen added this pull request to the merge queue Jan 31, 2025
Merged via the queue into master with commit 1893847 Jan 31, 2025
7 checks passed
@JoerivanEngelen JoerivanEngelen deleted the issue_#1395_fill_layers_with_better_missing_value branch January 31, 2025 13:42
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.

Incorrect value used to fill missing layers in angle1 array in NPF package
3 participants