/[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.8 by jmc, Mon Aug 11 22:28:06 2008 UTC revision 1.9 by dfer, Mon Sep 29 13:47:23 2008 UTC
# Line 180  C     compute Prandtl number (always gre Line 180  C     compute Prandtl number (always gre
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),
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22