I agree that since Tp >> Ts for all known bodies, that jarmonik's equation 7 (post 39) should be used.
Also, I tested it in the case of an Orbiter precession model for Triton, which has a very large angle between all three planes (ref, rel, ecliptic), and was able to create an Orbiter model...