--- src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C	2017-08-31 13:20:55.297970475 +0200
+++ /tmp/ThermoSurfaceFilm.C	2017-08-31 13:33:47.229357887 +0200
@@ -400,7 +400,7 @@
     }
 
     // Incident kinetic energy [J]
-    const scalar EKIn = 0.5*m*magSqr(Urel);
+    const scalar EKIn = 0.5*m*magSqr(Un);
 
     // Incident surface energy [J]
     const scalar ESigmaIn = np*sigma*p.areaS(d);