diff --git a/tests/test_solver.py b/tests/test_solver.py
index ab8ec4b0fc8100026b5becb3783e135a5afb6f49..05688129faa7cd2ed3b3de9445589f5425281772 100644
--- a/tests/test_solver.py
+++ b/tests/test_solver.py
@@ -347,3 +347,23 @@ def test_soc_analytic(T, n, alpha_soc):
         Js.append(J)
         Jans.append(Jy_an)
 
+
+def test_selfcons_h():
+    T = 0.4
+    h = 0.5
+
+    solver = Solver(nx=5, ny=2)
+    solver.mask[...] = MASK_NONE
+    solver.Lx = 5
+    solver.Ly = 5
+    solver.eta = 0
+    solver.Delta[...] = 1.0
+    solver.Omega[...] += h * np.kron(S_z, S_y)
+
+    T_c0 = np.exp(np.euler_gamma) / np.pi
+    Delta, I0, _, success = solver.self_consistency(T=T, T_c0=T_c0, workers=1)
+
+    Delta0 = BCS_Delta(T, T_c0, h=h)
+    Deltam = np.broadcast_to((Delta0 * S_0)[None, None, :, :], Delta.shape)
+
+    assert_allclose(Delta.copy(), Deltam, rtol=1e-3)