with shan chen we need to take into account
the full hydrodynamic velocity, otherwise
parasite currents at the interface border are
forcing particles in an unphysical  way.
Merge branch 'shanchen' of https://github.com/Marcello-Sega/espresso
Owen's fix for SHANCHEN