174 |
C compute Prandtl number (always greater than 0) |
C compute Prandtl number (always greater than 0) |
175 |
prTemp = 1. _d 0 |
prTemp = 1. _d 0 |
176 |
IF ( RiNumber .GE. 0.2 ) prTemp = 5.0 * RiNumber |
IF ( RiNumber .GE. 0.2 ) prTemp = 5.0 * RiNumber |
177 |
TKEPrandtlNumber(I,J,K) = MIN(10.,prTemp) |
TKEPrandtlNumber(I,J,K) = MIN(10.0 _d 0,prTemp) |
178 |
C mixing length |
C mixing length |
179 |
GGL90mixingLength(I,J,K) = |
GGL90mixingLength(I,J,K) = |
180 |
& SQRTTKE/SQRT( MAX(Nsquare,GGL90eps) ) |
& SQRTTKE/SQRT( MAX(Nsquare,GGL90eps) ) |