This are elemental_abund plot and iso_abund with decayed=True option.that should be possible to do with mass fraction average over a mass range within a radial profile from a multi-zone output such as mppnp or mesa se
Example:
pt.iso_abund(47740,amass_range=[6,32],mass_range=list(xxl),decayed=True)
Results in:
Decay option not yet implemented for mppnp - but it is easy do! Consider investing the time!