[Mesa-users] Heavy element profile

[Francis Timmes]

one could try accreting material of a specified composition

http://mesa.sourceforge.net/controls_defaults.html#composition_controls

fxt




> On Jan 11, 2019, at 3:06 AM, valletta <valletta at physik.uzh.ch> wrote:
> 
> 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
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>