Partial implementation viscosity/pressure.
wes, fixed r vec3, rdist. Viscosity not working.
Merges with andy_part2