/[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.2 by mlosch, Mon Sep 27 08:02:04 2004 UTC revision 1.5 by ce107, Fri Aug 5 23:46:09 2005 UTC
# Line 110  CEOP Line 110  CEOP
110        jMax = sNy+OLy-1        jMax = sNy+OLy-1
111    
112  C     set separate time step (should be deltaTtracer)  C     set separate time step (should be deltaTtracer)
113        deltaTggl90 = deltaTtracer        deltaTggl90 = dTtracerLev(1)
114  C      C    
115        kSurf = 1        kSurf = 1
116  C     implicit timestepping weights for dissipation  C     implicit timestepping weights for dissipation
# 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) )
# Line 365  C     Eq. (11), (18) and (21) Line 365  C     Eq. (11), (18) and (21)
365           KappaH = KappaM/TKEPrandtlNumber(I,J,K)           KappaH = KappaM/TKEPrandtlNumber(I,J,K)
366  C     Set a minium (= background) value  C     Set a minium (= background) value
367           KappaM = MAX(KappaM,viscAr)           KappaM = MAX(KappaM,viscAr)
368           KappaH = MAX(KappaH,diffKrT)           KappaH = MAX(KappaH,diffKrNrT(k))
369  C     Set a maximum and mask land point  C     Set a maximum and mask land point
370           GGL90viscAr(I,J,K,bi,bj) = MIN(KappaM,GGL90viscMax)           GGL90viscAr(I,J,K,bi,bj) = MIN(KappaM,GGL90viscMax)
371       &        * maskC(I,J,K,bi,bj)       &        * maskC(I,J,K,bi,bj)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22