Basic
Before working on a simulation, it is recommanded to be familiar with dependencies of AstroNbodySim.jl
The whole project is well documented. Use help? in REPL is the fastest way to get help if you forget the arguments or keywords of interfaces. For example,
julia> using AstroIO
help?> read_gadget2
search: read_gadget2 read_gadget2_pos read_gadget2_jld
read_gadget2(filename::AbstractString, units, fileunits = uGadget2; kw...)
Return a Tuple of header and particle data in snapshot file.
units is supported by PhysicalParticles: uSI, uCGS, uAstro,
uGadget2, nothing. fileunits is the internal units in the file,
and will be converted to units while reading the file.
Keywords
≡≡≡≡≡≡≡≡≡≡
• acc::Bool = false : read acceleration data if exist
• pot::Bool = false : read potential data if existPhysical vectors and particles
julia> # PhysicalParticles defines the basic physical vector and particle type using PhysicalParticlesjulia> # Set the default units to astrophysics. See `Unitful.preferunits` for more info astro()julia> # using UnitfulAstro is necessary, if you are using astrophysical units using UnitfulAstrojulia> # Define a physical vector a = PVector(3.0u"kpc", 4.0u"kpc", 12.0u"kpc")PVector(3.0 kpc, 4.0 kpc, 12.0 kpc)julia> b = PVector(1.0, 1.0, 1.0, u"kpc")PVector(1.0 kpc, 1.0 kpc, 1.0 kpc)julia> c = PVector2D()PVector2D{Float64}(0.0, 0.0)julia> d = PVector2D(0.0, 1.0)PVector2D{Float64}(0.0, 1.0)julia> # Basic algebra a * b19.0 kpc^2julia> c + dPVector2D{Float64}(0.0, 1.0)julia> norm(a)13.0 kpcjulia> normalize(a)PVector(0.23076923076923078 kpc, 0.3076923076923077 kpc, 0.9230769230769231 kpc)julia> # generate an array of physical vector points = rand(PVector{Float64}, 5) * u"kpc"5-element Vector{PVector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}}}: PVector(0.28510586229767176 kpc, 0.3624402810242151 kpc, 0.7529396250098489 kpc) PVector(0.759691455791893 kpc, 0.5165040074051774 kpc, 0.12978654761275865 kpc) PVector(0.11065640353676676 kpc, 0.5364186708379797 kpc, 0.2876987304007538 kpc) PVector(0.3288181631401986 kpc, 0.9813039245486699 kpc, 0.6515687063302311 kpc) PVector(0.8948097445167615 kpc, 0.4349660509797524 kpc, 0.8911796746901306 kpc)julia> p = randn_pvector2d(5)5-element Vector{PVector2D{Float64}}: PVector2D{Float64}(1.2479896129026782, 1.9284468356312068) PVector2D{Float64}(0.6401530366732445, 1.5797647736810008) PVector2D{Float64}(-0.050551960083584244, 1.1574890997858243) PVector2D{Float64}(0.18326045465612875, 0.4187339632534888) PVector2D{Float64}(-0.3541846848379587, -1.4300220355545683)julia> mean(p)PVector2D{Float64}(0.33333329186210175, 0.7308825273593905)julia> PhysicalParticles.center(p)PVector2D{Float64}(0.4469024640323598, 0.24921240003831924)julia> median(p)PVector2D{Float64}(0.18326045465612875, 1.1574890997858243)julia> # Define particles particles = [Star(uAstro, id = i) for i in 1:5]5-element Vector{Star{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Int64}}: Star 1 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 2 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 3 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 4 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 5 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2julia> # Assign particle positions with array of points assign_particles(particles, :Pos, points)julia> # average position average(particles, :Pos)PVector(0.4758163258566584 kpc, 0.5663265869591589 kpc, 0.5426346568087446 kpc)julia> # average position by mass assign_particles(particles, :Mass, rand(5) * u"Msun")julia> averagebymass(particles, :Pos)PVector(0.5854194302792788 kpc, 0.6043950993662042 kpc, 0.5654778105215538 kpc)julia> # StructArray is more efficient to manipulate field of struct in an array StructArray(particles)5-element StructArray(::Vector{PVector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}}, ::Vector{Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}}, ::Vector{Int64}, ::Vector{Collection}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}) with eltype Star{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Int64}: Star 1 STAR: Pos = PVector(0.28510586229767176 kpc, 0.3624402810242151 kpc, 0.7529396250098489 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.29898521664003597 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 2 STAR: Pos = PVector(0.759691455791893 kpc, 0.5165040074051774 kpc, 0.12978654761275865 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.834798065738929 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 3 STAR: Pos = PVector(0.11065640353676676 kpc, 0.5364186708379797 kpc, 0.2876987304007538 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.26901077366512294 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 4 STAR: Pos = PVector(0.3288181631401986 kpc, 0.9813039245486699 kpc, 0.6515687063302311 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.8606199174993377 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 5 STAR: Pos = PVector(0.8948097445167615 kpc, 0.4349660509797524 kpc, 0.8911796746901306 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.9465709958678674 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2julia> # Or construct directly s = StructArray(Star(uAstro, id = i) for i in 1:5)5-element StructArray(::Vector{PVector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}}, ::Vector{Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}}, ::Vector{Int64}, ::Vector{Collection}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}) with eltype Star{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Int64}: Star 1 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 2 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 3 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 4 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 5 STAR: Pos = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2julia> # StructArray is still supported by assign_particles, averagebymass, etc. # It is much more convenient to use dot operations s.Pos .= points5-element Vector{PVector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}}}: PVector(0.28510586229767176 kpc, 0.3624402810242151 kpc, 0.7529396250098489 kpc) PVector(0.759691455791893 kpc, 0.5165040074051774 kpc, 0.12978654761275865 kpc) PVector(0.11065640353676676 kpc, 0.5364186708379797 kpc, 0.2876987304007538 kpc) PVector(0.3288181631401986 kpc, 0.9813039245486699 kpc, 0.6515687063302311 kpc) PVector(0.8948097445167615 kpc, 0.4349660509797524 kpc, 0.8911796746901306 kpc)julia> s.Mass .= rand(5) * u"Msun"5-element Vector{Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}}: 0.4462748845220884 M⊙ 0.734937309128467 M⊙ 0.6185303745059065 M⊙ 0.11487296247070244 M⊙ 0.006332378682999851 M⊙julia> s5-element StructArray(::Vector{PVector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}}, ::Vector{Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}}, ::Vector{Int64}, ::Vector{Collection}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}) with eltype Star{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Int64}: Star 1 STAR: Pos = PVector(0.28510586229767176 kpc, 0.3624402810242151 kpc, 0.7529396250098489 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.4462748845220884 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 2 STAR: Pos = PVector(0.759691455791893 kpc, 0.5165040074051774 kpc, 0.12978654761275865 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.734937309128467 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 3 STAR: Pos = PVector(0.11065640353676676 kpc, 0.5364186708379797 kpc, 0.2876987304007538 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.6185303745059065 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 4 STAR: Pos = PVector(0.3288181631401986 kpc, 0.9813039245486699 kpc, 0.6515687063302311 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.11487296247070244 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 5 STAR: Pos = PVector(0.8948097445167615 kpc, 0.4349660509797524 kpc, 0.8911796746901306 kpc), Vel = PVector(0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1, 0.0 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 0.006332378682999851 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2
Generate initial conditions
julia> using AstroICjulia> config = PlummerStarCluster()Config of Plummer Star Cluster Initial Conditions: Gravity Model: AstroIC.Newton Particle Collection: STAR Number of Samples: 1000 Virial Radius: 0.01 kpc Total Mass: 100000.0 M⊙julia> particles = generate(config)1000-element StructArray(::Vector{PVector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}}, ::Vector{Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}}, ::Vector{Int64}, ::Vector{Collection}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}) with eltype Star{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Int64}: Star 1 STAR: Pos = PVector(-0.0016653158086158196 kpc, 0.00563252531227182 kpc, -0.008075736676002691 kpc), Vel = PVector(-2.3033025011029595 kpc Gyr^-1, -0.8554547657742699 kpc Gyr^-1, -1.0360885986522148 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 2 STAR: Pos = PVector(0.0026139980081741417 kpc, -0.004439589874897115 kpc, 0.0028663568610049215 kpc), Vel = PVector(-1.1430551434050997 kpc Gyr^-1, 1.2590390179017503 kpc Gyr^-1, -4.48983911205135 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 3 STAR: Pos = PVector(-0.00011606994079237597 kpc, 0.0027060738288232668 kpc, -0.003591062500196339 kpc), Vel = PVector(4.2293131123496455 kpc Gyr^-1, -1.8059400337137372 kpc Gyr^-1, -3.66907025861748 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 4 STAR: Pos = PVector(0.011658586688654934 kpc, -0.028686199375439828 kpc, 0.020743292513330524 kpc), Vel = PVector(-0.4158265743331983 kpc Gyr^-1, -1.2819538808873696 kpc Gyr^-1, 0.9338038899230501 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 5 STAR: Pos = PVector(-0.01240280603124357 kpc, -0.0024693452197793045 kpc, 0.0038751238284525295 kpc), Vel = PVector(3.6821510857717312 kpc Gyr^-1, -4.4829745349478545 kpc Gyr^-1, -0.024869726401783785 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 6 STAR: Pos = PVector(0.0018861869518809814 kpc, -0.004388466271410916 kpc, 0.001973000054921772 kpc), Vel = PVector(-0.22016998625353207 kpc Gyr^-1, -3.011759902011347 kpc Gyr^-1, -1.0935402498754094 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 7 STAR: Pos = PVector(0.0028895755059052 kpc, 0.01771750243226564 kpc, -0.0035142116413609216 kpc), Vel = PVector(-0.5080239240211444 kpc Gyr^-1, -1.568232399063682 kpc Gyr^-1, -2.3763439519155005 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 8 STAR: Pos = PVector(0.005450603837314128 kpc, -0.001787910149787971 kpc, -0.007419379764904896 kpc), Vel = PVector(2.0487069806460654 kpc Gyr^-1, -4.303919271904533 kpc Gyr^-1, 2.160346888250885 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 9 STAR: Pos = PVector(-0.0007041065887441831 kpc, -0.004472728331805591 kpc, -0.002363530776846375 kpc), Vel = PVector(4.692387193252354 kpc Gyr^-1, -1.936212927191336 kpc Gyr^-1, 2.4934461437767648 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 10 STAR: Pos = PVector(-0.003613067589961985 kpc, -0.010368351032841458 kpc, -0.008593581237645274 kpc), Vel = PVector(1.5499837959903162 kpc Gyr^-1, 2.0987084298977052 kpc Gyr^-1, 1.796589679491565 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 ⋮ Star 992 STAR: Pos = PVector(0.006990124380139745 kpc, -0.005206294985838018 kpc, -0.002206311890571247 kpc), Vel = PVector(-3.5604279883685304 kpc Gyr^-1, 1.3003661418054593 kpc Gyr^-1, 0.7123912608909204 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 993 STAR: Pos = PVector(0.006474836565285428 kpc, -0.010527451658388786 kpc, 0.007684641145100218 kpc), Vel = PVector(1.2158855275839824 kpc Gyr^-1, 1.106630448492298 kpc Gyr^-1, 0.01668133441734234 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 994 STAR: Pos = PVector(0.006443052503175608 kpc, -0.02206819754103854 kpc, -0.001289049509441723 kpc), Vel = PVector(0.6375651065525745 kpc Gyr^-1, -0.36866442600692 kpc Gyr^-1, 0.8531716110628953 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 995 STAR: Pos = PVector(-0.009907369804914282 kpc, -0.0032678595672352385 kpc, 0.005976296373120448 kpc), Vel = PVector(-2.2141064037830156 kpc Gyr^-1, -0.5043596314009984 kpc Gyr^-1, 1.2928382885031473 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 996 STAR: Pos = PVector(-0.0004141275976366911 kpc, 0.007617504851048121 kpc, -0.004019833892218843 kpc), Vel = PVector(2.9389660695472575 kpc Gyr^-1, 2.110054100369787 kpc Gyr^-1, 0.09877031047561456 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 997 STAR: Pos = PVector(-0.002940508169955408 kpc, -0.004923748195212407 kpc, 0.002436002591508488 kpc), Vel = PVector(0.8204676468015072 kpc Gyr^-1, -3.568177508689484 kpc Gyr^-1, -3.7255961767649217 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 998 STAR: Pos = PVector(0.0015214009049604894 kpc, -0.000644519982359962 kpc, 0.0022503130047443944 kpc), Vel = PVector(0.401775826200449 kpc Gyr^-1, -2.089103368006141 kpc Gyr^-1, 0.4709865659830987 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 999 STAR: Pos = PVector(-0.008167684132682597 kpc, 0.017623712287212556 kpc, -0.021503119301842642 kpc), Vel = PVector(1.6899848263085542 kpc Gyr^-1, -1.2215393402854389 kpc Gyr^-1, 1.2770351952914285 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 1000 STAR: Pos = PVector(0.0054037785918452895 kpc, 0.010081646800908155 kpc, 0.0010623079385243107 kpc), Vel = PVector(-2.446974157833299 kpc Gyr^-1, 2.343289028693592 kpc Gyr^-1, 1.0535010716453388 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2
Visualization
using AstroPlot
fig = plot_makie(particles)
fig
Snapshot File I/O
julia> using AstroIOjulia> if !isdir("output/") mkpath("output/") end"output"julia> write_csv("output/basic.csv", particles)truejulia> write_jld("output/basic.jld2", particles)truejulia> write_gadget2("output/basic.gadget2", particles) # This would generate a header automaticallytruejulia> # generate a header and modify it header = HeaderGadget2(particles)Gadget2 Header: Particle Info : | Type | Amount, Mass |-----------|------------------------- | Gas | 0, 0.0 M⊙ | Halo | 0, 0.0 M⊙ | Disk | 0, 0.0 M⊙ | Bulge | 0, 0.0 M⊙ | Star | 1000, 0.0 M⊙ | BlackHole | 0, 0.0 M⊙ (If zero, mass would be read from file) Start time: 0.0 Redshift: 0.0julia> header.time = 0.1 # Gyr0.1julia> write_gadget2("output/basicwithheader.gadget2", header, particles, uGadget2) # write in Gadget2 units (default)truejulia> # now load the files we had just written h, d = read_gadget2("output/basic.gadget2", uAstro)(Gadget2 Header: Particle Info : | Type | Amount, Mass |-----------|------------------------- | Gas | 0, 0.0 M⊙ | Halo | 0, 0.0 M⊙ | Disk | 0, 0.0 M⊙ | Bulge | 0, 0.0 M⊙ | Star | 1000, 0.0 M⊙ | BlackHole | 0, 0.0 M⊙ (If zero, mass would be read from file) Start time: 0.0 Redshift: 0.0 , AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}[AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.0016653158f0 kpc, 0.0056325253f0 kpc, -0.008075736f0 kpc), PVector(-2.3033025f0 kpc Gyr^-1, -0.85545474f0 kpc Gyr^-1, -1.0360886f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 1, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.002613998f0 kpc, -0.00443959f0 kpc, 0.0028663569f0 kpc), PVector(-1.1430551f0 kpc Gyr^-1, 1.2590389f0 kpc Gyr^-1, -4.489839f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 2, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.00011606994f0 kpc, 0.002706074f0 kpc, -0.0035910625f0 kpc), PVector(4.229313f0 kpc Gyr^-1, -1.80594f0 kpc Gyr^-1, -3.6690702f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 3, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.011658587f0 kpc, -0.0286862f0 kpc, 0.020743292f0 kpc), PVector(-0.41582656f0 kpc Gyr^-1, -1.2819538f0 kpc Gyr^-1, 0.9338039f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 4, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.012402806f0 kpc, -0.0024693452f0 kpc, 0.0038751238f0 kpc), PVector(3.682151f0 kpc Gyr^-1, -4.4829745f0 kpc Gyr^-1, -0.024869727f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 5, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.001886187f0 kpc, -0.0043884665f0 kpc, 0.0019730001f0 kpc), PVector(-0.22016999f0 kpc Gyr^-1, -3.01176f0 kpc Gyr^-1, -1.0935402f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 6, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.0028895754f0 kpc, 0.017717503f0 kpc, -0.0035142116f0 kpc), PVector(-0.5080239f0 kpc Gyr^-1, -1.5682324f0 kpc Gyr^-1, -2.376344f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 7, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.005450604f0 kpc, -0.0017879101f0 kpc, -0.00741938f0 kpc), PVector(2.048707f0 kpc Gyr^-1, -4.3039193f0 kpc Gyr^-1, 2.160347f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 8, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.0007041066f0 kpc, -0.0044727284f0 kpc, -0.0023635307f0 kpc), PVector(4.692387f0 kpc Gyr^-1, -1.936213f0 kpc Gyr^-1, 2.493446f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 9, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.0036130676f0 kpc, -0.010368351f0 kpc, -0.008593582f0 kpc), PVector(1.5499837f0 kpc Gyr^-1, 2.0987086f0 kpc Gyr^-1, 1.7965897f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 10, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K) … AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.002362137f0 kpc, 0.0017144242f0 kpc, -0.0024365047f0 kpc), PVector(1.7473718f0 kpc Gyr^-1, -9.259498f0 kpc Gyr^-1, 0.33293977f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 991, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.0069901245f0 kpc, -0.005206295f0 kpc, -0.0022063118f0 kpc), PVector(-3.560428f0 kpc Gyr^-1, 1.3003662f0 kpc Gyr^-1, 0.71239126f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 992, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.0064748367f0 kpc, -0.0105274515f0 kpc, 0.007684641f0 kpc), PVector(1.2158855f0 kpc Gyr^-1, 1.1066304f0 kpc Gyr^-1, 0.016681336f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 993, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.0064430526f0 kpc, -0.022068197f0 kpc, -0.0012890495f0 kpc), PVector(0.63756514f0 kpc Gyr^-1, -0.3686644f0 kpc Gyr^-1, 0.85317165f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 994, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.0099073695f0 kpc, -0.0032678596f0 kpc, 0.0059762965f0 kpc), PVector(-2.2141066f0 kpc Gyr^-1, -0.50435966f0 kpc Gyr^-1, 1.2928383f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 995, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.0004141276f0 kpc, 0.007617505f0 kpc, -0.004019834f0 kpc), PVector(2.938966f0 kpc Gyr^-1, 2.110054f0 kpc Gyr^-1, 0.098770306f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 996, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.002940508f0 kpc, -0.0049237483f0 kpc, 0.0024360025f0 kpc), PVector(0.82046765f0 kpc Gyr^-1, -3.5681775f0 kpc Gyr^-1, -3.7255962f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 997, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.0015214009f0 kpc, -0.00064452f0 kpc, 0.002250313f0 kpc), PVector(0.40177584f0 kpc Gyr^-1, -2.0891035f0 kpc Gyr^-1, 0.47098657f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 998, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(-0.008167684f0 kpc, 0.017623711f0 kpc, -0.021503119f0 kpc), PVector(1.6899848f0 kpc Gyr^-1, -1.2215393f0 kpc Gyr^-1, 1.2770352f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 999, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K), AstroIO.Gadget2Particle{Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float32, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Float32, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-3, Unitful.FreeUnits{(kpc^-3, M⊙), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float32, 𝚯, Unitful.FreeUnits{(K,), 𝚯, nothing}}, Unitful.Quantity{Float32, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float32, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(K^-1, kpc^2, M⊙, Gyr^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float32, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kpc^-1, M⊙, Gyr^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float32, 𝐓^-1, Unitful.FreeUnits{(Gyr^-1,), 𝐓^-1, nothing}}, Int32}(PVector(0.0054037785f0 kpc, 0.010081647f0 kpc, 0.001062308f0 kpc), PVector(-2.446974f0 kpc Gyr^-1, 2.343289f0 kpc Gyr^-1, 1.053501f0 kpc Gyr^-1), PVector(0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc Gyr^-2), 100.0f0 M⊙, 1000, STAR, 0, 0, 0, 0.0f0 kpc^2 Gyr^-2, 0.0f0 kpc Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-2, 0.0f0 M⊙ kpc^-3, 0.0f0 kpc, 0.0f0, 0.0f0, 0, PVector(0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1, 0.0f0 kpc Gyr^-1), 0.0f0 Gyr^-1, 0.0f0 Gyr^-1, 0.0f0 kpc, 0.0f0 M⊙ kpc^-1 Gyr^-2, 0.0f0 kpc^2 M⊙ K^-1 Gyr^-3, 0.0f0 kpc Gyr^-1, 0.0f0 kpc^2 Gyr^-2, 0.0f0 K)])julia> d = read_jld("output/basic.jld2")1000-element StructArray(::Vector{PVector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}}, ::Vector{Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}}, ::Vector{Int64}, ::Vector{Collection}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}}) with eltype Star{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(kpc,), 𝐋, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(kpc, Gyr^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(kpc, Gyr^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(M⊙,), 𝐌, nothing}}, Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(kpc^2, Gyr^-2), 𝐋^2 𝐓^-2, nothing}}, Int64}: Star 1 STAR: Pos = PVector(-0.0016653158086158196 kpc, 0.00563252531227182 kpc, -0.008075736676002691 kpc), Vel = PVector(-2.3033025011029595 kpc Gyr^-1, -0.8554547657742699 kpc Gyr^-1, -1.0360885986522148 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 2 STAR: Pos = PVector(0.0026139980081741417 kpc, -0.004439589874897115 kpc, 0.0028663568610049215 kpc), Vel = PVector(-1.1430551434050997 kpc Gyr^-1, 1.2590390179017503 kpc Gyr^-1, -4.48983911205135 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 3 STAR: Pos = PVector(-0.00011606994079237597 kpc, 0.0027060738288232668 kpc, -0.003591062500196339 kpc), Vel = PVector(4.2293131123496455 kpc Gyr^-1, -1.8059400337137372 kpc Gyr^-1, -3.66907025861748 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 4 STAR: Pos = PVector(0.011658586688654934 kpc, -0.028686199375439828 kpc, 0.020743292513330524 kpc), Vel = PVector(-0.4158265743331983 kpc Gyr^-1, -1.2819538808873696 kpc Gyr^-1, 0.9338038899230501 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 5 STAR: Pos = PVector(-0.01240280603124357 kpc, -0.0024693452197793045 kpc, 0.0038751238284525295 kpc), Vel = PVector(3.6821510857717312 kpc Gyr^-1, -4.4829745349478545 kpc Gyr^-1, -0.024869726401783785 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 6 STAR: Pos = PVector(0.0018861869518809814 kpc, -0.004388466271410916 kpc, 0.001973000054921772 kpc), Vel = PVector(-0.22016998625353207 kpc Gyr^-1, -3.011759902011347 kpc Gyr^-1, -1.0935402498754094 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 7 STAR: Pos = PVector(0.0028895755059052 kpc, 0.01771750243226564 kpc, -0.0035142116413609216 kpc), Vel = PVector(-0.5080239240211444 kpc Gyr^-1, -1.568232399063682 kpc Gyr^-1, -2.3763439519155005 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 8 STAR: Pos = PVector(0.005450603837314128 kpc, -0.001787910149787971 kpc, -0.007419379764904896 kpc), Vel = PVector(2.0487069806460654 kpc Gyr^-1, -4.303919271904533 kpc Gyr^-1, 2.160346888250885 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 9 STAR: Pos = PVector(-0.0007041065887441831 kpc, -0.004472728331805591 kpc, -0.002363530776846375 kpc), Vel = PVector(4.692387193252354 kpc Gyr^-1, -1.936212927191336 kpc Gyr^-1, 2.4934461437767648 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 10 STAR: Pos = PVector(-0.003613067589961985 kpc, -0.010368351032841458 kpc, -0.008593581237645274 kpc), Vel = PVector(1.5499837959903162 kpc Gyr^-1, 2.0987084298977052 kpc Gyr^-1, 1.796589679491565 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 ⋮ Star 992 STAR: Pos = PVector(0.006990124380139745 kpc, -0.005206294985838018 kpc, -0.002206311890571247 kpc), Vel = PVector(-3.5604279883685304 kpc Gyr^-1, 1.3003661418054593 kpc Gyr^-1, 0.7123912608909204 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 993 STAR: Pos = PVector(0.006474836565285428 kpc, -0.010527451658388786 kpc, 0.007684641145100218 kpc), Vel = PVector(1.2158855275839824 kpc Gyr^-1, 1.106630448492298 kpc Gyr^-1, 0.01668133441734234 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 994 STAR: Pos = PVector(0.006443052503175608 kpc, -0.02206819754103854 kpc, -0.001289049509441723 kpc), Vel = PVector(0.6375651065525745 kpc Gyr^-1, -0.36866442600692 kpc Gyr^-1, 0.8531716110628953 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 995 STAR: Pos = PVector(-0.009907369804914282 kpc, -0.0032678595672352385 kpc, 0.005976296373120448 kpc), Vel = PVector(-2.2141064037830156 kpc Gyr^-1, -0.5043596314009984 kpc Gyr^-1, 1.2928382885031473 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 996 STAR: Pos = PVector(-0.0004141275976366911 kpc, 0.007617504851048121 kpc, -0.004019833892218843 kpc), Vel = PVector(2.9389660695472575 kpc Gyr^-1, 2.110054100369787 kpc Gyr^-1, 0.09877031047561456 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 997 STAR: Pos = PVector(-0.002940508169955408 kpc, -0.004923748195212407 kpc, 0.002436002591508488 kpc), Vel = PVector(0.8204676468015072 kpc Gyr^-1, -3.568177508689484 kpc Gyr^-1, -3.7255961767649217 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 998 STAR: Pos = PVector(0.0015214009049604894 kpc, -0.000644519982359962 kpc, 0.0022503130047443944 kpc), Vel = PVector(0.401775826200449 kpc Gyr^-1, -2.089103368006141 kpc Gyr^-1, 0.4709865659830987 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 999 STAR: Pos = PVector(-0.008167684132682597 kpc, 0.017623712287212556 kpc, -0.021503119301842642 kpc), Vel = PVector(1.6899848263085542 kpc Gyr^-1, -1.2215393402854389 kpc Gyr^-1, 1.2770351952914285 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2 Star 1000 STAR: Pos = PVector(0.0054037785918452895 kpc, 0.010081646800908155 kpc, 0.0010623079385243107 kpc), Vel = PVector(-2.446974157833299 kpc Gyr^-1, 2.343289028693592 kpc Gyr^-1, 1.0535010716453388 kpc Gyr^-1), Acc = PVector(0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2, 0.0 kpc Gyr^-2), Mass = 100.0 M⊙, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 kpc^2 Gyr^-2, OldAcc = 0.0 kpc Gyr^-2