divergence
 

Hi,
I want to create a volumeScalar field with the value for the divergence
of a quantity S*V where S is a scalar and V a vector (velocity).
For this I do

phiS = linearInterpolate(S*V)&mesh.Sf();
forAll(mesh.cells()cellI){
netFluxVol[cellI] = 0.0;
forAll(mesh.cells()[cellI],iface)
{
netFluxVol[cellI]+=phiS[iface]
}
}
Is this right?
many thanks,
Radu