[mesa-users] MESA3 paper 7_binary_with_rot
Pablo Marchant
pamarca at gmail.com
Tue Jun 16 03:28:36 EDT 2015
No, but the use_other_tsync option is there for you to play with it ;)
On Tue, Jun 16, 2015 at 9:10 AM, Roni Waldman <roni181066 at gmail.com> wrote:
> 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
>>
>
--
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/e3fcfc82/attachment.html>
More information about the Mesa-users
mailing list