Mave backwards a bit to test second order method on simple problems

There was a major typo in the plm_gradient function switching the
xs (coordinates) and the fs (solution values). It's unclear why this
appeared to give good results with the wind problem in spherical
coordinates.
Show that a cylindrical MHD explosion works with second order