The evolution defined in the following source code has one of the star's envelope locking and unlocking properly, but for the other star, the envelope seems to remain locked in a 2:1 spin-orbit ratio even when the eccentricity is negligible, which should not happen:
The evolution defined in the following source code has one of the star's envelope locking and unlocking properly, but for the other star, the envelope seems to remain locked in a 2:1 spin-orbit ratio even when the eccentricity is negligible, which should not happen:
test_evol.cpp.txt