Check that periodic BC's work in 1d,2d,3d; hydro and serial only
Carefully outline evolution scheme in m2sim_runge_kutta_substep function