Examples
pos = [[7000.e3; 0; 0], [0; 7000.e3; 0], [0; 0; 7000.e3]]; CL_fo_zonHarmPot(pos) // all harmonics CL_fo_zonHarmPot(pos, 2:6) CL_fo_zonHarmPot(pos, [2,4,6]) |  |  |
Potential due to gravity (zonal harmonics)
[pot] = CL_fo_zonHarmPot(pos, [nz, er, mu, j1jn])
Potential due to zonal harmonics.
By definition, the acceleration derived from the potential is +grad(potential).
Notes:
- nz is an array of the zonal term numbers to consider.
- The origin for the position vector must be the central body.
- The coordinates frame can be any frame whose Z axis is the revolution axis of the perturbation field.
- If nz contains numbers in the range [2,6], explicit formulas are used (faster).
- If nz is equal to -1 or "all", all the harmonics are used.
See Force models for more details.
Position vector from the central body [m]. (3xN)
(optional, integer) Zonal term numbers to considerer (1xP) or "all" or -1
(optional) Equatorial radius [m]. Default value is %CL_eqRad.
(optional) Gravitational constant [m^3/s^2]. Default value is %CL_mu.
(optional) Zonal harmonics. Default value is %CL_j1jn.
Potential [m^2/s^2]. (1xN)
CNES - DCT/SB
D. Vallado - Fundamentals of Astrodynamics and Applications, 2nd edition
pos = [[7000.e3; 0; 0], [0; 7000.e3; 0], [0; 0; 7000.e3]]; CL_fo_zonHarmPot(pos) // all harmonics CL_fo_zonHarmPot(pos, 2:6) CL_fo_zonHarmPot(pos, [2,4,6]) | | |