ENH: RADModel.C: clip log input

This commit is contained in:
mattijs
2011-06-14 10:52:48 +01:00
parent a47d5f8113
commit 27ab6fcf94

View File

@ -165,7 +165,7 @@ scalar RASModel::yPlusLam(const scalar kappa, const scalar E) const
for (int i=0; i<10; i++)
{
ypl = log(E*ypl)/kappa;
ypl = log(max(E*ypl, 1))/kappa;
}
return ypl;