<< CL_op_gauss Orbit properties CL_op_latSwath >>

celestlab >> Orbit properties > CL_op_incdotSsoCirc

CL_op_incdotSsoCirc

Long-term time derivative of inclination (thirdbody perturbation)

Calling Sequence

[incdot] = CL_op_incdotSsoCirc(mltan, sma [[, inc, mu_b, r_b, inc_b, er, mu, j2, rotr_pla_sun]])

Description

Parameters

mltan:

Mean local time of the ascending node. [h] (1xN or 1x1)

sma:

Mean semi-major axis. [m] (1xN or 1x1)

inc:

(optional) Mean inclination. Default: []. [rad] (1xN or 1x1)

mu_b:

(optional) Gravitational constant of perturbing body. Default: %CL_body.Sun.mu. [m^3/s^2] (1x1)

r_b:

(optional) Distance from central body to perturbing body. Default: %CL_au. [m] (1x1 or 1xN)

inc_b:

(optional) Inclination of orbit of perturbing body. Default: 23.44 deg. [rad] (1x1 or 1xN)

er:

(optional) Equatorial radius. Default: %CL_eqRad. [m] (1x1).

mu:

(optional) Gravitational constant. Default: %CL_mu. [m^3/s^2] (1x1)

j2:

(optional) Second zonal harmonic (unnormalized). Default: %CL_j1jn(2). [-] (1x1)

rotr_pla_sun:

(optional) Mean apparent rotation rate of the Sun around the body. Default: %CL_rotrBodySun. [rad/s] (1x1)

incdot:

Time derivative of inclination. [rad/s] (1xN)

Examples

mltan = 0 : 0.1 : 24; // [hours]
sma = 6378.e3 + 700.e3; // [m]

// Inclination average time derivative [rad/s]
incdot = CL_op_incdotSsoCirc(mltan, sma);

// Plot result [deg/year]
scf();
plot(mltan, CL_rad2deg(incdot) * 86400 * 365.25);

Report an issue
<< CL_op_gauss Orbit properties CL_op_latSwath >>