Meshblock averaged values

[leeyp2009]

Hi everyone,
Is there a simply method to calculate the azimuthal averaged values (e.g., velocity or density averaged over azimuth) in user-defined source term. This needs the summation over the different meshblocks. Is there a simply example to implement this. Many thanks!