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

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

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

revision 1.1 by jmc, Fri Oct 22 01:39:58 2004 UTC revision 1.2 by dimitri, Mon Apr 23 20:46:49 2007 UTC
# Line 47  C     i,j,k     :: Loop counters Line 47  C     i,j,k     :: Loop counters
47  C--   Add Vertical diffusivity contribution from KPP  C--   Add Vertical diffusivity contribution from KPP
48        IF ( kArg .EQ. 0 ) THEN        IF ( kArg .EQ. 0 ) THEN
49  C-    do all levels :  C-    do all levels :
50          DO k=1,MIN(Nr,kSize)         DO k=1,MIN(Nr,kSize)
51           DO j=jMin,jMax          DO j=jMin,jMax
52            DO i=iMin,iMax           DO i=iMin,iMax
53              KappaRT(i,j,k) = KappaRT(i,j,k)            KappaRT(i,j,k) = KPPdiffKzT(i,j,k,bi,bj)
      &                 + ( KPPdiffKzT(i,j,k,bi,bj)  
 #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))  
      &                     - diffKr(i,j,k,bi,bj) )  
 #else  
      &                     - diffKrNrT(k) )  
 #endif  
           ENDDO  
54           ENDDO           ENDDO
55          ENDDO          ENDDO
56           ENDDO
57        ELSE        ELSE
58  C-    do level k=kArg only :  C-    do level k=kArg only :
59           k = MIN(kArg,kSize)         k = MIN(kArg,kSize)
60           DO j=jMin,jMax         DO j=jMin,jMax
61            DO i=iMin,iMax          DO i=iMin,iMax
62              KappaRT(i,j,k) = KappaRT(i,j,k)           KappaRT(i,j,k) = KPPdiffKzT(i,j,kArg,bi,bj)
63       &                 + ( KPPdiffKzT(i,j,kArg,bi,bj)          ENDDO
64  #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))         ENDDO
      &                     - diffKr(i,j,kArg,bi,bj) )  
 #else  
      &                     - diffKrNrT(kArg) )  
 #endif  
           ENDDO  
          ENDDO  
65        ENDIF        ENDIF
66    
67  #endif /* ALLOW_KPP */  #endif /* ALLOW_KPP */
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22