Add 2010 TK7 trojan analysis example #126
Conversation
dahlend
commented
Sep 27, 2024
dahlend
commented
There was a problem hiding this comment.
I'm not sure if the state propagation is quite right. I see now it's two_body, not n_body, but the frac*period_? part only will keep the right spacing if the periods are basically the same. If they are different, you can get drift between the JDs of state_a and state_b. I think you need to just clock the MA and then convert from elements to a state, then evaluate, rather than propagate
|
|
||
| total = 0 | ||
| for frac in np.linspace(0, 1, m_steps + 1)[:-1]: | ||
| state_a = kete.propagate_two_body([state_a], jd + (frac + shift) * period_a)[0] |
There was a problem hiding this comment.
this is not quite right. I don't believe you want to do propagation here. I think rather you want to sweep mean anomaly for the instantaneous orbits, then evaluate the state at each of those MA clockings. Time would be a third axis to consider; e.g. make one plot for time=jd, then integrate the jd of the system a step forward and make another plot.
There was a problem hiding this comment.
i.e., I'm not sure that this and the following line maintain the same delta-longitude for all steps in m_steps when period_a and period_b aren't the same, and I think you need to maintain that spacing to get the energy correct
