function iglo = ij2n(i,j,N)
   iglo = i + (j-1)*N;
end