The resolution of 192 grid cells along the slab of their calculations, from which they derived their results, is comparable with the one used in our high resolution case. We use the same method as described in Sect. 4.3 to produce the initial conditions. The SPHI calculation slightly underestimates and , as is also observed for the grid code. The units of the x-axis are normalized to the position of the shock front.
Ionizing radiation comes from the left.