From 21a9076441f0fc81968e7be6e8121eb1d0f5350c Mon Sep 17 00:00:00 2001
From: sergio <s.ferraris@opencfd.co.uk>
Date: Fri, 22 Apr 2022 10:38:43 -0700
Subject: [PATCH] ENH: Correcting mass rate output

---
 .../humidityTemperatureCoupledMixedFvPatchScalarField.C         | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.C b/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.C
index d64f100256..1baf7faece 100644
--- a/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.C
+++ b/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.C
@@ -747,7 +747,7 @@ void Foam::humidityTemperatureCoupledMixedFvPatchScalarField::updateCoeffs()
 
     if (debug && fluid_)
     {
-        scalar Qdm = gSum(dm);
+        scalar Qdm = gSum(dm*magSf);
         scalar QMass = gSum(mass_);
         scalar Qt = gSum(myKDelta_*(Tp - Tin)*magSf);
         scalar QtSolid = gSum(KDeltaNbr*(Tp - nbrIntFld)*magSf);