A verification benchmark has been carried out between the M3D-C1 and NIMROD extendedmagnetohydrodynamic codes for simulations of impurity-induced disruption mitigation. Disruptions are a significant concern for future tokamaks and high-fidelity simulations are required in order to ensure the success of disruption mitigation techniques (e.g. shattered-pellet injection) in large-scale fusion reactors. Both magnetohydrodynamic (MHD) codes have been coupled to the Killer Pellet RADiation code for impurity dynamics. The codes show excellent agreement in four axisymmetric, nonlinear simulations, particularly during the thermal quench. This agreement is seen in the time histories of global plasma quantities such as thermal energy, radiated power, and total number of electrons, as well as 2D contours of temperature and current density. The simulations predict that, given the same number of atoms injected, argon quenches the plasma two-to-three times as fast as neon. Furthermore, the inclusion of temperaturedependent Spitzer resistivity causes the current to diffuse and to decay, inducing axisymmetric MHD instabilities that result in a current quench. This work represents an important verification of the coupled impurity and MHD models implemented in M3D-C1 and NIMROD, giving greater confidence in the ability of both codes to perform more sophisticated disruption mitigation simulations.