36 |
#include "GMREDI.h" |
#include "GMREDI.h" |
37 |
|
|
38 |
C == Local variables == |
C == Local variables == |
39 |
C I, J, K - Loop counters |
C I, J - Loop counters |
40 |
INTEGER i,j |
INTEGER i,j |
41 |
|
|
42 |
DO j=jMin,jMax |
DO j=jMin,jMax |
43 |
DO i=iMin,iMax |
DO i=iMin,iMax |
44 |
KappaRT(i,j,k) = KappaRT(i,j,k)+maskUp(i,j)* |
KappaRT(i,j,k) = KappaRT(i,j,k)+maskUp(i,j)*Kwz(i,j,k,bi,bj) |
|
& ( GM_background_K |
|
|
#ifdef GM_VISBECK_VARIABLE_K |
|
|
& + VisbeckK(i,j,bi,bj) |
|
|
#endif |
|
|
& )*Kwz(i,j,k,bi,bj) |
|
45 |
ENDDO |
ENDDO |
46 |
ENDDO |
ENDDO |
47 |
|
|
48 |
DO j=jMin,jMax |
DO j=jMin,jMax |
49 |
DO i=iMin,iMax |
DO i=iMin,iMax |
50 |
KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)* |
KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)*Kwz(i,j,k,bi,bj) |
|
& ( GM_background_K |
|
|
#ifdef GM_VISBECK_VARIABLE_K |
|
|
& + VisbeckK(i,j,bi,bj) |
|
|
#endif |
|
|
& )*Kwz(i,j,k,bi,bj) |
|
51 |
ENDDO |
ENDDO |
52 |
ENDDO |
ENDDO |
53 |
|
|