Visualizing mean average value for each body

[vencels]

Hi,

I want to calculate the mean/average value for a field (e.g., joule loss) within each body separately and save it as a field.

I intend to use it to validate the homogenization model.

Potential solutions I identified:
1) Using SaveScalars with many different masks, applying mean operation on them.
2) Writing a solver that just does that.

Are there easy and clever ways to do it?