Integration of a value over cross-section

Hello everyone,

I have a model of molecular diffusion in a cylinder (axisymmetric geometry), and after solving (stationary solver) I get the concentration field c(r,z).

I would like to get rid of the r dependency and have access to the average concentration per cross-section

c_average(z)=2*pi*int(r*c(r,z),0,R)/(pi*R^2)

where R is the radius of the cylinder.

I managed to get the value for any fixed value of z using Results > Derived Value > Line integration, but I would like to define this variable for any value of z to be able to plot c_average(z) vs z.

Is there a way to do this?

Thanks a lot for your help!

Fabien