Reduced data products and models from "Limits on the post-eclipse spectrum of HD 80606 b from high-resolution spectroscopy"

This repository contains reduced data products, planetary and stellar models, and posteriors from Finnerty et al. 2026a. The repository contains all of the planet models plotted in that paper and/or used to compute Kp-vsys maps in conjunction with the included reduced data products. HD80606_jan21_NIRSPEC.h5 contains the reduced data products. The archive keys are "baryrv" (barycentric radial velocity for each frame in km/s), "time" (midpoint time of each frame, in JD), "flux" (extracted fluxes, shape nframes x norders x npixels, in counts), and "wave" (wavelength, shape norders x npixels, in microns). Raw data products are available via the Keck Observatory Archive. The files ending in "post_equal_weights.dat" are the equally weighted posterior samples produced by the MultiNest runs. Each row is a sample. Columns are the parameters, in the order listed in Table 2 of the paper. The final column is the computed log-likelihood of the sample. The files ending in "_emission_petit.dat" are planet models produced with petitRADTRANs. The first column is the wavenumber, in cm^-1, and the second column is the flux, in erg/s/cm^3. HD80606_T5500_logg4.5_Z0.5_vsini1.8_orders.npy is the stellar model used, which is based on the PHOENIX spectral library. This file must be opened with allow_pickel=True option in numpy.load. The 9 entries correspond to the nine spectral orders. For each order, there is an array with the full-resolution stellar model, with rotational broadening applied. The first column of the array is the wavelength (in micron), and the second is the flux in erg/s/cm^3.