sunkit-spex
sunkit-spex copied to clipboard
Compute the thermal emission data in Python
For the thermal plasma emission function, we currently rely on IDL .sav files generated externally.
I hope this issue starts a conversation around generating the data in Python.
I made a github gist which shows how to save out all of the emission lines on a temperature grid. This uses the Fiasco package made by @wtbarnes. Those lines can get binned into an (E, T) matrix much like the IDL version, and summed with the continuum. This is one part of the problem.
I stopped working on this after the lines, so I haven't touched the continuum.
Things that need to happen:
- [ ] make a continuum emission function analogous to the lines one
- [ ] figure out how to integrate the two
- [ ] decide how we want to bin things
- [ ] ... ?