[Mesa-users] Heavy element profile

[valletta]

Dear MESA users

Hello. I am modelling the evolution of a giant planet, I start the model 
with a given composition of X, Y and Z, constant through the envelope. 
Now I need to modify these compositions at every time-step, depending on 
the accretion of heavy element by the planet, and the amount of X,Y and 
Z must depend on the position of the layer in the envelope (so, a 
non-constant profile through the envelope). My idea to do that was to go 
in the run_star_extras file and write something like

for(i = 1, i< layer number)

    s% x(i) = new x(i)
    s% y(i) = new y(i)
    s% z(i) = 1 - new x - new y

the problem that I have is that these values are not saved by mesa, but 
at every time step the code restart from the initial profile of X,Y and 
Z.

Thank you in advance for every help

Best,
Claudio