An integrated model for the ionization, radiation, and advection of impurities in the extended-magnetohydrodynamic code M3D-C1 is described. This implementation makes use of the KPRAD model, which calculates bremsstrahlung radiation and impurity ionization, recombination, and radiation rates using a model in which the density of each charge state is advanced separately. The integrated model presented here allows the independent evolution of electron and ion temperatures, which is necessary to accurately model cases where the electron temperature drops more quickly than the electron–ion thermal equilibration time. This model is used to simulate the disruption of a model NSTX discharge caused by the introduction of argon impurities, using physically realistic resistivity. Despite well-mixed impurities, contraction of the current channel is found to lead to magnetohydrodynamic instabilities that result in stochastization of the magnetic field, a fast thermal quench, and localized parallel electric fields that can exceed the axisymmetric values by a factor of five for brief periods.