[Zhiheng Wang]

Hi my createField.H has defination of T as

const volScalarField& T = thermo.T();
I want to impliment my code on boundaryCondition correction as

forAll(T.boundaryField()[patchID],i)
{

T.boundaryField)([patchID][i] = T.boundaryField()[patchID][i]-DT;
------------
-------------------------- (some code for spoces) ........
T.write();

}

I am getting following error

error: assignment of read-only location ‘(&(&(& T)->Foam::GeometricField<Type, PatchField, GeoMesh>::boundaryField<double,.....

How to resolve this condition I want to correct the temperature T with some differet function based on spices flux. please help

[chriss85]

As of OpenFOAM 4.0 you need to use boundaryFieldRef() instead of boundaryField() if you want a non-const access to the boundaries. Same goes for internalField() and internalFieldRef() (and possibly other functions as well).

[Zeppo]

Originally Posted by chriss85

As of OpenFOAM 4.0 you need to use boundaryFieldRef() instead of boundaryField() if you want a non-const access to the boundaries. Same goes for internalField() and internalFieldRef() (and possibly other functions as well).

And also:
instead of

const volScalarField& T = thermo.T();

use

volScalarField& T = thermo.T();

otherwise it won't even allow you to call boundaryFieldRef().

[Jerryfan]

Hi wang,


As others have already mentioned above, you shouldn't use that

const

qualifier. In C++, a variable defined as const-qualifier is ready-only so that you can't change the value of the variable. So in your case, const should be removed.