Open ezapartas opened 1 month ago
The idea is to use M4 and m4 for estimates of Ekin of each SN event. WIth @dimsour94 we already have some results with this method, but we relied up to now only to nearest neighbour interpolation, where M4 and m4 we recalculated. We want them to be added in the post-processin of the interpolated quantities
We should discuss this in the next Dev meeting, whether we need those values, which are specific to a SN prescription or whether the new values at He depletion are sufficient to simply recalculate M4 and mu4 in case they are needed.
Keeping this discussion of whether to add these values in the interpolated quantities too.
avg_c_in_c_core_at_He_depletion and co_core_mass_at_He_depletion should in general also be stored to be used in post-processing to get info about the collapse. In principle we have functions to recalculate them after a POSYDON run, but the functions are working on a singlestar POSYDON object not the data for the h5d5 file, so it is much much more straightforward to keep these values during the calculation of step_SN (and the interpolator).
The values at He depletion are ready to be used in the star object, which takes them from the ones from the hdf5 file of the grids, see https://github.com/POSYDON-code/POSYDON/pull/281/files#diff-e4b192b1c4e79f131a84d9c341737303dc7932b4a3c7136ea75c33592320bf86
…inglestar object as scalars