The integrated autocorrelation time for the MCMC chains of each physical parameter from our atmospheric model analysis, computed with different chain lengths based on the Goodman & Weare (2010) method (GW2010; solid squares) and the revised version suggested by Foreman-Mackey (DFM 2017; open squares), both of which produce consistent estimates. As suggested by the emcee documentation, the chains are supposed to converge once their lengths exceeds 50 times the estimated autocorrelation time from both (dashed line). While the chains for different parameters require different chain lengths (purple shadows), all of our chains converge after 2.5 × 104 samples (purple vertical lines).