[Aaron_L]

Hi, foamers,
I want to get mean temperature field from T field if z coordinate is equal , this is my test code,

//get T mean along z axis
    forAll(mesh.C(),celli) 
    {
        scalar zi = mesh.C()[celli].component(vector::Z);
        scalar& Tmeani = T_mean[celli];
        scalar& Ti = T[celli];
        scalar Vi = mesh.V()[celli];
        scalar Vtotal = Vi;
        scalar Ttotal = Ti*Vi;
        forAll(mesh.C(),cellj)
        {
            scalar zj = mesh.C()[cellj].component(vector::Z);
            if(zi == zj)
            {
                scalar& Tj = T[cellj];
                scalar Vj = mesh.V()[cellj];
                Vtotal = Vtotal + Vj;
                Ttotal = Ttotal + Tj*Vj;
            }
        }
        Tmeani = Ttotal/Vtotal ; 
    }

but it was run too slow because I use for two nested loop, how can I handle for it for one loop?
any suggestions? or someone can give me some brief solutions?

[Aaron_L]

attachment picture is T_mean equation, I want to get T_mean along z direction using xy plane average

[miha23]

Dear Aaron!

Have you made any process with it ever since?

[Aaron_L]

no, it is not very important, so I didn't complete it.