Semi major axis and inclination corresponding to (N,P,Q)
[sma, inc] = CL_op_repeat2smaInc(N, P, Q, ecc, sso [[, incInput, er, mu, j2, rotr_pla, rotr_pla_sun]])
Computes the semi major axis (and inclination for sunsynchronous orbits) of repeat (also called phased) orbits with given repeat parameters N, P, Q.
The 3 integers N, P, Q are such that:
Tc (duration of repeat cycle) = (N * Q + P) * T = Q * Tr, where T is the orbit period (more exactly the mean nodal period), and Tr is the revolution period of the planet with respect to the orbit's (ascending or descending) node.
Then:
Q (= Tc/Tr) is the number of planet revolutions per repeat cycle.
N+P/Q (= Tr/T) is the number of orbits per planet revolutions with respect to the orbit's (ascending or descending) node.
N*Q+P (= Tc/T) is the (whole) number of orbits per repeat cycle.
Note that Tr depends on the orbit parameters because the node drifts under the effect of J2.
If sso = 1, the orbits are sun-synchronous (the inclination is then computed)
If sso = 0, the function uses incInput as the inclination to be used.
Values of N. [-] (1x1 or 1xNt)
Values of P. [-] (1x1 or 1xNt)
Values of Q. [-] (1x1 or 1xNt)
Eccentricity. [-] (1x1 or 1xNt)
1 = sun synchronous orbits, 0 = fixed inclination. (1x1)
(mandatory in case sso=0, unused otherwise) Inclinations. [rad] (1x1 or 1xNt)
(optional) Equatorial radius. Default: %CL_eqRad. [m] (1x1)
(optional) Gravitational constant. Default: %CL_mu. [m^3/s^2] (1x1)
(optional) Second zonal harmonic (unnormalized). Default: %CL_j1jn(2). [-] (1x1)
(optional) Body rotation rate. Default: %CL_rotrBody. [rad/s] (1x1)
(optional) Mean apparent rotation rate of the Sun around the body. Default: %CL_rotrBodySun. [rad/s] (1x1)
Semi-major axis. 'sma' is '%nan' if value cannot be computed. [m] (1xNt)
Inclination. 'inc' is '%nan' if value cannot be computed. [rad] (1xNt)
CNES - DCT/SB