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

Throw error for insufficient initial heating #93

Merged
merged 4 commits into from
Jul 13, 2024

Conversation

jwreep
Copy link
Member

@jwreep jwreep commented Jul 10, 2024

Fixes #32

Throws an error for small initial heating values, such that it is approaching the range where ebtel++ would fail to produce equilibrium conditions. The threshold is slightly higher than the actual point of failure, but the initial temperatures produced right at the threshold are on the order of a few x $10^{4}$ K, which is perhaps a questionable initial condition for ebtel++ anyway.

@wtbarnes
Copy link
Member

Thanks for taking this on! This looks good to me. Would you mind including a test for this as well?

@jwreep
Copy link
Member Author

jwreep commented Jul 12, 2024

Right, good point. I've added a test to make sure it returns a runtime error. I had to make an edit to util because it assumes that the data file is output when ebtel++ is run, which is not the case if it crashes.

examples/util.py Outdated Show resolved Hide resolved
@wtbarnes wtbarnes merged commit 824f76d into rice-solar-physics:main Jul 13, 2024
2 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.

Warn or error if insufficient initial heating
2 participants