[AbhishekTAMU]

Hi,

I am exactly trying to do same thing (printing out cellSkewness values when the simulation is running and assigning these values to a volScalarField so that I can use it as refinement parameter. I am using openfoam-6 version.

But I end up getting an error. This is what I write in my solver (which is .C file)-

Here. the 'createCellQualityFields.H' is the same file Vassili Kitsios has posted in this thread.

//code lines

# include "createCellQualityFields.H"
cellQuality qualityCheck(mesh);
cellSkewness.ref() = qualityCheck.skewness();
cellNonOrthogonality.ref() = qualityCheck.nonOrthogonality();

error that I get-

error: no operator "=" matches these operands
operand types are: Foam:: DimensionedField<Foam::scalar, Foam::volMesh> = Foam::tmp<Foam::Field<Foam::scalar>>
cellSkewness.ref() = qualityCheck.skewness() ;
^

error: no operator "=" matches these operands
operand types are: Foam:: DimensionedField<Foam::scalar, Foam::volMesh> = Foam::tmp<Foam::Field<Foam::scalar>>
cellNonOrthogonality.ref() = qualityCheck.nonOrthogonality() ;

any help would be highly appreciated.
Thanks.