diff --git a/stdpopsim/catalog/HomSap/dfes.py b/stdpopsim/catalog/HomSap/dfes.py
index ec690d440..7c1a6ce6d 100644
--- a/stdpopsim/catalog/HomSap/dfes.py
+++ b/stdpopsim/catalog/HomSap/dfes.py
@@ -38,7 +38,6 @@ def _KimDFE():
     # in how fitness is defined
     # (1+s for homozygote in SLiM versus 1+2s in dadi)
     gamma_mean = (-gamma_shape * gamma_scale * 2) / (2 * Na)  # expected value
-    print(gamma_mean)
     h = 0.5  # dominance coefficient
     negative = stdpopsim.MutationType(
         dominance_coeff=h,