diff --git a/examples/MHDTests/DecayingTurbulence_MHD/plotSpectrum.py b/examples/MHDTests/DecayingTurbulence_MHD/plotSpectrum.py
index 5b4666411b279cc799505fbe3db3e08fa0754af2..9d0c5e23e9c8ec2b29a9daaac1de750db67b8420 100644
--- a/examples/MHDTests/DecayingTurbulence_MHD/plotSpectrum.py
+++ b/examples/MHDTests/DecayingTurbulence_MHD/plotSpectrum.py
@@ -47,7 +47,7 @@ divB = data.gas.magnetic_divergences
 minh = np.min(h.value)
 
 # Renormalize quantities
-mu0 = 1.25663706127e-1 * unyt.g * unyt.cm / (unyt.s ** 2 * unyt.A ** 2)
+mu0 = 1.25663706127e-1 * unyt.g * unyt.cm / (unyt.s ** 2 * unyt.statA ** 2)
 v = v * (rho[:, None]/2)**0.5
 B = B / np.sqrt(2*mu0)