diff --git a/Mesh/meshRefine.cpp b/Mesh/meshRefine.cpp
index 899613cbf43a04461f89ec425d6858845052c3af..368f937056648ed9dde1753eed63b1e343cdecb2 100644
--- a/Mesh/meshRefine.cpp
+++ b/Mesh/meshRefine.cpp
@@ -197,6 +197,7 @@ static void Subdivide(GRegion *gr, bool splitIntoHexas, faceContainer &faceVerti
           else{
             SPoint3 pc = face.barycenter();
             newv.push_back(new MVertex(pc.x(), pc.y(), pc.z(), gr));
+            faceVertices[face].push_back(newv.back());
             gr->mesh_vertices.push_back(newv.back());
           }
         }
@@ -233,6 +234,7 @@ static void Subdivide(GRegion *gr, bool splitIntoHexas, faceContainer &faceVerti
           else{
             SPoint3 pc = face.barycenter();
             newv.push_back(new MVertex(pc.x(), pc.y(), pc.z(), gr));
+            faceVertices[face].push_back(newv.back());
             gr->mesh_vertices.push_back(newv.back());
           }
         }