[saba_saeb]

Hi,

I need to take the average of a field on a boundary during the simulation. I am doing that using the normal approach:

scalar tmp= average(frictionvelocity);

where frictionvelocity is a scalar field. Everything is running like a charm in serial, but the simulation would stop if I run it in parallel. I tried to tackle the problem using this approach:

scalar tmp= average(frictionVelocity);
scalar tauw = returnReduce(tmp, sumOp<scalar>());

But the output is not correct since the calculated averages on every processors are being added together which make tauw that seems incorrect. Any idea would be appreciated.

Tnx,
Saba