============================================================= ======================== Solar-Sail ========================= ============================================================= Each of the subdirectories, Sun-Jupiter, Sun-Uranus, and Sun-Earth, has the following further subdirectories: Libration: Determines libration points, loci of libration points, unstable manifolds of libration points, and connecting orbits from libration points to (typically) tori. Orbits: Computes families of periodic orbits Manifolds: Computes the unstable manifold of periodic orbits The solar sail equations defined in directory "Equations" are used in all subdirectories. NOTE: This demo is subject to further verification. Comments are welcome. ============================================================= ==================== AUTO-Parameter use ===================== ============================================================= PAR(1) unfolding parameter lambda PAR(2) mass-ratio mu PAR(3) solar sail parameter alpha PAR(4) solar sail parameter beta PAR(5) solar sail parameter delta PAR(6) Jacobi constant (when alpha and delta are zero) PAR(9) eps: the initial point of the orbits that make up the manifold lies along the line X + eps V in 6D PAR(10) log(mu): the natural logarithm of mu, used as continuation parameter instead of mu itself PAR(11) period of periodic orbits PAR(12) length of manifold orbits PAR(13) integration time of manifold orbits PAR(14) norm of Floquet eigenfunction: used in Manifolds/flq.f90 PAR(15) largest real positive eigenvalue or Floquet multiplier (zero if there is no positive eigenvalue or multiplier) PAR(16) next largest real positive eigenvalue or multiplier PAR(17-19) period(s) T1, T2, T3, corresponding to purely imaginary eigenvalues (if any) PAR(20) log_10(fmx), where fmx is the magnitude of the largest Floquet multiplier: used in Libration/Connection1/man.f90 PAR(21) x-coordinate of libration point (and max |x(t)| in "Orbits") PAR(22) y-coordinate of libration point (but max |y(t)| in "Orbits") PAR(23) z-coordinate of libration point (but max |t(t)| in "Orbits") PAR(24) max |x'(t)| in "Orbits" PAR(25) max |y'(t)| in "Orbits" PAR(26) max |z'(t)| in "Orbits" PAR(27) x-coordinate of the end point of manifold orbits PAR(28) y-coordinate of the end point of manifold orbits PAR(29) z-coordinate of the end point of manifold orbits Used in Manifolds/man.f90 and Manifolds/data.py: PAR(31-36) starting point X on the base periodic orbit PAR(41-46) direction vector V at the starting point X =============================================================