sum_square_vel_tile (bi,bj) = 0. _d 0
sum_square_vel_tile (bi,bj) =
& sum_square_vel_tile (bi,bj) +
& sum_square_vel_tile, sum_square_vel, myThid )
_RL sum_square_vel_tile (nSx,nSy)