Top: similar to Figure 3 but for 10 realizations with different initial perturbation seeds. Bottom: ratios between the model with and without thermal conduction for hot gas density (upper left), hot gas pressure (middle left) hot gas temperature (lower left), total cooling rate (upper right), hot gas cooling rate (middle right), and hot gas volume (lower right) as a function of the radius at different times in the 10 realizations. The lines show the median value, and the shaded regions show the range of all 10 pairs of models. Despite a large variance by a factor of ∼2–3, the SNRs with thermal conduction show systematically higher density, lower pressure, and thus lower temperature. The total cooling rate is overall similar. The variance is slightly larger at the early stage in the central region due to the distribution of the cold clouds.