major bugfix in nebular.h, Omega has to be calculated in
ion::transition_levels()