J_max (top panel) and H_max (bottom panel) residuals, between simulations with coeval coeval J- and H-band epoch and reference value. The weighted mean and uncertainty on the weighted mean in bins of two days are shown in blue. A “correction snake” and its uncertainty are calculated by fitting the residuals with GPs (red line and shaded region), while N is the total number of simulations.