This package is a subproject of GadJet.jl and provides some basic IO functionality to work with the SPH code "Gadget" by Volker Springel (doi:10.1111/j.1365-2966.2005.09655.x).
Any help and contribution is greatly appreciated, as this is still a work in progress.
If you want to read a simulation snapshot into memory with GadJet.jl, it's as easy as this:
data = read_snap(filename)
This will return a dictionary with the header information in
data["Header"] and the blocks sorted by particle type.
As an example, this is how you would access the positions of the gas particles:
If you only want to read a specific block for a single particle type (e.g. positions of gas particles) you can use the function with a specified blockname and particle type like so:
pos = read_snap(filename, "POS", 0)
This will return an array of the datatype of your simulation, usually Float32.