[scfd]

Hi,

I'm trying to get average value of temperature with additional statements, like if the Velocity magnitude is zero take that cell to account.

I've just 2 conditions.

1. First one is the positional, check the cell z-position if it is between P1 and P2.
2. Second one is the velocity magnitude.

If the cell satisfies the conditions, add the temperature value of that cell by Tm. Divide Tm by counter after the for loop.

	
temperatureAverage
{
    type coded;
    libs (utilityFunctionObjects);
    writeControl writeTime;
    name temp;
	
    codeWrite
    #{
	scalar P1 = 1.0;
	scalar P2 = 2.0;
	scalar Tm = 0;
	scalar TposAve = 0;
	scalar counter = 0;
        const scalarField& zPos= mesh().C().component(2);
        const volVectorField& U = mesh().lookupObject<volVectorField>("U");
	const volScalarField& T = mesh().lookupObject<volScalarField>("T");
	forAll(zPos,i)
	{
		if (zPos[i] >= P1 && zPos[i] <= P2)
		{
			if (U[i] > 0)
			{
				Tm += T[i];
				counter += 1;
			}
		}
	}
	TposAve = Tm/counter;
#};
}

ERROR:

/aveTest/system/controlDict.functions.temperatureAverage: In member function ‘virtual bool Foam::tempFunctionObject::write()’:
/aveTest/system/controlDict.functions.temperatureAverage:96:14: error: no match for ‘operator>’ (operand types are ‘const Foam::Vector<double>’ and ‘int’)

I attached the log file for further reading.

regards,
scfd