Fix a minus sign in the magnetic_pressure_gradient_force function
Replace the force-free field option, using Chandrasekhar's prescription