Version v2.2.1bn22 (September. 2014)

Photoionization integration routine now uses thresholds calculated on energy bin boundaries. This allows for better evaluation of the Milne integral, though it may affect photoionization rates in the case of very coarse energy bins. Removed redundant subroutine phint53new.f Added code to print ion column densities as part of lpri=2 output Increased buffer size in subroutine fstepr2 which writes to xo01_detal2.fits such that table of lines is not artificially truncated. Add fine structure to He-like ion level and radiative decay data.