Copy IP viscosity scheme from IPCS branch

Not sure this is working correctly yet.