Instructions for reading in the FLASH data (on small scales, 1k^3 and less)

The simulation setup is here: https://www.nature.com/articles/s41550-025-02551-5

For the small datasets, probably the easiest way to directly read them is using my PLASMAtools repo. 

You can find it here:

https://github.com/AstroJames/PLASMAtools

I have provided a reader.py file, that shows an example of how to read the data in. 

FLASH is stored in its native block decomposition. This reader will read it into (x,y,z) coords. 

I would highly recommend making sure you have numba when you use the libray. 

Vector fields are read in such that:

(3, x, y, z)

Scalar fields are read in such that:

(1, x, y, z)

Let me know if you have any questions.