Would be worth benchmarking against https://github.com/johnson-c/ColRadPy and https://github.com/fsciortino/Aurora
Would be worth benchmarking against https://github.com/johnson-c/ColRadPy and https://github.com/fsciortino/Aurora