[ ] Zhdanov closure: need this to ensure our equations are properly valid for multiple species. We currently have something Braginskii-like which is only valid in the trace approximation.
[ ] 2D velocity fluid neutral model: Full momentum equation for neutrals in cartesian directions. This would likely much improve our fluid neutral model transport, and may have an additional benefit for molecules. See PR https://github.com/bendudson/hermes-3/pull/105
[ ] Molecular reactions: These are very important for detachment. Not started yet.
[ ] Support for more impurities: Need to implement Argon, Nitrogen and Carbon ADAS reactions.
Development
[ ] Balance tests: We need to have CI tests that check for energy and particle balance in a 2D case.
Input/Output
[ ] Flow diagnostics: It would be really useful to have all of the flow operators output diagnostics so that we can post-process fluxes. Some of these are already done.
[ ] LaTeX output: It would be good to have a way to check which equations the code was solving. One way would be to output LaTeX equations.
The terms included may depend on the quantities available in the state, so this should be taken into account. This would help catch errors caused by putting components in the wrong order.
Core physics
Development
Input/Output