Implementation of externally-applied magnetic fields to a non-ideal MHD solver for MPD thruster simulation