[cfd@kgp]

Dear Foam experts',

here is a coded function to calculate a derived variable using a existiong variable IO object variable sMS1:alpha1. The following quoted text was added in the controlDict,

TIntegral
{
functionObjectLibs ("libutilityFunctionObjects.so");
type coded;
redirectType integral;
outputControl timeStep;
code
#{
const volScalarField& fl = mesh().lookupObject<volScalarField>("sMS1:alpha1") ;
// volScalarField magU(mag(U));

scalar volIntegral = 0;
scalar volIntegralV = 0;

forAll (fl, cellI)
{
volIntegral += fl[cellI]*mesh().V()[cellI];
volIntegralV += 1.0*mesh().V()[cellI];
}
volIntegral=volIntegral/volIntegralV;
Info<<"ANSWER: " << volIntegral << endl;
//foo.write();

#};
}

In serial running of the solver it performs and does the job. But when executed in parallel the output is always zero.

Somebody kindly help.

Thanks and regards,
cfd@kgp