Reformat the outputs from pydisort into an "Atmosphere" Python object / class constructed around the optical properties of the atmosphere (more specifically the eigenpairs).
A Python class would enable easier data queries, e.g. getting the optical depths tau_arr. In addition, this would enable the addition of Fourier modes of the intensity and the changing of sources without having to solve the entire problem again (at the very least the eigendecomposition would not need to be repeated).