Viscous stress (split bulk from deviator)

File viscsplit.h

viscosity:

mu = f(rho,T)
stress = mu . grad u
type name
split_viscous
parameters
   mubulk(rho,T)
   mudev(rho,T)
where both are function2d<double,double>.