A Multiscale Central Difference Scheme Applied to Magnetohydrodynamic Simulations of Cometary Atmospheres