Solar System

julia> using AstroIC, Dates
julia> solarsystem(now())9-element StructArray(::Vector{PVector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(m,), 𝐋, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(m, s^-1), 𝐋 𝐓^-1, nothing}}}}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(m, s^-2), 𝐋 𝐓^-2, nothing}}}}, ::Vector{Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(kg,), 𝐌, nothing}}}, ::Vector{Int64}, ::Vector{Collection}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(m^2, s^-2), 𝐋^2 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(m, s^-2), 𝐋 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(kg, K^-1, m^2, s^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐌 𝐋^-3, Unitful.FreeUnits{(kg, m^-3), 𝐌 𝐋^-3, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(m,), 𝐋, nothing}}}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Int64}, ::Vector{PVector{Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(m, s^-1), 𝐋 𝐓^-1, nothing}}}}, ::Vector{Unitful.Quantity{Float64, 𝐓^-1, Unitful.FreeUnits{(s^-1,), 𝐓^-1, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐓^-1, Unitful.FreeUnits{(s^-1,), 𝐓^-1, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(m,), 𝐋, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kg, m^-1, s^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(kg, K^-1, m^2, s^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}}, ::Vector{Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(m, s^-1), 𝐋 𝐓^-1, nothing}}}) with eltype Star{Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(m,), 𝐋, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-1, Unitful.FreeUnits{(m, s^-1), 𝐋 𝐓^-1, nothing}}, Unitful.Quantity{Float64, 𝐋 𝐓^-2, Unitful.FreeUnits{(m, s^-2), 𝐋 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐌, Unitful.FreeUnits{(kg,), 𝐌, nothing}}, Unitful.Quantity{Float64, 𝐋^2 𝐓^-2, Unitful.FreeUnits{(m^2, s^-2), 𝐋^2 𝐓^-2, nothing}}, Float64, Unitful.Quantity{Float64, 𝐋^2 𝐌 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(kg, K^-1, m^2, s^-2), 𝐋^2 𝐌 𝚯^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐌 𝐋^-3, Unitful.FreeUnits{(kg, m^-3), 𝐌 𝐋^-3, nothing}}, Unitful.Quantity{Float64, 𝐋, Unitful.FreeUnits{(m,), 𝐋, nothing}}, Unitful.Quantity{Float64, 𝐋^2 𝐌 𝚯^-1 𝐓^-3, Unitful.FreeUnits{(kg, K^-1, m^2, s^-3), 𝐋^2 𝐌 𝚯^-1 𝐓^-3, nothing}}, Unitful.Quantity{Float64, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(kg, m^-1, s^-2), 𝐌 𝐋^-1 𝐓^-2, nothing}}, Unitful.Quantity{Float64, 𝐓^-1, Unitful.FreeUnits{(s^-1,), 𝐓^-1, nothing}}, Int64}:
 Star 1 STAR: Pos = PVector(-5.892357845115905e10 m, 1.313274332518838e9 m, -3.564815354332039e10 m), Vel = PVector(-23042.31785583496 m s^-1, 4555.733298540115 m s^-1, 32036.04567718506 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 3.3011e23 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1
 Star 2 STAR: Pos = PVector(6.642650536070359e10 m, -2.4628036214706607e9 m, 8.54399719443391e10 m), Vel = PVector(27501.92010498047 m s^-1, -25.482051849365234 m s^-1, -21682.1609954834 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 4.8675e24 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1
 Star 3 STAR: Pos = PVector(1.3260680319233925e11 m, -5.926262433968833e6 m, 6.4914145433069626e10 m), Vel = PVector(12985.274475097656 m s^-1, -1.1817566640675068 m s^-1, -27239.87921142578 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 5.97237e24 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1
 Star 4 STAR: Pos = PVector(-1.450800353600163e11 m, -9.440379738459054e8 m, -1.8667592990111627e11 m), Vel = PVector(-17067.772186279297 m s^-1, 337.87816858291626 m s^-1, 15787.40185546875 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 3.3011e23 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1
 Star 5 STAR: Pos = PVector(-3.082499282740924e11 m, 5.758093723800179e9 m, 6.819805537561738e11 m), Vel = PVector(12528.296569824219 m s^-1, -161.32886123657227 m s^-1, 5226.3546142578125 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 1.89819e27 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1
 Star 6 STAR: Pos = PVector(-1.0746573085398428e12 m, 1.590045570125274e10 m, 1.0212825548058986e12 m), Vel = PVector(6634.746337890625 m s^-1, 175.95999145507812 m s^-1, 6483.4405517578125 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 5.6834e25 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1
 Star 7 STAR: Pos = PVector(2.0025120003788352e12 m, -1.4013012005645533e10 m, 2.166965696846139e12 m), Vel = PVector(4684.336669921875 m s^-1, 18.90427017211914 m s^-1, -4672.373046875 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 8.6813e25 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1
 Star 8 STAR: Pos = PVector(-6.270106679921497e11 m, 1.2493571225714912e10 m, 4.43071784391448e12 m), Vel = PVector(5413.958251953125 m s^-1, -89.79300689697266 m s^-1, 728.6328125 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 1.02413e26 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1
 Star 0 STAR: Pos = PVector(0.0 m, 0.0 m, 0.0 m), Vel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), Acc = PVector(0.0 m s^-2, 0.0 m s^-2, 0.0 m s^-2), Mass = 1.988409870698051e30 kg, Ti_endstep = 0, Ti_begstep = 0, Potential = 0.0 m^2 s^-2, OldAcc = 0.0 m s^-2, Entropy = 0.0 kg m^2 K^-1 s^-2, Density = 0.0 kg m^-3, Hsml = 0.0 m, Left = 0.0, Right = 0.0, NumNgbFound = 0, RotVel = PVector(0.0 m s^-1, 0.0 m s^-1, 0.0 m s^-1), DivVel = 0.0 s^-1, CurlVel = 0.0 s^-1, dHsmlRho = 0.0 m, Pressure = 0.0 kg m^-1 s^-2, DtEntropy = 0.0 kg m^2 K^-1 s^-3, MaxSignalVel = 0.0 m s^-1