diff --git a/src/probabilist/evaluator_quantiles.cpp b/src/probabilist/evaluator_quantiles.cpp
index a92c76a5dd4d60dd252ed0fbd59e4e973bd2cb6d..b0dc6f609068d5943a88b3be928033732b77364e 100644
--- a/src/probabilist/evaluator_quantiles.cpp
+++ b/src/probabilist/evaluator_quantiles.cpp
@@ -90,7 +90,7 @@ namespace evalhyd
             // integrate with trapezoidal rule
             crps = xt::view(
                     // xt::trapz(y, dx=1/(n+1), axis=0)
-                    xt::trapz(qs, 1./(double(n_mbr) + .1), 0),
+                    xt::trapz(qs, 1./(double(n_mbr) + 1.), 0),
                     xt::newaxis(), xt::all()
             );
         }
diff --git a/tests/test_probabilist.cpp b/tests/test_probabilist.cpp
index ec65961be1d1f03094f03a5c9c6e8e1ecb6f4025..5906cbc9847a85307a5428d32b2cccd6c6059fb2 100644
--- a/tests/test_probabilist.cpp
+++ b/tests/test_probabilist.cpp
@@ -115,7 +115,7 @@ TEST(ProbabilistTests, TestQuantiles)
 
     // Continuous ranked probability scores
     xt::xtensor<double, 3> crps =
-            {{{ 257.412129}}};
+            {{{252.956919}}};
     EXPECT_TRUE(xt::allclose(metrics[1], crps));
 }