[woody]

Dear All,

I am actually implementing a boundary condition that accesses some Field values like p, U, rho from cutting planes inside the mesh. As only the average of the cut Field is of interest, I have to create some kind of weighting factor for each processor part.

So far the Code looks like this:

   // Create the cutting Planes
   point planeLocation = planeLocations[planeI];
   plane extPlane(planeLocation,orientPatch);
   cuttingPlane cuttedPlane(extPlane,this->dimensionedInternalField().mesh());
    
   vector averageUPRho(0.0,0.0,0.0);

   scalar weigthPatch;

   // Compute the the cut of each processor
   if (cuttedPlane.cut()== 1)
   {
    
    averageUPRho.component(0)= (average(cuttedPlane.sample(Ufield)) & orientPatch); // The mean Velocity in normal Direction
    averageUPRho.component(1)= average(cuttedPlane.sample(pfield));// The mean Pressure
    averageUPRho.component(2)=average(cuttedPlane.sample(rhofield));// The mean Density 
    weigthPatch= cuttedPlane.size(); // The weighting factor
    }
    else
    {
    weigthPatch=0;
    }

Now I forward the averageUPRho and the weigthPatch to the Master and do the averaging with

rhoMean=sum(averageUPRho*weigthPatch)/sum(weigthPatch);

unfortunately I weight by the number of faces and not by their Area (magSf for patches works ... not for planes?).

How can I acces them??

Any hints?