fixed: wrong calculation of H1

8883e2b8 · Andrew Heather · f97e559b · 8883e2b8
Commit 8883e2b8 authored 15 years ago by Andrew Heather
--- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C
@@ -110,7 +110,7 @@ void Foam::ThermoParcel<ParcelType>::calc
    scalar T1 = calcHeatTransfer(td, dt, cellI, d0, U0, rho0, T0, cp0, Sh);

    // Calculate new enthalpy state
-    scalar H1 = cp0 - T1;
+    scalar H1 = cp0*T1;


    // Motion