diff --git a/util/neml2_verification/README.md b/util/neml2_verification/README.md
index 20757a1a..541bf54f 100644
--- a/util/neml2_verification/README.md
+++ b/util/neml2_verification/README.md
@@ -19,7 +19,7 @@ viscoplastic.xml viscoplastic_linear_isotropic xx xx no_temperature 1e-5 1e-8
 Compare Perzyna viscoplasticity with linear isotropic hardening
 # Test input as a space-separated list of time strain stress and 
 # (optionally) temperature values
-# time strain_xx strain_yy strain_zz strain_yz strain_xz strain_xy stress_xx stress_yy stress_zz stress_yz stress_xz stress_xy (temperature)
+time strain_xx strain_yy strain_zz strain_yz strain_xz strain_xy stress_xx stress_yy stress_zz stress_yz stress_xz stress_xy (temperature)
 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
 1.0 0.001 -0.0005 -0.0005 0.0 0.0 0.0 100 0 0 0 0 0 0
 ...
diff --git a/util/neml2_verification/run_neml2_test.py b/util/neml2_verification/run_neml2_test.py
index 9d9e4731..2cf18038 100755
--- a/util/neml2_verification/run_neml2_test.py
+++ b/util/neml2_verification/run_neml2_test.py
@@ -94,7 +94,7 @@ def drive_all_strain(model, emax, erate, nsteps, T):
         f.write("# Comment line\n")
         f.write("\n")
         f.write("# The following uses the Mandel convention to store the strain and stress\n")
-        data_header = "# time strain_xx strain_yy strain_zz strain_yz strain_xz strain_xy stress_xx stress_yy stress_zz stress_yz stress_xz stress_xy"
+        data_header = "time strain_xx strain_yy strain_zz strain_yz strain_xz strain_xy stress_xx stress_yy stress_zz stress_yz stress_xz stress_xy"
         if use_T:
             data_header += " temperature"
         f.write(data_header + "\n")