research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

atmBoundaryLayer: Correct for SP compilation

Browse source 
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1536
This commit is contained in:
Henry 2015-02-17 13:38:00 +00:00
parent 9af2d6e22e
commit 449c461a29
1 changed files with 2 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C
View file

@ -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;