setup shear_heating with fenics solution and it's working successfully
next steps
--to migrate to the gryphon implicit runge-kutta adaptive time-stepping
--to migrate to a 2d solution
--to add in dynamics with the pair of advection equations

I'm separately integrating the shear-heating from the main diffusion
equation, because the round-off error is harmful when they are added
together.