###################################################################### # # Sample HNBody Integration Input File. # # This file sets up a simple 3-body integration. Most of the file consists # of explanatory comments describing the most important options. Options that # have been commented out have the default values listed and can be omitted if # the default is satisfactory; they are included here to orient the user as # to what is available. The file does not include every possible option; the # complete list of options is described in the options document. # # # # INTEGRATION PROPERTIES # #Integrator: Symplectic # choices are: Symplectic [Default], Bulirsh-Stoer, Runge-Kutta # for B-S and R-K there is an optional second value, the accuracy to use # (default: 1e-12) # # The following three options are meaningful only with Symplectic integrations; # they specify the canonical coordinates and splitting to use, whether to # utilize a symplectic corrector, and whether to follow close encounters: # #IntegCoord: Jacobi Order2 Kick-Drift # primary values are: Jacobi (or Barycentric) [Default], # RegJacobi (or RegBarycentric), Bodycentric # secondary values are: Order2 [Default], Order4 # tertiary values are: Kick-Drift [Default], Drift-Kick, Shift-Drift # #Corrector: False # choices are: false/no [Default], true/yes, Order2, Order4, Order6 # (true or yes means to use the highest available order) # # PROPERTIES OF INITIAL CONDITIONS # The following options are all required; units can contain numbers # and metric prefixes: AngleUnit: deg # base unit choices: deg, rad LengthUnit: 5.2*AU # base unit choices: m, AU (or ua), pc MassUnit: 2e27 kg # base unit choices: g, Msun TimeUnit: yr # base unit choices: s, h, d, yr StepSize: 0.01 # integration time step (units of TimeUnit) M = 1 # mass of the dominant object (units of MassUnit) N = 3 # number of objects (including the dominant mass) #OblateJ2 = 0 # J2 parameter of dominant mass (default is 0). #OblateRadius = 0 # J2 radius of dominant mass (default is 0). #PostNewtonian: No # Include leading-order post-Newtonian corrections? # (default is false/no) # INITIAL CONDITIONS # Specifying the initial conditions for each particle requires seven # numbers, its mass and six independent phase space coordinates. The # ordering and meaning of the seven required values is specified by # the InputOrder option, below. Units must be consistent with the units # declared above. If the Encounters option is enabled (see above), an # encounter radius (tag: EncRadius) must also be included in the initial # conditions. There are two basic input coordinate choices: # # Elements are body-centered on the dominant mass; the dominant mass must # NOT be specified since elements are relative to its local rest frame. # The osculating orbital elements can be defined in numerous equivalent # representations; any six can be specified in InputOrder, so long as # they are independent. # # Cartesian coordinates require that particle positions and velocities be # specified in some inertial reference frame. If the dominant mass is # not specified, it is assumed to be at rest at the origin; if it is # specified, its mass must match that given by the M option, above. # InputOrder: Mass SemiMajorAxis Eccentricity Inclination \ LongAscendNode ArgPeriapse MeanAnomaly # # Tags ignore case and can be abbreviated (if unique); # for example, you could write the above as: # #INPUTORD= mass semi ecc incl longasc argperi meananom # # alternate choices for elements: # PeriDistance, LongPeriapse, TrueLongitude, MeanLongitude, # TrueLatitude, MeanLatitude, TrueAnomaly, (Epoch, TimePeriapse)** # Peri means periapse, Long means longitude, and Arg means argument. # ** ( , ) means both are required if either is specified. # # The following would be a typical Cartesian example: # #InputOrder: x1 x2 x3 v1 v2 v3 mass # # InputOrder must also include EncRadius when Encounters are enabled. # # Here we specify the initial conditions: # #ParticleType: Default # choices are: HWPs [=Default], LWPs, ZWPs # (sets implied particle type for subsequent particle specifications) 0.001 5.2 0.0 0.0 0.0 0.0 0.0 0.001 6.0 0.0 0.0 0.0 0.0 0.0 # # The above could also be written as: HWP = 0.001 5.2 ... etc. # # To include initial conditions (or anything else) from a file or files, use: # #Include: thisfile.hnb thatfile.hnb ... # # Output options are not required for initialization; the include mechanism # makes it easy to maintain them separately: # include: output.hnb # # $Id: sample.hnb,v 1.2 2003/08/04 15:53:00 rauch Exp $ #