diff --git a/Project.toml b/Project.toml
index 667edae..269fd25 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "ITensorMPS"
 uuid = "0d1a4710-d33b-49a5-8f18-73bdf49b47e2"
 authors = ["Matthew Fishman <mfishman@flatironinstitute.org>", "Miles Stoudenmire <mstoudenmire@flatironinstitute.org>"]
-version = "0.3.23"
+version = "0.3.24"
 
 [deps]
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
diff --git a/examples/finite_temperature/metts.jl b/examples/finite_temperature/metts.jl
index 56a7f37..43ff74d 100644
--- a/examples/finite_temperature/metts.jl
+++ b/examples/finite_temperature/metts.jl
@@ -50,8 +50,9 @@ function main(; N = 10, cutoff = 1.0e-8, δτ = 0.1, beta = 2.0, NMETTS = 3000,
     # Make y-rotation gates to use in METTS collapses
     Ry_gates = ops([("Ry", n, (θ = π / 2,)) for n in 1:N], s)
 
-    # Arbitrary initial state
-    psi = random_mps(s)
+    # Arbitrary initial Z-product state
+    init_state = rand(["Z+", "Z-"], N)
+    psi = MPS(s, init_state)
 
     # Make H for measuring the energy
     terms = OpSum()
@@ -108,6 +109,7 @@ function main(; N = 10, cutoff = 1.0e-8, δτ = 0.1, beta = 2.0, NMETTS = 3000,
             samp = sample!(psi)
             new_state = [samp[j] == 1 ? "Z+" : "Z-" for j in 1:N]
         end
+        @printf("  Next Product State = %s\n", new_state)
         psi = MPS(s, new_state)
     end