[Mojtaba.a]

Hello,

I am stuck in adding a source term to the turbulence kinetic energy dissipation Eqn in kEpsilon turbulence model. the source term consists of cell coordinate, for example:

S_eps (z) = C*z

I am unable to get cell coordinate. I believe I have to load the mesh, but I don't know how. I tried:

    fvMesh & mesh  =  mesh();
    volScalarField z = mesh.C().component(vector::Z);

Buy I get error in compilation:

error: no match for call to ‘(Foam::fvMesh) ()’
     fvMesh & mesh  =  mesh();

I have no idea how I can do that.
It would be great help if you could teach me how.

Thanks.