diff --git a/acc/mhd_solver/stencil_process.sps b/acc/mhd_solver/stencil_process.sps
index 7074a4b..440ffe3 100644
--- a/acc/mhd_solver/stencil_process.sps
+++ b/acc/mhd_solver/stencil_process.sps
@@ -118,6 +118,7 @@ accretion_profile(int3 globalVertexIdx, in Scalar lnrho){
      
     // multiplying the truelove density by a wave function to avoid step-function like accretion profile.   
     const Scalar weight = exp(-(accretion_distance/profile_range));    
+    const Scalar sink_mass = DCONST_REAL(AC_M_sink);
 //    const Scalar rate = truelove_density(lnrho);
     const Scalar B = Scalar(3.0);
     const Scalar k = Scalar(2.0);