*shorten lines for gfortran
*fixed issue in one of the scattering calls
*added support for inelastic scattering
*for tables, the inelastic kernels are stored in a 2D format: logT and logeta