method
NeumannZ.
set_iload
Compute the indices that are needed (symmertic velocities and space indices).