[mesa-users] Extra diffusion

[Corentin Cadiou]

Hi MESA users,

I am trying to add a custom routine accounting for diffusion of chemical
elements that is related to the rotation rate (Omega) of the star and other
parameters.

To do that, I'm using a custom 'other_D_mix' routine in the
'src/run_star_extras.f' file that is mostly doing:
> s%D_mix = s%D_mix + factor*diffusion_coefficient
Where factor is 0.8 and the diffusion coefficient is ~15 (cm²/s).

If I then launch MESA, I get the error:

> mixing_type conv_vel           1           0    5.9868160194578010D-07
> STOP set mixing info

Looking at star/private/mix_info.f90:364 (where the error is thrown), I
see that MESA is stopping if the mixing_type is 0 and the convection
velocity positive. What should be mixing_type?

I also saw there was an 'other_diffusion' hook. What's the difference
with 'other_D_mix'? And which one should I use?

Thank you,
Corentin Cadiou



-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 819 bytes
Desc: OpenPGP digital signature
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160809/36f92403/attachment.sig>