[mesa-users] Custom Irradiation

[Bill Paxton]

Hi,

The simplest way to get started would be to use the "other_wind" hook to define your special mass loss formula.

Take a look at mesa/star/other/other_wind.f90 for a few comments and a sample "null" version.

Using the null version as a template, implement your actual routine in your "run_star_extras",

make it accessible to the star mass loss code by setting the pointer
	s% other_wind => my_wind_routine
and tell star to call it on each timestep by setting 
	use_other_wind = .true.

Keep in mind that big jumps in mass loss rates from step to step can cause convergence problems.
You may need to limit timesteps by setting max_years_for_timesteps.

good luck,
b




On Mar 16, 2015, at 10:43 AM, etheridge d.j. (dje1g11) wrote:

> Hi,
> 
> I'm aware that the binary module in MESA can be customised using run_star_extras.f. How could I implement a different prescription for the mass lost from a donor star due to pulsar irradiation? I would be using equation 3 from Stevens et al. 1992, MNRAS, 254, 19P.
> Any help would be greatly appreciated.
> 
> Many thanks,
> 
> Daniel Etheridge
> ------------------------------------------------------------------------------
> Dive into the World of Parallel Programming The Go Parallel Website, sponsored
> by Intel and developed in partnership with Slashdot Media, is your hub for all
> things parallel software development, from weekly thought leadership blogs to
> news, videos, case studies, tutorials and more. Take a look and join the 
> conversation now.http://goparallel.sourceforge.net/_______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150316/84f8c672/attachment.html>