412 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
413 |
af(i,j) = af(i,j) |
af(i,j) = af(i,j) |
414 |
& + (maskC(i,j,k,bi,bj)-maskC(i,j,k-1,bi,bj))* |
& + (maskC(i,j,k,bi,bj)-maskC(i,j,k-1,bi,bj))* |
415 |
& rTrans(i,j)*localTijk(i,j,k) |
& rTrans(i,j)*tracer(i,j,k,bi,bj) |
416 |
|
c & rTrans(i,j)*localTijk(i,j,k) |
417 |
ENDDO |
ENDDO |
418 |
ENDDO |
ENDDO |
419 |
ELSE |
ELSE |
420 |
C- Surface "correction" term at k=1 : |
C- Surface "correction" term at k=1 : |
421 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly,sNy+Oly |
422 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
423 |
af(i,j) = rTrans(i,j)*localTijk(i,j,k) |
af(i,j) = rTrans(i,j)*tracer(i,j,k,bi,bj) |
424 |
|
c af(i,j) = rTrans(i,j)*localTijk(i,j,k) |
425 |
ENDDO |
ENDDO |
426 |
ENDDO |
ENDDO |
427 |
ENDIF |
ENDIF |