function pos = reduz_posicao(pos,NC)
%reduz a posicao {pos} de uma celula ao intervalo {0..NC-1}, supondo
%topologia toroidal.

pos = mod(pos,NC);