The heating rate should be a scalar function.

The heating rate is currently defined as a vector-valued function. This makes no sense.