& ptracer_half , ptracer_full
_RL ptracer_half(nSx,nSy), ptracer_full(nSx,nSy)
ptracer_full(bi,bj) = 0.
ptracer_full(bi,bj) = ptracer_full(bi,bj)+deltaTclock
& ptracer_full, 1, bi, bj, myThid )