diff --git a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
index e480d6e102..03b6de267a 100644
--- a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
+++ b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
@@ -31,6 +31,7 @@ License
 #include "zeroGradientFvPatchFields.H"
 #include "mappedFieldFvPatchField.H"
 #include "mapDistribute.H"
+#include "constants.H"
 
 // Sub-models
 #include "filmThermoModel.H"
@@ -255,19 +256,13 @@ void thermoSingleLayer::updateSubmodels()
 
 tmp<fvScalarMatrix> thermoSingleLayer::q(volScalarField& hs) const
 {
-    dimensionedScalar Tstd("Tstd", dimTemperature, 298.15);
-
-    volScalarField boundedAlpha(max(alpha_, ROOTVSMALL));
-    volScalarField htcst(htcs_->h()*boundedAlpha);
-    volScalarField htcwt(htcw_->h()*boundedAlpha);
-
-    htcst.correctBoundaryConditions();
-    htcwt.correctBoundaryConditions();
-
     return
     (
-      - fvm::Sp(htcst/Cp_, hs) - htcst*(Tstd - TPrimary_)
-      - fvm::Sp(htcwt/Cp_, hs) - htcwt*(Tstd - Tw_)
+      - fvm::Sp(htcs_->h()/Cp_, hs)
+      - htcs_->h()*(constant::standard::Tstd - TPrimary_)
+
+      - fvm::Sp(htcw_->h()/Cp_, hs)
+      - htcw_->h()*(constant::standard::Tstd - Tw_)
     );
 }
 
diff --git a/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/constant/reactingCloud1Properties b/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/constant/reactingCloud1Properties
index e2ea29a4f8..c9d32f4468 100644
--- a/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/constant/reactingCloud1Properties
+++ b/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/constant/reactingCloud1Properties
@@ -83,12 +83,12 @@ subModels
                 ((0.6 0.35 1.45) (0 0 -1))
             );
 
-            massTotal       5;
+            massTotal       10;
             parcelsPerInjector 20000;
             parcelsPerSecond 500;
             parcelBasisType mass;
-            flowRateProfile constant 1;
-            Umag            constant 1;
+            flowRateProfile constant 0.1;
+            Umag            constant 3.0;
             thetaInner      constant 0;
             thetaOuter      constant 45;