35 |
C tracer :: tracer field |
C tracer :: tracer field |
36 |
C myThid :: thread number |
C myThid :: thread number |
37 |
INTEGER bi,bj,k |
INTEGER bi,bj,k |
38 |
_RL KappaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
_RL KappaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
39 |
_RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
_RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
40 |
INTEGER myThid |
INTEGER myThid |
41 |
|
|
62 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
63 |
dfr(i,j) = |
dfr(i,j) = |
64 |
& -_rA(i,j,bi,bj) |
& -_rA(i,j,bi,bj) |
65 |
& *KappaR(i,j,k)*recip_drC(k) |
& *KappaR(i,j)*recip_drC(k) |
66 |
& *(Tracer(i,j,km1,bi,bj)-Tracer(i,j,k,bi,bj))*rkFac |
& *(Tracer(i,j,km1,bi,bj)-Tracer(i,j,k,bi,bj))*rkFac |
67 |
ENDDO |
ENDDO |
68 |
ENDDO |
ENDDO |