CL_op_driftJ2

Description

• Computes:

  - the secular drifts (pomdot=d(pom)/dt, gomdot=d(gom)/dt, anmdot=d(anm)/dt) on the keplerian elements, where pom=argument of periapsis, gom=right ascension of ascending node,anm=mean anomaly, considering the effects of the central (Keplerian) force and 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

  

Parameters

sma:

Semi-major axis [m] (1x1 or 1xN)

ecc:

Eccentricity (1x1 or 1xN)

inc:

Inclination [rad] (1x1 or 1xN)

er:

(optional) Equatorial radius [m] (default is %CL_eqRad)

mu:

(optional) Gravitational constant [m^3/s^2] (default value is %CL_mu)

j2:

(optional) Second zonal harmonic) (default is %CL_j1jn(2))

pomdot:

Mean time derivative of argument of periapsis [rad/s] (1xN)

gomdot:

Mean time derivative of right ascension of ascending node [rad/s] (1xN)

anmdot:

Mean time derivative of mean anomaly [rad/s] (1xN)

dpomdotdaei:

Partial derivative of pomdot with respect to sma, ecc and inc (3xN)

dgomdotdaei:

Partial derivative of gomdot with respect to sma, ecc and inc (3xN)

danmdotdaei:

Partial derivative of anmdot with respect to sma, ecc and inc (3xN)

Bibliography

• 1) CNES - MSLIB FORTRAN 90, Volume E (me_deriv_secul_J2)

Authors

• CNES - DCT/SB

Examples

sma = [42000.e3, 7000.e3]; ecc = [0.001, 0.001]; inc = CL_deg2rad([51.6, 51.6]); [pomdot,gomdot,anmdot] = CL_op_driftJ2(sma,ecc,inc)