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) |
|
#ifdef ALLOW_GGL90_SMOOTH |
|
|
& +( GGL90diffKrS(i,j,k,bi,bj) |
|
|
#else |
|
54 |
& +( GGL90diffKr(i,j,k,bi,bj) |
& +( GGL90diffKr(i,j,k,bi,bj) |
|
#endif |
|
|
#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL)) |
|
|
& - diffKr(i,j,k,bi,bj) ) |
|
|
#else |
|
55 |
& - diffKrNrT(k) ) |
& - diffKrNrT(k) ) |
|
#endif |
|
56 |
ENDDO |
ENDDO |
57 |
ENDDO |
ENDDO |
58 |
ENDDO |
ENDDO |
62 |
DO j=jMin,jMax |
DO j=jMin,jMax |
63 |
DO i=iMin,iMax |
DO i=iMin,iMax |
64 |
KappaRx(i,j,k) = KappaRx(i,j,k) |
KappaRx(i,j,k) = KappaRx(i,j,k) |
|
#ifdef ALLOW_GGL90_SMOOTH |
|
|
& +( GGL90diffKrS(i,j,kArg,bi,bj) |
|
|
#else |
|
65 |
& +( GGL90diffKr(i,j,kArg,bi,bj) |
& +( GGL90diffKr(i,j,kArg,bi,bj) |
|
#endif |
|
|
#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL)) |
|
|
& - diffKr(i,j,kArg,bi,bj) ) |
|
|
#else |
|
66 |
& - diffKrNrT(kArg) ) |
& - diffKrNrT(kArg) ) |
|
#endif |
|
67 |
ENDDO |
ENDDO |
68 |
ENDDO |
ENDDO |
69 |
ENDIF |
ENDIF |