[mesa-users] Strange boundary mass reporting?

[Josh Shiode]

Users,

I'm looking at a series of what become wolf-rayet stars (Helium-burning, no hydrogen envelope), and I'm seeing that the he4_boundary_mass and c12_boundary_mass are being reported as identical (and equal to the total mass) once the Hydrogen envelope is removed. 

This is confusing because c12 > c12_boundary_limit (= 1e-4, the default) for the outermost few solar masses in most cases. This doesn't seem to happen when the hydrogen layer is not removed (but other masses are similar..).

Looking through mesa/star/private I find that this seems to be set in report.f, I see that there is a pre-setting of the variables equal to one another, but then get_info_at_mass_bdy( ) is called and should reset values to their correct ones. Anyone know what I'm missing here?

Thanks,
Josh