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

Bug Fix: Diagnostics and Second Radial Derivatives when Using Finite-Difference Mode #566

Merged
merged 4 commits into from
Jul 16, 2024

Conversation

feathern
Copy link
Contributor

This fixes a bug related to finite-difference mode. When diagnostic output quantities are requested that require a second radial derivative, that derivative is calculated at output time. The subroutine that handles the derivative calculations at output time was lacking any logic that would call the appropriate finite-difference derivative routine. Instead, a Chebyshev derivative routine was called, which caused a crash because the Chebyshev grid object had not been initialized. This PR fixes this bug by adding appropriate finite-difference logic to Diagnostics_Second_Derivatives.F90.

@feathern
Copy link
Contributor Author

I forgot to update the dependencies file originally. I just did, and now the tests pass.

@tukss
Copy link
Contributor

tukss commented Jul 2, 2024

@rpvin Maybe that fixes your crashes?

Copy link
Contributor

@tukss tukss left a comment

Choose a reason for hiding this comment

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

Looks good to me. Thanks for fixing that.

@feathern feathern merged commit 65a1ca1 into geodynamics:main Jul 16, 2024
7 checks passed
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.

2 participants