diff --git a/usadelndsoc/solver.py b/usadelndsoc/solver.py
index 04c81ea6803ac082d16dd810f20a1326fc80a56e..96b33e2058dfc0f9159968d3b03e533effebe1eb 100644
--- a/usadelndsoc/solver.py
+++ b/usadelndsoc/solver.py
@@ -1317,9 +1317,9 @@ def cpr(
             raise RuntimeError("Failed to converge")
 
         ds_cur = get_ds(Delta, state.Deltas[-1], phi, state.phis[-1])
-        state.dss.append(float(ds_cur))
+        state.dss.append(float(np.squeeze(ds_cur)))
         state.Deltas.append(Delta.copy())
-        state.phis.append(float(phi))
+        state.phis.append(float(np.squeeze(phi)))
         state.Js.append(solver.J_tot.copy())
 
         t3 = np.diag([1, 1, -1, -1])