[celestial]

Hi everybody

Either this will be obvious to most of you or else it will be baffling to all.
I have no idea and so I have to give it a shot.

Here is a short excerpt from my createFields.H :

.........


Info<< "Reading field Fa, Fs \n" << endl;

volScalarField Fa
(
IOobject
(
"Fa",
runTime.timeName(),
airMesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
airMesh
);

volScalarField Fs
(
IOobject
(
"Fs",
runTime.timeName(),
solidMesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
solidMesh
);

..........

and what follows is a relevant excerpt from my C code:



...............

// patch airsolid is defined in the case dictionaries polyMesh for both meshes //airMesh and solidMesh
// the two meshes were created independently of each other. They only share the common interface airsolid in the geometry


const label id1= airMesh.boundaryMesh().findPatchID"airsolid");
const fvPatch& patch1 = airMesh.boundary()[id1];


const label id2= solidMesh.boundaryMesh().findPatchID"airsolid");
const fvPatch& patch2 = solidMesh.boundary()[id2];




scalarField dummy(patch1.size(), scalar(0) );


forAll(patch1, i)
{

dummy[i] = Fa[patch1.faceCells()[i] ] - ????;


// the above line is my problem.
// by ???? I mean some kind of interpolation of the values on the solid side:

// Fs[patch2.faceCells()[j] ] where the values of j will depend on the // current value of i
// (of course, j lies in the range 0 to patch2.size)


}

............

The question is whether any such mapping actually exist in OF.


At any rate, best regards to the community regardless of what come out of this inquiry.