# # $Id: qtd.hnb,v 1.2 2002/05/08 20:27:55 rauch Exp $ # # Duplicate the orbital integration of the nine planets performed by # Quinn, Tremaine, & Duncan (1991); requires a driver implementing an # extra_kick() representing the Earth-Moon quadrupole as described # in the reference. Note that the mass of Pluto used by QTD (inherited from # DE 102) is 45x higher than its correct value! # [See Null & Owen (1996) for an HST derived value.] # # The values below produce results somewhat more accurate than those of QTD. # IntegCoord: Jacobi Order2 Drift-Kick Corrector: True AngleUnit: rad LengthUnit: AU MassUnit: Msun TimeUnit: d StepSize: -6.6644 Tinit = 2433280.5 M = 1 N = 10 PostNewtonian: True Tfinal = -1113787075.5 # final QTD integration time StateFiles: qtd%d.dat StateInterval: 1116220356 # output only the final values StateOrder: x1 x2 x3 v1 v2 v3 StateCoord: Bodycentric StateDigits: 15 EnergyFile: energy.dat EnergyInterval: 1e5 steps # Planet state vectors for JD 2433280.5 in DE 102 reference frame: # InputOrder: mass x1 x2 x3 v1 v2 v3 # Mercury 1.66013679527193E-7 \ 0.343926450169642 0.0456154799533135 -0.0109240372119325 \ -0.00846653204500986 0.0256146053072818 0.0145868453362396 # Venus 2.44783959796682E-6 \ 0.142965184343246 0.647005066033887 0.28248240066038 \ -0.0198938122793425 0.00311311946611859 0.0026594458057458 # Earth-Moon barycenter 3.04043292359547E-6 \ -0.136364695954795 0.893397922857 0.387458344639337 \ -0.0173199988485296 -0.00224430473176756 -0.000973361115758044 # Mars 3.22714936215393E-7 \ -1.36983397618342 0.843135248017904 0.42383290661143 \ -0.00738456383127117 -0.0094773586392742 -0.00415165513666213 # Jupiter 9.54786104043042E-4 \ 3.34936422369601 -3.47376144901258 -1.5721496863938 \ 0.00558564314958231 0.00496226113722251 0.00199227692673937 # Saturn 2.85836787194512E-4 \ -8.97250506828211 2.27968200813286 1.33033860971146 \ -0.00185825195699671 -0.00498385858141744 -0.00198025741280725 # Uranus 4.37273164545892E-5 \ -1.00300399532732 17.3235084732718 7.60482504641591 \ -0.00395525416301772 -0.000375913785112941 -0.000108849991287794 # Neptune 5.17759138448794E-5 \ -29.1945807386112 -7.71928519199847 -2.42724537828877 \ 0.000820748057818085 -0.00277209825958023 -0.00115611603043592 # Pluto 3.33333333333333E-7 \ -26.2336507820155 20.5619754200559 14.4445571277807 \ -0.00131588869828116 -0.00262012820549352 -0.000427083355026298