offSet = 0.
& offSet = -psiVel(iPsi0(bi,bj),jPsi0(bi,bj),bi,bj)
_GLOBAL_SUM_RL( offSet, myThid )
psiLoc(1) = psiLoc(1) + offSet
psiLoc(2) = psiLoc(2) + offSet
psiVel(i,j,bi,bj) = psiVel(i,j,bi,bj) + offSet
_RL offSet