[pardoa]

Hello,

I wish to write out to time directories the pressure field and Z component from one of the patches of my domain to later calculate the hydraulic head (p/(rho*g)+Z). I managed to write both values to the screen but fail to do the same to a file. Can anyone help me with this? I'm pretty sure that this is a very simple thing but I don't find the key to make it work.

This is the coded functionObject that I have written so far:

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

label patchID = mesh().boundaryMesh().findPatchID("riverbed");
//Info << "patchID = " << patchID << endl;

//Look up "p" field at the riverbed patch

const volScalarField& p = mesh().lookupObject<volScalarField>("p");
//Info << "p= " << p << endl;

const fvPatchScalarField& p_river = p.boundaryField()[patchID];
//Info << "p_river = " << p_river << endl;

//Look up "XYZ coordinates" of the riverbed patch

const fvPatchVectorField& faceCentres_river = mesh().C().boundaryField()[patchID];
//Info << "(" << faceCentres_river << ") " << p_river << endl;

fvPatchScalarField p_river
(
IOobject
(
"p_river",
mesh().time().timeName(),
mesh().boundaryMesh(),
IOobject::NO_READ
),
p.boundaryField()[patchID]
);

p_river.write();

// Create the hydraulic head field
/*fvPatchScalarField head
(
IOobject
(
"head",
mesh().time().timeName(),
boundaryMesh(),
IOobject::NO_READ
),
p_river/(9.81*1000)
);

// Write data to the screen
Info<< "head": << head.value() << endl;

// Write values to time directories
head.write();*/
#};
}

Many thanks in advance!
Álvaro