438 |
|
|
439 |
C viscosity and diffusivity |
C viscosity and diffusivity |
440 |
KappaM = GGL90ck*GGL90mixingLength(i,j,k)*SQRTTKE(i,j,k) |
KappaM = GGL90ck*GGL90mixingLength(i,j,k)*SQRTTKE(i,j,k) |
441 |
GGL90visctmp(i,j,k) = MAX(KappaM,diffKrNrT(k)) |
GGL90visctmp(i,j,k) = MAX(KappaM,diffKrNrS(k)) |
442 |
& * maskC(i,j,k,bi,bj) |
& * maskC(i,j,k,bi,bj) |
443 |
C note: storing GGL90visctmp like this, and using it later to compute |
C note: storing GGL90visctmp like this, and using it later to compute |
444 |
C GGL9rdiffKr etc. is robust in case of smoothing (e.g. see OPA) |
C GGL9rdiffKr etc. is robust in case of smoothing (e.g. see OPA) |
652 |
tmpVisc = GGL90visctmp(i,j,k) |
tmpVisc = GGL90visctmp(i,j,k) |
653 |
#endif |
#endif |
654 |
tmpVisc = MIN(tmpVisc/TKEPrandtlNumber(i,j,k),GGL90diffMax) |
tmpVisc = MIN(tmpVisc/TKEPrandtlNumber(i,j,k),GGL90diffMax) |
655 |
GGL90diffKr(i,j,k,bi,bj)= MAX( tmpVisc , diffKrNrT(k) ) |
GGL90diffKr(i,j,k,bi,bj)= MAX( tmpVisc , diffKrNrS(k) ) |
656 |
ENDDO |
ENDDO |
657 |
ENDDO |
ENDDO |
658 |
ENDDO |
ENDDO |