diff --git a/cdn/qpw_to_nmt.f90 b/cdn/qpw_to_nmt.f90
index a008e09b6a29706fd3f3f684003793912589cbf6..382c825855cb797aeea396b748d0acfe36f4697f 100644
--- a/cdn/qpw_to_nmt.f90
+++ b/cdn/qpw_to_nmt.f90
@@ -133,6 +133,8 @@ CONTAINS
        n1 = 1
        DO in = 1 , nrm
           n2 = mr(in)
+          bsl_r = 0.0
+          bsl_i = 0.0
           DO j = 1,atoms%jri(n1)
              cprr2 = cp*atoms%rmsh(j,n1)*atoms%rmsh(j,n1)
              gr = stars%sk3(k)*atoms%rmsh(j,n1)