[Mesa-users] Stationary model

[g.mirouh at surrey.ac.uk]

Dear MESA users,

I am new at using the code and to test some specific physics, I would 
like to stop the nuclear reactions in the core while keeping the 
luminosities to the value they had before stopping the reactions. 
Basically I want to create a stationary model and let time go on with no 
changes in the composition or the structure.

I tried by setting

s%rate_factors=0d0
s%weak_rate_factor = 0d0

which keeps my composition constant, but obviously my "star" collapses 
and cools slowly. How can I "force" the luminosity/radiation 
pressure/etc. to keep my structure invariant with time?

Thanks a lot for your help, and apologies to whom is not concerned.
Best.
-- 
Giovanni Mirouh.