Added printouts of level opacitites, and level populations to final printout if the print switch is set to 2.
Added a column to the printout of the file xout_detail.fits for the upp level index of each line.
Repaired and streamlined the printing of the file xout_detail.fits.
Added more informative statement when the code stops because the rate matrix overflows (ipmat too large).
Added rate type 42: Auger decay
Fixed arithmetic error which affected recombination rate calculation when kT E
Added new data type (85) for photoionization resonances below threshold, along with new subroutine to calculate cross subsection (PEXS.f).
Streamlined the photoionization rate calculation (phint53)
More accurate treatment of line damping, Voigt profiles.