CL_mod_precessionAngles — Precession angles (Lieske 1976)
[K[,KP[,KPP]]]=CL_mod_precessionAngles(jj_tai,sec_tai[,flag])
if P0 is the pole of J2000 equator, and P is the pole of the equator at date t; if G0 is the equinox of J2000 and G is the equinox at date t:
zetaa : angle between the meridians (pole P0) P0-P and P0-G0
thetaa: angle between the equateur J2000 and the equateur t
za : angle between the meridians (pole P) P-P0 and P-G
CNES Julian day (TAI time) (1xN)
seconds (TAI time) (1xN)
(optional) "n" -> K; "p" -> K and Kp; "s" -> K, Kp, Kpp (default depends on number of outputs)
precession angles K = [zetaa;tetaa;za] [rad] (3xN)
(optional) first time derivative of K [rad/s] (3xN)
(optional) second time derivative of K [rad/s^2] (3xN)