Function neighbours 

fn neighbours(
    nei: &mut [Option<usize>; 6],
    ind: usize,
    nx: usize,
    ny: usize,
    nz: usize,
)

get_neighbors: x−1, x+1, y−1, y+1, z−1, z+1, None past the array’s edge (the C’s −1).