Load Data from Existing Outputs

using Mera

Load data from a sequence of snapshots

# mw_L10 has one complete output (300); 301 is an incomplete output
for i in [300]
    info = getinfo(output=i, "/Volumes/FASTStorage/Simulations/Mera-Tests/mw_L10", verbose=false)
    #...gethydro(info)...getparticles(info)... etc.
end

Load data from existing simulations in a given folder

List the content of a given folder:

path = "/Volumes/FASTStorage/Simulations/Mera-Tests/mw_L10"
readdir(path)

3-element Vector{String}:
 ".DS_Store"
 "output_00300"
 "output_00301"

Get the relevant simulation output-numbers:

N = checkoutputs(path);

Outputs - existing: 1 betw. 300:300 - missing: 1

N.outputs

1-element Vector{Int64}:
 300

List of empty simulation folders:

N.miss

1-element Vector{Int64}:
 301

Load the data:

for i in N.outputs
    println("Output: $i")
    info = getinfo(output=i, path, verbose=false)
    #...gethydro(info)...getparticles(info)... etc.
end

Output: 300

Get the physical time of all existing outputs:

gettime.(N.outputs, path, :Myr)

1-element Vector{Float64}:
 445.8861174695

One call: timeseries

The whole pattern above — discover the outputs, load each one, read its physical time, collect a quantity — is what timeseries automates into a single call. You give it a reducer (data -> scalar | NamedTuple); it loads one snapshot at a time (RAM-safe) and returns one table with an output column and a physical time column in Myr (the same gettime(:Myr) shown above), plus redshift/aexp columns for a cosmological run:

ts = timeseries(path, d -> (mass = msum(d, :Msol), rho_max = maximum(getvar(d, :rho))))
# columns: output | time [Myr] | mass | rho_max   (one row per output)

Use the manual loop when you want full control per snapshot; reach for timeseries when you just want X(t) as a table. See Time Series for output selection, memory control, mera-file and cosmological runs, projections-over-time, and plotting.