function icel = encontra_celula_pela_posicao(posi,malha)
%dada uma malha possivelmente incompleta e a posicao dessa malha encontra o
%indice {icel} da celula correspondente a essa posicao. Devolve zero se a
%celula com posicao posi nao estah na malha

NC = malha.NC;
assert((posi(1) >= 0) &&(posi(1) < NC(1)))
assert((posi(2) >= 0) &&(posi(2) < NC(2)))
icel = malha.ixc(posi(1)+1,posi(2)+1);