Viscous stress (split bulk from deviator)

viscosity:
mu = f(rho,T)
stress = mu . grad u
type name
split_viscous
parameters
   mubulk(rho,T)
   mudev(rho,T)
where both are Realfunction2ds.