[Mesa-users] Angular momentum mixing error

[Pablo Marchant]

Hi Andriy,

I can't manage to compile your code, apparently because your
run_star_extras includes

s% other_am_mixing => my_mixing

but there's no my_mixing subroutine. It would also be useful if you specify
more clearly where the equations you use come from (specific references to
equation numbers in papers are ideal). For instance, I can't find the
source of this equation

jdot = s% x_ctrl(3) * 1.1265d47 * s%omega(1) * s%omega(1) * s%omega(1)

Cheers.

On Tue, May 26, 2020 at 9:11 PM Razinkov, Andriy Vasylhovych (1313xxx) <
andriy.razinkov at edu.uni-graz.at> wrote:

> ​​Dear MESA community
>
> I am writing with a question regarding the implementation of rotation in
> MESA (version 12115). I am trying to reproduce the spin-down of
> main-sequence stars, following Gallet et al (2013).
>
> So far I have implemented a disk-locking phase and wind loss following
> Kawaler (1988). The model can reproduce the spin down if solid body
> rotation is assumed and starting from a 1Myr 1Mo model, with up to 10 times
> the current solar rotation rate.
>
> However, when in this setup the star is allowed to rotate non-uniformly
> using "set_uniform_am_nu_non_rot = .false.”, MESA stops with the error
> message "failed in mixing angular momentum“.
>
> I would be grateful if someone could give me some suggestions on how to
> better understand this error message and how solve this problem in my
> models.
>
> For further reference, I attach the inlists and run_star_extra for my
> project
>
> Best Regards
> Andriy​
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>
>

-- 
Pablo Marchant Campos
M.Sc on Astrophysics, Universidad Católica de Chile
PhD on Astrophysics, Argelander-Institut für Astronomie, Universität Bonn
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200603/015fe6d4/attachment.htm>