[Harris_Daniel]

Hi,
I am using OpenFOAM-v1612+ and modified InterFoam solver by introducing the Fourier's Heat Equation from laplacian Foam. Now I need to calculate the snGrad for all the faces of all the cells by looping over the cells and faces. I wrote a code but there seems to be a problem with it.

I wrote the following code:

forAll(C.internalField(), cellI)
    {
        const labelList& faces = mesh.cells()[cellI];
        double Tin = Tf.ref()[cellI];
        scalarField Ti = Tin.sndGrad();
        forAll(faces,I)
        {
           const scalar Tsngrad(T.ref()[faces[I]].snGrad());
        }
    }

And got the following a error message:

error: request for member ‘sndGrad’ in ‘Tin’, which is of non-class type ‘double’
scalarField Ti = Tin.sndGrad();

error: request for member ‘snGrad’ in ‘(& T.Foam::GeometricField<Type, PatchField, GeoMesh>::ref<double, Foam::fvPatchField, Foam::volMesh>())->Foam:imensionedField<double, Foam::volMesh>::<anonymous>.Foam::Field<double>::< anonymous>.Foam::List<double>::<anonymous>.Foam::U List<T>:perator[]<double>(((int)(& faces)->Foam::List<int>::<anonymous>.Foam::UList<T>:per ator[]<int>(CellI)))’, which is of non-class type ‘double’
const scalar Tsngrad(T.ref()[faces[CellI]].snGrad());

Please help me clear this error. I have had a look through the snGradScheme.C file but in vain.

What does sndGrad scheme return? It says that it return the sndGrad for a given cell field. Can anyone elaborate on this.

More importantly what should be the inputs to a snGrad and sndGrad scheme(I am asking about the data types).

Regards,
S. Harris Daniel.