From 67ba5a6ac6a052fd287cd5ff6850c6ba034057c0 Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 15:47:55 -0400 Subject: [PATCH] Added new methods for adding rotation rates of planets in initial conditions files --- examples/symba_chambers_2013/.idea/.gitignore | 3 + examples/symba_chambers_2013/init_cond.py | 56 ++ examples/symba_chambers_2013/param.in | 69 +- .../symba_chambers_2013/pl_chambers_2013.in | 944 +----------------- examples/symba_chambers_2013/sun_MsunAUYR.in | 18 +- python/swiftest/swiftest/constants.py | 1 + python/swiftest/swiftest/init_cond.py | 69 +- python/swiftest/swiftest/io.py | 34 +- python/swiftest/swiftest/simulation_class.py | 10 +- 9 files changed, 215 insertions(+), 989 deletions(-) create mode 100644 examples/symba_chambers_2013/.idea/.gitignore create mode 100755 examples/symba_chambers_2013/init_cond.py diff --git a/examples/symba_chambers_2013/.idea/.gitignore b/examples/symba_chambers_2013/.idea/.gitignore new file mode 100644 index 000000000..26d33521a --- /dev/null +++ b/examples/symba_chambers_2013/.idea/.gitignore @@ -0,0 +1,3 @@ +# Default ignored files +/shelf/ +/workspace.xml diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py new file mode 100755 index 000000000..5cda5056e --- /dev/null +++ b/examples/symba_chambers_2013/init_cond.py @@ -0,0 +1,56 @@ +#!/usr/bin/env python3 +import swiftest +import numpy as np +import os + +# Initialize simulation object +sim = swiftest.Simulation() + +# Set unit conversion factors +MU2KG = swiftest.MSun +TU2S = swiftest.YR2S +DU2M = swiftest.AU2M +sim.param['MU2KG'] = MU2KG +sim.param['TU2S'] = TU2S +sim.param['DU2M'] = DU2M + +# Simulation time parameters +sim.param['T0'] = 0.0 +sim.param['TSTOP'] = 300e6 +sim.param['DT'] = 6 * swiftest.JD2S / sim.param['TU2S'] +t_print = 1000.0 +iout = int(np.ceil(t_print / sim.param['DT'])) +sim.param['ISTEP_OUT'] = iout +sim.param['ISTEP_DUMP'] = iout + +# Optional output file names +sim.param['PARTICLE_OUT'] = "particle.dat" +sim.param['ENERGY'] = "YES" +sim.param['ENERGY_OUT'] = "energy.dat" +sim.param['PL_IN'] = "pl_chambers_2013.in" +sim.param['CB_IN'] = "sun_MsunAUYR.in" + +# Simulation parameters +sim.param['FRAGMENTATION'] = "YES" +sim.param['ROTATION'] = "YES" +sim.param['CHK_RMAX'] = 1000.0 +sim.param['CHK_EJECT'] = 1000.0 + +# Add central body +sim.add("Sun") +sim.add("Earth") + +# Add bodies described in Chambers (2013) Sec. 2.1, with the uniform spatial distribution and two bodies sizes (big and small) +Nb = 14 +Ns = 140 +Mb = 2.8e-7 +Ms = 2.8e-8 +dens = 3000.0 / (MU2KG / DU2M**3) +Rb = (3 * Mb / (4 * np.pi * dens) )**(1.0 / 3.0) +Rs = (3 * Ms / (4 * np.pi * dens) )**(1.0 / 3.0) + + +sim.save("param.in") + + + diff --git a/examples/symba_chambers_2013/param.in b/examples/symba_chambers_2013/param.in index bf59cfb97..d096d13fe 100644 --- a/examples/symba_chambers_2013/param.in +++ b/examples/symba_chambers_2013/param.in @@ -1,32 +1,37 @@ -! -! Parameter file for Chambers 2013 in units of Msun, AU, year -! -T0 0.0e0 -TSTOP 1e8 ! simulation length in years -DT 0.016 ! stepsize in years -CB_IN sun_MsunAUYR.in -PL_IN pl_chambers_2013.in -TP_IN tp.in -IN_TYPE ASCII -ISTEP_OUT 6250 ! output cadence -ISTEP_DUMP 6250 ! system dump cadence -BIN_OUT bin.dat -PARTICLE_OUT particle.dat -OUT_TYPE REAL8 ! double precision real output -OUT_FORM EL ! osculating element output -OUT_STAT REPLACE -CHK_CLOSE yes ! check for planetary close encounters -CHK_RMAX 100000.0 ! discard outside of -EXTRA_FORCE no ! no extra user-defined forces -BIG_DISCARD no ! output all planets if anything discarded -RHILL_PRESENT yes ! Hill's sphere radii in input file -MU2KG 1.98847e30 ! (M_sun-> kg) -DU2M 1.495979e11 ! distance unit to meters (AU --> m) -TU2S 3.15569259747e7 ! time unit to seconds (years --> seconds) -GMTINY 1e-10 -ENERGY yes -ENERGY_OUT energy.dat -ROTATION yes -FRAGMENTATION yes -DISCARD_OUT discard.out -SEED 8 12261555 871132 92734722 21132443 36344777 4334443 219291656 3848566 +! VERSION Swiftest parameter input +T0 0.0 +TSTOP 300000000.0 +DT 0.01642710472279261 +ISTEP_OUT 60876 +ISTEP_DUMP 60876 +OUT_FORM EL +OUT_TYPE REAL8 +OUT_STAT REPLACE +IN_TYPE ASCII +PL_IN pl_chambers_2013.in +TP_IN tp.in +CB_IN sun_MsunAUYR.in +BIN_OUT bin.dat +CHK_QMIN -1.0 +CHK_RMIN -1.0 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE -1.0 -1.0 +MU2KG 1.988409870698051e+30 +TU2S 31557600.0 +DU2M 149597870700.0 +IN_FORM XV +ENC_OUT enc.dat +EXTRA_FORCE NO +DISCARD_OUT discard.out +PARTICLE_OUT particle.dat +BIG_DISCARD NO +CHK_CLOSE YES +RHILL_PRESENT YES +FRAGMENTATION YES +ROTATION YES +TIDES NO +ENERGY YES +GR NO +ENERGY_OUT energy.dat diff --git a/examples/symba_chambers_2013/pl_chambers_2013.in b/examples/symba_chambers_2013/pl_chambers_2013.in index 1ba4c845e..56d3f87e4 100644 --- a/examples/symba_chambers_2013/pl_chambers_2013.in +++ b/examples/symba_chambers_2013/pl_chambers_2013.in @@ -1,937 +1,7 @@ -156 ! Solar System in unit system AU, M_sun, and years -500 3.76722155e-02 1.08899907e-02 ! ID / G*Mass / Rhill Jupiter -4.67326079e-04 ! Radius -3.46016026e+00 -3.69444209e+00 -6.20863032e-02 ! x y z -1.97779378e+00 2.01435697e+00 -5.26023669e-02 ! vx vy vz -0.0 0.0 0.2756 ! ip --81.01800465 -2388.16786145 5008.65350675 ! rot -600 1.12798314e-02 1.99069165e-02 ! ID / G*Mass / Rhill Saturn -3.89256801e-04 ! Radius -5.80057555e+00 -8.09777935e+00 -9.01449584e-02 ! x y z -1.54008846e+00 1.17938756e+00 -8.18949214e-02 ! vx vy vz -0.0 0.0 0.22 ! ip -445.81155409 381.95827345 5182.3159189 ! rot -4 1.10539569e-05 2.76930511e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius -2.72988779e-01 1.29092883e+00 1.33913146e-01 ! x y z --5.30347096e+00 1.10519810e+00 4.11733597e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -5 1.10539569e-05 2.47028309e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius --7.52767237e-01 -7.89903716e-01 4.41780109e-01 ! x y z -3.97724375e+00 -4.13307219e+00 -7.22805540e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -6 1.10539569e-05 1.32328966e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius -6.20377608e-01 -7.87272103e-02 7.12067297e-02 ! x y z -1.02503121e+00 7.84536983e+00 -1.81237075e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -7 1.10539569e-05 9.44674456e-04 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius -6.98175547e-02 -4.40604282e-01 -6.26515066e-02 ! x y z -9.16286233e+00 1.29526313e+00 1.28292031e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -8 1.10539569e-05 2.10455183e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius --7.37344424e-01 -6.59393259e-01 -3.39334297e-02 ! x y z -4.23617128e+00 -4.72560240e+00 1.90202709e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -9 1.10539569e-05 7.93127033e-04 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius --3.72282885e-01 1.00137160e-02 4.74495877e-02 ! x y z --5.53222723e-02 -1.01517333e+01 1.56417032e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -10 1.10539569e-05 3.72890834e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius --1.71148151e+00 4.03705975e-02 -4.41629664e-01 ! x y z --5.30815676e-01 -4.44829638e+00 1.51595205e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -11 1.10539569e-05 1.24950838e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius -2.86212864e-01 5.18828229e-01 -9.65137375e-03 ! x y z --7.11468376e+00 3.99045778e+00 4.08951506e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -12 1.10539569e-05 1.29201140e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius -6.01579717e-01 -1.24985784e-01 -6.68579475e-02 ! x y z -1.03137581e+00 6.97294287e+00 -3.70526296e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -13 1.10539569e-05 8.45210453e-04 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius -1.66504352e-01 3.65136786e-01 -2.37270588e-02 ! x y z --9.01075623e+00 4.10307061e+00 -1.73991242e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -14 1.10539569e-05 1.80973111e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius --1.00013418e-01 8.46049108e-01 1.29409189e-01 ! x y z --6.69397403e+00 -6.88216193e-01 -6.43089317e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -15 1.10539569e-05 6.77162299e-04 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius --6.64074759e-02 -2.85943101e-01 -1.33132689e-01 ! x y z -1.04794195e+01 -3.18219029e+00 1.50287735e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -16 1.10539569e-05 2.22462802e-03 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius -6.93683497e-01 -7.97060245e-01 -1.56701341e-03 ! x y z -4.27957357e+00 3.72233196e+00 -2.27779518e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -17 1.10539569e-05 9.07412459e-04 ! ID / G*Mass / Rhill -2.36555201e-05 ! Radius --4.30243754e-01 3.92163409e-02 1.61581635e-02 ! x y z --8.96241375e-01 -9.13377010e+00 -2.60722121e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -18 1.10539569e-06 2.16228637e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.49171166e-01 9.89231336e-01 9.47431817e-02 ! x y z --5.97909060e+00 -1.61462871e+00 5.17955232e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -19 1.10539569e-06 2.44552162e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.13055691e+00 -2.73349719e-01 -8.64446048e-04 ! x y z -1.35697138e+00 5.66139385e+00 -8.67279260e-03 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -20 1.10539569e-06 2.69931846e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --7.95952391e-01 8.54202216e-01 5.12437562e-01 ! x y z --4.18257122e+00 -3.66903213e+00 -4.16936339e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -21 1.10539569e-06 3.64709598e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.50222087e+00 2.48717046e-01 -8.19081229e-01 ! x y z --8.04975551e-01 -4.71387035e+00 7.84362176e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -22 1.10539569e-06 1.17536449e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.44682948e-01 -4.96503483e-01 -8.53414160e-02 ! x y z -7.40695872e+00 -3.28819700e+00 -2.14066023e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -23 1.10539569e-06 1.50460742e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.36805798e-01 -5.86582522e-01 3.17351611e-01 ! x y z -6.83507157e+00 3.05184627e+00 5.35881505e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -24 1.10539569e-06 2.39048166e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.80090550e-01 1.00147849e+00 -3.41420172e-01 ! x y z --5.59868356e+00 -1.78046309e+00 9.70115472e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -25 1.10539569e-06 2.01027356e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --7.39485525e-01 6.02196012e-01 2.65945142e-02 ! x y z --4.01389113e+00 -4.93502846e+00 -9.72482631e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -26 1.10539569e-06 3.44715330e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --4.30456826e-01 1.57570303e+00 -1.09595941e-01 ! x y z --4.64706737e+00 -1.19846831e+00 1.04119454e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -27 1.10539569e-06 1.76967616e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --7.35819630e-01 -4.18940068e-01 -5.26182368e-02 ! x y z -3.13485099e+00 -5.73938205e+00 1.82527292e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -28 1.10539569e-06 4.06342312e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --7.54229262e-01 -1.76044238e+00 -1.80128626e-01 ! x y z -4.18274269e+00 -1.75246993e+00 -1.48893698e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -29 1.10539569e-06 6.51392241e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --6.97001957e-02 -2.83866631e-01 1.02094271e-01 ! x y z -1.06859367e+01 -3.16734146e+00 -1.78495005e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -30 1.10539569e-06 1.50723858e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.50547968e-01 5.85483214e-01 -2.02242303e-01 ! x y z --6.51912702e+00 -3.31902976e+00 1.50673282e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -31 1.10539569e-06 9.91191586e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -5.34148347e-02 -4.67706829e-01 -4.50777179e-04 ! x y z -9.07741521e+00 1.03526094e+00 6.27189957e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -32 1.10539569e-06 1.27434773e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.76081124e-01 -1.72561939e-01 2.32546498e-02 ! x y z -2.26640367e+00 -7.48739107e+00 2.18424337e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -33 1.10539569e-06 1.88899665e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.21441400e-01 8.67808697e-01 6.09525334e-02 ! x y z --6.38752841e+00 1.69289574e+00 -5.34224727e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -34 1.10539569e-06 2.51153352e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -6.27618072e-01 -9.85796209e-01 -2.48455975e-01 ! x y z -4.57388053e+00 3.25533389e+00 -1.21335726e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -35 1.10539569e-06 1.31898109e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.04370616e-01 5.13132374e-01 1.83051491e-01 ! x y z --6.94545994e+00 3.58830246e+00 1.54416575e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -36 1.10539569e-06 1.31747151e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.72209409e-01 5.58609144e-01 -5.24341572e-02 ! x y z --7.11229047e+00 -3.50089370e+00 -8.23500441e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -37 1.10539569e-06 8.04475879e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.84218244e-01 3.32444349e-01 -3.37468460e-02 ! x y z --8.87199586e+00 4.71912663e+00 -1.61815818e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -38 1.10539569e-06 8.53415964e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.71574947e-01 2.99863851e-01 2.77485703e-02 ! x y z --7.27373701e+00 6.59527467e+00 -9.50013137e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -39 1.10539569e-06 1.54874159e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.42177998e-01 4.64346195e-01 -1.77247174e-01 ! x y z --4.32593022e+00 -5.71371010e+00 -1.52200642e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -40 1.10539569e-06 1.11811225e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.37513378e-01 4.78960335e-01 -4.18799485e-02 ! x y z --7.64314149e+00 3.80464264e+00 -4.40507044e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -41 1.10539569e-06 6.46344589e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.17502253e-01 2.74368057e-01 6.52591009e-02 ! x y z --1.04985621e+01 -4.38785797e+00 -5.95356835e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -42 1.10539569e-06 3.14895117e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.08577487e+00 -1.03445606e+00 -1.55590020e-01 ! x y z -3.53014787e+00 3.59960362e+00 7.43100049e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -43 1.10539569e-06 3.30212225e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.85350014e-01 -1.34034680e+00 7.46972274e-01 ! x y z -4.81219105e+00 1.19508883e+00 -3.85665415e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -44 1.10539569e-06 7.23534141e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.13136821e-02 3.42286865e-01 -3.05898632e-02 ! x y z --1.01368906e+01 -6.23201390e-01 -3.40648593e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -45 1.10539569e-06 2.39770005e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.04671093e+00 1.11212978e-01 -4.33163084e-01 ! x y z --3.60007968e-01 -5.84725364e+00 -6.15822401e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -46 1.10539569e-06 3.11101805e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.28031558e+00 6.71820187e-01 3.13521220e-01 ! x y z --2.10656104e+00 -4.54924436e+00 1.23122282e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -47 1.10539569e-06 9.97006532e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.03695384e-01 -3.60967985e-01 5.51306924e-02 ! x y z -6.85056674e+00 -5.91778788e+00 -9.60485047e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -48 1.10539569e-06 7.74544738e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.54176241e-01 6.02822917e-02 7.63020359e-02 ! x y z --1.61729423e+00 -1.02381855e+01 4.69477979e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -49 1.10539569e-06 2.60533414e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -4.70647812e-01 -1.14442870e+00 -2.16624680e-04 ! x y z -5.22089775e+00 2.15396431e+00 6.97355593e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -50 1.10539569e-06 1.73280111e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --8.47538292e-02 -7.76819569e-01 -2.57912415e-01 ! x y z -6.69507282e+00 -1.17185683e+00 1.33518772e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -51 1.10539569e-06 3.17790663e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -6.65614001e-01 -1.36216721e+00 3.04841771e-02 ! x y z -4.55635977e+00 2.26593667e+00 1.30678452e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -52 1.10539569e-06 7.23862851e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.31702394e-01 -9.11930665e-02 2.41707377e-02 ! x y z -2.70205684e+00 -1.02322038e+01 -1.45694907e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -53 1.10539569e-06 4.09271910e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.43010053e-01 1.89928569e+00 1.63915030e-01 ! x y z --4.44471826e+00 8.12964188e-01 1.91799605e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -54 1.10539569e-06 9.00526678e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.76279203e-01 3.69578224e-01 -1.37212342e-01 ! x y z --8.58447413e+00 -3.12914501e+00 2.65447958e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -55 1.10539569e-06 1.66875324e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --7.45616929e-01 -2.47231702e-01 -1.19477629e-01 ! x y z -2.19161600e+00 -6.68909269e+00 -1.07794911e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -56 1.10539569e-06 2.61296171e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.21771549e+00 8.22263396e-02 -2.28049506e-01 ! x y z --7.77955309e-01 -5.04420236e+00 2.39397114e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -57 1.10539569e-06 7.84583650e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.35031870e-01 -1.49752491e-01 7.21876006e-02 ! x y z -3.53139295e+00 -9.14063147e+00 -3.02347388e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -58 1.10539569e-06 8.15353611e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.16127731e-01 2.19457827e-01 -2.86640647e-02 ! x y z --5.75118301e+00 -8.30744977e+00 -7.52286441e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -59 1.10539569e-06 9.46502491e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -4.14567897e-01 -1.70158031e-01 -2.06385746e-02 ! x y z -3.52462324e+00 8.64965730e+00 9.71699819e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -60 1.10539569e-06 1.06023870e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.68326520e-01 3.43200640e-01 -3.18358134e-03 ! x y z --5.99833054e+00 6.50593164e+00 -3.52076060e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -61 1.10539569e-06 2.85747597e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.22743812e+00 4.54544048e-01 -3.53405167e-01 ! x y z --2.13063294e+00 4.75488906e+00 -1.41396521e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -62 1.10539569e-06 1.11156710e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --4.90230366e-01 8.52524184e-02 1.64539092e-01 ! x y z --1.33697442e+00 -8.59490437e+00 4.64492986e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -63 1.10539569e-06 7.47979433e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.26924952e-01 1.39246422e-01 -4.37240131e-02 ! x y z --4.02802590e+00 9.60769401e+00 9.33600395e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -64 1.10539569e-06 2.36256862e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.11787086e+00 1.27230662e-01 -2.97194116e-02 ! x y z --6.51392127e-01 -5.87014579e+00 -2.95793081e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -65 1.10539569e-06 1.36083964e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.68540327e-01 1.32853908e-01 -2.74248856e-01 ! x y z --2.28636154e+00 -7.40032660e+00 1.15230861e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -66 1.10539569e-06 2.12644561e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --9.02814420e-01 4.51623081e-01 3.31973730e-02 ! x y z --2.80126142e+00 -5.58560982e+00 -1.92746316e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -67 1.10539569e-06 8.14605985e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.85382182e-01 1.43901530e-03 1.15478430e-02 ! x y z --2.59844799e-02 -1.01271155e+01 4.30703613e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -68 1.10539569e-06 9.27538405e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.61232954e-01 2.20844429e-01 -1.22567174e-01 ! x y z --4.75444340e+00 8.15211658e+00 6.72406103e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -69 1.10539569e-06 1.70813096e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -7.62285263e-01 2.77535790e-01 7.50460161e-03 ! x y z --2.38561781e+00 6.55492046e+00 -5.72763423e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -70 1.10539569e-06 2.71135229e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -6.33966979e-01 1.05145872e+00 -3.79026357e-01 ! x y z --4.80055687e+00 2.37239115e+00 -1.45444659e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -71 1.10539569e-06 9.29779721e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.09086355e-01 -4.25927399e-01 -5.92609052e-02 ! x y z -8.83426715e+00 -1.97333711e+00 -2.57514440e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -72 1.10539569e-06 2.48562722e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.52243615e-01 -1.03472811e+00 -9.20470885e-03 ! x y z -5.14745205e+00 -2.71638935e+00 1.85243774e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -73 1.10539569e-06 1.57390173e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.37564413e-01 -7.02794064e-01 -1.30853933e-01 ! x y z -6.84384400e+00 -2.13579218e+00 -7.74318491e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -74 1.10539569e-06 1.33469063e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -5.24670248e-01 -3.60722160e-01 -5.89339466e-03 ! x y z -4.43972950e+00 6.48102200e+00 9.11331663e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -75 1.10539569e-06 6.47192566e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -5.61797829e-02 -3.04486498e-01 4.89626567e-03 ! x y z -1.10606484e+01 2.04593839e+00 -1.21551157e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -76 1.10539569e-06 1.47162297e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.43061104e-01 -6.84452631e-01 -2.55174546e-02 ! x y z -7.11009943e+00 1.43236802e+00 1.93726736e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -77 1.10539569e-06 8.63668982e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.85992031e-02 -4.08122100e-01 2.78731319e-02 ! x y z -9.76823729e+00 -6.46270469e-01 6.80833278e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -78 1.10539569e-06 2.36804818e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -6.91157068e-01 -7.91937234e-01 3.82540684e-01 ! x y z -4.02153165e+00 4.19828096e+00 1.29880450e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -79 1.10539569e-06 9.14394348e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.51098607e-01 2.59366529e-01 -1.55386862e-02 ! x y z --5.06657998e+00 -7.18889701e+00 -3.53716670e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -80 1.10539569e-06 1.26910472e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -4.23152493e-01 3.88117899e-01 -2.00809903e-01 ! x y z --4.53578552e+00 6.18184215e+00 2.34837457e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -81 1.10539569e-06 6.33003056e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.28512111e-01 2.64697808e-01 -7.13482769e-02 ! x y z --9.33728759e+00 5.41684777e+00 3.59675928e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -82 1.10539569e-06 2.22271119e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.71802215e-01 -9.92089651e-01 5.32987899e-02 ! x y z -5.68468560e+00 2.17000210e+00 -9.37982171e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -83 1.10539569e-06 1.74051421e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -7.86568140e-01 2.32350855e-01 6.38082505e-02 ! x y z --1.88303694e+00 6.63245612e+00 -8.28072975e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -84 1.10539569e-06 9.31517858e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.72567042e-01 -2.35531924e-01 4.96465701e-02 ! x y z -4.70939503e+00 -7.96544948e+00 -1.77358479e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -85 1.10539569e-06 2.46953975e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.20593248e-01 1.08358397e+00 3.52755169e-01 ! x y z --5.53042972e+00 1.53963953e+00 2.87126462e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -86 1.10539569e-06 1.50876657e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.56037714e-01 -4.53117788e-01 -4.40686710e-02 ! x y z -4.32334289e+00 -5.53881012e+00 2.32535969e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -87 1.10539569e-06 6.79814743e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.71412212e-01 -2.72520775e-01 2.76595629e-02 ! x y z -9.35632724e+00 -5.87401170e+00 1.71700376e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -88 1.10539569e-06 1.37728604e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -4.52644996e-01 4.52246640e-01 1.36128602e-01 ! x y z --5.47133577e+00 5.50145396e+00 -3.76988545e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -89 1.10539569e-06 1.75150424e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.17549909e-01 -7.69243942e-01 2.30477642e-01 ! x y z -6.02966869e+00 -2.44216734e+00 -2.26348882e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -90 1.10539569e-06 1.82243275e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -8.53300892e-01 1.05420241e-01 6.10357233e-03 ! x y z --7.23180710e-01 6.14709501e+00 -2.81238589e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -91 1.10539569e-06 1.04051592e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -4.76117727e-02 -4.92462500e-01 -7.89640418e-03 ! x y z -8.70186553e+00 7.97359443e-01 -1.82242405e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -92 1.10539569e-06 9.04059360e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.13222540e-01 -3.73863988e-01 -4.13231963e-04 ! x y z -8.30795424e+00 -4.74104957e+00 -9.20828866e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -93 1.10539569e-06 2.08546876e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --4.00253331e-02 9.90845414e-01 3.30248466e-02 ! x y z --6.29625004e+00 -2.72856242e-01 4.06580221e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -94 1.10539569e-06 6.61548636e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.95782619e-01 1.07222597e-01 -2.12686228e-02 ! x y z --3.48575733e+00 1.03058724e+01 2.52560893e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -95 1.10539569e-06 1.26032554e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.86024380e-01 -5.25659190e-01 -1.93449627e-02 ! x y z -6.91385395e+00 -3.87796179e+00 1.75371157e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -96 1.10539569e-06 2.28076998e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -8.65990862e-01 6.01183979e-01 2.43471322e-01 ! x y z --3.24575804e+00 5.02670056e+00 -8.53001170e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -97 1.10539569e-06 7.05221824e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.39637261e-01 -3.01107928e-01 -4.31558543e-02 ! x y z -9.00360957e+00 -4.71459755e+00 3.83918395e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -98 1.10539569e-06 6.39269462e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.44254383e-01 -1.79341582e-01 -1.04974267e-03 ! x y z -6.75937463e+00 9.18023255e+00 7.52738258e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -99 1.10539569e-06 1.79085174e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --8.13843608e-01 1.99633270e-01 -1.80207407e-01 ! x y z --1.29939496e+00 -6.46953788e+00 -1.46932255e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -100 1.10539569e-06 1.32458326e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --6.32126102e-01 -2.47196136e-02 -2.42353871e-02 ! x y z -2.92190575e-01 -7.70428120e+00 1.58876111e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -101 1.10539569e-06 1.52095833e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --6.91088291e-01 -1.16411821e-01 -1.71283865e-01 ! x y z -1.53252063e+00 -7.10600791e+00 -1.39483264e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -102 1.10539569e-06 1.58389732e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.73673516e-01 -4.82303150e-01 -6.92080251e-02 ! x y z -4.56925372e+00 -5.55908510e+00 8.00568861e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -103 1.10539569e-06 6.78999949e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.70326377e-01 1.68246296e-01 4.59682324e-02 ! x y z --6.01982148e+00 9.02506097e+00 2.30638849e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -104 1.10539569e-06 1.01600407e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --4.27129315e-01 -1.60689947e-01 1.56748833e-01 ! x y z -2.47314908e+00 -8.47972366e+00 -1.94836513e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -105 1.10539569e-06 2.39114531e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -9.42557812e-01 5.29448416e-01 3.49342967e-01 ! x y z --3.26429676e+00 4.46194553e+00 2.04186289e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -106 1.10539569e-06 7.00793819e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.89445458e-01 1.51184189e-01 6.04053839e-02 ! x y z --4.58420182e+00 9.69365003e+00 -2.04005101e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -107 1.10539569e-06 1.90283015e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -8.56335053e-01 -2.54662167e-01 1.41334047e-01 ! x y z -2.10228445e+00 5.79525323e+00 -2.36721595e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -108 1.10539569e-06 2.96231967e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.36313774e+00 -1.39484993e-03 -3.33851852e-01 ! x y z --1.80152731e-01 -5.22233557e+00 9.45811316e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -109 1.10539569e-06 1.11806853e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.46999479e-01 4.69183085e-01 1.48273567e-02 ! x y z --7.62737255e+00 4.03088377e+00 -2.87208744e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -110 1.10539569e-06 3.12227105e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -8.92685552e-01 -1.19778080e+00 7.96108948e-02 ! x y z -4.05065890e+00 2.97803027e+00 -9.38454309e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -111 1.10539569e-06 2.24585100e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --8.16543398e-01 6.44376993e-01 2.28097541e-01 ! x y z --3.89719354e+00 -4.47060943e+00 -1.39962365e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -112 1.10539569e-06 3.76823843e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.19421123e+00 1.32084250e+00 -7.72400022e-02 ! x y z --3.08133686e+00 2.88161171e+00 2.10766368e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -113 1.10539569e-06 1.85829977e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -5.64160987e-01 6.79541760e-01 1.90515744e-02 ! x y z --5.13431676e+00 4.27566051e+00 6.03635853e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -114 1.10539569e-06 3.49744654e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.56537224e+00 5.82910904e-01 -8.86453767e-02 ! x y z --1.64314439e+00 4.51635617e+00 5.81379936e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -115 1.10539569e-06 2.44627752e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.21489220e-01 -9.80146980e-01 3.22082629e-01 ! x y z -4.69945017e+00 -3.07448039e+00 -1.66849603e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -116 1.10539569e-06 1.48946902e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.60792560e-01 -6.87644911e-01 3.66413788e-02 ! x y z -6.37303278e+00 -1.67465930e+00 -3.52540571e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -117 1.10539569e-06 9.48472009e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -4.15102464e-01 -1.75149973e-01 -2.85549434e-02 ! x y z -3.17699950e+00 8.02980191e+00 -3.56193766e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -118 1.10539569e-06 1.46132866e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.45225828e-01 5.98258841e-01 -7.45003268e-02 ! x y z --6.51260143e+00 -3.60690443e+00 1.16008799e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -119 1.10539569e-06 9.12611171e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.84725575e-01 -3.85457639e-01 -5.09774260e-02 ! x y z -8.51718303e+00 -4.26677573e+00 1.26149786e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -120 1.10539569e-06 1.22021916e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.44841297e-01 -5.39901332e-01 1.54755735e-01 ! x y z -7.76971331e+00 1.39084378e+00 -2.38817543e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -121 1.10539569e-06 6.51657309e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.95970426e-02 -3.07995454e-01 -2.51102574e-02 ! x y z -1.12576282e+01 -6.78821641e-01 -4.95031697e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -122 1.10539569e-06 8.35698784e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -3.58295595e-01 -4.07805583e-02 -1.62744429e-01 ! x y z -5.65842623e-02 9.77741717e+00 -2.12236760e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -123 1.10539569e-06 1.08014635e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.57112225e-01 3.26038545e-01 -1.68902691e-01 ! x y z --5.23071639e+00 -6.75806936e+00 -2.04006316e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -124 1.10539569e-06 1.02550770e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --4.30953626e-01 -1.21419747e-01 -1.92524368e-01 ! x y z -2.35863321e+00 -8.68071542e+00 1.79033554e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -125 1.10539569e-06 6.95814419e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.10223043e-02 3.29698672e-01 -2.33216443e-02 ! x y z --1.05036416e+01 -1.18257626e-01 2.99156731e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -126 1.10539569e-06 1.21782856e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.37540070e-01 -4.68051579e-01 3.12017441e-02 ! x y z -6.67938807e+00 -4.87057933e+00 -1.80758676e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -127 1.10539569e-06 1.27648491e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -4.70107622e-01 3.81394049e-01 2.94976958e-02 ! x y z --5.04123374e+00 6.27657958e+00 5.88259014e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -128 1.10539569e-06 1.62877006e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --4.55138385e-01 -6.23311530e-01 5.02376028e-02 ! x y z -5.28132760e+00 -4.06425983e+00 -2.57712837e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -129 1.10539569e-06 2.66986578e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --7.13493035e-02 -1.22710114e+00 -3.17757228e-01 ! x y z -5.28592709e+00 -7.54459447e-01 1.59593511e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -130 1.10539569e-06 6.56076786e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.44014453e-01 -1.91481938e-01 3.64286209e-02 ! x y z -6.89410930e+00 -7.97798070e+00 3.86568598e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -131 1.10539569e-06 1.47791991e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --8.58860487e-02 6.98105661e-01 -6.05263656e-02 ! x y z --7.39418369e+00 -8.88083586e-01 3.68318815e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -132 1.10539569e-06 1.18316891e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.30463297e-01 5.39537186e-01 8.94135251e-02 ! x y z --7.91933257e+00 2.19758719e+00 -1.62232353e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -133 1.10539569e-06 1.53219837e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.10956489e-01 -7.21494396e-01 5.19355079e-02 ! x y z -7.24010486e+00 -1.09301700e+00 -1.60852706e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -134 1.10539569e-06 6.43300411e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.88041757e-02 -3.04722867e-01 -3.14812865e-02 ! x y z -1.07993595e+01 -1.40429215e+00 2.96552456e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -135 1.10539569e-06 8.51935749e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.86591612e-01 -4.29771812e-02 -1.05863039e-01 ! x y z -1.97234762e+00 -9.03918811e+00 -3.56317371e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -136 1.10539569e-06 7.63507559e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.88157405e-01 2.88666278e-01 1.10921677e-01 ! x y z --8.87614786e+00 -5.45415368e+00 -8.46055713e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -137 1.10539569e-06 8.42313899e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.55471151e-01 3.65848047e-01 -8.97493453e-03 ! x y z --9.20044180e+00 -3.87411273e+00 -4.84829815e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -138 1.10539569e-06 6.72600246e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.87375637e-01 -1.36801833e-01 -1.72137474e-02 ! x y z -4.76338219e+00 1.00717170e+01 2.97994928e-03 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -139 1.10539569e-06 1.02461717e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.70166209e-01 -4.01549931e-01 6.36091815e-02 ! x y z -7.45419134e+00 4.96228856e+00 -6.59900898e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -140 1.10539569e-06 6.87840289e-04 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.61307437e-01 -2.72164290e-01 8.71883296e-02 ! x y z -8.60206869e+00 6.06917549e+00 2.98663214e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -141 1.10539569e-06 3.92140957e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.78539217e+00 5.34215902e-01 -1.20626554e-01 ! x y z --1.35880150e+00 -4.33539223e+00 6.63369806e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -142 1.10539569e-06 2.70800787e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --2.31436788e-01 5.48162971e-01 1.64734173e-02 ! x y z --7.51357402e+00 -3.16907615e+00 -1.32982328e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -143 1.10539569e-06 2.04544912e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.94940717e-01 -3.37625386e-01 -4.75480276e-02 ! x y z -6.93658760e+00 6.20070236e+00 -1.01079366e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -144 1.10539569e-06 6.37587428e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -9.71735603e-01 -9.61566886e-01 3.24968572e-01 ! x y z -3.81159593e+00 3.54522020e+00 -1.00291449e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -145 1.10539569e-06 2.08933536e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.75759753e-01 4.21380582e-01 6.28003104e-02 ! x y z --8.52008102e+00 3.47130251e+00 9.83757295e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -146 1.10539569e-06 1.50397186e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --1.39588743e-02 -3.30180335e-01 -1.05230287e-02 ! x y z -1.07866352e+01 -4.97967898e-01 1.76757426e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -147 1.10539569e-06 3.38562140e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --7.29920730e-01 5.89061070e-02 -1.20920053e-01 ! x y z --7.74235710e-01 -7.10977778e+00 1.52865443e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -148 1.10539569e-06 2.32741739e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -4.96020721e-01 7.20447253e-02 -1.09717020e-01 ! x y z --2.02019828e+00 7.69560863e+00 -3.69254498e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -149 1.10539569e-06 2.12332437e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --4.28791847e-01 -1.88187724e-01 1.60873382e-02 ! x y z -3.57426531e+00 -8.26404401e+00 -1.76157927e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -150 1.10539569e-06 1.76485902e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.55182279e-01 1.52620292e-01 -6.22833990e-02 ! x y z --3.88412930e+00 -9.19453465e+00 -7.37459129e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -151 1.10539569e-06 1.57284431e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --3.47785192e-01 9.69744008e-03 -5.97342340e-03 ! x y z --3.78784685e-01 -1.06242178e+01 1.68422812e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -152 1.10539569e-06 5.30159037e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.29951992e-01 1.09612807e+00 -3.27460884e-01 ! x y z --5.22124087e+00 1.65375984e+00 1.97215321e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -153 1.10539569e-06 2.55861101e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.21007349e-02 5.32952897e-01 -1.70855331e-01 ! x y z --8.35827288e+00 -7.84724150e-01 2.96312578e-02 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -154 1.10539569e-06 7.10383113e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius --5.72959788e-01 -1.44639128e+00 -1.28014693e-01 ! x y z -4.55403004e+00 -1.72466270e+00 -1.28754288e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -155 1.10539569e-06 1.36743334e-03 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -2.14648375e-01 2.10911522e-01 -9.00355256e-03 ! x y z --7.51711446e+00 7.37660713e+00 -4.51391831e+00 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -156 1.10539569e-06 2.65261394e-02 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -5.50228741e-02 3.86702728e-01 2.53976723e-02 ! x y z --9.89354846e+00 1.39418207e+00 5.42336207e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot -157 1.10539569e-06 4.92529404e-02 ! ID / G*Mass / Rhill -1.09799198e-05 ! Radius -1.00902512e+00 -1.88090001e-01 3.31281617e-01 ! x y z -1.00573195e+00 5.95735671e+00 3.15059243e-01 ! vx vy vz -4.00000000e-01 4.00000000e-01 4.00000000e-01 ! ip -0.00000000e+00 0.00000000e+00 0.00000000e+00 ! rot +1 +3 0.000120026935827952456416 0.010044948626436044442 +4.25875607065040958e-05 +0.87405601066282501144 -0.50863684909036299686 2.1125685010193420466e-05 +3.058408673972023793 5.407722806790176349 -0.00026983464502404899522 +0.0 0.0 0.33069999999999999396 +4.774949887540911433 0.034051061217069346033 2301.2115366660092022 diff --git a/examples/symba_chambers_2013/sun_MsunAUYR.in b/examples/symba_chambers_2013/sun_MsunAUYR.in index ceb04eea1..45e4aa7ef 100644 --- a/examples/symba_chambers_2013/sun_MsunAUYR.in +++ b/examples/symba_chambers_2013/sun_MsunAUYR.in @@ -1,7 +1,11 @@ -0 ! id -39.476926408897626 ! G*Mass -0.004650467260962157 ! Radius -0.0 !4.7535806948127355e-12 ! J2 -0.0 !-2.2473967953572827e-18 ! J4 -0.0 0.0 0.07 ! Principle axes moments of inertia -11.2093063 -38.75937204 82.25088158 ! Rotation vector (rad/TU) +0 +39.476926408897626 +0.004650467260962157 +4.7535806948127355e-12 +-2.2473967953572827e-18 +0.0 +0.0 +0.07 +11.209306302144773 +-38.759372036774764 +82.25088158389266 diff --git a/python/swiftest/swiftest/constants.py b/python/swiftest/swiftest/constants.py index 7b031b4b7..caf8f8fe7 100644 --- a/python/swiftest/swiftest/constants.py +++ b/python/swiftest/swiftest/constants.py @@ -13,3 +13,4 @@ # Solar oblatenes values: From Mecheri et al. (2004), using Corbard (b) 2002 values (Table II) J2Sun = np.longdouble(2.198e-7) J4Sun = np.longdouble(-4.805e-9) + diff --git a/python/swiftest/swiftest/init_cond.py b/python/swiftest/swiftest/init_cond.py index 76c973533..297b9df68 100644 --- a/python/swiftest/swiftest/init_cond.py +++ b/python/swiftest/swiftest/init_cond.py @@ -1,6 +1,8 @@ import swiftest import numpy as np from astroquery.jplhorizons import Horizons +import astropy.units as u +from astropy.coordinates import SkyCoord import datetime from datetime import date import xarray as xr @@ -74,6 +76,33 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): 'Pluto': np.longdouble(1188.3e3 / swiftest.AU2M) } + planetrot = { + 'Sun' : np.longdouble(2*np.pi / 25.05) / swiftest.JD2S, # Approximate + 'Mercury': np.longdouble(2*np.pi / 58.646) / swiftest.JD2S, + 'Venus': np.longdouble(2*np.pi / 243.0226 ) / swiftest.JD2S, + 'Earth': np.longdouble(2*np.pi / 0.99726968) / swiftest.JD2S, + 'Mars': np.longdouble(2*np.pi / 1.025957) / swiftest.JD2S, + 'Jupiter': np.longdouble(2*np.pi / (9.9250 / 24.0) ) / swiftest.JD2S, + 'Saturn': np.longdouble(2*np.pi / (10.656 / 24.0) ) / swiftest.JD2S, + 'Uranus': np.longdouble(2*np.pi / 0.71833) / swiftest.JD2S, + 'Neptune': np.longdouble(2*np.pi / 0.6713) / swiftest.JD2S, + 'Pluto': np.longdouble(2*np.pi / 6.387230) / swiftest.JD2S + } + + planetIpz = { # Only the polar moments of inertia are used currently. Where the quantity is unkown, we just use the value of a sphere = 0.4 + 'Sun' : np.longdouble(0.070), + 'Mercury' : np.longdouble(0.346), + 'Venus': np.longdouble(0.4), + 'Earth': np.longdouble(0.3307), + 'Mars': np.longdouble(0.3644), + 'Jupiter': np.longdouble(0.2756), + 'Saturn': np.longdouble(0.22), + 'Uranus': np.longdouble(0.23), + 'Neptune': np.longdouble(0.23), + 'Pluto': np.longdouble(0.4) + } + + # Unit conversion factors DCONV = swiftest.AU2M / param['DU2M'] VCONV = (swiftest.AU2M / swiftest.JD2S) / (param['DU2M'] / param['TU2S']) @@ -84,8 +113,14 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): Rcb = np.array([swiftest.RSun / param['DU2M']]) J2RP2 = np.array([swiftest.J2Sun * (swiftest.RSun / param['DU2M']) ** 2]) J4RP4 = np.array([swiftest.J4Sun * (swiftest.RSun / param['DU2M']) ** 4]) + + solarpole = SkyCoord(ra=286.13 * u.degree, dec=63.87 * u.degree) + solarrot = planetrot['Sun'] * param['TU2S'] + rotcb = solarpole.cartesian * solarrot + Ipsun = np.array([0.0, 0.0, planetIpz['Sun']]) + cbid = np.array([0]) - cvec = np.vstack([GMcb, Rcb, J2RP2, J4RP4]) + cvec = np.vstack([GMcb, Rcb, J2RP2, J4RP4, Ipsun[0], Ipsun[1], Ipsun[2], rotcb.x, rotcb.y, rotcb.z ]) # Horizons date time internal variables tstart = datetime.date.fromisoformat(ephemerides_start_date) @@ -155,6 +190,12 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): Rhill = [] Rpl = [] GMpl = [] + Ip_x = [] + Ip_y = [] + Ip_z = [] + rot_x = [] + rot_y = [] + rot_z = [] pldata = {} if ispl: @@ -191,19 +232,35 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): p10.append(pldata[key].vectors()['vx'][0] * VCONV) p11.append(pldata[key].vectors()['vy'][0] * VCONV) p12.append(pldata[key].vectors()['vz'][0] * VCONV) + pvec = np.vstack([p1, p2, p3, p4, p5, p6]) if ispl: Rpl.append(planetradius[key] * DCONV) GMpl.append(GMcb[0] / MSun_over_Mpl[key]) + pvec = np.vstack([pvec, GMpl, Rpl]) + # Generate planet value vectors if (param['RHILL_PRESENT'] == 'YES'): Rhill.append(pldata[key].elements()['a'][0] * DCONV * (3 * MSun_over_Mpl[key]) ** (-THIRDLONG)) - pvec = np.vstack([p1, p2, p3, p4, p5, p6, GMpl, Rpl, Rhill, p7, p8, p9, p10, p11, p12]) - else: - pvec = np.vstack([p1, p2, p3, p4, p5, p6, GMpl, Rpl, p7, p8, p9, p10, p11, p12]) + pvec = np.vstack([pvec, Rhill]) + if (param['ROTATION'] == 'YES'): + RA = pldata[key].ephemerides()['NPole_RA'][0] + DEC = pldata[key].ephemerides()['NPole_DEC'][0] + + rotpole = SkyCoord(ra=RA * u.degree, dec=DEC * u.degree) + rotrate = planetrot[key] * param['TU2S'] + rot = rotpole.cartesian * rotrate + Ip = np.array([0.0, 0.0, planetIpz[key]]) + Ip_x.append(Ip[0]) + Ip_y.append(Ip[1]) + Ip_z.append(Ip[2]) + rot_x.append(rot.x) + rot_y.append(rot.y) + rot_z.append(rot.z) + pvec = np.vstack([pvec, Ip_x, Ip_y, Ip_z, rot_x, rot_y, rot_z]) else: - pvec = np.vstack([p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12]) plab = tlab.copy() - + pvec = np.vstack([pvec, p7, p8, p9, p10, p11, p12]) + if idval is None: plid = np.array([planetid[key]], dtype=int) else: diff --git a/python/swiftest/swiftest/io.py b/python/swiftest/swiftest/io.py index ecc5051b0..1116f6509 100644 --- a/python/swiftest/swiftest/io.py +++ b/python/swiftest/swiftest/io.py @@ -741,6 +741,12 @@ def swiftest_xr2infile(ds, param, framenum=-1): RSun = np.double(cb['Radius']) J2 = np.double(cb['J_2']) J4 = np.double(cb['J_4']) + Ip_xcb = np.double(cb['Ip_x']) + Ip_ycb = np.double(cb['Ip_y']) + Ip_zcb = np.double(cb['Ip_z']) + rot_xcb = np.double(cb['rot_x']) + rot_ycb = np.double(cb['rot_y']) + rot_zcb = np.double(cb['rot_z']) cbid = int(0) if param['IN_TYPE'] == 'ASCII': @@ -751,6 +757,13 @@ def swiftest_xr2infile(ds, param, framenum=-1): print(RSun, file=cbfile) print(J2, file=cbfile) print(J4, file=cbfile) + if param['ROTATION'] == 'YES': + print(Ip_xcb, file=cbfile) + print(Ip_ycb, file=cbfile) + print(Ip_zcb, file=cbfile) + print(rot_xcb, file=cbfile) + print(rot_ycb, file=cbfile) + print(rot_zcb, file=cbfile) cbfile.close() plfile = open(param['PL_IN'], 'w') @@ -770,6 +783,9 @@ def swiftest_xr2infile(ds, param, framenum=-1): print(pli['capom'].values, pli['omega'].values, pli['capm'].values, file=plfile) else: print(f"{param['IN_FORM']} is not a valid input format type.") + if param['ROTATION'] == 'YES': + print(pli['Ip_x'].values, pli['Ip_y'].values, pli['Ip_z'].values, file=plfile) + print(pli['rot_x'].values, pli['rot_y'].values, pli['rot_z'].values, file=plfile) plfile.close() # TP file @@ -795,6 +811,14 @@ def swiftest_xr2infile(ds, param, framenum=-1): cbfile.write_record(np.double(RSun)) cbfile.write_record(np.double(J2)) cbfile.write_record(np.double(J4)) + if param['ROTATION'] == 'YES': + cbfile.write_record(np.double(Ip_xcb)) + cbfile.write_record(np.double(Ip_ycb)) + cbfile.write_record(np.double(Ip_zcb)) + cbfile.write_record(np.double(rot_xcb)) + cbfile.write_record(np.double(rot_ycb)) + cbfile.write_record(np.double(rot_zcb)) + cbfile.close() plfile = FortranFile(param['PL_IN'], 'w') @@ -829,9 +853,15 @@ def swiftest_xr2infile(ds, param, framenum=-1): plfile.write_record(v6) plfile.write_record(Gmass) if param['RHILL_PRESENT'] == 'YES': - rhill = pl['Rhill'].values - plfile.write_record(rhill) + plfile.write_record(pl['Rhill'].values) plfile.write_record(radius) + if param['ROTATION'] == 'YES': + plfile.write_record(pl['Ip_x'].values) + plfile.write_record(pl['Ip_y'].values) + plfile.write_record(pl['Ip_z'].values) + plfile.write_record(pl['rot_x'].values) + plfile.write_record(pl['rot_y'].values) + plfile.write_record(pl['rot_z'].values) plfile.close() tpfile = FortranFile(param['TP_IN'], 'w') ntp = tp.id.count().values diff --git a/python/swiftest/swiftest/simulation_class.py b/python/swiftest/swiftest/simulation_class.py index 78f262df8..e7b8c560b 100644 --- a/python/swiftest/swiftest/simulation_class.py +++ b/python/swiftest/swiftest/simulation_class.py @@ -27,12 +27,12 @@ def __init__(self, codename="Swiftest", param_file=""): 'OUT_TYPE': 'REAL8', 'OUT_FORM': "EL", 'OUT_STAT': "REPLACE", - 'CHK_RMAX': "1000.0", - 'CHK_EJECT': "1000.0", - 'CHK_RMIN': f"{constants.RSun / constants.AU2M}", - 'CHK_QMIN': f"{constants.RSun / constants.AU2M}", + 'CHK_RMAX': "-1.0", + 'CHK_EJECT': "-1.0", + 'CHK_RMIN': "-1.0", + 'CHK_QMIN': "-1.0", 'CHK_QMIN_COORD': "HELIO", - 'CHK_QMIN_RANGE': f"{constants.RSun / constants.AU2M} 1000.0", + 'CHK_QMIN_RANGE': "-1.0 -1.0", 'ENC_OUT': "enc.dat", 'MU2KG': constants.MSun, 'TU2S': constants.JD2S,