diff --git a/src/USER-TALLY/compute_heat_flux_tally.cpp b/src/USER-TALLY/compute_heat_flux_tally.cpp
index 65f57b7678096e5a066f3710a54f2c4f70faedd3..b366b92be39442cb98208171e64494425eea9fe5 100644
--- a/src/USER-TALLY/compute_heat_flux_tally.cpp
+++ b/src/USER-TALLY/compute_heat_flux_tally.cpp
@@ -246,7 +246,7 @@ void ComputeHeatFluxTally::compute_vector()
       double ke_i;
 
       if (rmass) ke_i = pfactor * rmass[i];
-      else ke_i *= pfactor * mass[type[i]];
+      else ke_i = pfactor * mass[type[i]];
       ke_i *= (vi[0]*vi[0] + vi[1]*vi[1] + vi[2]*vi[2]);
       ke_i += eatom[i];