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

can we get rid of rotation on import of imod5 models? #1015

Open
luitjansl opened this issue May 7, 2024 · 0 comments
Open

can we get rid of rotation on import of imod5 models? #1015

luitjansl opened this issue May 7, 2024 · 0 comments

Comments

@luitjansl
Copy link
Contributor

luitjansl commented May 7, 2024

Currently when we import imod5 models we always get an ANGLE defined in the npf package, even when the major and minor axes of the conductivity ellipsoid are aligned with the grid axes. But when that's the case, we can swap the names of K and K22 and then we don't need the angle anymore. This simplifies calculations on the modflow side.

This issue comes from a review comment by @Huite:

Just a note, maybe you're aware already: iMOD5 chooses the north (or 90 degrees with x-axis) as the principal direction of the horizontal conductivity. MODFLOW 6 chooses the x-axis (0 degrees) as the principal direction (as you'd expect from the unit circle). If I recall correctly, in the iMOD5 conversion, this results in there always being an 90 degree angle defined even if there's no anisotropy. However, this makes the formulation step of MODFLOW 6 more expensive, because it requires trigonometry evaluations. It's worth considering, I think, to consistently transform the iMOD5 conductivity to an x-axis aligned one for k11.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
Status: 📯 New
Development

No branches or pull requests

1 participant