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