Well, as I needed to apply free-slip conditions to a flat surface normal to the z axis with the d3q19 descriptor, I managed to copy the BounceBack class and write an ad-hoc class with the right reflection rules. However, isn’t there a more elegant way to do the job?