need to add max_atom_lj or already included ? to check
correct the normalisation of rdf with different types of solutes
check if resid 3d / 2d with sdamm is possible
At present calculation not parallelised
Could add option to skip certain number of trajectory frames
Virial