[mesa-users] Extending MESA for mass-loss

[Bill Paxton]

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>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160728/9a0123cf/attachment.html>