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

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

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

revision 1.7 by cnh, Sun Feb 4 14:38:50 2001 UTC revision 1.8 by heimbach, Mon Aug 13 18:05:26 2001 UTC
# Line 6  C $Name$ Line 6  C $Name$
6        SUBROUTINE KPP_CALC_DIFF(        SUBROUTINE KPP_CALC_DIFF(
7       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
8       I        maskUp,       I        maskUp,
9       U        KappaRT,KappaRS,KappaRU,KappaRV,       U        KappaRT,KappaRS,
10       I        myThid)       I        myThid)
11    
12  C     /==========================================================\  C     /==========================================================\
# Line 19  C     == GLobal variables == Line 19  C     == GLobal variables ==
19  #include "SIZE.h"  #include "SIZE.h"
20  #include "EEPARAMS.h"  #include "EEPARAMS.h"
21  #include "PARAMS.h"  #include "PARAMS.h"
22    #include "DYNVARS.h"
23  #include "GRID.h"  #include "GRID.h"
24    
25  C     == Routine arguments ==  C     == Routine arguments ==
# Line 30  C Line 31  C
31        _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
32        _RL KappaRT(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRT(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
33        _RL KappaRS(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRS(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
       _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)  
       _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)  
34        INTEGER myThid        INTEGER myThid
35    
36  #ifdef ALLOW_KPP  #ifdef ALLOW_KPP
# Line 44  C     I, J, K - Loop counters Line 43  C     I, J, K - Loop counters
43        DO j=jMin,jMax        DO j=jMin,jMax
44         DO i=iMin,iMax         DO i=iMin,iMax
45          KappaRT(i,j,k) = KappaRT(i,j,k)+maskUp(i,j)*          KappaRT(i,j,k) = KappaRT(i,j,k)+maskUp(i,j)*
46       &                   (KPPdiffKzT(i,j,k,bi,bj)-diffKrT)       &                   ( KPPdiffKzT(i,j,k,bi,bj)
47    #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
48         &                     - diffKr(i,j,k,bi,bj) )
49    #else
50         &                     - diffKrT )
51    #endif
52         ENDDO         ENDDO
53        ENDDO        ENDDO
54    
55        DO j=jMin,jMax        DO j=jMin,jMax
56         DO i=iMin,iMax         DO i=iMin,iMax
57          KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)*          KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)*
58       &                   (KPPdiffKzS(i,j,k,bi,bj)-diffKrS)       &                   ( KPPdiffKzS(i,j,k,bi,bj)
59         ENDDO  #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
60        ENDDO       &                     - diffKr(i,j,k,bi,bj) )
61    #else
62        DO j=jMin,jMax       &                     - diffKrS )
63         DO i=iMin,iMax  #endif
         KappaRU(i,j,k) = max ( KappaRU(i,j,k),  
      &         KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) *  
      &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i-1,j,k,bi,bj)) )  
        ENDDO  
       ENDDO  
   
       DO j=jMin,jMax  
        DO i=iMin,iMax  
         KappaRV(i,j,k) = max ( KappaRV(i,j,k),  
      &         KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) *  
      &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i,j-1,k,bi,bj)) )  
64         ENDDO         ENDDO
65        ENDDO        ENDDO
66    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22