[mesa-users] Convergence issues when simulating massive eruptions by modifying .mod file

[Josiah Schwab]

> I do this by artificially modifying the saved model file and restarting the evolution with loading that file.

I tried manipulating model files once and in my experience (and as you seem to be finding) that way lies madness :)

Is there something about MESA's built-in routines for changing the mass of model that are insufficient/inappropriate for what you want to explore?  What other (if any) approaches have you tried?

My personal recommendation is to use an existing relaxation routine or to add something you write yourself to run_star_extras.f to do the relaxation.  I think it will save you a lot of grief.

Perhaps others will share their experiences.

Josiah