180 |
IF ( RiNumber .GE. 0.2 ) prTemp = 5.0 * RiNumber |
IF ( RiNumber .GE. 0.2 ) prTemp = 5.0 * RiNumber |
181 |
TKEPrandtlNumber(I,J,K) = MIN(10.0 _d 0,prTemp) |
TKEPrandtlNumber(I,J,K) = MIN(10.0 _d 0,prTemp) |
182 |
C mixing length |
C mixing length |
183 |
GGL90mixingLength(I,J,K) = |
GGL90mixingLength(I,J,K) = SQRTTWO * |
184 |
& SQRTTKE/SQRT( MAX(Nsquare,GGL90eps) ) |
& SQRTTKE/SQRT( MAX(Nsquare,GGL90eps) ) |
185 |
C impose upper bound for mixing length (total depth) |
C impose upper bound for mixing length (total depth) |
186 |
GGL90mixingLength(I,J,K) = MIN(GGL90mixingLength(I,J,K), |
GGL90mixingLength(I,J,K) = MIN(GGL90mixingLength(I,J,K), |