Added material model based on rheology in Morency and Doin, 2004

you can also make 'z' const if you write this as
  const double z = (dim==2 ? ... : ...);