[mesa-users] adding code and variables to mesa to experiment with new capabilities

[Bill Paxton]

Please send whatever is necessary for other folks to reproduce the problem.  inlists, starting models, run_star_extras......

I've found that i'm unable to guess what is going wrong inside mesa!   i have to do it the hard way and debug.  ;D

-B




On Apr 30, 2014, at 9:16 AM, Tomer Shacham wrote:

> Hi Bill,
> Inlists will be useless since I'm using my other routines...
> On Apr 30, 2014 6:39 PM, "Bill Paxton" <paxton at kitp.ucsb.edu> wrote:
> Hi Tomer,
> 
> I think we can get some answers from other folks on mesa-users, especially concerning your #1.
> 
> And please send inlists, starting models, etc. so we can reproduce the behavior of your #2.
> 
> Cheers,
> Bill
> 
> On Apr 30, 2014, at 2:17 AM, Tomer Shacham wrote:
> 
>> Hi Bill,
>> I am continuing to work on my routines, and have encountered some weird behaviour:
>> 1. L_rad/L_edd can become larger than 1. (up to 2 or 3...) how is this possible
>> 2. My star is rapidly loosing mass  (msun/year...), and this loss is unseen by  s% mstar_dot. 
>> 3. (how) are 1. and 2. related in MESA?
>> Thanks for your help,
>> Tomer
>> 
>> --
>> Tomer Shacham
>> 
>> 
>> High Energy Physics Group
>> Institute     of     Physics
>> The  Hebrew  University
>> Jerusalem,   Israel  91904
>> 
>> 
>> 
>> On 19 March 2014 19:03, Bill Paxton <paxton at kitp.ucsb.edu> wrote:
>> 
>> On Mar 19, 2014, at 8:31 AM, Tomer Shacham wrote:
>> 
>>> Hi Bill,
>>> okay.  write routines for me that decide if an iteration is necessary and, if so, which parameters are changed and how for the next iteration so that it might get closer to being an okay solution.
>>> 
>>> I'm on it. There is one last thing though - I need to be able to change the boundary conditions of the star. I can change the pressure using 
>>> 's% surface_extra_Pgas' as you mentioned, but I also need to change rho. (what I really need is boundary conditions on T,and Rho, but can use the EOS to set P with it). 
>>> these boundary conditions need to be iterated along with my wind routine and opacity_factor routine (which is iterated already).
>> 
>> You'll want to use the hook "other_atm".  It sets the required T and P for the outer edge of the outer most cell.  That location is usually the photosphere, but it is not required to be at any particular optical depth.  The routine also returns the value of Teff -- although in your case that may not be well defined?   In normal operation, the surface T and P are functions of M, R, L, and kap at the surface; the routine returns partials wrt these.
>> 
>> -B
>> 
>> 
>> 
>> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140430/5b52199c/attachment.html>