The heating rate should be a scalar function.
Merge pull request #176 from bangerth/fix-default-function

The heating rate should be a scalar function.