[sreehahaha]

Dear Foamers,

I'm new to use OF and I'm developing my own solver for my problem. I have a simple cubic geometry. In it, I have to define a 'surfaceScalarField payal' and I want to give pre-calculated values to it at every time step.

It is easy to give values to the internal faces, so I'll skip that here. I give values to the boundary patches with:

// to loop over all boundary faces
forAll(payal.boundaryFieldRef(), patchID)
{
	forAll(payal.boundaryFieldRef()[patchID],facei) 
   {
   	payal.boundaryFieldRef()[patchID][facei]=<aPre-calculatedValue>;
   }
}

But I use cyclicBC, so I have to do something like below for two neighbourPatches with patchIDs 0,1:

// set values on patchID=1 same as patchID=0, because cyclicBC between patches 0,1
forAll(payal.boundaryFieldRef()[1],facei)
{
	payal.boundaryFieldRef()[1][facei]=payal.boundaryFieldRef()[0][facei];
}

and likewise for other pairs of neighbourPatches. This works fine in a serial run. But, I imagine this would be problematic in a parallel run, as both patches may be dealt with by different processes. Can somebody please suggest an elegant way which works for both serial and parallel runs? It would be very helpful.

Many thanks and have a great day/weekend!