[pardoa]

Hi everyone,

I have defined a new function object to calculate the hydraulic head from the values of water pressure at the bottom of my problem domain. Since the area that I'm simulating is too big, I have split it into several pieces so as to run it in parallel using several cores.

The computation of the heads works very well at each processor but once I reconstruct the domain, the values of head are not merged as it happens to the rest of calculated fields.

Can anyone give me a hand on this? Thanks!

Alvaro

[pardoa]

This is the object function I wrote:

head
{
type coded;
functionObjectLibs ("libutilityFunctionObjects.so");
name head;
writeControl writeTime;
codeExecute
#{
// Find the riverbed patch ID

label patchID = mesh().boundaryMesh().findPatchID("riverbed");

// Select pressure field at the riverbed patch

const volScalarField& p = mesh().lookupObject<volScalarField>("p");
const fvPatchScalarField& p_river = p.boundaryField()[patchID];

// Select face center coordinates of the riverbed patch

const fvPatchVectorField& faceCentres_river = mesh().C().boundaryField()[patchID];

// Create the hydraulic head field
IOField<scalar> head
(
IOobject
(
"head",
mesh().time().timeName(),
mesh(),
IOobject::NO_READ
),
(p_river/(-9.81*1000))+(faceCentres_river & vector(0,0,1))
);

// Write out only at output times selected for the simulation
if(mesh().time().outputTime())
{
head.write();
}
#};
}

As I said in my previous post, the resulting values are recorded into a file in the time sub-directories included in each of the processor directories.

I have checked the reconstructPar code and apparently it is not possible to reconstruct patchScalarFields but only volumes. Does anyone face this issue before? Many thanks in advance.

Álvaro