As suggested by @Jeppembo, would be good to create a feature in the mesh intersection in "InterfaceSurface.map_power", where we output the total power deposited on the surfaces intersected with. This would be a good way to check for bugs in the implementation and make it easier to check for power balance.
The normalised power profile that is used for power sampling could be moved to a hidden attribute.