Implement new kick prescription

[TomWagg]

tldr;

This PR replaces the default COSMIC prescription for handling kick velocities. Instead of using Kiel & Hurley 2009, it now follows Pfahl+2002. The former produces strange results for disrupted systems, in which the secondary ejection velocity showed a dependence on the natal kick strength, which goes against physical intuition (see #675).

Changes

Breaking code changes:

• Kick info table is now (2, 18) not (2, 17), where the extra column is included so every Euler angle can be recorded
• kickflag is now an option that can be 1, 2, 3, or 4 to replace the old options. Additionally, if users want to use the old K&H prescription, you use the negative of any of those values.

Other code fixes (note: in the interest of reproducibility I haven't fixed these in the old Kiel prescription)

• A coordinate transformation for secondary stars reaching supernova in an already disrupted binary was wrong (assumed co-lateral coordinates erroneously), I've fixed this
• Not sure if this is a fix per-se but the criteria for collisions were weird in the old version, I now do proper vector calculations to find their closest separation after the SN and check if it is within the secondary radius

Physical evolution changes

• Ejection velocities for secondaries are now much more tightly distributed around the pre-SN orbital velocity

Comparison of a disrupted system

For unbound systems, the velocity of the companion has been changed a fair bit. Here's 5000 random draws for the same binary (m1=20, m1=15, porb=200, Z=0.02, ecc=0)

Comparison of a bound system

And for bound systems, the final separations are pretty much unchanged (these are just different random draws so you see a little scatter).

This will close #675.

[katiebreivik]

This looks great! if you can add a test to make sure that we have a way to test whether or not a future development for the kick subroutine breaks what you have now, then we are good!

[codecov]

Codecov Report

Attention: Patch coverage is 39.90826% with 393 lines in your changes missing coverage. Please review.

Project coverage is 79.63%. Comparing base (8772c07) to head (00e2c2d).
Report is 87 commits behind head on develop.

Files with missing lines	Patch %	Lines
src/cosmic/src/kick.f	2.73%	392 Missing ⚠️
src/cosmic/utils.py	0.00%	1 Missing ⚠️

Additional details and impacted files

@@             Coverage Diff             @@
##           develop     #678      +/-   ##
===========================================
- Coverage    86.91%   79.63%   -7.28%     
===========================================
  Files           40       48       +8     
  Lines        25542    26992    +1450     
===========================================
- Hits         22198    21494     -704     
- Misses        3344     5498    +2154     

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.