OD Code is a program that allows the user to determine the orbital elements of an asteroid and uncertanties of the calculation, given the right ascension, declination, time, and sun vectors at 3 points on the orbit.
- Calculate sun to asteroid position and velocity vectors using OD_final
- Generate the orbital elements using the given vectors through orb_elements and print using orb_main
- Generate the ephemerus using the orbital elements through eph_gen and print using eph_main
- Complete error correction and Monte Carlo accuracy simulation
OD_final(ra1, dec1, ra2, dec2, ra3, dec3, t1, t2, t3, R1, R2, R3) Calculates the position and velocity of the asteroid at the second observation and returns as r2, r2_dot
read_csv(name) Reads a file seperated by spaces and newlines and returns as t1, ra1, dec1, R1, t2, ra2, dec2, R2, t3, ra3, dec3, R3
f_g_series(tau, r2, r2dot) Calculates the f and g series
d_vals(ra1, dec1, ra2, dec2, ra3, dec3, R) Calculates the d values used in OD_final and returns as D_0, D_1[0], D_1[1], D_1[2], D_2[0], D_2[1], D_2[2], D_3[0], D_3[1], D_3[2]
eph_main(expected, orbit) Outputs ephemerus and errors given the expected and actual right ascension and declination
eph_gen(a, e, i, O, w, T, t) Calculates and outputs the right ascension and declination of an asteroid given its orbital elements
error(calculated, actual) Calculates percent error in a given calculation
orb_main(expected, orbit) Outputs orbital elements and errors given the expected and actual orbital elements
orb_elements(r, r_dot, deg, t) Calculates and outputs orbital elements given position and velocity of second observation, whether it should be in degrees, and the time of observation
Written by Ellie Zhang with help from Ilisha Gupta, Tejus Deo-Dixit, Vincent Stone, Andres Rodriguez, Michael Rodruck, Michael Dubson, Alan Tondryk, Emily Zhang, Kevin Kim, Kendra Nguyen et al.
2024-07-17