#Copyright (C) 2008 Bernhard W. Adams #This file is free software: you can redistribute it and/or modify #it under the terms of the GNU General Public License as published by #the Free Software Foundation, either version 3 of the License, or #(at your option) any later version. #This file is distributed in the hope that it will be useful, #but WITHOUT ANY WARRANTY; without even the implied warranty of #MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the #GNU General Public License for more details. #You should have received a copy of the GNU General Public License #along with this program. If not, see . #depth of nested defines: 3 #define degrees Pi/180 #define Pi 2*asin(1) #define AU 1.49597870691e8*km #define day 24*3600*sec # this program uses km and seconds as units #define km 1 #define km^2 1 #define km^3 1 #define sec 1 #define sec^2 1 #define Cy (3.15576e9*sec) #define J2000 2451545.0 --orbits[]: "Earth", "Apophis", "Venus", "Mars", "Mercury" orbits[]: "Earth", "Apophis", "E_A0" --orbits[]: "Earth", "Mars", "2007WD5", "Apophis" --orbits[]: "Earth" --todo: "plot trajectories, plot relative trajectories" --todo: "orbital data" todo: "transfer trajectory, plot trajectories, plot relative trajectories" --todo: "scan Delta v, plot trajectories" # #http://ssd.jpl.nasa.gov/txt/p_elem_t1.txt #http://ssd.jpl.nasa.gov/txt/aprx_pos_planets.pdf output unit of time: day output unit of length: km output unit of velocity: km/sec output dt of velocity vectors: 30*day output, draw velocity vectors: "yes" grav. param. mu: 1.32712440018e11*km^3/sec^2 maximum Delta v: 100*km/sec #################################### orbit E_A0, name: "trans_E_A_0" # the parameters "from" and "to" refer to orbit names "orbit xxx, name" orbit E_A0, from: "Earth" orbit E_A0, to: "Apophis" orbit E_A0, no. of loops: 0 orbit E_A0, start posn file: "E_A_pos0_s" orbit E_A0, end posn file: "E_A_pos0_e" orbit E_A0, traj file: "E_A_traj0" orbit E_A0, orbit file: "E_A_orb0" orbit E_A0, rel start posn file: "E_A_rpos0_s" orbit E_A0, rel end posn file: "E_A_rpos0_e" orbit E_A0, rel traj file: "E_A_rtraj0" orbit E_A0, rel orbit file: "E_A_rorb0" orbit E_A0, Delta v file: "E_A_Dv0" orbit E_A0, tt log file: "se0.log" #################################### orbit E_A-1, name: "trans_E_A_-1" # the parameters "from" and "to" refer to orbit names "orbit xxx, name" orbit E_A-1, from: "Earth" orbit E_A-1, to: "Apophis" orbit E_A-1, no. of loops: -1 orbit E_A-1, start posn file: "E_A_pos-1_s" orbit E_A-1, end posn file: "E_A_pos-1_e" orbit E_A-1, traj file: "E_A_traj-1" orbit E_A-1, orbit file: "E_A_orb-1" orbit E_A-1, rel start posn file:"E_A_rpos-1_s" orbit E_A-1, rel end posn file: "E_A_rpos-1_e" orbit E_A-1, rel traj file: "E_A_rtraj-1" orbit E_A-1, rel orbit file: "E_A_rorb-1" orbit E_A-1, Delta v file: "E_A_Dv-1" orbit E_A-1, tt log file: "se-1.log" ##################################### # # 365.256 363 051 # orbit Earth, name: "Earth" orbit Earth, a: 1.0000001124*AU orbit Earth, da/dt: 5.62e-6*AU/Cy orbit Earth, e: 0.01671123 orbit Earth, de/dt: -4.392e-5/Cy orbit Earth, i: -0.00001531*degrees orbit Earth, di/dt: -0.01294668*degrees/Cy orbit Earth, asc. node: 0.0*degrees orbit Earth, d_an/dt: 0.0*degrees/Cy orbit Earth, epoch: J2000*day --orbit Earth, arg. of periapsis: 102.93768193*degrees orbit Earth, long. of periapsis: 102.93768193*degrees orbit Earth, d_om/dt: 0.32327364*degrees/Cy orbit Earth, L_0 at ref. epoch: 100.46457166*degrees --orbit Earth, M at ref. epoch: 0*degrees --orbit Earth, time of periapsis: 0 orbit Earth, start posn file: "Earth_pos_s" orbit Earth, end posn file: "Earth_pos_e" orbit Earth, traj file: "Earth_traj" orbit Earth, orbit file: "Earth_orb" orbit Earth, rel start posn file:"Earth_rpos0_s" orbit Earth, rel end posn file: "Earth_rpos0_e" orbit Earth, rel traj file: "Earth_rtraj0" orbit Earth, rel orbit file: "Earth_rorb0" ########################################### orbit Venus, name: "Venus" orbit Venus, a: 0.72333566*AU orbit Venus, da/dt: 3.90e-6*AU/Cy orbit Venus, e: 0.00677672 orbit Venus, de/dt: -4.107e-5/Cy orbit Venus, i: 3.39467605*degrees orbit Venus, di/dt: -7.889e-4*degrees/Cy orbit Venus, asc. node: 76.67984255*degrees orbit Venus, d_an/dt: -0.27769418*degrees/Cy orbit Venus, epoch: J2000*day orbit Venus, long. of periapsis: 131.60246718*degrees orbit Venus, d_om/dt: 0.00268329*degrees/Cy orbit Venus, L_0 at ref. epoch: 181.9790995*degrees --orbit Venus, M at ref. epoch: 0*degrees --orbit Venus, time of periapsis: 0 orbit Venus, start posn file: "Venus_pos_s" orbit Venus, end posn file: "Venus_pos_e" orbit Venus, traj file: "Venus_traj" orbit Venus, orbit file: "Venus_orb" ########################################### orbit Mars, name: "Mars" orbit Mars, a: 1.52371034*AU orbit Mars, da/dt: 1.847e-5*AU/Cy orbit Mars, e: 0.09339410 orbit Mars, de/dt: 7.882e-5/Cy orbit Mars, i: 1.84969142*degrees orbit Mars, di/dt: -0.00813131*degrees/Cy orbit Mars, asc. node: 49.55953891*degrees orbit Mars, d_an/dt: -0.29257343*degrees/Cy orbit Mars, epoch: J2000*day orbit Mars, long. of periapsis: -23.94362959*degrees orbit Mars, d_om/dt: 0.44441088*degrees/Cy orbit Mars, L_0 at ref. epoch: -4.55343205*degrees --orbit Mars, M at ref. epoch: 357.52688973*degrees --orbit Mars, time of periapsis: 0 orbit Mars, start posn file: "Mars_pos_s" orbit Mars, end posn file: "Mars_pos_e" orbit Mars, traj file: "Mars_traj" orbit Mars, orbit file: "Mars_orb" # orb. per: 686.98 days ########################################### orbit Mercury, name: "Mercury" orbit Mercury, a: 0.38709927*AU orbit Mercury, da/dt: 3.7e-7*AU/Cy orbit Mercury, e: 0.20563593 orbit Mercury, de/dt: 1.906e-5/Cy orbit Mercury, i: 7.00497902*degrees orbit Mercury, di/dt: -0.00594749*degrees/Cy orbit Mercury, asc. node: 48.33076593*degrees orbit Mercury, d_an/dt: -0.12534081*degrees/Cy orbit Mercury, epoch: J2000*day orbit Mercury, long. of periapsis: 77.45779628*degrees orbit Mercury, d_om/dt: 0.16047689*degrees/Cy orbit Mercury, L_0 at ref. epoch: 252.2503235*degrees --orbit Mercury, M at ref. epoch: 357.52688973*degrees --orbit Mercury, time of periapsis: 0 orbit Mercury, start posn file: "Mercury_pos_s" orbit Mercury, end posn file: "Mercury_pos_e" orbit Mercury, traj file: "Mercury_traj" orbit Mercury, orbit file: "Mercury_orb" ########################################### # # JPL data on apophis, solution of 2006-Sep-01 # Orbital Elements at Epoch 2454200.5 (2007-Apr-10.0) TDB # see also http://newton.dm.unipi.it/cgi-bin/neodys/neoibo?objects:Apophis;main # orbit Apophis, name: "Apophis" orbit Apophis, a: 0.92226141513913*AU orbit Apophis, da/dt: 0*AU/Cy orbit Apophis, e: 0.191059415070228 orbit Apophis, de/dt: 0/Cy orbit Apophis, i: 3.33131464186276*degrees orbit Apophis, di/dt: 0*degrees/Cy orbit Apophis, asc. node: 204.4591523029436*degrees orbit Apophis, d_an/dt: degrees/Cy orbit Apophis, epoch: 2454200.5*day orbit Apophis, arg. of periapsis: 126.3855713119676*degrees orbit Apophis, d_om/dt: 0*degrees/Cy o--rbit Apophis, L_0 at ref. epoch: 307.3630785264554*degrees orbit Apophis, M at ref. epoch: 307.3630785264554*degrees o--rbit Apophis, time of periapsis: 2454247.800688934931*sec # 2007-May-27.30068894 JED orbit Apophis, start posn file: "Apophis_pos_s" orbit Apophis, end posn file: "Apophis_pos_e" orbit Apophis, traj file: "Apophis_traj" orbit Apophis, orbit file: "Apophis_orb" # 2nd orbit, period: 323.503874084533*day # mean motion 2nd orbit, n: 1.112815112396244*degrees/day # perihelion distance 2nd orbit, q: 0.746054688620808*AU # aphelion distance 2nd orbit, Q: 1.098468141657453*AU # #################################### # http://ssd.jpl.nasa.gov/sbdb.cgi?sstr=2007%20WD5&orb=1 orbit 2007WD5, name: "2007WD5" orbit 2007WD5, a: 2.54137995816249*AU orbit 2007WD5, da/dt: 0*AU/Cy orbit 2007WD5, e: 0.602518686382147 orbit 2007WD5, de/dt: 0/Cy orbit 2007WD5, i: 2.37134953892162*degrees orbit 2007WD5, di/dt: 0*degrees/Cy orbit 2007WD5, asc. node: 67.53998738637327*degrees orbit 2007WD5, d_an/dt: 0*degrees/Cy orbit 2007WD5, epoch: 2454200.5*day orbit 2007WD5, arg. of periapsis: 312.6459818181322*degrees orbit 2007WD5, d_om/dt: 0*degrees/Cy --orbit 2007WD5, L_0 at ref. epoch: 0*degrees orbit 2007WD5, M at ref. epoch: 313.3896582797448*degrees --orbit 2007WD5, time of periapsis: "2007/10/18" # tp = 2454392.094313025830 = 2007-Oct-18.59431303 # 2007-May-27.30068894 JED orbit 2007WD5, start posn file: "2007WD5_pos_s" orbit 2007WD5, end posn file: "2007WD5_pos_e" orbit 2007WD5, traj file: "2007WD5_traj" orbit 2007WD5, orbit file: "2007WD5_orb" ########################### --start time: "2012/7/15" --start time: "2012/7/23" start time: 2456084.81999356*day --end time: "2014/2/15" --end time: "2013/4/28" end time: 2456371.07635661*day # 23.07.2012 - 20.02.2013 scan, start time from: "2011/7/1" scan, start time to: "2013/6/30" scan, start time step: 10*day scan, end time from: "2012/7/1" scan, end time to: "2014/6/30" scan, end time step: 10*day no. of points for orbits: 1000 no. of points for trajectories: 1000 ## diagnostics: no #