/[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.4 by jmc, Sat Dec 4 00:16:49 2004 UTC revision 1.5 by ce107, Fri Aug 5 23:46:09 2005 UTC
# Line 174  C     vertical shear term (dU/dz)^2+(dV/ Line 174  C     vertical shear term (dU/dz)^2+(dV/
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) )
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22