[mesa-users] Extending MESA for mass-loss

[Francis Timmes]

additional guidance is also available at 

http://mesa.sourceforge.net/run_star_extras.html#toc-3

and the 2015 summer school tutorial, "Part 2c: Changing input physics"

http://jschwab.github.io/mesa-2015/

fxt






> On Jul 28, 2016, at 10:41 AM, Bill Paxton <paxton at kitp.ucsb.edu> wrote:
> 
> Hi,
> 
> To tell mesa/star to use your routine, set 'use_other_wind' = .true. in the &controls section of your inlist, or set the wind scheme to 'other'.
> 
> During your initial debugging you will of course want to put write's in your other_wind routine to ensure that it is being called and to sanity check the values.
> 
> -B 
> 
> 
> 
> 
> 
> On Jul 28, 2016, at 10:30 AM, Farag, Ebraheem Khaled wrote:
> 
>> Hello,
>> 
>>     Thank you for your response. I've tried implementing the mass loss equation into MESA using the example_solar_model. The Problem I'm encountering is that my wind routine doesn't seem to be operating, since the mass is staying constant throughout the evolution. I've attached my inlist files and run_star_extras.f . Have I misused the other_wind routine, or am I not calling it correctly into the code?
>> 
>> -I also have attached the final_profile.data file showing full convergence of the initial model at 1M(when initial_mass is commented out and no mass-loss) with z/x = 0.0181 for AGSS09.
>> 
>> Thanks,
>> 
>> Ebraheem 
>> From: Bill Paxton <paxton at kitp.ucsb.edu>
>> Sent: Wednesday, July 27, 2016 11:58 AM
>> To: Farag, Ebraheem Khaled
>> Cc: mesa-users
>> Subject: Re: [mesa-users] Extending MESA for mass-loss
>>  
>> Hi, 
>> 
>> The arguments to other_wind include the current mass (in grams) along with L, R, and T.   
>> 
>> I suggest using s% time for age in seconds (the s% star_age variable is not updated until after the call on other_wind).
>> 
>> Good luck,
>> 
>> Bill
>> 
>> 
>> 
>> 
>> 
>> On Jul 27, 2016, at 10:38 AM, Farag, Ebraheem Khaled wrote:
>> 
>>> Hello all,
>>> 
>>>     I'm currently working to implement an exponential mass loss function as  implemented in the LANL Stellar Pulsation Codes.( The function can be seen in the attached image.)
>>> 
>>> I believe this can be done through extending MESA. My current idea is to modify the other_winds subroutine in run_star_extras.f, as to create function for the value of 'w'. One of my problems is that I don't know what variable in star_def to use for time and Solar mass at that given time. Do you think this is the best approach toward implementing mass loss in this form for a 1-1.15M solar evolution model of the sun, or is there a better approach? 
>>> 
>>> Thanks,
>>> 
>>> Ebraheem
>>> <mass_loss.png>------------------------------------------------------------------------------
>>> What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic
>>> patterns at an interface-level. Reveals which users, apps, and protocols are 
>>> consuming the most bandwidth. Provides multi-vendor support for NetFlow, 
>>> J-Flow, sFlow and other flows. Make informed decisions using capacity planning
>>> reports.http://sdm.link/zohodev2dev_______________________________________________
>>> mesa-users mailing list
>>> mesa-users at lists.sourceforge.net
>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>> 
>> <run_star_extras.f><inlist><inlist_solar_model><final_profile.data>
> 
> ------------------------------------------------------------------------------
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users