Check that the code and ICs have the same dimensionality

Does what it says on the tin. :)

Implementation of #203 following @bvandenbroucke's request. When reading ICs we check whether the ICs have the same dimensionality as the code was compiled with. We assume d=3 by default to retain compatibility with Gadget-[234].

Branch can be removed.

See merge request !221