Above the star formation threshold, Arepo uses the HI fraction
on the eEOS, while we want it for the cold fraction.
Therefore for gas with nH > 0.1 cm^-3, use the Rahmati fitting formula.

get_reproc_rhoHI is now the correct neutral hydrogen density.