Upper panel: HRD showing model luminosities and effective temperatures. The colour of the markers represents the necessary dust temperature required to match the remnant colours. For reference, we overplot single-star evolutionary tracks from the BPASS models (Eldridge et al. 2017; Stanway & Eldridge 2018). We highlight the progenitor appearance in early 2016 from Brennan et al. (2022b). Lower panel: matching spectra from DUSTY. Markers with error bars represent the extinction-corrected flux from the 2021 data given in Table 1. The colour of each spectral energy distribution represents the luminosity of the scaled model. We treat the F275W measurement from March 2021 and the F350LP data from December 2021 as upper limits.