diff --git a/Geo/MElementCut.cpp b/Geo/MElementCut.cpp
index 7331cb5862dac454567213d6973e58920ffd85cc..7c80f70d1b47704b1e52e2e1d59f9fa0ab30fe2d 100644
--- a/Geo/MElementCut.cpp
+++ b/Geo/MElementCut.cpp
@@ -761,6 +761,7 @@ GModel *buildCutMesh(GModel *gm, gLevelset *ls,
   for(unsigned int i = 0; i < gmEntities.size(); i++) {
     for(unsigned int j = 0; j < gmEntities[i]->getNumMeshElements(); j++) {
       MElement *e = gmEntities[i]->getMeshElement(j);
+      e->setVolumePositive();
       elementCutMesh (e, ls, gmEntities[i], gm, numEle,
                       vertexMap, newVertices, elements, border, physicals, entityCut);
       cutGM->getMeshPartitions().insert(e->getPartition());