function CD_c = centro_da_celula(pos,CD,LD,NC)
%determina o centro da celula de posicao pos de uma malha de celulas
%regular completa com NC(1) x NC(2) celulas no dominio CD, LD.

LD_c = LD ./ NC;
CD_c = CD - 0.5*LD + (pos .* LD_c) + 0.5*LD_c;