Tisserand invariant with respect to Jupiter

[mkelley]

For isolating objects in cometary orbits and discriminating between cometary dynamical classes, the Tisserand invariant with respect to Jupiter is very helpful:

T_J = a_J / a + 2 * cos(i) * sqrt(a / a_J * (1 - e**2))

Having this in the mpcorb table, or even coded into a table view, would make quick work for the above use case.

What to do when e=1? It might be T_J = 0, but need to verify that some other time.