CL_gm_intersectPlanes

Description

• Computes the argument of latitude at the intersection of 2 orbit planes.

  The orbit planes are defined by their inclinations and right ascensions of ascending node

• Notes:

  - 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.

  - 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.

Parameters

inc1:

Inclination of orbit 1 [rad] (1xN or 1x1)

raan1:

Right ascension of ascending node of orbit 1 [rad] (1xN or 1x1)

inc2:

Inclination of orbit 2 [rad] (1xN or 1x1)

raan2:

Right ascension of ascending node of orbit 2 [rad] (1xN or 1x1)

pso1:

Argument of latitude in orbit 1 where the 2 planes intersect [rad] (1xN)

pso2:

Argument of latitude in orbit 2 where the 2 planes intersect [rad] (1xN)

inters:

Flag indicating if the planes intersect (inters = 1) or is there exists an infinity of intersections (inters=0) (1xN)

Authors

• CNES - DCT/SB

Examples

inc1 = 1; // rad raan1 = 2; inc2 = 0.1; raan2 = 1.9; [pso1,pso2] = CL_gm_intersectPlanes(inc1,raan1,inc2,raan2)