TeamCOMPAS / COMPAS

COMPAS rapid binary population synthesis code
http://compas.science
MIT License
64 stars 67 forks source link

Units on post-CEE semi-major axis #876

Closed nrsegovia closed 1 year ago

nrsegovia commented 1 year ago

Describe the bug A comparison between the separation given by the RLOF logs against the CEE ones showed me that separation post CEE was incredibly small, implying that most stars were touching post-CEE.

Investigating this issue, I found that the ratio of semiMajorAxisPostRLOF/semiMajorAxisPostCEE is equal to 46238.77350477 in all cases, which is consistent with ~ 10000000/215 where 215 is just ~1AU in Rsol units. So I believe this is just a misplaced/missing unit conversion.

Note that, as far as I have been able to test, this issue only affects the log files. The flag for stellar merger is not triggered by this unusual post-event separation value.

Label the issue

urgency_low - This issue is not urgent

severity_minor - This is a minor bug with minimal impact

To Reproduce

Expected behavior The separation post event shown in the CEE and RLOF logs should be similar, as both should be expressed in Rsol units (according to the documentation).

Versioning (please complete the following information):

nrsegovia commented 1 year ago

Line 1659 in BaseBinaryStar.cpp seems to be the issue. The / operator should be *, I will submit a PR ~tomorrow.

nrsegovia commented 1 year ago

Fixed in #877