[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