With the discussion that took place on answer(s) to How to calculate cone angle between two satellites given their look angles? I would like to ask if my proposal solving that question works like intended:

sun = ephemeris['sun']
earth = ephemeris['earth']
vecObserverSat = (s1 - loc).at(tEv)
vecObsserverSun = (sun - (earth + loc)).at(tEv)
phase = vecObserverSat.separation_from(vecObserverSun)


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.