[Voulet]

Hello everybody.
I would like to print out some scalarFields defined in boundary conditions files in paraview. I use OF8.

I'm coding my own boundary condition for a volScalarField called pa.

In the void Foam::myOwnBC::updateCoeffs() member function, I've defined some intermediate scalarField and vectorField, for example :

 scalarField Un = phip / patch().magSf();

Then I update my pa BC field with the usual

operate==myCustomField1 + Un + myCustomField2 + ... ;

For a debbuging purpose I would like to see in paraview the intermediate fields Un or myCustomField1.

My first attempt was to create a volScalarField in CreateFields.H :

volScalarField testField
(
   IOobject
   (
     "testField",
    runTime.timeName(),
    mesh,
    IOobject::READ_IF_PRESENT,
    IOobject::AUTO_WRITE
  ),
  mesh,
  dimensionedScalar("test", dimless, 0.0)
);

and access it in my BC file inside the ::updateCoeffs() :

 const fvPatchField<scalar> &testFieldp =
                patch().lookupPatchField<volScalarField, scalar>("testField");

and then I try to update his value (still inside ::updateCoeffs() ) :

testFieldp = Un ;

but it doesn't work :

error: passing ‘const Foam::fvPatchField<double>’ as ‘this’ argument discards qualifiers [-fpermissive]
      testFieldp = Un ;

I don't understand this error. Does anyone can advice me in order to see these fields in paraview ?