Generalize implementation of Beavers-Joseph(-Saffmann) slip condition

So far ths BJS condition only considers du/dy. We should extend this to (du/dy + dv/dx)