_RL maskm1, maskp1
maskm1 = 1.
IF (k.EQ.1) maskm1 = 0.
& ( maskm1*_maskS(i,j, k ,bi,bj)*tauyEddy(i,j, k ,bi,bj)
& ( maskm1*_maskW(i,j, k ,bi,bj)*tauxEddy(i,j, k ,bi,bj)