diff --git a/Resample/Coherence/SumDWBA.cpp b/Resample/Coherence/SumDWBA.cpp
index c98290044c793f69c2bffdca09c0f2cfc52f96a8..58e4089bee24fe7c76de78b3884e04eb87da3621 100644
--- a/Resample/Coherence/SumDWBA.cpp
+++ b/Resample/Coherence/SumDWBA.cpp
@@ -87,10 +87,14 @@ complex_t SumDWBA::coherentFF(const DiffuseElement& ele) const
     // Note that the order of multiplication matters:
     // If a prefactor is 0, then theFF() won't be called.
 
-    const complex_t term_S = T_in * T_out * m_ff->theFF(q_TT);
+    std::cout << "DEBUG TERM q=" << q_TT.getQ() << std::endl;
+    complex_t ff = m_ff->theFF(q_TT);
+    std::cout << "DEBUG TERM ->" << ff << std::endl << std::endl;
+    const complex_t term_S = T_in * T_out * ff;
     const complex_t term_RS = R_in * T_out * m_ff->theFF(q_RT);
     const complex_t term_SR = T_in * R_out * m_ff->theFF(q_TR);
     const complex_t term_RSR = R_in * R_out * m_ff->theFF(q_RR);
+    std::cout << "DEBUG TERM done with remaining terms" << ff << std::endl << std::endl;
 
     return term_S + term_RS + term_SR + term_RSR;
 }