Added a convergence test for the global ionized fraction in order
to stop the iteration loop over all the sources. This should lead
to faster convergence (3 - 5 times?) with minimal loss of accuracy.