2.1.0

• 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.