diff --git a/src/finiteVolume/fvMesh/fvMeshGeometry.C b/src/finiteVolume/fvMesh/fvMeshGeometry.C index 638d244530..7844ebade1 100644 --- a/src/finiteVolume/fvMesh/fvMeshGeometry.C +++ b/src/finiteVolume/fvMesh/fvMeshGeometry.C @@ -31,6 +31,7 @@ License #include "slicedVolFields.H" #include "slicedSurfaceFields.H" #include "SubField.H" +#include "cyclicFvPatchFields.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -146,6 +147,27 @@ void fvMesh::makeC() const cellCentres(), faceCentres() ); + + + // Need to correct for cyclics transformation since absolute quantity. + // Ok on processor patches since hold opposite cell centre (no + // transformation) + slicedVolVectorField& C = *CPtr_; + + forAll(C.boundaryField(), patchi) + { + if (isA(C.boundaryField()[patchi])) + { + // Note: cyclic is not slice but proper field + C.boundaryField()[patchi] == static_cast + ( + static_cast&> + ( + boundary_[patchi].patchSlice(faceCentres()) + ) + ); + } + } }