[Mesa-users] How to use custom (or make modifications) mass loss recipe?
Evan Bauer
ebauer at physics.ucsb.edu
Sat Jun 30 16:33:57 EDT 2018
Hi Matheus,
For applications like this, making direct changes to the source code in star/private is not a good practice.
Instead, the files in star/other provide examples of hooks that you can include in your run_star_extras.f to modify input physics. The file other_wind.f90 should be useful for you. You can find a tutorial on using these hooks in Josiah’s lecture from MESA summer school: https://jschwab.github.io/mesa-2017/ <https://jschwab.github.io/mesa-2017/>
Cheers,
Evan
> On Jun 30, 2018, at 8:23 AM, Matheus Bernini via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>
> Hello!
>
> I'm an undergraduate student. I'm using MESA (r-9575) to study the evolution of massive stars.
>
> I'd like to know how to write on MESA a custom mass loss recipe (actually, change an already existent). I tried change one of the recipes on the $MESADIR/star/private/winds.f90:
>
> subroutine eval_van_Loon_wind(w)
> ! van Loon et al. 2005, A&A, 438, 273
> real(dp), intent(out) :: w
> real(dp) :: log10w
> include 'formats'
> !log10w = -5.65d0 + 1.05*log10_cr(L1/(1d4*Lsun)) - 6.3d0*log10_cr(T1/35d2) !original
> log10w = -50 + 70*log10_cr(L1/(Lsun))
> w = exp10_cr(log10w) !changed (just for testing)
> end subroutine eval_van_Loon_wind
>
> but when I run a model after making this change on the winds.f90 nothing different happens. It's like I have done nothing. What could I do to make my model 'notice' the change?
>
> Thank you!
>
> Matheus
>
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180630/cc423e2e/attachment.html>
More information about the Mesa-users
mailing list