Function circ_next 

fn circ_next(l: usize, axis: usize, n: [usize; 3]) -> usize

Circular neighbour l + e_axis on the grid n.