Drift of orbital elements due to J2
[drift_pom,drift_gom,drift_M[,dpomdaei[,dgomdaei[,dMdaei]]]] = CL_op_driftJ2(sma,ecc,inc[,er[,mu[,j2]]])
Computes:
- the secular drifts (drift_pom=d(pom)/dt, drift_gom=d(gom)/dt, drift_M=dM/dt) on the keplerian elements, where pom=argument of periapsis, gom=right ascension of ascending node,M=mean anomaly, considering the effect of the first zonal harmonic J2,
- the partial derivatives of these drifts with respect to semi major axis, eccentricity and inclination.
(See the formulas, where: R = equatorial radius, n = keplerian mean motion
Semi-major axis [m] (1x1 or 1xN)
Eccentricity (1x1 or 1xN)
Inclination [rad] (1x1 or 1xN)
(optional) Equatorial radius [m] (default is %CL_eqRad)
(optional) Gravitational constant [m^3/s^2] (default value is %CL_mu)
(optional) Second zonal harmonic) (default is %CL_j2)
drift of the argument of periapsis [rad/s] (1xN)
drift of the right ascension of ascending node [rad/s] (1xN)
drift of the mean anomaly [rad/s] (1xN)
partial derivative of drift_pom with respect to sma, ecc and inc (3xN)
partial derivative of drift_gom with respect to sma, ecc and inc (3xN)
partial derivative of drift_M with respect to sma, ecc and inc (3xN)
1) CNES - MSLIB FORTRAN 90, Volume E (me_deriv_secul_J2)
CNES - DCT/SB