[maksjood]

Hi. I'm in need of summing up the values of a specific part of a scalar field. I've tried a for loop. The statement is as follows:

sum+=field[celli];

but I need the sum to be dimensioned and the above statement doesn't do that.

is there a way I can use sum() or gSum for this purpose?

[zhangyan]

just like this:

template<class Type, class GeoMesh>
dimensioned<Type> DimensionedField<Type, GeoMesh>::weightedAverage
(
    const DimensionedField<scalar, GeoMesh>& weightField
) const
{
    return
    (
        dimensioned<Type>
        (
            this->name() + ".weightedAverage(weights)",
            this->dimensions(),
            gSum(weightField*field())/gSum(weightField)
        )
    );
}

[fumiya]

You can use the volFieldValue function object:

    volFieldValue1
    {
        type            volFieldValue;
        libs            ("libfieldFunctionObjects.so");

        log             true;
        writeControl    writeTime;
        writeFields     true;

        regionType      cellZone; // Create your cellZone that the sum needs to be calculated
        name            c0;       // and set its name
        operation       sum;

        fields
        (
            p // List of fields to operate on
        );
    }

[maksjood]

Thanks everybody. I've accomplished my goal.