[Mesa-users] Trying to inject energy into a star

[Bill Wolf]

And for a more guided tour of using run_star_extras.f, I highly recommend working through Josiah Schwab's tutorial from this year's MESA summer school:

https://jschwab.github.io/mesa-2018/

which also links to a mini-guide on fortran syntax. 

Regards,

Bill


From: Aaron Dotter via Mesa-users <mesa-users at lists.mesastar.org>
Reply: Aaron Dotter <aaron.dotter at gmail.com>
Date: November 8, 2018 at 9:28:07 AM
To: yotam.cohen at yale.edu <yotam.cohen at yale.edu>
Cc: mesa-users <mesa-users at lists.mesastar.org>
Subject:  Re: [Mesa-users] Trying to inject energy into a star  

Hi Yotam,

Have a look at mesa/star/test_suite/envelope_inflation/  

Try running the test case and reading src/run_star_extras.f.

It uses the "other_energy" hook provided by MESA and the array x_ctrl(:) which is accessible from the &controls inlist to input the parameters.

Cheers,
Aaron

On Thu, Nov 8, 2018 at 11:03 AM Yotam Cohen via Mesa-users <mesa-users at lists.mesastar.org> wrote:
Hello,

My name is Yotam Cohen and I’m an astronomy graduate student at Yale University.

This isn’t exactly troubleshooting, but I was hoping someone could help me set up the following situation in MESA. Basically, I want to set up a normal stellar model but I want to inject an anomalous source of energy into the star. More specifically, I would like to be able to inject the energy at a specific zone or zones, say at each zone between R_i and R_j for this example, and I want to do it at specific timesteps, say between T_i and T_j, and I want to inject a certain fraction of the overall stars luminosity, say 1%, into each of those zones during each of those timesteps.

Now, I understand how to setup the vanilla stellar model, and I understand that I can achieve what I’m trying to do using the run_star_extras module, but I’m having trouble figuring out how to exactly go about it. I’m a MESA novice and a Fortran novice, so my main problems are: 1) I don’t know the names of or how to properly access the arrays/quantities that I’m trying to use and 2) I don’t know how/where exactly within the run_star_extras module to set up this code.

Any help would be greatly appreciated. If I can provide any more information, please let me know.

Thanks,
Yotam
————————
Yotam Cohen
Graduate Student
Gruber Science Fellow
Dept. of Astronomy
Yale University

_______________________________________________
mesa-users at lists.mesastar.org
https://lists.mesastar.org/mailman/listinfo/mesa-users

_______________________________________________  
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/20181108/84f3fc08/attachment.html>