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

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

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

revision 1.3 by jmc, Fri Oct 22 01:42:06 2004 UTC revision 1.5 by jmc, Sun Aug 2 19:30:19 2009 UTC
# Line 35  C     KappaRx   :: vertical diffusivity Line 35  C     KappaRx   :: vertical diffusivity
35  C     myThid    :: Instance number for this innvocation of GGL90_CALC_DIFF  C     myThid    :: Instance number for this innvocation of GGL90_CALC_DIFF
36    
37        INTEGER bi,bj,iMin,iMax,jMin,jMax,kArg,kSize        INTEGER bi,bj,iMin,iMax,jMin,jMax,kArg,kSize
38        _RL KappaRx(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRx(1-Olx:sNx+Olx,1-Oly:sNy+Oly,kSize)
39        INTEGER myThid        INTEGER myThid
40    
41  #ifdef ALLOW_GGL90  #ifdef ALLOW_GGL90
# Line 51  C-    do all levels : Line 51  C-    do all levels :
51           DO j=jMin,jMax           DO j=jMin,jMax
52            DO i=iMin,iMax            DO i=iMin,iMax
53              KappaRx(i,j,k) = KappaRx(i,j,k)              KappaRx(i,j,k) = KappaRx(i,j,k)
54    #ifdef ALLOW_GGL90_SMOOTH
55         &                  +( GGL90diffKrS(i,j,k,bi,bj)
56    #else
57       &                  +( GGL90diffKr(i,j,k,bi,bj)       &                  +( GGL90diffKr(i,j,k,bi,bj)
58    #endif
59  #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))  #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
60       &                     - diffKr(i,j,k,bi,bj) )       &                     - diffKr(i,j,k,bi,bj) )
61  #else  #else
# Line 66  C-    do level k=kArg only : Line 70  C-    do level k=kArg only :
70           DO j=jMin,jMax           DO j=jMin,jMax
71            DO i=iMin,iMax            DO i=iMin,iMax
72              KappaRx(i,j,k) = KappaRx(i,j,k)              KappaRx(i,j,k) = KappaRx(i,j,k)
73    #ifdef ALLOW_GGL90_SMOOTH
74         &                  +( GGL90diffKrS(i,j,kArg,bi,bj)
75    #else
76       &                  +( GGL90diffKr(i,j,kArg,bi,bj)       &                  +( GGL90diffKr(i,j,kArg,bi,bj)
77    #endif
78  #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))  #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
79       &                     - diffKr(i,j,kArg,bi,bj) )       &                     - diffKr(i,j,kArg,bi,bj) )
80  #else  #else
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22