Rates and fractions (saturated and decreased ISRF and CR rates)

Similar to Ploeckinger & Schaye (2020) but with the following changes: - turbulence of 6km/s is added in the Cloudy calculations - turbulent pressure of 6 km/s is included in the calculation of the Jeans column density (pressure = max of turbulent and thermal pressure) - ISRF scales with the new Jeans column density (incl. turbulent pressure) but saturates at a column density of 1.3e22 cm-2 (corresponds to 100 pc turbulent Jeans length) - cosmic ray rate scales with the new Jeans column density according to the Lacki et al. (2010) scaling (CR rate proportional to the column density) and saturates at a Milky Way value of 3.63e20 cm-2 (this is the log N_0 = 20.56 value used in PS20) - both the CR rate and the ISRF are scaled down by 2 dex (1 dex more than in the fiducial model from PS20 for a given column density but note that the column densities are higher in the ISM due to the turbulent pressure contribution).