93 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly,sNy+Oly |
94 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
95 |
Rjp=(tracer(i,j,kp1,bi,bj)-tracer(i,j,k,bi,bj)) |
Rjp=(tracer(i,j,kp1,bi,bj)-tracer(i,j,k,bi,bj)) |
96 |
& *maskC(i,j,kp1,bi,bj) |
& *maskC(i,j,kp1,bi_arg,bj_arg) |
97 |
Rj=(tracer(i,j,k,bi,bj)-tracer(i,j,kM1,bi,bj)) |
Rj=(tracer(i,j,k,bi,bj)-tracer(i,j,kM1,bi,bj)) |
98 |
Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,kM2,bi,bj)) |
Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,kM2,bi,bj)) |
99 |
& *maskC(i,j,km2,bi,bj) |
& *maskC(i,j,km2,bi_arg,bj_arg) |
100 |
IF (Rj.NE.0.) THEN |
IF (Rj.NE.0.) THEN |
101 |
IF (rTrans(i,j).LT.0.) THEN |
IF (rTrans(i,j).LT.0.) THEN |
102 |
Cr=Rjm/Rj |
Cr=Rjm/Rj |