/[MITgcm]/MITgcm/pkg/kpp/kpp_routines.F
ViewVC logotype

Diff of /MITgcm/pkg/kpp/kpp_routines.F

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

revision 1.42 by heimbach, Fri Feb 13 21:58:36 2009 UTC revision 1.43 by dfer, Thu May 28 20:20:59 2009 UTC
# Line 848  c     fRi, fcon             function of Line 848  c     fRi, fcon             function of
848        _RL Rig        _RL Rig
849        _RL fRi, fcon        _RL fRi, fcon
850        _RL ratio        _RL ratio
851        integer i, ki        integer i, ki, kp1
852        _RL c1, c0        _RL c1, c0
853    
854  #ifdef ALLOW_KPP_VERTICALLY_SMOOTH  #ifdef ALLOW_KPP_VERTICALLY_SMOOTH
# Line 934  c -------------------------------------- Line 934  c --------------------------------------
934  c            evaluate diffusivities and viscosity  c            evaluate diffusivities and viscosity
935  c    mixing due to internal waves, and shear and static instability  c    mixing due to internal waves, and shear and static instability
936    
937                kp1 = MIN(ki+1,Nr)
938  #ifndef EXCLUDE_KPP_SHEAR_MIX  #ifndef EXCLUDE_KPP_SHEAR_MIX
939              if ( .NOT. inAdMode ) then              if ( .NOT. inAdMode ) then
940                 diffus(i,ki,1) = viscAr + fcon * difmcon + fRi * difm0                 diffus(i,ki,1) = viscAr + fcon * difmcon + fRi * difm0
941                 diffus(i,ki,2) = diffusKzS(i,ki)+fcon*difscon+fRi*difs0                 diffus(i,ki,2) = diffusKzS(i,kp1)+fcon*difscon+fRi*difs0
942                 diffus(i,ki,3) = diffusKzT(i,ki)+fcon*diftcon+fRi*dift0                 diffus(i,ki,3) = diffusKzT(i,kp1)+fcon*diftcon+fRi*dift0
943              else              else
944                 diffus(i,ki,1) = viscAr                 diffus(i,ki,1) = viscAr
945                 diffus(i,ki,2) = diffusKzS(i,ki)                 diffus(i,ki,2) = diffusKzS(i,kp1)
946                 diffus(i,ki,3) = diffusKzT(i,ki)                 diffus(i,ki,3) = diffusKzT(i,kp1)
947              endif              endif
948  #else  #else
949              diffus(i,ki,1) = viscAr              diffus(i,ki,1) = viscAr
950              diffus(i,ki,2) = diffusKzS(i,ki)              diffus(i,ki,2) = diffusKzS(i,kp1)
951              diffus(i,ki,3) = diffusKzT(i,ki)              diffus(i,ki,3) = diffusKzT(i,kp1)
952  #endif  #endif
953    
954           end do           end do
Legend:
Removed from v.1.42  
changed lines
  Added in v.1.43
  ViewVC Help
Powered by ViewVC 1.1.22