MPI parallel solver for a set of 3D Magnetohydrodynamics equations.  Currently solves Boussinesq equations, but this will be extended to be more flexible.