Start putting in Newtonian MHD

Solution is still wrong, because MHD wavespeeds are not written yet.