# # $Id: jovian.hnb,v 1.1 2002/05/08 22:51:29 rauch Exp $ # # A simple example file integrating the jovian planets, with Pluto # added as a lightweight (and essentially massless) particle. # # The parameters have been chosen to closely match the "example" calculation # in the SWIFT 2.0 package, enabling an easy comparison. # Integrator = Symplectic IntegCoord = Jacobi Order2 Drift-Kick Corrector = True AngleUnit = rad LengthUnit = AU MassUnit = Msun TimeUnit = d StepSize = 365.25 Tinitial = 0 M = 1.00000597682 N = 6 NLWP = 1 InputOrder = Mass x1 x2 x3 v1 v2 v3 # Jupiter 9.547861040430417E-4 \ 3.409530427945 3.635870038323 3.424028779974999E-02 \ -5.604670182013000E-03 5.524493219597000E-03 -2.663981907246999E-06 # Saturn 2.855837331505598E-4 \ 6.612079829705 6.386934883415 -.1361443021015 \ -4.180230691976999E-03 4.003576742277000E-03 1.672376407859000E-05 # Uranus 4.372731645458918E-5 \ 11.16769742622999 16.04343604329 .3617849409932999 \ -3.265538550416999E-03 2.070723353855000E-03 -2.176677219660999E-05 # Neptune 5.177591384487935E-5 \ -30.17366005485 1.917641287545 -.1538859339981 \ -2.241622739519000E-04 -3.107271885462999E-03 3.583760257057000E-05 # Pluto LWP = 1e-30 \ -21.381833467487 32.07799861155299 2.4924585571843 \ -1.776056261431200E-03 -2.060870159021399E-03 6.580950635152800E-04 # # Output options. # Tfinal = 365.25e6 # 1 Myr OutputInterval = 365.25e5 # 0.1 Myr OutputFiles = plan%d.dat OutputOrder = time semi ecc incl longasc argperi meananom OutputCoord = Bodycentric OutputDigits = 16 SaveInterval = 365.25e5 # 0.1 Myr SaveFiles = save.dat EnergyInterval = 365.25e5 # 0.1 Myr EnergyFile = energy.dat