/[MITgcm]/MITgcm/pkg/ggl90/ggl90_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/ggl90/ggl90_calc.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.30 by jmc, Sat Feb 21 20:11:20 2015 UTC revision 1.31 by jmc, Mon Feb 23 21:20:15 2015 UTC
# Line 438  C     vertical shear term (dU/dz)^2+(dV/ Line 438  C     vertical shear term (dU/dz)^2+(dV/
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)
# Line 652  C     =============================== Line 652  C     ===============================
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
Legend:
Removed from v.1.30  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.22