When fitting with the new evolutionary tracks, the pdf1d output files do not all have the same shape.
Namely, the hdus for M_ini and log(M_ini) vary in shape (N_stars, N_mass_gridpoints), where N_mass_gridpoints is often missing gridpoints.
Specifically, I have 11 pdf1d files corresponding to 11 source density bins, and the M_ini and log(M_ini) hdus have either 30, 28, or 27 grid points included.
It seems like the BEAST is lopping off the smallest mass bins where they aren't being used. It would be helpful if these unused bins were still included in the pdf1d files, even if the associated values are NaNs or zeros.
When fitting with the new evolutionary tracks, the pdf1d output files do not all have the same shape.
Namely, the hdus for M_ini and log(M_ini) vary in shape (N_stars, N_mass_gridpoints), where N_mass_gridpoints is often missing gridpoints.
Specifically, I have 11 pdf1d files corresponding to 11 source density bins, and the M_ini and log(M_ini) hdus have either 30, 28, or 27 grid points included.
It seems like the BEAST is lopping off the smallest mass bins where they aren't being used. It would be helpful if these unused bins were still included in the pdf1d files, even if the associated values are NaNs or zeros.