/[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.16 by gforget, Fri Aug 6 18:37:05 2010 UTC revision 1.17 by gforget, Mon Aug 9 20:34:02 2010 UTC
# Line 137  C     Initialize local fields Line 137  C     Initialize local fields
137           KappaE(I,J,K)            = 0. _d 0           KappaE(I,J,K)            = 0. _d 0
138           TKEPrandtlNumber(I,J,K)  = 1. _d 0           TKEPrandtlNumber(I,J,K)  = 1. _d 0
139           GGL90mixingLength(I,J,K) = GGL90mixingLengthMin           GGL90mixingLength(I,J,K) = GGL90mixingLengthMin
140             GGL90visctmp(I,J,K)      = 0. _d 0
141          ENDDO          ENDDO
142         ENDDO         ENDDO
143        ENDDO        ENDDO
# Line 538  C     impose minimum TKE to avoid numeri Line 539  C     impose minimum TKE to avoid numeri
539  C     end of time step  C     end of time step
540  C     ===============================  C     ===============================
541    
542        DO K=1,Nr        DO K=2,Nr
543         DO J=jMin,jMax         DO J=1,sNy
544          DO I=iMin,iMax          DO I=1,sNx
545  #ifdef ALLOW_GGL90_SMOOTH  #ifdef ALLOW_GGL90_SMOOTH
546           tmpVisc=           tmpVisc=
547       &  (       &  (
# Line 575  C     =============================== Line 576  C     ===============================
576    
577    
578    
579        DO K=1,Nr        DO K=2,Nr
580         DO J=jMin,jMax         DO J=1,sNy
581          DO I=iMin,iMax          DO I=1,sNx+1
582  #ifdef ALLOW_GGL90_SMOOTH  #ifdef ALLOW_GGL90_SMOOTH
583          tmpVisc =          tmpVisc =
584       & (       & (
# Line 609  C     =============================== Line 610  C     ===============================
610        ENDDO        ENDDO
611    
612    
613        DO K=1,Nr        DO K=2,Nr
614         DO J=jMin,jMax         DO J=1,sNy+1
615          DO I=iMin,iMax          DO I=1,sNx
616  #ifdef ALLOW_GGL90_SMOOTH  #ifdef ALLOW_GGL90_SMOOTH
617          tmpVisc =          tmpVisc =
618       & (       & (
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22