Top panel: continuum-subtracted iSHELL data (black) with the best-fitting model (red) based on the various slab model fits to the rotational diagrams over the 4.90–5.25 µm wavelength region. The middle panel shows the best fitting model binned to the JWST-MIRI resolution together with the JWST-MIRI observations, whereas the bottom panel shows the model convolved with the pseudo-Voigt profile shown in Fig. D.6.