method
NeumannZ.
prepare_rhs
Compute the distribution function at the equilibrium with the value on the border.
simulation class