[mesa-users] MESA3 paper 7_binary_with_rot

[Roni Waldman]

Thanks Pablo.
Do you happen to have a subroutine for the Tassoul formula?
I'm trying to compare to Yoon et al 2010.

Roni

‫בתאריך יום ג׳, 16 ביוני 2015 ב-10:04 מאת ‪Pablo Marchant‬‏ <‪
pamarca at gmail.com‬‏>:‬

> My bad, here it is. Will update the ones on mesastar
>
> On Tue, Jun 16, 2015 at 8:37 AM, Roni Waldman <roni181066 at gmail.com>
> wrote:
>
>> Hi,
>> More on this one.
>> The inlist_project has:
>>    ! use_other_tsync will use the Zahn timescale as defined in
>> run_binary_extras.f
>>    !use_other_tsync = .true.
>>
>> I dont't see any other_tsync in the default run_binary_extras.f, and
>> there's no special run_binary_extras.f in the 7_binary_with_rot
>> directory.
>>
>> Roni
>>
>>
>> ‫בתאריך יום ב׳, 15 ביוני 2015 ב-21:11 מאת ‪Roni Waldman‬‏ <‪
>> roni181066 at gmail.com‬‏>:‬
>>
>>> Same here.
>>>
>>> ‫בתאריך יום ב׳, 15 ביוני 2015 ב-19:37 מאת ‪Willie Strickland‬‏ <‪
>>> cwskas at gmail.com‬‏>:‬
>>>
>>>> Thanks for the confirmation.  I thought they were probably for
>>>> debugging, but didnt see obvious way to turn them off.
>>>>
>>>> Willie
>>>>
>>>>
>>>> On Jun 15, 2015, at 11:13 AM, Pablo Marchant <pamarca at gmail.com> wrote:
>>>>
>>>> The w > w_crit messages come from the implicit method that increases
>>>> the wind strength so the star does not evolve avobe critical rotation.
>>>> Maybe they could be turned off, I mainly used that info for debugging. And
>>>> your results are the same as mine.
>>>>
>>>> On Mon, Jun 15, 2015 at 5:59 PM, Willie Strickland <cwskas at gmail.com>
>>>> wrote:
>>>>
>>>>> Did you get a lot of w > w_crit lines in the rotating one?  I got the
>>>>> same results as you it sounds like, but had all these redo lines.
>>>>>
>>>>> 8,821 lines fix w > w_crit: bisect mdots and redo
>>>>>      89 lines w > w_crit: revise mdot and redo
>>>>>
>>>>> My run finished at step 6429, age 1.254e7, 2.45 mass1, 18.7 mass2.
>>>>>
>>>>> accretor photosphere has reached its Roche Lobe
>>>>>                                      accretor (r-rl)/rl        6429
>>>>>  5.0151227593021562D-02
>>>>>
>>>>> Willie
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On Jun 15, 2015, at 7:49 AM, Roni Waldman <roni181066 at gmail.com>
>>>>> wrote:
>>>>>
>>>>> Hi Pablo,
>>>>>
>>>>> Oops... right, I put the file in the wrong place.
>>>>> Rerunning, and let's see what happens.
>>>>>
>>>>> I also ran the rotating binary example.
>>>>> Ended up with donor mass 2.45 Msun and accretor reaching Roche-lobe.
>>>>> Does that sound right?
>>>>>
>>>>> Thanks,
>>>>> Roni
>>>>>
>>>>> ‫בתאריך יום ב׳, 15 ביוני 2015 ב-14:11 מאת ‪Pablo Marchant‬‏ <‪
>>>>> pamarca at gmail.com‬‏>:‬
>>>>>
>>>>>> Roni, together with the inlists there is a run_star_extras.f provided
>>>>>> (and in the inlists it is said that you need it). Are you using it? From
>>>>>> you output it seems not. This implements the switch from a Kudritski wind
>>>>>> to a Hamann one as described in Wellstein & Langer 1999.
>>>>>>
>>>>>> I guess changing this will make it consistent with that of the paper.
>>>>>> Your model is doing some crazy things though.
>>>>>>
>>>>>> On Mon, Jun 15, 2015 at 9:06 AM, Roni Waldman <roni181066 at gmail.com>
>>>>>> wrote:
>>>>>>
>>>>>>> Hi Pablo,
>>>>>>> Too big for one email, will send several.
>>>>>>> Thanks,
>>>>>>> Roni
>>>>>>>
>>>>>>> ‫בתאריך יום א׳, 14 ביוני 2015 ב-22:55 מאת ‪Pablo Marchant‬‏ <‪
>>>>>>> pamarca at gmail.com‬‏>:‬
>>>>>>>
>>>>>>>> Roni, can you send me your history.data for each star and
>>>>>>>> binary_history.data to compare?
>>>>>>>>
>>>>>>>> On Sun, Jun 14, 2015 at 8:35 PM, Roni Waldman <roni181066 at gmail.com
>>>>>>>> > wrote:
>>>>>>>>
>>>>>>>>> Hi,
>>>>>>>>>
>>>>>>>>> I downloaded the inlists of the MESA3 paper, and ran the example
>>>>>>>>> 6_massive_binary, without changing anything in the inlists or the code,
>>>>>>>>> version 7624.
>>>>>>>>> The run ends with the donor having 2Msun, the accretor filling its
>>>>>>>>> Roche-lobe, and the time steps being super small (1e-12).
>>>>>>>>> This not what seems to happen in the paper.
>>>>>>>>> Am I doing anything wrong?
>>>>>>>>>
>>>>>>>>> Roni
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> ------------------------------------------------------------------------------
>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> mesa-users mailing list
>>>>>>>>> mesa-users at lists.sourceforge.net
>>>>>>>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> --
>>>>>>>> Pablo Marchant Campos
>>>>>>>> M.Sc on Astrophysics, Universidad Católica de Chile
>>>>>>>> PhD student, Argelander-Institut für Astronomie
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>>
>>>>>> --
>>>>>> Pablo Marchant Campos
>>>>>> M.Sc on Astrophysics, Universidad Católica de Chile
>>>>>> PhD student, Argelander-Institut für Astronomie
>>>>>>
>>>>>
>>>>> ------------------------------------------------------------------------------
>>>>> _______________________________________________
>>>>> mesa-users mailing list
>>>>> mesa-users at lists.sourceforge.net
>>>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>>>>
>>>>>
>>>>>
>>>>
>>>>
>>>> --
>>>> Pablo Marchant Campos
>>>> M.Sc on Astrophysics, Universidad Católica de Chile
>>>> PhD student, Argelander-Institut für Astronomie
>>>>
>>>>
>>>>
>>>> ------------------------------------------------------------------------------
>>>> _______________________________________________
>>>> mesa-users mailing list
>>>> mesa-users at lists.sourceforge.net
>>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>>>
>>>
>
>
> --
> Pablo Marchant Campos
> M.Sc on Astrophysics, Universidad Católica de Chile
> PhD student, Argelander-Institut für Astronomie
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150616/fea302c2/attachment.html>