# [Mesa-users] some problems in using Mesa to calculate wind transfer

Lingrui Lin llrw at 163.com
Thu Dec 12 04:25:19 EST 2019

oh, I am so sorry. I am new to mesa and new to do research
here is my inlist_project:
&binary_job

inlist_names(1) = 'inlist1'
inlist_names(2) = 'inlist2'

evolve_both_stars = .false.
change_ignore_rlof_flag = .true.
/ ! end of binary_job namelist

&binary_controls

m1 = 1.0d0  ! donor mass in Msun
m2 = 1.2d0 ! companion mass in Msun
initial_period_in_days = 1d0

fr = 0.05
fr_limit = 1d-2

!do_enhance_wind_1 = .true.
!do_enhance_wind_2 = .true.

do_wind_mass_transfer_1= .true.
limit_retention_by_mdot_edd = .true.
use_other_binary_wind_transfer = .true.
/ ! end of binary_controls namelist

and here is my other_wind subroutine:
and of cource I have set  b% other_binary_wind_transfer => wind_transfer_routine
The print* in the subroutine work all properly
subroutine wind_transfer_routine(binary_id,s_i,ierr)
integer,intent(in)::binary_id,s_i
integer,intent(out)::ierr
real(dp)::etas,etaa,phi,etahydr,etahe,lim_mtransdot,C1,loghe
type(binary_info),pointer::b
ierr=0
call binary_ptr(binary_id,b,ierr)
if(ierr/=0)then
write(*,*)'failed in binary_ptr'
return
end if
etahydr=0.15
phi = 0.5
etas=0.9
etaa=1
if (b%m(b%d_i).lt.b%m(b%a_i)) then
b%wind_xfer_fraction(b%d_i) = &
-b%s_donor%mstar_dot*secyer*1d7/3.5**2* &
b%m(b%d_i)**(-5/3)* &
(b%m(b%a_i)+b%m(b%d_i))**(2/3) &
/(etas*etaa)/phi**2
else
b%wind_xfer_fraction(b%d_i) = &
-b%s_donor%mstar_dot*secyer*1d7/3.5**2* &
b%m(b%d_i)**(-1.9)* &
(b%m(b%a_i)+b%m(b%d_i))**(2/3) &
*b%m(b%a_i)**0.24/(etas*etaa)/phi**2
end if
print*,'b%m(b%d_i)/Msun',b%m(b%d_i)/Msun
print*,'b%m(b%a_i)/Msun',b%m(b%a_i)/Msun
print*,'lim_mtransdot000',lim_mtransdot
lim_mtransdot=7.2d-6*(b%m(b%a_i)/Msun-0.6)
print*,'lim_mtransdot',lim_mtransdot
!C1=3
if (-b%mtransfer_rate*secyer/Msun &
.gt.lim_mtransdot) then
print*,'critical accretion'
b%mdot_wind_transfer(b%d_i)=-lim_mtransdot* &
Msun/secyer
end if

loghe=log10(-b%mtransfer_rate*secyer/Msun* &
b%s_donor%center_he4)
!print*,'centerhe4',b%s_donor%center_he4
print*,'loghe',loghe
print*,'b%mdot_wind_transfer(b%d_i)', &
b%mdot_wind_transfer(b%d_i)
print*, 'b%mtransfer_rate',b%mtransfer_rate
print*,'b%s_donor%mstar_dot', &
b%s_donor%mstar_dot
print*,'b%wind_xfer_fraction(b%d_i)', &
b%wind_xfer_fraction(b%d_i)
if(loghe.gt.-7.06d0.and.loghe.lt.-5.95d0) then
etahe=0.54*loghe+4.16
else if (loghe.gt.-5.95d0.and.loghe.lt.-5.76d0)then
etahe=0.54*(loghe+5.6)**2+1.01
else
etahe=1
end if
end subroutine wind_transfer_routine

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20191212/b1a2304f/attachment-0001.htm>