Remember if MPI has already started.
Add unittest for advection.
Add group functions and constants to mpi header.
Add documentation and unittests for velocity and density calculations.