Move bed shear stress arrays into roughness module
Refactor VINTERP function