[Mesa-users] Failure of getting the effective temperature in run_binary_extras

[M Sun]

Hello Koushik,

Yes, it should work. As you may see, Rob already helped me how to grab the
surface quantities in the binary_wind_mass_transfer subroutine (e.g.
luminosity s% xh(s% i_lum, 1) and radius s% xh(s% i_lnR, 1) at the 1st
zone). And thank you so much for pointing me to the right further step! :)

cheers,
Meng

On Mon, Mar 15, 2021 at 3:51 AM Koushik Sen <senkoushik1995 at gmail.com>
wrote:

> Hi Meng,
>
> It's not a direct answer to your question but you can always calculate the
> effective temperature using the Stefan-Boltzmann equation if the pointers
> to radius and luminosity are working for you.
>
> Cheers,
> Koushik
>
> PhD student
> Argelander Institute for Astronomy
>
>
> On Sun, 14 Mar 2021, 1:34 am M Sun via Mesa-users, <
> mesa-users at lists.mesastar.org> wrote:
>
>> Dear MESA users,
>>
>> I am trying to add a wind mass accretion mechanism in MESA 12115. I
>> failed to get the effective temperature of the donor star in the subroutine
>> my_other_binary_wind_transfer from src/run_binary_extras.f. I guess that I
>> misuse/misunderstand the pointers. Here is the detail ---
>>
>> The declaration and association of the star and binary pointers are:
>> subroutine my_other_binary_wind_transfer(binary_id, s_i, ierr)
>>       implicit none
>>       integer, intent(in) :: binary_id, s_i ! s_i is index of the wind
>> mass losing star
>>       integer, intent(out) :: ierr
>>       type (binary_info), pointer :: b
>>       type (star_info), pointer :: s
>>       ierr = 0
>>       call binary_ptr(binary_id, b, ierr)
>>       if (ierr /= 0) then
>>          write(*,*) 'failed in binary_ptr'
>>          return
>>       end if
>>
>>       if (s_i == 1) then
>>          s => b% s1
>>       else
>>          s => b% s2
>>       end if
>>       print*, b% s1% T(1)
>>
>> I call the effective temperature by "b% s1% T(1)". The output looks good
>> at the very first step, then it becomes "NaN" at the following steps.
>>
>> I realize T(1) is surface temperature, not the effective temperature,
>> just use it as an example of calling the star information from the binary
>> module. The run_binary_extras.f and the binary inlist are attached
>> (starting at line 99 in run_binary_extras is my code). Please let me know
>> if you need more information. Many thanks!!
>>
>> cheers,
>> --------------
>> Sun, Meng
>> Research Associate
>> Department of Astronomy
>> University of Wisconsin - Madison
>> https://sunmeng1118.wixsite.com/mysite
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210315/ce1051c7/attachment.htm>