Intersection of 2 orbit planes
[pso1, pso2, inters] = CL_gm_intersectPlanes(inc1, raan1, inc2, raan2)
Computes the (true) argument of latitude at the intersection of 2 orbit planes.
The orbit planes are defined by their inclinations and right ascensions of the ascending node.
By definition the argument of latitude that is computed is the angle between the line from the (implicit) frame center to the ascending node, and the line defined by the intersection of the 2 planes.
The two output arguments are:
- pso1: true argument of latitude in plane 1 (from ascending node 1).
- pso2: true argument of latitude in plane 2 (from ascending node 2).
Only one solution is computed such that pso1 and pso2 belong to [0, pi]. The arguments of latitude for the second solution are pso1+pi and pso2+pi respectively.
Notes:
- If the orbits are circular, then pso1 and pso2 define the positions in the orbits where the orbit paths intersect.
- If the planes are identical, the number of intersections is infinite. The output argument (inters) is then set to 0, otherwise it is equal to 1.
Inclination of orbit 1 [rad] (1xN or 1x1)
Right ascension of the ascending node of orbit 1 [rad] (1xN or 1x1)
Inclination of orbit 2 [rad] (1xN or 1x1)
Right ascension of the ascending node of orbit 2 [rad] (1xN or 1x1)
True argument of latitude in orbit 1 where the 2 planes intersect [rad] (1xN)
True Argument of latitude in orbit 2 where the 2 planes intersect [rad] (1xN)
Flag indicating if the planes intersect (inters = 1) or is there exists an infinity of intersections (inters=0) (1xN)
CNES - DCT/SB