From ce539bcd13c17ee630da38a3a608de9786bc2993 Mon Sep 17 00:00:00 2001 From: Henry <Henry> Date: Fri, 23 Jan 2015 14:06:02 +0000 Subject: [PATCH] atmBoundaryLayer: updated for SP compilation --- .../derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/TurbulenceModels/turbulenceModels/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C b/src/TurbulenceModels/turbulenceModels/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C index 0ec40850460..bc666d79cef 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_, 0.001))); + Ustar_ = kappa_*Uref_/(log((Zref_ + z0_)/max(z0_, scalar(0.001)))); } @@ -134,7 +134,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; -- GitLab