From 67ba5a6ac6a052fd287cd5ff6850c6ba034057c0 Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 15:47:55 -0400 Subject: [PATCH 1/8] 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, From 7203c8149b3935f8082c16e29d17972dd3d8c81f Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 16:40:36 -0400 Subject: [PATCH 2/8] Added initial conditions generator for Chambers run --- examples/symba_chambers_2013/init_cond.py | 38 +- examples/symba_chambers_2013/param.in | 3 +- .../symba_chambers_2013/pl_chambers_2013.in | 944 +++++++++++++++++- python/swiftest/swiftest/init_cond.py | 35 +- python/swiftest/swiftest/simulation_class.py | 12 +- 5 files changed, 1007 insertions(+), 25 deletions(-) diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py index 5cda5056e..5f6782491 100755 --- a/examples/symba_chambers_2013/init_cond.py +++ b/examples/symba_chambers_2013/init_cond.py @@ -1,7 +1,7 @@ #!/usr/bin/env python3 import swiftest import numpy as np -import os +from numpy.random import default_rng # Initialize simulation object sim = swiftest.Simulation() @@ -10,9 +10,11 @@ MU2KG = swiftest.MSun TU2S = swiftest.YR2S DU2M = swiftest.AU2M +GU = swiftest.GC / (DU2M**3 / (MU2KG * TU2S**2)) sim.param['MU2KG'] = MU2KG sim.param['TU2S'] = TU2S sim.param['DU2M'] = DU2M +sim.param['GU'] = GU # Simulation time parameters sim.param['T0'] = 0.0 @@ -35,10 +37,14 @@ sim.param['ROTATION'] = "YES" sim.param['CHK_RMAX'] = 1000.0 sim.param['CHK_EJECT'] = 1000.0 +sim.param['IN_FORM'] = 'EL' +sim.param['OUT_FORM'] = 'EL' # Add central body sim.add("Sun") -sim.add("Earth") +GMcb = sim.ds['GMass'].values[0] +sim.add("Jupiter") +sim.add("Saturn") # Add bodies described in Chambers (2013) Sec. 2.1, with the uniform spatial distribution and two bodies sizes (big and small) Nb = 14 @@ -49,7 +55,35 @@ Rb = (3 * Mb / (4 * np.pi * dens) )**(1.0 / 3.0) Rs = (3 * Ms / (4 * np.pi * dens) )**(1.0 / 3.0) +# Define the initial orbital elements of the big and small bodies +avalb = default_rng().uniform(0.3, 2.0, Nb) +avals = default_rng().uniform(0.3, 2.0, Ns) +evalb = default_rng().uniform(0.0, 0.001, Nb) +evals = default_rng().uniform(0.0, 0.001, Ns) +incvalb = default_rng().uniform(0.0, 0.0005 * 180 / np.pi, Nb) +incvals = default_rng().uniform(0.0, 0.0005 * 180 / np.pi, Ns) +capomvalb = default_rng().uniform(0.0, 360.0, Nb) +capomvals = default_rng().uniform(0.0, 360.0, Ns) +omegavalb = default_rng().uniform(0.0, 360.0, Nb) +omegavals = default_rng().uniform(0.0, 360.0, Ns) +capmvalb = default_rng().uniform(0.0, 360.0, Nb) +capmvals = default_rng().uniform(0.0, 360.0, Ns) +GMvalb = np.full(Nb, Mb * GU) +GMvals = np.full(Ns, Ms * GU) +Rvalb = np.full(Nb, Rb) +Rvals = np.full(Ns, Rs) +Rhb = avalb * (GMvalb / (3 * GMcb))**(1.0/3.0) +Rhs = avals * (GMvals / (3 * GMcb))**(1.0/3.0) +# Give the bodies unique ids +idb = np.arange(100, 100 + Nb) +ids = np.arange(100 + Nb, 100 + Nb + Ns) + +# Populate the simulation object with the two types of bodies +sim.addp(idb, avalb, evalb, incvalb, capomvalb, omegavalb, capmvalb, GMpl=GMvalb, Rpl=Rvalb, Rhill=Rhb) +sim.addp(ids, avals, evals, incvals, capomvals, omegavals, capmvals, GMpl=GMvals, Rpl=Rvals, Rhill=Rhs) + +# Save everything to a set of initial conditions files sim.save("param.in") diff --git a/examples/symba_chambers_2013/param.in b/examples/symba_chambers_2013/param.in index d096d13fe..19e98e87a 100644 --- a/examples/symba_chambers_2013/param.in +++ b/examples/symba_chambers_2013/param.in @@ -21,7 +21,7 @@ CHK_QMIN_RANGE -1.0 -1.0 MU2KG 1.988409870698051e+30 TU2S 31557600.0 DU2M 149597870700.0 -IN_FORM XV +IN_FORM EL ENC_OUT enc.dat EXTRA_FORCE NO DISCARD_OUT discard.out @@ -34,4 +34,5 @@ ROTATION YES TIDES NO ENERGY YES GR NO +GU 39.476926408897626 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 56d3f87e4..67d6989d2 100644 --- a/examples/symba_chambers_2013/pl_chambers_2013.in +++ b/examples/symba_chambers_2013/pl_chambers_2013.in @@ -1,7 +1,937 @@ -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 +156 +5 0.03769225108898567778 0.35527124629249239585 +0.00046732617030490929307 +5.2035240151214816606 0.048518677336151047608 1.303568872392780964 +100.51672855018520636 273.38350380059858935 317.24733983021741324 +0.0 0.0 0.27560000000000001164 +-80.96882015576837074 -2387.9972166758053216 5008.7356633679251408 +6 0.01128589982009127331 0.43766275227962878915 +0.00038925687730393611812 +9.581843040035108672 0.052193289679895830957 2.486259885154824989 +113.5952699440524043 335.70375572759667193 225.29914557362599226 +0.0 0.0 0.22000000000000000111 +441.93102714484842597 378.52937229654817222 5135.9113502360652523 +100 1.105353939449133635e-05 0.0076985326356758955327 +2.365334758890906669e-05 +1.6971843938499941995 0.0005816364466293669124 0.026187138610341612749 +199.66295678003845637 349.1838313868397563 239.5201319664316486 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +101 1.105353939449133635e-05 0.0033890547724920599881 +2.365334758890906669e-05 +0.7471359987644425882 0.0001484221862068907502 0.015912119531244738047 +7.72758829713570794 98.74578965353445881 292.48405218591960875 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +102 1.105353939449133635e-05 0.0032403861730961017119 +2.365334758890906669e-05 +0.71436117806925114593 0.0006667054143762592925 0.02623486192401118175 +42.268600860678539277 152.49698993249327827 151.45780278371265126 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +103 1.105353939449133635e-05 0.00848716742472066395 +2.365334758890906669e-05 +1.8710433251236489038 0.0009470568874946718764 0.022705939228401819663 +38.80730095910011812 143.34784343975428555 218.97185967315348876 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +104 1.105353939449133635e-05 0.004200287321695052973 +2.365334758890906669e-05 +0.9259767321153884456 0.00016780856600745607692 0.008111571611109064106 +275.52852024631755512 255.91035494922905968 283.72241775520802776 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +105 1.105353939449133635e-05 0.008200656366547455989 +2.365334758890906669e-05 +1.8078803667251066578 0.0005806820778782817395 0.0043050823832896005852 +262.76594876192444872 197.55827100651990236 359.12597183160295344 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +106 1.105353939449133635e-05 0.0031948304428970366643 +2.365334758890906669e-05 +0.70431816363998178687 0.0006166870106254423523 0.028312766571177344632 +100.38580777052690962 247.08615435134802851 28.274821850602428697 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +107 1.105353939449133635e-05 0.00645602527104512983 +2.365334758890906669e-05 +1.4232667256020528068 0.0003621838900236009598 0.012197930522699866954 +189.06372569967658137 282.1051532665543391 175.25556542499867874 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +108 1.105353939449133635e-05 0.0041177716184015826697 +2.365334758890906669e-05 +0.90778568578166862757 0.00088526942168412907277 0.0052077153285733724167 +355.61137098230727815 88.33816057346382422 35.585620637422806567 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +109 1.105353939449133635e-05 0.00862613046390682299 +2.365334758890906669e-05 +1.9016785010188288485 0.00044486325805035754379 0.0175194201364592822 +74.14924815238387623 142.82151728160823723 178.34909057749027284 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +110 1.105353939449133635e-05 0.004510164863223143653 +2.365334758890906669e-05 +0.9942909619010631417 0.00081444817122918432024 0.00707240791737578315 +68.31335729815678803 294.6246189285234891 122.92252049539324332 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +111 1.105353939449133635e-05 0.004651281910070013023 +2.365334758890906669e-05 +1.0254010007810474558 0.0004312636045460260673 0.01674019444109412641 +337.31604521335424351 171.74369728670052382 266.0027373696372024 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +112 1.105353939449133635e-05 0.0038244126458011504717 +2.365334758890906669e-05 +0.84311306651055351047 2.7875407573919153787e-05 0.017986416975313480882 +225.7229782864748131 308.33509788803951324 112.41779338822905743 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +113 1.105353939449133635e-05 0.0029643452153518679503 +2.365334758890906669e-05 +0.6535064116198798789 0.00037425733872322575592 0.0041403864217484038016 +325.7386399373389736 300.94512487157606984 333.7954171286211249 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +114 1.1053539394491335305e-06 0.0038455497484053390705 +1.0978911404624207699e-05 +1.826471255610046418 0.00087780859820767002683 0.009722642750510790757 +120.71268570369510087 352.6610169497902234 335.05498321080250435 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +115 1.1053539394491335305e-06 0.0017443709230869269861 +1.0978911404624207699e-05 +0.82850139995235116963 0.00031729137906426232822 0.01699303328109756997 +182.6703622068744437 213.3557405525360764 250.28249688543755269 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +116 1.1053539394491335305e-06 0.0020702769104938987164 +1.0978911404624207699e-05 +0.98329277101103618186 0.00058251749032120785 0.011797523288670572553 +312.4959308001525642 359.33855098907400816 162.46970124565356741 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +117 1.1053539394491335305e-06 0.00080530024674720453575 +1.0978911404624207699e-05 +0.38248309059826302025 0.0005969943570699224548 0.01380590580233940004 +48.708972362298950998 330.19391063831329802 102.38386339213280962 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +118 1.1053539394491335305e-06 0.002285280534438147877 +1.0978911404624207699e-05 +1.085410274275427911 0.00025130462779814554623 0.006036721527856849477 +45.867981468353910657 64.60926497341390018 253.24516188255682891 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +119 1.1053539394491335305e-06 0.002898627766350067431 +1.0978911404624207699e-05 +1.3767239126595514076 0.0002530591745570577502 0.006701891775872563864 +256.03674703821604908 136.96778530539529584 171.0933268407595449 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +120 1.1053539394491335305e-06 0.0014542058638660382582 +1.0978911404624207699e-05 +0.6906854374182271883 0.0006003363904860686438 0.004833197414767716714 +56.14180364632849063 273.3261312610867435 199.6250014185157795 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +121 1.1053539394491335305e-06 0.0024981790226714743672 +1.0978911404624207699e-05 +1.1865279283331482407 0.00082336466808019992576 0.018265656692325986377 +132.47659416127314103 48.61694507533226073 207.43492631988763719 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +122 1.1053539394491335305e-06 0.0028828564897507973385 +1.0978911404624207699e-05 +1.3692332324558205148 0.00040127628444019838947 0.007585985018489528435 +343.91044751539715207 196.83240087837353371 62.629566549905092643 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +123 1.1053539394491335305e-06 0.0038943677929657182489 +1.0978911404624207699e-05 +1.8496577337414494835 0.00043859595304759446763 0.011023337424484119101 +72.91657450189477174 344.12080798821472172 100.72699073828245275 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +124 1.1053539394491335305e-06 0.001675555086167977485 +1.0978911404624207699e-05 +0.7958168278400463702 0.00034926098253153968964 0.0099212607809325507185 +135.68666178212228601 296.83136368103197356 107.90664998734122548 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +125 1.1053539394491335305e-06 0.0023748127187634624592 +1.0978911404624207699e-05 +1.1279342232088616083 0.0007746490841826023793 0.017642382003092794973 +268.43319860669271293 252.98510234252850637 120.5395807841271818 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +126 1.1053539394491335305e-06 0.0012043684437486966761 +1.0978911404624207699e-05 +0.5720233744428824263 0.00034480510921434492748 0.017251867284702383787 +154.90093862136961889 356.5658098586184792 57.388355737939328094 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +127 1.1053539394491335305e-06 0.0007890385972653567634 +1.0978911404624207699e-05 +0.37475950429965448896 0.00031644523056148786714 0.027616371909429979908 +171.3578942788073789 299.01930445458822305 225.76712007928702519 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +128 1.1053539394491335305e-06 0.001777100120316001991 +1.0978911404624207699e-05 +0.84404636539789934524 0.00087953695475614056525 0.0006577653108544906282 +70.992958768969344874 130.55955257913041123 40.84173864503155471 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +129 1.1053539394491335305e-06 0.0040678503543026272943 +1.0978911404624207699e-05 +1.9320545124498424983 0.0005834355620338948882 0.0002608548084170523012 +270.88252512263818517 195.45419716249062958 140.3411838625577559 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +130 1.1053539394491335305e-06 0.0026829522072340253071 +1.0978911404624207699e-05 +1.2742872689972433164 0.00060055603766298343727 0.0064920160946518451484 +162.55438763596669105 178.06818403372230364 68.91417398625686985 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +131 1.1053539394491335305e-06 0.00071194132454311543046 +1.0978911404624207699e-05 +0.33814160524074998726 0.00091270005510230359236 0.0011893404294311912296 +250.26053020749216671 223.06302666485407826 192.24977287221634015 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +132 1.1053539394491335305e-06 0.002812359973542913556 +1.0978911404624207699e-05 +1.3357504097390577869 0.00026337894924635706691 0.027464957742126525075 +199.22420478899164209 171.93696705940942593 94.48037873747355775 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +133 1.1053539394491335305e-06 0.0028963111872931898259 +1.0978911404624207699e-05 +1.3756236369290164578 0.0006573629715563674607 0.021333593911504113161 +87.075565489431383526 153.68132922246732619 173.6496245283211124 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +134 1.1053539394491335305e-06 0.002785993277843155817 +1.0978911404624207699e-05 +1.3232273597327499637 0.00073604309841295738944 0.011001193704390446401 +339.40562860328299166 83.8138037375181284 177.19280309751334812 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +135 1.1053539394491335305e-06 0.0019038535805104330423 +1.0978911404624207699e-05 +0.9042488245366071542 0.0009361201833287777522 0.011953628791727587133 +153.92620389616467946 44.723571716565317047 318.60665806354387541 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +136 1.1053539394491335305e-06 0.0010869785583050439427 +1.0978911404624207699e-05 +0.5162682118549849486 0.00022103242271550004863 0.017892470794599633216 +219.23592171773313453 303.7464473661763691 122.639027547649121175 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +137 1.1053539394491335305e-06 0.0032939926784115201725 +1.0978911404624207699e-05 +1.5645052949330442438 0.000991170480938109201 0.019190182292000351055 +346.07086810942291777 54.79631567078042309 335.26549358631774567 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +138 1.1053539394491335305e-06 0.0010896488129816246724 +1.0978911404624207699e-05 +0.5175364683413183897 0.0002088585855651833357 0.025066384881598170048 +68.72020084550582908 232.28654212667152024 35.431071301701685172 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +139 1.1053539394491335305e-06 0.0008099901996317258652 +1.0978911404624207699e-05 +0.38471061714041626667 0.0005249334614017682433 0.023288096472191086028 +254.41357781115257808 117.993033232943290045 239.40340264257432068 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +140 1.1053539394491335305e-06 0.0026223707179148065143 +1.0978911404624207699e-05 +1.2455136589537159075 0.00084338241955415730895 0.0023232936840040861563 +269.3199834748557464 357.25116637136363806 20.298945789952838936 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +141 1.1053539394491335305e-06 0.0026159634140416132347 +1.0978911404624207699e-05 +1.2424704643219992573 0.00044638238521322806543 0.0012966745044005083872 +154.64047440242666198 236.78631050296556282 350.38453034745816694 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +142 1.1053539394491335305e-06 0.003456693657444609839 +1.0978911404624207699e-05 +1.6417812842988768196 0.0006303353332627267564 0.016812028430699579862 +144.4420533567427185 153.14985299997016455 271.69228258365166084 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +143 1.1053539394491335305e-06 0.0007301401336089313031 +1.0978911404624207699e-05 +0.34678525928756914753 0.00045501001474390856362 0.02594919197139722647 +116.332998235860372915 256.36327622207414834 161.64230688982479478 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +144 1.1053539394491335305e-06 0.0034349976066789206724 +1.0978911404624207699e-05 +1.6314765903860701801 0.000856330141764382333 0.014786839646800843326 +42.332000671474226294 138.1316800890192269 216.26130379036857221 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +145 1.1053539394491335305e-06 0.0025439979715687957412 +1.0978911404624207699e-05 +1.2082899646084366285 0.0004070939489201934802 0.004503954508818640984 +48.07928625757347163 276.96063567133671768 187.66327227579773762 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +146 1.1053539394491335305e-06 0.003859828194739848691 +1.0978911404624207699e-05 +1.8332529054419246517 0.0003598790557617636025 0.014970896975960788777 +271.0689558389924514 52.94488824830970941 1.9029339518958510169 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +147 1.1053539394491335305e-06 0.0040866242766629467682 +1.0978911404624207699e-05 +1.9409713206540264441 0.00011426903199890425854 0.011184164105792333138 +67.28912042774895497 289.72221855724853867 184.15766072550803756 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +148 1.1053539394491335305e-06 0.00316122359322290065 +1.0978911404624207699e-05 +1.5014456718372271027 0.00086962582925130050334 0.01914542568965043115 +250.74751474204416013 182.77230460121793953 167.04339861860248106 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +149 1.1053539394491335305e-06 0.00097271678791404052593 +1.0978911404624207699e-05 +0.46199876980165455187 0.00034328565960485866526 0.021012970927227814233 +297.8391700642537785 90.09431163085170624 45.254280687068330735 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +150 1.1053539394491335305e-06 0.0033482251435429574658 +1.0978911404624207699e-05 +1.5902633906967305055 0.00019727479448159425065 0.025485920919319147299 +39.29235960954184037 250.38592799994987104 142.50585121059737048 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +151 1.1053539394491335305e-06 0.0025466847865406780281 +1.0978911404624207699e-05 +1.2095660865250288385 0.00016751618714245874271 0.0016713720493250405416 +333.37472688335139992 145.95098571742994409 234.96494950660738255 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +152 1.1053539394491335305e-06 0.000719807963238783997 +1.0978911404624207699e-05 +0.3418779213453242427 0.00074674271745437542386 0.0011890436742735912591 +195.44114509251284062 283.06165279842355176 315.09159268499945483 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +153 1.1053539394491335305e-06 0.00157025465995761385 +1.0978911404624207699e-05 +0.74580364006202581884 0.0002094126871885914243 0.022074802063619617526 +153.41729982837850343 124.788130215813424684 122.51830102672674627 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +154 1.1053539394491335305e-06 0.001874785507921179697 +1.0978911404624207699e-05 +0.89044273632717141354 0.00085481426099495228573 0.021979234537361322377 +202.10509721048100573 320.04294945394877914 309.06206051991284767 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +155 1.1053539394491335305e-06 0.0019122234532341110998 +1.0978911404624207699e-05 +0.9082241552287277919 0.0004597575061389033409 0.020148347940879648976 +163.72826672182648622 305.80294187113088356 283.64822483177448476 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +156 1.1053539394491335305e-06 0.0039778936569127347484 +1.0978911404624207699e-05 +1.8893289380113937881 0.00049610582898170353487 0.0014385636633325129576 +157.23049053490680649 323.571089698826313 102.65627395667121391 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +157 1.1053539394491335305e-06 0.0019708158646589233432 +1.0978911404624207699e-05 +0.93605303855254251566 0.00019680356564931200973 0.020455589191385341319 +278.47203022847156717 255.0643652208907497 258.27522081614978333 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +158 1.1053539394491335305e-06 0.003967327173706975909 +1.0978911404624207699e-05 +1.8843103115182096019 0.0009097170816114348319 0.0018044172991811243192 +117.34832961121314554 234.13353289800707557 259.1454183769003521 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +159 1.1053539394491335305e-06 0.0013842445575880614169 +1.0978911404624207699e-05 +0.6574568164714708063 1.4540080482081997083e-05 0.009892527806868892798 +133.30970394955991765 89.751557170898067284 193.05936528552655318 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +160 1.1053539394491335305e-06 0.0035162555550013580372 +1.0978911404624207699e-05 +1.670070631969415409 0.00091422718885765909876 0.028436761433306107644 +353.37301866523387162 51.552218623340351655 332.64071506725503014 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +161 1.1053539394491335305e-06 0.0019396124736002073748 +1.0978911404624207699e-05 +0.9212327656201906301 0.00057552402003210683355 0.00950969025161889879 +283.93076060868486366 267.5628271717484381 143.52189652919906848 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +162 1.1053539394491335305e-06 0.0031661554024814644338 +1.0978911404624207699e-05 +1.503788069787651116 0.00018115407388301875608 0.02283324079244578031 +177.01583573865022458 234.49703162631655573 336.92462545405476249 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +163 1.1053539394491335305e-06 0.0033722481017907932155 +1.0978911404624207699e-05 +1.6016732658992518523 0.00053832824323583156837 0.017356963073014760462 +210.48701348074004613 252.36586170021701037 241.77353994833666206 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +164 1.1053539394491335305e-06 0.0018549151304582611073 +1.0978911404624207699e-05 +0.8810051589588920873 0.00039455774939218580828 0.014186890964832121059 +268.78131437644134394 234.86755048529201417 135.82402939345988102 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +165 1.1053539394491335305e-06 0.0011704042480178788963 +1.0978911404624207699e-05 +0.5558918376585835208 5.564504104615730423e-05 0.023585563294025682635 +48.050891171653226763 199.47056209951909977 219.93699201694707313 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +166 1.1053539394491335305e-06 0.0040699423683521412583 +1.0978911404624207699e-05 +1.9330481294299255435 0.00052600080391025812437 0.021919559654086234746 +39.096865156286462195 128.8303250981957433 331.23155076245132022 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +167 1.1053539394491335305e-06 0.0028820213749619458248 +1.0978911404624207699e-05 +1.3688365887359974504 0.00026969964684439974665 0.017275229592356205338 +330.69633084561246505 147.5243863767787218 153.22005163183507648 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +168 1.1053539394491335305e-06 0.003980615057015220282 +1.0978911404624207699e-05 +1.8906214863822126127 0.00036100510956251820067 0.01521335462518577826 +340.23428559641661195 130.00201841928748081 169.2953836810935968 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +169 1.1053539394491335305e-06 0.0012767502410615831949 +1.0978911404624207699e-05 +0.60640162485459569996 0.00020696978447336191083 0.0056857512796010872477 +181.57921819815450704 115.817629846161409546 320.0491270864600324 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +170 1.1053539394491335305e-06 0.002399637727624715802 +1.0978911404624207699e-05 +1.1397250380654750934 0.0006442965482497394836 0.002941317764456898541 +333.45417971179358574 185.18080690840253055 269.38711004629510626 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +171 1.1053539394491335305e-06 0.0035439123990565288466 +1.0978911404624207699e-05 +1.6832064471304613029 0.00072650143891821498086 0.02076627535231315505 +82.57443269249242235 19.640854808506546192 343.42020390750877823 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +172 1.1053539394491335305e-06 0.0021893049972247140684 +1.0978911404624207699e-05 +1.0398260089737605405 0.00095820846394839131327 0.0047667113246031159526 +231.34174007276843099 346.25149299156703364 359.3927866655562866 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +173 1.1053539394491335305e-06 0.003808558580359920666 +1.0978911404624207699e-05 +1.8089020367553392266 0.00080486608391464400675 0.026061626097182894607 +333.73026949841113264 290.72940981027193175 353.58605983201931622 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +174 1.1053539394491335305e-06 0.0036046989829116192528 +1.0978911404624207699e-05 +1.7120774682852628601 0.000408327401359263474 0.023141693393884565227 +190.15243963894701551 10.9019960190977514 206.20261856852459914 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +175 1.1053539394491335305e-06 0.0025289888373368178023 +1.0978911404624207699e-05 +1.201161269352922556 0.0008119307695628591708 0.02330748815658333778 +293.74129631246364625 228.8273707230899845 50.496213420126203175 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +176 1.1053539394491335305e-06 0.0035813927117425375896 +1.0978911404624207699e-05 +1.7010079887177620428 0.0008261188698801552469 0.0033817239341808955078 +270.3371253433707011 303.21090622901044753 55.41742303760840116 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +177 1.1053539394491335305e-06 0.0029864597633055024009 +1.0978911404624207699e-05 +1.4184403454864715766 0.0006100544512734605053 0.0031745002208599282408 +128.80488878502711714 46.8536787103431962 190.83918585108236243 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +178 1.1053539394491335305e-06 0.0023567486788243735888 +1.0978911404624207699e-05 +1.119354578719118809 0.00053509985437350090966 0.020556636907334471104 +93.86909946810600047 329.53011231398716063 68.53400707162742833 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +179 1.1053539394491335305e-06 0.0034140907141169699521 +1.0978911404624207699e-05 +1.6215467127855098806 0.0007155502442420363365 0.021636466572293039379 +15.484181145735732343 319.8815184503466753 126.01422430110272899 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +180 1.1053539394491335305e-06 0.0025675289816375823053 +1.0978911404624207699e-05 +1.2194661855177959531 0.00048216964694683286767 0.017262139448479693615 +335.34868702893140835 186.94351621903143723 26.64334186676326155 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +181 1.1053539394491335305e-06 0.0034248039486638383743 +1.0978911404624207699e-05 +1.6266350398738742555 0.00081883719141492166717 0.003694641256859914659 +218.92219023702784853 218.29951433993917931 259.11630406545350525 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +182 1.1053539394491335305e-06 0.0028977454974775277425 +1.0978911404624207699e-05 +1.3763048727716702402 0.00040844368595388948033 0.020281148678402180496 +289.91356665087704414 165.61033347014003425 288.23615258185753873 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +183 1.1053539394491335305e-06 0.0011783105169288593343 +1.0978911404624207699e-05 +0.5596469764163169991 0.00026931995881909574905 0.0010082118613559110359 +154.06847791408287662 39.901773134600148296 187.47166921424076236 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +184 1.1053539394491335305e-06 0.002491583136722370224 +1.0978911404624207699e-05 +1.1833951653006802296 0.00048583037529452476142 0.027591397048547783993 +149.2419352144119955 63.397457239281493457 295.88510332284619153 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +185 1.1053539394491335305e-06 0.0039004676796874942429 +1.0978911404624207699e-05 +1.8525549184065084329 0.0005185930845717148926 0.019438436095256288877 +292.5205667956682305 105.44538380282074286 123.70099114339642199 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +186 1.1053539394491335305e-06 0.0034369089438783143244 +1.0978911404624207699e-05 +1.6323843936087221174 0.00010509062853767581647 0.0032505177502269231964 +143.85173028989797217 189.31567203063778493 223.19791686246551876 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +187 1.1053539394491335305e-06 0.004128326569038383845 +1.0978911404624207699e-05 +1.9607781215797892571 0.0007065471846638792094 0.021817391680966469192 +267.38223511629621498 257.29519900708868363 77.6961655587950446 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +188 1.1053539394491335305e-06 0.0020499009381892692037 +1.0978911404624207699e-05 +0.97361505776992030015 0.00093181583258145170724 0.0006097968157406988537 +145.87863557509786006 138.78820490355309403 73.09644258927524163 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +189 1.1053539394491335305e-06 0.0013791485719345369476 +1.0978911404624207699e-05 +0.65503644177237263513 0.00014046622778450534863 0.015442926715237644333 +113.70608882665497674 243.3631092677804304 81.13225672390240106 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +190 1.1053539394491335305e-06 0.0024853665935430262675 +1.0978911404624207699e-05 +1.1804425738198291995 0.00096027971308601327814 0.023639897027866934448 +334.99323146372387328 76.84265628244520485 155.49317295086675017 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +191 1.1053539394491335305e-06 0.003654973676990583536 +1.0978911404624207699e-05 +1.7359557924853059596 0.00054611453698988137676 0.026616017411125021674 +173.99931107991119461 359.07752431447528352 78.019305599684969366 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +192 1.1053539394491335305e-06 0.003360253875484283402 +1.0978911404624207699e-05 +1.5959765226465592036 0.00089913403289640294203 0.014826170842651739679 +204.39386065798206005 17.787041065446423715 47.062099086162277217 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +193 1.1053539394491335305e-06 0.0036411616313956547123 +1.0978911404624207699e-05 +1.7293956630081683823 0.0004269520101434239978 0.011446361145646259325 +225.78458330355539374 323.21819225965253963 348.78801121716753642 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +194 1.1053539394491335305e-06 0.0026716566036685272943 +1.0978911404624207699e-05 +1.2689223415936377037 0.00043093138180629422026 0.007198866767697694423 +49.75485744273438371 283.14995135249375835 251.76161931871857291 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +195 1.1053539394491335305e-06 0.0021962918766307708317 +1.0978911404624207699e-05 +1.0431444771347480049 0.00086121939260353213164 0.004059279588003008102 +286.63158878609289104 179.04178066528515956 146.49836003344253754 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +196 1.1053539394491335305e-06 0.0029447439089289198142 +1.0978911404624207699e-05 +1.3986271031916246432 0.00018272629779208073367 0.018596268807723595473 +131.73550112256867806 122.63667230710116485 118.00875344636077102 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +197 1.1053539394491335305e-06 0.0036073925872425995972 +1.0978911404624207699e-05 +1.7133568148563382838 0.00075563849669159832353 0.0037783795914888887967 +128.81574846279312396 26.084496990571235386 338.71113698392537117 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +198 1.1053539394491335305e-06 0.0035771279460317551398 +1.0978911404624207699e-05 +1.6989824078535733509 0.0007581912267808143613 0.009405275951210755395 +94.72942154944794879 234.18898767453580945 120.281591054773343785 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +199 1.1053539394491335305e-06 0.0035858601416377149508 +1.0978911404624207699e-05 +1.7031298263804728332 0.00095813031245370230537 0.0036113709533022904374 +146.16269587061592006 115.18323006911683137 268.9501223857537866 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +200 1.1053539394491335305e-06 0.0037229358585092736563 +1.0978911404624207699e-05 +1.7682349148822829044 0.0009838787345463682714 0.0041291372775860849323 +1.7590124360204661969 28.267489335889418811 140.55607459951542637 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +201 1.1053539394491335305e-06 0.001651848468190447193 +1.0978911404624207699e-05 +0.78455720070308299086 0.00013641617523719662343 0.019427991831481188079 +130.8007266845419565 110.181352575198033605 106.00139496728822053 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +202 1.1053539394491335305e-06 0.0027063636542107174076 +1.0978911404624207699e-05 +1.2854067025640301125 0.00080079146808039282975 0.0108589657051855435765 +46.771264060493685122 32.497203081429894667 53.366093926082044163 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +203 1.1053539394491335305e-06 0.0017526730483667882186 +1.0978911404624207699e-05 +0.8324445535132765084 9.953863767625981303e-05 0.027893078159839147095 +116.997421293990342406 10.0440500665721010876 63.77086305348634454 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +204 1.1053539394491335305e-06 0.0013663552841492629148 +1.0978911404624207699e-05 +0.6489601785763914332 0.0007547835818466555732 0.028400265298084318571 +183.67148924061424964 352.2190315881923084 274.3720362883656776 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +205 1.1053539394491335305e-06 0.0016416918432609031937 +1.0978911404624207699e-05 +0.7797332393187537125 0.00076392586396015573803 0.020535102200658142402 +221.9135945165130579 194.21835875557150075 137.1052917217631375 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +206 1.1053539394491335305e-06 0.0041503320593118873817 +1.0978911404624207699e-05 +1.9712297859918352927 0.00025071053464458946825 0.014913864673346191939 +319.79234195814888153 349.28387326164579463 138.32956542869692385 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +207 1.1053539394491335305e-06 0.0017463766922070864303 +1.0978911404624207699e-05 +0.8294540542886734702 0.00021197645574354129889 0.018247549715991118885 +105.00262688872396666 232.80424007039977141 84.949104918041584256 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +208 1.1053539394491335305e-06 0.0037525194240774044083 +1.0978911404624207699e-05 +1.7822858401553345331 0.00080294190095561387456 0.023444883786946393162 +356.27577921605848132 11.694296440931752912 96.8991254541982272 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +209 1.1053539394491335305e-06 0.0033673867988775121108 +1.0978911404624207699e-05 +1.5993643554399337159 0.00045136667731674606644 0.027694134961144371149 +122.09817852172396613 331.56919146778068352 23.29494948461385917 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +210 1.1053539394491335305e-06 0.002601544338314537602 +1.0978911404624207699e-05 +1.235622021558025363 0.00052748956446136828925 0.014390768296716713662 +47.688625110281321895 218.81863989673610149 71.091061616253483635 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +211 1.1053539394491335305e-06 0.0034572364291313065697 +1.0978911404624207699e-05 +1.6420390775791540161 0.00012278039477614099018 0.013364601572277041985 +37.099447743349671214 25.711537395277808571 46.771237173257247832 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +212 1.1053539394491335305e-06 0.0014692126295102979798 +1.0978911404624207699e-05 +0.69781300769612641943 0.00020831235893638435199 0.0105463257356791569735 +45.282702057561735387 269.00333228724332457 141.05705036356778237 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +213 1.1053539394491335305e-06 0.0012506227963414040885 +1.0978911404624207699e-05 +0.5939922088058925409 0.0005111276050050357622 0.006922367726696107386 +193.68387748392143521 15.642417562054848901 17.908891439446314564 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +214 1.1053539394491335305e-06 0.0041476312180326669947 +1.0978911404624207699e-05 +1.9699470021806242404 0.00058751583022024968673 0.02496086898048604022 +247.93088919368153711 124.53408226858653052 258.30383247130868085 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +215 1.1053539394491335305e-06 0.0029875489394385165544 +1.0978911404624207699e-05 +1.4189576574520941232 0.0004909660898418238828 0.0023047035739488575846 +238.65714156827573333 164.51271027675758774 223.70978205956649276 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +216 1.1053539394491335305e-06 0.0023921676837597391677 +1.0978911404624207699e-05 +1.1361770875017920268 0.00023970389209557418876 0.026923977514779456888 +58.044978520057981086 79.883336583186263624 94.38325642624329248 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +217 1.1053539394491335305e-06 0.0026816001073412998205 +1.0978911404624207699e-05 +1.2736450795184046658 0.0007910064207286296922 0.028329584831009391177 +35.296218138003794706 187.99468430047886613 356.03440184202310093 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +218 1.1053539394491335305e-06 0.001331364345185033964 +1.0978911404624207699e-05 +0.6323409827770205549 0.00028739168058106059726 0.009019690678395121675 +338.27263574838821114 230.96954437756605216 126.72414489683163197 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +219 1.1053539394491335305e-06 0.0032121580786483427635 +1.0978911404624207699e-05 +1.5256373686388187139 0.0002481645483046961251 0.0066040479374564345708 +60.055442506513337264 346.99690826233506868 179.83160228793084912 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +220 1.1053539394491335305e-06 0.0041836951887061120445 +1.0978911404624207699e-05 +1.9870758420364929009 0.00038625231628182366255 0.014342679397268678562 +335.15357535816582413 299.7702497912664512 317.6921811746384492 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +221 1.1053539394491335305e-06 0.0007589345479948146076 +1.0978911404624207699e-05 +0.3604613715832816201 0.0008130249542857412217 0.0132530748254426049976 +199.40136765025368959 44.661857657232609142 303.65462581497462224 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +222 1.1053539394491335305e-06 0.0011646339137693681694 +1.0978911404624207699e-05 +0.55315117628902576996 0.00058320518586504638704 0.025565557231859135229 +88.51316300040842577 64.3025295238989969 307.91222823830349853 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +223 1.1053539394491335305e-06 0.004180113624911285161 +1.0978911404624207699e-05 +1.9853747527906451253 0.00094918793041903390915 0.01898810058981277582 +339.11739337814219653 326.5281938064337055 348.13403023570242567 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +224 1.1053539394491335305e-06 0.002907903028903174131 +1.0978911404624207699e-05 +1.3811292647027828551 0.00017854471163964626287 0.01757285626967412806 +324.68564254821177428 173.27665949324585881 212.81811478340242161 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +225 1.1053539394491335305e-06 0.0030424513946792162559 +1.0978911404624207699e-05 +1.4450339697923877402 0.00083102791963415834596 0.016743369452382143164 +72.37540677365963404 25.69124406916907688 36.560572732805965757 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +226 1.1053539394491335305e-06 0.0014152066797417541456 +1.0978911404624207699e-05 +0.6721624970181487857 0.0004194325219653943738 0.02466884060344504137 +312.10194954625785613 75.06433419730545609 213.15995579142392558 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +227 1.1053539394491335305e-06 0.00074807664212413927804 +1.0978911404624207699e-05 +0.35530433181877674187 0.0006196274885235493945 0.0057712733317203191233 +112.28915034524422367 308.71008190145914796 250.26404455810629202 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +228 1.1053539394491335305e-06 0.0015240975037274374223 +1.0978911404624207699e-05 +0.7238809698040012952 0.00038945125604979613423 0.0034764685260635502795 +100.18804344114703042 289.83313211239982365 254.06542835983233886 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +229 1.1053539394491335305e-06 0.0013613878149108384452 +1.0978911404624207699e-05 +0.64660084366442682757 0.0009530460362832646236 0.00080068641185022341056 +290.2519853457544059 126.69300624511200226 129.6226360067927601 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +230 1.1053539394491335305e-06 0.0007497560275515179048 +1.0978911404624207699e-05 +0.3561019678944690181 6.2453564361704131855e-05 0.008646200757580674995 +51.77354471090342969 304.1036845586841082 280.42024752572388024 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +231 1.1053539394491335305e-06 0.0020677692895160400756 +1.0978911404624207699e-05 +0.9821017585588054377 0.00092676675850360832123 0.006735246833832509404 +263.97158005173145057 104.45971633042557869 156.07773457283062157 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +232 1.1053539394491335305e-06 0.00095787667212301207604 +1.0978911404624207699e-05 +0.4549503510590604538 0.00052515340653565799334 0.019061642884057490577 +284.31553463421840888 213.15965776812532795 3.9051805628102265189 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +233 1.1053539394491335305e-06 0.0037735814895644029364 +1.0978911404624207699e-05 +1.7922894182423769749 0.00017129142152418842175 0.025294321601488489043 +70.01841692477869117 111.37830955236358932 314.1754665499188377 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +234 1.1053539394491335305e-06 0.0007847310095064422164 +1.0978911404624207699e-05 +0.37271358479856425694 0.0003231942299594923215 0.01628512097921978874 +224.60279807470539026 103.680766227009058866 197.02880532870963748 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +235 1.1053539394491335305e-06 0.0015320077326663605468 +1.0978911404624207699e-05 +0.72763798940522450565 0.00094252324011187993847 0.01692070785353958634 +73.3265875193739447 291.12276980910894508 97.10035705120415628 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +236 1.1053539394491335305e-06 0.0017944863754330241527 +1.0978911404624207699e-05 +0.85230409115664484965 3.7676763751743336793e-05 0.021560198984918479276 +216.79211114566814445 74.30371365646033155 295.38458279934900474 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +237 1.1053539394491335305e-06 0.0022077289096507573257 +1.0978911404624207699e-05 +1.0485765774655604332 0.00092835655556198343457 0.002473122352105824439 +207.85994871205122081 242.10322219988140091 164.79329124450524091 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +238 1.1053539394491335305e-06 0.0017780507509334266778 +1.0978911404624207699e-05 +0.84449787418364552494 0.00012149108796338048226 0.022688746545291040707 +105.00259291531052952 236.90787466426584729 258.88778385063903897 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +239 1.1053539394491335305e-06 0.0009633437960757620414 +1.0978911404624207699e-05 +0.4575470005380317673 0.00021936857586009873258 0.024584544078017318142 +141.46008032785664454 156.81449724386195044 299.68311915392922629 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +240 1.1053539394491335305e-06 0.0035009114125315640493 +1.0978911404624207699e-05 +1.6627828221641500761 0.00083380110219331913645 0.016595030400960369271 +68.173689361921276486 154.75502404200838669 311.70470139259805364 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +241 1.1053539394491335305e-06 0.0006470573938361051598 +1.0978911404624207699e-05 +0.30732452000176924134 0.00021338219268733039035 0.006499688048810576234 +359.85986519138003814 244.23351243313086911 98.01858239568616682 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +242 1.1053539394491335305e-06 0.00065901961137865321196 +1.0978911404624207699e-05 +0.31300606046393031345 0.00084676255106667984527 0.00028332565690039374258 +44.286246649765288907 20.245198329963866257 276.88008765639568765 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +243 1.1053539394491335305e-06 0.003806998050224602309 +1.0978911404624207699e-05 +1.8081608518475495995 0.00027745564680026148143 0.021943110426273559899 +266.24251634239237774 252.3739085379230005 283.488252705583875 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +244 1.1053539394491335305e-06 0.0029050019908176934909 +1.0978911404624207699e-05 +1.3797513959918767945 0.0008840097714271473667 0.011119175742524006062 +218.83470399188448141 307.5328062303744332 53.607438100582484708 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +245 1.1053539394491335305e-06 0.0041506206755627920196 +1.0978911404624207699e-05 +1.9713668663368231204 0.00072844746599102529454 0.014679484165514528754 +216.30736619477607974 127.43423126496449527 280.14881244016112305 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +246 1.1053539394491335305e-06 0.0007277677857464638177 +1.0978911404624207699e-05 +0.34565849576542212462 0.00059173314782810218417 0.007831481826486607831 +243.4161800124451247 113.710960693445571223 254.21623664849960278 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +247 1.1053539394491335305e-06 0.0017717047625134947208 +1.0978911404624207699e-05 +0.84148380176337678193 0.00052751195938902041473 0.014082285132383451365 +194.45249461995499018 318.31107915541781495 157.54667905036907882 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +248 1.1053539394491335305e-06 0.004169267887052278984 +1.0978911404624207699e-05 +1.980223492309989286 0.00067309606461422974606 0.017396151114924100073 +184.39865388000126245 210.62388541157199029 167.14339720153236613 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +249 1.1053539394491335305e-06 0.0019522337566574509366 +1.0978911404624207699e-05 +0.92722733394492307113 0.00058907724893334585843 0.021547757024762807104 +190.98028164050376176 93.657838091120083845 154.55544413858271469 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +250 1.1053539394491335305e-06 0.004103199565563551535 +1.0978911404624207699e-05 +1.9488438722295042815 0.00023535689304344910105 0.020747164931323787418 +235.4258950080400723 210.09500044028584398 329.76890939133602387 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +251 1.1053539394491335305e-06 0.00096993472194017628167 +1.0978911404624207699e-05 +0.46067740774293175487 6.722535556689789992e-05 0.0099933940728082100435 +339.95817378663616637 181.27958061831904502 34.322050877546928405 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +252 1.1053539394491335305e-06 0.0009661516009396302622 +1.0978911404624207699e-05 +0.45888058746597204518 0.0005025170114595179727 0.0024236270953698216615 +57.831023059303113598 16.636326610996761133 83.793613608447557795 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +253 1.1053539394491335305e-06 0.0032298402438196435668 +1.0978911404624207699e-05 +1.5340356389864386877 0.00062884613640746395407 0.013163328547706066446 +348.2310414025072305 241.09038942507976344 294.8270165591735008 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 diff --git a/python/swiftest/swiftest/init_cond.py b/python/swiftest/swiftest/init_cond.py index 297b9df68..cf8f109cc 100644 --- a/python/swiftest/swiftest/init_cond.py +++ b/python/swiftest/swiftest/init_cond.py @@ -102,7 +102,6 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): 'Pluto': np.longdouble(0.4) } - # Unit conversion factors DCONV = swiftest.AU2M / param['DU2M'] VCONV = (swiftest.AU2M / swiftest.JD2S) / (param['DU2M'] / param['TU2S']) @@ -274,8 +273,22 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): return ds -def vec2xr(param, idvals, v1, v2, v3, v4, v5, v6, GMpl=None, Rpl=None, Rhill=None, t=0.0): - +def vec2xr(param, idvals, v1, v2, v3, v4, v5, v6, GMpl=None, Rpl=None, Rhill=None, Ip_x=None, Ip_y=None, Ip_z=None, rot_x=None, rot_y=None, rot_z=None, t=0.0): + + if param['ROTATION'] == 'YES': + if Ip_x is None: + Ip_x = np.full_like(v1, 0.4) + if Ip_y is None: + Ip_y = np.full_like(v1, 0.4) + if Ip_z is None: + Ip_z = np.full_like(v1, 0.4) + if rot_x is None: + rot_x = np.full_like(v1, 0.0) + if rot_y is None: + rot_y = np.full_like(v1, 0.0) + if rot_z is None: + rot_z = np.full_like(v1, 0.0) + dims = ['time', 'id', 'vec'] if GMpl is not None: ispl = True @@ -290,18 +303,18 @@ def vec2xr(param, idvals, v1, v2, v3, v4, v5, v6, GMpl=None, Rpl=None, Rhill=Non return None clab, plab, tlab = swiftest.io.make_swiftest_labels(param) + vec = np.vstack([v1, v2, v3, v4, v5, v6]) if ispl: + vec = np.vstack([vec, GMpl, Rpl]) if param['RHILL_PRESENT'] == 'YES': - vec = np.vstack([v1, v2, v3, v4, v5, v6, GMpl, Rpl, Rhill]).T - else: - vec = np.vstack([v1, v2, v3, v4, v5, v6, GMpl, Rpl]).T - else: - vec = np.vstack([v1, v2, v3, v4, v5, v6]).T - bodyframe = np.expand_dims(vec, axis=0) + vec = np.vstack([vec, Rhill]) + if param['ROTATION'] == 'YES': + vec = np.vstack([vec, Ip_x, Ip_y, Ip_z, rot_x, rot_y, rot_z]) + bodyframe = np.expand_dims(vec.T, axis=0) if ispl: - bodyxr = xr.DataArray(bodyframe, dims=dims, coords={'time': [t], 'id': tpid, 'vec': plab}) + bodyxr = xr.DataArray(bodyframe, dims=dims, coords={'time': [t], 'id': idvals, 'vec': plab}) else: - bodyxr = xr.DataArray(bodyframe, dims=dims, coords={'time': [t], 'id': tpid, 'vec': tlab}) + bodyxr = xr.DataArray(bodyframe, dims=dims, coords={'time': [t], 'id': idvals, 'vec': tlab}) ds = bodyxr.to_dataset(dim='vec') return ds \ No newline at end of file diff --git a/python/swiftest/swiftest/simulation_class.py b/python/swiftest/swiftest/simulation_class.py index e7b8c560b..e0ec698a2 100644 --- a/python/swiftest/swiftest/simulation_class.py +++ b/python/swiftest/swiftest/simulation_class.py @@ -4,6 +4,7 @@ from swiftest import constants from datetime import date import xarray as xr +import numpy as np class Simulation: """ @@ -52,6 +53,7 @@ def __init__(self, codename="Swiftest", param_file=""): self.codename = codename if param_file != "" : self.read_param(param_file, codename) + self.param['GU'] = constants.GC / (param['DU2M'])**3 / (param['MU2KG'] * param['TU2S']**2) return @@ -73,7 +75,7 @@ def add(self, plname, date=date.today().isoformat(), idval=None): return - def addp(self, idvals, t1, t2, t3, t4, t5, t6, Gmass=None, radius=None, Rhill=None, t=0.0): + def addp(self, idvals, t1, t2, t3, t4, t5, t6, GMpl=None, Rpl=None, Rhill=None, Ip_x=None, Ip_y=None, Ip_z=None, rot_x=None, rot_y=None, rot_z=None): """ Adds a body (test particle or massive body) to the internal DataSet given a set up 6 vectors (orbital elements or cartesian state vectors, depending on the value of self.param). Input all angles in degress @@ -90,16 +92,18 @@ def addp(self, idvals, t1, t2, t3, t4, t5, t6, Gmass=None, radius=None, Rhill=No Gmass : Optional: Array of G*mass values if these are massive bodies radius : Optional: Array radius values if these are massive bodies Rhill : Optional: Array Rhill values if these are massive bodies + Ip_x,y,z : Optional: Principal axes moments of inertia + rot_x,y,z: Optional: Rotation rate vector components Returns ------- self.ds : xarray dataset """ - dsnew = init_cond.vec2xr(self.param, idvals, t1, t2, t3, t4, t5, t6, Gmass, radius, Rhill) + t = self.param['T0'] + + dsnew = init_cond.vec2xr(self.param, idvals, t1, t2, t3, t4, t5, t6, GMpl, Rpl, Rhill, Ip_x, Ip_y, Ip_z, rot_x, rot_y, rot_z, t) if dsnew is not None: self.ds = xr.combine_by_coords([self.ds, dsnew]) return - - def read_param(self, param_file, codename="Swiftest"): From f82d6480e27938f193f9d33644e4916d8c407f92 Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 16:43:24 -0400 Subject: [PATCH 3/8] Added missing GMTINY parameter --- examples/symba_chambers_2013/init_cond.py | 2 ++ 1 file changed, 2 insertions(+) diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py index 5f6782491..38b023ed1 100755 --- a/examples/symba_chambers_2013/init_cond.py +++ b/examples/symba_chambers_2013/init_cond.py @@ -54,6 +54,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) +MTINY = 1e-2 * Ms +sim.param['GMTINY'] = MTINY * GU # Define the initial orbital elements of the big and small bodies avalb = default_rng().uniform(0.3, 2.0, Nb) From 291184508afd676b52876bb045f3199f329af8db Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 16:47:23 -0400 Subject: [PATCH 4/8] Fixed central body rotation parameter formatting for input file generation. --- python/swiftest/swiftest/io.py | 8 ++------ 1 file changed, 2 insertions(+), 6 deletions(-) diff --git a/python/swiftest/swiftest/io.py b/python/swiftest/swiftest/io.py index 1116f6509..dfecd5811 100644 --- a/python/swiftest/swiftest/io.py +++ b/python/swiftest/swiftest/io.py @@ -758,12 +758,8 @@ def swiftest_xr2infile(ds, param, framenum=-1): 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) + print(Ip_xcb, Ip_ycb, Ip_zcb, file=cbfile) + print(rot_xcb, rot_ycb, rot_zcb, file=cbfile) cbfile.close() plfile = open(param['PL_IN'], 'w') From b8d31d0e54cc3c2da01c15e6176e93aabd3e0e98 Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 16:51:59 -0400 Subject: [PATCH 5/8] Removed unneeded GU parameter --- python/swiftest/swiftest/simulation_class.py | 1 - 1 file changed, 1 deletion(-) diff --git a/python/swiftest/swiftest/simulation_class.py b/python/swiftest/swiftest/simulation_class.py index e0ec698a2..af9b36475 100644 --- a/python/swiftest/swiftest/simulation_class.py +++ b/python/swiftest/swiftest/simulation_class.py @@ -53,7 +53,6 @@ def __init__(self, codename="Swiftest", param_file=""): self.codename = codename if param_file != "" : self.read_param(param_file, codename) - self.param['GU'] = constants.GC / (param['DU2M'])**3 / (param['MU2KG'] * param['TU2S']**2) return From 0456a568ef90b384611f06d2e730f757431c3a32 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Mon, 23 Aug 2021 17:04:20 -0400 Subject: [PATCH 6/8] Improved error reporting in the param reader --- src/io/io.f90 | 2 +- src/symba/symba_collision.f90 | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/io/io.f90 b/src/io/io.f90 index a4b619847..58e233dd1 100644 --- a/src/io/io.f90 +++ b/src/io/io.f90 @@ -555,7 +555,7 @@ module subroutine io_param_reader(self, unit, iotype, v_list, iostat, iomsg) read(param_value, *, err = 667, iomsg = iomsg) param%maxid case ("NPLMAX", "NTPMAX", "GMTINY", "MIN_GMFRAG", "PARTICLE_OUT", "FRAGMENTATION", "SEED", "YARKOVSKY", "YORP") ! Ignore SyMBA-specific, not-yet-implemented, or obsolete input parameters case default - write(iomsg,*) "Unknown parameter -> ",param_name + write(*,*) "Unknown parameter -> ",param_name iostat = -1 return end select diff --git a/src/symba/symba_collision.f90 b/src/symba/symba_collision.f90 index 3672f959b..fd4eae517 100644 --- a/src/symba/symba_collision.f90 +++ b/src/symba/symba_collision.f90 @@ -283,7 +283,7 @@ module function symba_collision_casemerge(system, param, family, x, v, mass, rad param%Ecollisions = param%Ecollisions + pe param%Euntracked = param%Euntracked - pe - ! Update any encounter lists that have the removedbodies in them so that they instead point to the new + ! Update any encounter lists that have the removed bodies in them so that they instead point to the new do k = 1, system%plplenc_list%nenc do j = 1, nfamily i = family(j) From ae44e29fe6459028e33dac4ec43ccc8b276a67d9 Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 17:05:12 -0400 Subject: [PATCH 7/8] Fixed bad parameter --- examples/symba_chambers_2013/init_cond.py | 1 - 1 file changed, 1 deletion(-) diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py index 38b023ed1..aa213a5ab 100755 --- a/examples/symba_chambers_2013/init_cond.py +++ b/examples/symba_chambers_2013/init_cond.py @@ -14,7 +14,6 @@ sim.param['MU2KG'] = MU2KG sim.param['TU2S'] = TU2S sim.param['DU2M'] = DU2M -sim.param['GU'] = GU # Simulation time parameters sim.param['T0'] = 0.0 From 56e745b90ce3efb8d033e22d7bf5039c10487990 Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 17:08:08 -0400 Subject: [PATCH 8/8] Added minimum fragment mass to input file --- examples/symba_chambers_2013/init_cond.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py index aa213a5ab..5ea4b1205 100755 --- a/examples/symba_chambers_2013/init_cond.py +++ b/examples/symba_chambers_2013/init_cond.py @@ -53,8 +53,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) -MTINY = 1e-2 * Ms -sim.param['GMTINY'] = MTINY * GU +sim.param['GMTINY'] = 1e-2 * GU * Ms +sim.param['MIN_GMFRAG'] = 1e-4 * GU * Ms # Define the initial orbital elements of the big and small bodies avalb = default_rng().uniform(0.3, 2.0, Nb)