UDF to calculate average heat flux
 

Hi,
I am running a steady axisymmetric simulation (turbomachinery flow). During steady simulations, the wall heat flux oscillated due to unsteady and 3D nature of the flow (and I am forcing to solve it using steady axisymmetric approach).

I want to cumulatively average the wall heat flux w.r.t iterations ie.,
Let q_i = heat flux on a wall face at 'i' th iteration

i=1
qsum=q_1
qav=qsum/1

i=2
qsum=qsum+q_2
qav=qsum/2

i=3
qsum=qsum+q_3
qav=qsum/3

and so on

After the computations, I am planning to output the average heat flux (qav) on all the walls using
file/write-profile [wall_name]

Which UDF I need to use and please give some leads so that I can write it myself.

Thanks in advance