diff --git a/dumux/boxmodels/common/boxfvelementgeometry.hh b/dumux/boxmodels/common/boxfvelementgeometry.hh
index bf2cc56c51395ece30cc25056c322000eca3bef0..7c39894bfec451ae0bcb3df182b7e513e1a8df6e 100644
--- a/dumux/boxmodels/common/boxfvelementgeometry.hh
+++ b/dumux/boxmodels/common/boxfvelementgeometry.hh
@@ -779,6 +779,7 @@ public:
                     {
                         jacInvT.mv(localJac[vert][0], boundaryFace[bfIdx].grad[vert]);
                         boundaryFace[bfIdx].shapeValue[vert] = Scalar(shapeVal[vert]);
+                        boundaryFace[bfIdx].fapIndices[vert] = vert;
                     }
 
                     //                    std::cout << "boundary face " << face << ", vert = " << vertInElement << ", ipLocal = "