tmpSumTile(bi,bj)=0. _d 0
tmpSumTile(bi,bj)=tmpSumTile(bi,bj)+
CALL GLOBAL_SUM_TILE_RL( tmpSumTile, tmpSumGlo, myThid )
_RL tmpSumTile(nSx,nSy),tmpVolTile(nSx,nSy)