Skip to content

Conversation

@awilson110
Copy link
Collaborator

Fixes #76

Added functionality to compute 1/a0 values for comets by backwards (or forwards, if necessary) integrating using an assist simulation (or rebound, if the simulation goes outside the assist timeframe). Includes a demo for layup comet.

Review Checklist for Source Code Changes

  • Does pip install still work?
  • Have you written a unit test for any new functions?
  • Do all the units tests run successfully?
  • Does Layup run successfully on a test set of input files/databases?
  • Have you used black on the files you have updated to confirm python programming style guide enforcement?

reads in the lpc file, converts to KEP to check a values, if they're too small (indicating a short period comet) remove that data.
First draft of the code, there are a few edge cases to work out still
Switching to pure rebound is still a work in progress, but if the comet remains inside the timeframe of the assist ephemeris then the program works as intended
Also cleaned up the code further, now converting to rebound works as expected.
Each test has its own file to read in and a file of expected values too
Changed the order of arguments in the functions so they're the same, made changes to for loops to make them more pythonic (also switched to storing output data as a dict for the same reason), made the tolerances on some unit tests smaller
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

calculate 1/ao for comets

1 participant