Orbits of Milky Way satellites in an ensemble of Galactic potentials including the LMC

This archive contains the samples from the MCMC analysis of Milky Way potential,
including the dynamical perturbation from the LMC, and corresponding orbits
of other Galactic satellites in each choice of the MW+LMC potential.

There are 1000 samples of the MW+LMC potential parameters from the chain,
and for each of them, one sample from the posterior distribution of present-day
position/velocity for each satellite (i.e., sampled from its measurement
uncertainties and weighted by the probability of finding this phase-space point
in the DF corresponding to the given potential in the chain).
The archive contains pre-computed trajectories for 63 objects,
stored as a 4d numpy array  orbits.npy  with shape
(63 objects, 1000 samples, 151 timesteps, 6 phase-space coordinates).
The timesteps are equally spaced between -3 Gyr and now (also stored in
orbit_times.npy).
Object names are listed in  names.npy;  LMC comes 0th.
The potential parameters for each of these 1000 samples are stored in
potential_params.npy - each row has 6 parameters, 5 for the MW halo
and the last one is the LMC mass.
The script  integrate_orbits.py  contains a routine for constructing
the MW potential with the given parameters (taken from the chain),
computing the past trajectories of MW+LMC and constructing the time-dependent
potential of both galaxies, which can then be used to integrate orbits of
test particles, such as other satellites. Doing this for all 1000 samples
and 62 objects would take some time, that's why they are provided in already
pre-computed form.