[PanPeter]

Hi,
I have the following problem: for an internal boundary the following class returns the velocity at every boundary face. That works fine. But I also want to write an interim value at every boundary face.
In the code snippet below it is the quantity 'jout' that I also want to have printed (which is calculated at every boundary face).
What I tried was to define another volScalarField registry object jout, with NO_READ and AUTO_WRITE options. But this will not compile, if I refer 'mesh' instead of some function in createFields.H.

So any ideas how to extract my jout directly in updateCoeffs?
Cheers!

void Foam::explicitFOmembraneVelocityFvPatchVectorField::updateCoeffs()
{
    if(updated())
    {
        return;
    }

    {
        // get the surface normals
        tmp<vectorField> tvfnf = patch().nf();
        const vectorField& vfnf = tvfnf();

            // get the pressure field
            const fvPatchField<scalar>& ppsf = patch().lookupPatchField<volScalarField, scalar>(pName_);

            scalar rho = rho0_.value();
            forAll(ppsf, facei)
            {
                // set the velocity
                operator[](facei) = vfnf[facei] * (A_*(ppsf[facei]-ppsf[fm_[facei]])) / rho;
                jout = rho*ppsf[facei];
            }
            
}