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 |
KappaRS(i,j,k) = KappaRS(i,j,k) |
KappaRS(i,j,k) = KPPdiffKzS(i,j,k,bi,bj) |
|
& + ( KPPdiffKzS(i,j,k,bi,bj) |
|
|
#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL)) |
|
|
& - diffKr(i,j,k,bi,bj) ) |
|
|
#else |
|
|
& - diffKrNrS(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 |
KappaRS(i,j,k) = KappaRS(i,j,k) |
KappaRS(i,j,k) = KPPdiffKzS(i,j,kArg,bi,bj) |
63 |
& + ( KPPdiffKzS(i,j,kArg,bi,bj) |
ENDDO |
64 |
#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL)) |
ENDDO |
|
& - diffKr(i,j,kArg,bi,bj) ) |
|
|
#else |
|
|
& - diffKrNrS(kArg) ) |
|
|
#endif |
|
|
ENDDO |
|
|
ENDDO |
|
65 |
ENDIF |
ENDIF |
66 |
|
|
67 |
#endif /* ALLOW_KPP */ |
#endif /* ALLOW_KPP */ |