diff --git a/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C b/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C
index 541e9831556c32be4f3b5d4a212f6f904e2f1b4b..cfa6b785ca525815ff562d1398e9dcf73940261e 100644
--- a/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C
+++ b/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C
@@ -467,7 +467,7 @@ void Foam::ThermoSurfaceFilm<CloudType>::splashInteraction
 
         // Apply correction to velocity for 2-D cases
         meshTools::constrainDirection(mesh, mesh.solutionD(), pPtr->U());
-Info<< "NEW PARTICLE: " << *pPtr << endl;
+
         // Add the new parcel
         this->owner().addParticle(pPtr);
 
diff --git a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
index e6ce90d3af1e3c520a50c2e6c54160305c7125aa..48c6762b39bb0dc35b47181888e52fee7a2592c2 100644
--- a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
+++ b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
@@ -273,15 +273,12 @@ void thermoSingleLayer::solveEnergy()
 
     updateSurfaceTemperatures();
 
-    dimensionedScalar hs0("SMALL", hs_.dimensions(), SMALL);
-
     solve
     (
         fvm::ddt(deltaRho_, hs_)
       + fvm::div(phi_, hs_)
      ==
-//      - hsSp_
-      - fvm::Sp(hsSp_/(hs_ + hs0), hs_)
+      - hsSp_
       + q(hs_)
       + radiation_->Shs()
       - fvm::SuSp(rhoSp_, hs_)