- Add `cpu_list` argument to load function to enable selective loading of cpu files: `data = osyris.load(71, scale="au", path="data", cpu_list=[5, 9, 20, 30, 60])`
- Add hilbert based loader, when selective loading is carried out on the x,y,z coordinates: `data = osyris.load(71, scale="au", select={"x": lambda x : x > 5500. * osyris.units('au')})`
- `lmax` has been removed in arguments of `load`, instead use `data = osyris.load(71, scale="au", select={"level": lambda l : l < 12})`
- New algorithm for plane plots: hybrid method with first histogramming and then filling the blanks by finding the actual cell in the 3d domain that contains the blank image pixel.