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

[Lingrui Lin]

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.


   !accretion_prowered_irradiation = .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>