Use GTH algorithm for all species with 0 source terms
This is a preliminary bug fix.
When the total sources are 0, atom_levelN computes negative populations. The fix here uses the GTH algorithm, which guarantees positivity. It is preliminary in the sense that the algorithm is designed to only work with 0 source terms, i.e., it cannot handle finite values.