diff --git a/dumux/implicit/mpnc/energy/mpnclocalresidualenergykinetic.hh b/dumux/implicit/mpnc/energy/mpnclocalresidualenergykinetic.hh
index cbda91eddd3ebd0ea1b8e83c0e72166ed4394ba9..bd5bac4e84beb881c36ee6f5c5b8ddc64b8737d1 100644
--- a/dumux/implicit/mpnc/energy/mpnclocalresidualenergykinetic.hh
+++ b/dumux/implicit/mpnc/energy/mpnclocalresidualenergykinetic.hh
@@ -755,7 +755,8 @@ public:
     	const Scalar deltaRho = fs.density(wPhaseIdx) - fs.density(nPhaseIdx) ;
 	const Scalar firstBracket = std::pow(g * deltaRho / gamma, 0.5);
 	const Scalar cp = FluidSystem::heatCapacity(fs, wPhaseIdx) ;
-	const Scalar TFluid		= volVars.temperature(temperatureFluidIdx);
+    // This use of Tsat is only justified if the fluid is always boiling (tsat equals boiling conditions)
+    // If a different state is to be simulated, please use the actual fluid temperature instead.
 	const Scalar Tsat = FluidSystem::vaporTemperature(fs, nPhaseIdx ) ;
 	const Scalar deltaT = TSolid - Tsat ;
 	const Scalar secondBracket = std::pow( (cp *deltaT / (0.006 * deltahv)  ) , 3.0 ) ;