From 70f1bad0887e7fd27f30eff23d1977e150fe20e3 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs@hunt.opencfd.co.uk>
Date: Wed, 4 Mar 2009 12:35:15 +0000
Subject: [PATCH] mesh.C on cyclics

---
 src/finiteVolume/fvMesh/fvMeshGeometry.C | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)

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<cyclicFvPatchVectorField>(C.boundaryField()[patchi]))
+        {
+            // Note: cyclic is not slice but proper field
+            C.boundaryField()[patchi] == static_cast<const vectorField&>
+            (
+                static_cast<const List<vector>&>
+                (
+                    boundary_[patchi].patchSlice(faceCentres())
+                )
+            );
+        }
+    }
 }