Mass-weighted n_H–T diagrams from our simulation suite. Columns represent different CR transport methods, while rows correspond to varying initial CR pressure fractions, X_{cr, 0}. The 1D probability densities of the respective variables are shown above each axis: columns reflect the resulting n_H distributions for a fixed CR transport method across different X_{cr, 0} and rows show the resulting T distributions for a fixed X_{cr, 0} across different CR transport methods. The phase diagrams are constructed using 70 bins of equal (logarithmic) width for both variables. All simulations are analysed at t=6 t_cool.