diff --git a/src/parameters.f90 b/src/parameters.f90 index 6dab424191e0786d80c413d9cc767f74fb568c35..8c9bcc49bd98e1ec13379611a020c316d936f2c2 100644 --- a/src/parameters.f90 +++ b/src/parameters.f90 @@ -255,7 +255,11 @@ subroutine parameter(input_i3d) dy2 = dy * dy dz2 = dz * dz - xnu=one/re + if (re.gt.0.001) then + xnu=one/re + else + xnu=zero + endif !! Constant pressure gradient, re = Re_tau -> use to compute Re_centerline if (cpg) then re_cent = (re/0.116_mytype)**(1.0_mytype/0.88_mytype)