[Mesa-users] NaN values when getting zones' opacity (in run_star_extras.f90)

[Matheus Bernini]

Thanks for the answer!

So, a few more questions:
(1) If the opacities are calculated after the mass adjustment, how can I
get (if there is a way to) the values used in previous steps?
(2) So, can I access the cells' opacities (even if they are from the
previous step) with the xh(:,:) ! (nvar, nz)? If so, where can I find the
index corresponding to the opacities?
(3) I called the subroutine you mentioned but it still giving NaN results.
Probably I didn't do this the right way, since my knowledge in Fortran is
not very advanced.



Em ter, 30 de jul de 2019 às 18:20, Pablo Marchant <pamarca at gmail.com>
escreveu:

> Hi Matheus,
>
> different variables are set at different points during a step in the
> evolution. You're trying to implement a new wind routine, there are
> guidelines on how to implement new wind physics under
> $MESA_DIR/star/other/other_wind.f90, and in particular there you'll find
> that
>
>            ! NOTE: surface is outermost cell. not necessarily at
> photosphere.
>            ! NOTE: don't assume that vars are set at this point.
>            ! so if you want values other than those given as args,
>            ! you should use values from s% xh(:,:) and s% xa(:,:) only.
>            ! rather than things like s% Teff or s% lnT(:) which have not
> been set yet.
>
> so you cannot directly access things like opacity at this point, it is
> computed further ahead in the code. In particular, all those variables are
> set after mass is adjusted. Perhaps one way you can work around this is by
> calling the function star_set_vars which is part of star_lib (look in
> $MESA_DIR/star/public/star_lib.f90 to see what arguments the function
> receives). That will set many structure variables including opacity.
>
> Cheers
>
> On Tue, Jul 30, 2019 at 2:24 PM Matheus Bernini via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
>> I've forgotten some files :)
>>
>> Em ter, 30 de jul de 2019 às 16:19, Matheus Bernini <
>> matheus.b.peron at gmail.com> escreveu:
>>
>>> Hello everyone.
>>>
>>> I'm having a problem to get the opacity in the star zones in
>>> run_star_extras, using the pointer star_info during the evolution of the
>>> model.
>>>
>>> Basically, when the model starts, before it relax the initial rotation,
>>> the opacities(:) have real values, but after the rotation relaxation, the
>>> returned opacities become NaN, and I can't figure out why this happens.
>>>
>>> Cheers
>>>
>>> [image: image.png]
>>>
>> _______________________________________________
>> 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/20190731/c904544f/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 120534 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190731/c904544f/attachment-0001.png>