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

Bad solar flux data corrupts MSIS #43

Closed
greglucas opened this issue Feb 21, 2024 · 1 comment
Closed

Bad solar flux data corrupts MSIS #43

greglucas opened this issue Feb 21, 2024 · 1 comment

Comments

@greglucas
Copy link
Member

Plotting with bad solar flux values causes all subsequent values to produce "nan" output (not just that single day).

>>> from pymsis import msis
>>> msis.run("2005-09-10", 0, 0, 150)
array([[          nan,           nan,           nan,           nan,
                  nan,           nan,           nan,           nan,
        3.0979247e-17,           nan, 9.8023016e+02]], dtype=float32)

Specifically, the date 2005/9/10 causes issues if you try to plot over that day.

  • Should we put limits on valid/realistic F10.7 values?
  • Should we automatically interpolate over those values if there is a radio burst or something else going on that day?

Related to: #28

@greglucas
Copy link
Member Author

This was fixed by #46

UserWarning: There is data that was either interpolated or predicted (not observed), use at your own risk.
  warnings.warn(
array([[1.5667554e-09, 2.2181889e+16, 2.3551942e+15, 1.5285736e+16,
        1.6321905e+13, 1.1822177e+12, 4.9340983e+13, 1.2026543e+13,
        3.0979247e-17, 4.1002158e+12, 6.0943359e+02]], dtype=float32)

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

No branches or pull requests

1 participant