Top. WFC3+GNIRS+MIRI consistency plot. The posterior predictive distribution $p(f(\theta )+ϵ|x_{\text{obs}})$ of noisy simulations spectra for the 99.7%, 95% and 68.7% quartiles (hues of blue), overlaid on top of the WFC3+GNIRS+MIRI observation x_obs (black line). Bottom. Residuals of the posterior predictive samples, normalized by the inflated standard deviation of the noise distribution for each spectral channel and a horizontal line at 0 for reference (in black).