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

[Koushik Sen]

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/d1198257/attachment.htm>