Add a HeatFlux model that wraps existing models

`model = HeatFlux.wrap(model)` will create a new model with the
same outputs, and (multiple) extra outputs for the heat_flux,
one for each energy variant.

Author: GardevoirX