[Mesa-users] Issues with modifying double_bh test suite to switch accretor and donor

[Levi]

Dear MESA-users group,



I am currently using MESA version 24.03.1 and have been attempting to modify the double_bh test suite to switch the accretor and donor stars based on certain conditions. To achieve this, I added the following code to the run_binary_extras.f90 file:
if (b% rl_relative_gap(b% d_i) > -0.05 .and. b% rl_relative_gap(b% a_i) > -0.05) then      write(*,*) "rl_relative_gap(b% a_i)", b% rl_relative_gap(b% a_i)      write(*,*) "s_accretor% w_div_w_crit_avg_surf", b% s_accretor% w_div_w_crit_avg_surf      if (b% s_accretor% w_div_w_crit_avg_surf > 0.95 .and. b% s_donor% w_div_w_crit_avg_surf < 0.95) then         switch_donor = .true.      else if (b% s_donor% w_div_w_crit_avg_surf > 0.95 .and. b% s_accretor% w_div_w_crit_avg_surf > 0.95) then         extras_binary_finish_step = terminate    end if   end if  
However, upon running the simulation, I encountered the following error:
get_xa_for_accretion: accretion species mass fractions do not add to 1.0                                       sum(xa(1:species))    0.0000000000000000D+00                                                       h1           1    0.0000000000000000D+00                                                      he3           2    0.0000000000000000D+00                                                      he4           3    0.0000000000000000D+00                                                      c12           4    0.0000000000000000D+00                                                      n14           5    0.0000000000000000D+00                                                      o16           6    0.0000000000000000D+00                                                     ne20           7    0.0000000000000000D+00                                                     mg24           8    0.0000000000000000D+00   retry: get_xa_for_accretion failed in adjust mass     176   retry: get_xa_for_accretion failed in adjust mass     176
I am unsure what might be causing this issue, and I have attached the terminal output in the file rn.out for your reference.
Could anyone please provide some insight into what might be causing this error and how I might resolve it?
Thank you very much for your time and assistance.


Best regards,
Levi



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240717/6092c36a/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rn.out
Type: application/octet-stream
Size: 282634 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240717/6092c36a/attachment.obj>