From 449c461a299d61ef1afc82976ee6b959b88c8d5e Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Tue, 17 Feb 2015 13:38:00 +0000
Subject: [PATCH] atmBoundaryLayer: Correct for SP compilation Resolves
 bug-report http://www.openfoam.org/mantisbt/view.php?id=1536

---
 .../derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C  | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C
index d01a0079..eeb30de7 100644
--- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C
+++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C
@@ -73,7 +73,7 @@ atmBoundaryLayer::atmBoundaryLayer(const vectorField& p, const dictionary& dict)
     flowDir_ /= mag(flowDir_);
     zDir_ /= mag(zDir_);
 
-    Ustar_ = kappa_*Uref_/(log((Zref_ + z0_)/max(z0_, 0.001)));
+    Ustar_ = kappa_*Uref_/(log((Zref_ + z0_)/max(z0_, scalar(0.001))));
 }
 
 
@@ -136,7 +136,7 @@ tmp<vectorField> atmBoundaryLayer::U(const vectorField& p) const
     scalarField Un
     (
         (Ustar_/kappa_)
-       *log(((zDir_ & p) - zGround_ + z0_)/max(z0_, 0.001))
+        *log(((zDir_ & p) - zGround_ + z0_)/max(z0_, scalar(0.001)))
     );
 
     return flowDir_*Un;