diff --git a/Geo/GEdgeLoop.cpp b/Geo/GEdgeLoop.cpp
index c29b288a8b586b0ab7747c19cd2835258101e7de..7025cf611d8671a334a600641363bfebb505c550 100644
--- a/Geo/GEdgeLoop.cpp
+++ b/Geo/GEdgeLoop.cpp
@@ -61,6 +61,7 @@ GEdgeSigned nextOne(GEdgeSigned *thisOne, std::list<GEdge*> &wire)
       if(v1 == gv) return GEdgeSigned(1, ge);   
       if(v2 == gv) return GEdgeSigned(-1, ge);   
       Msg::Error("Something wrong in edge loop");
+      thisOne->print();
     }
     ++it;
   }
@@ -77,6 +78,7 @@ GEdgeSigned nextOne(GEdgeSigned *thisOne, std::list<GEdge*> &wire)
       if(v1 == gv) return GEdgeSigned(1, ge);   
       if(v2 == gv) return GEdgeSigned(-1, ge);
       Msg::Error("Something wrong in edge loop");
+      thisOne->print();
     }   
     ++it;
   }
@@ -107,7 +109,7 @@ GEdgeLoop::GEdgeLoop(const std::list<GEdge*> &cwire)
   while(wire.size()){
     ges = nextOne(prevOne, wire);
     if(ges.getSign() == 0){ // oops
-      Msg::Error("Something wrong in edge loop");
+      Msg::Error("Something wrong in edge loop, no sign !");
       break;
     }
     prevOne = &ges;
diff --git a/Geo/GModel.cpp b/Geo/GModel.cpp
index bc96d1eab25b04e7a2c1bd077e3c2c246841fd4e..45bff1cc8633e5c0b6eb06e1120f6c48f2464b0f 100644
--- a/Geo/GModel.cpp
+++ b/Geo/GModel.cpp
@@ -28,6 +28,7 @@
 #include "SmoothData.h"
 #include "Context.h"
 #include "OS.h"
+#include "GEdgeLoop.h"
 
 #include "OpenFile.h"
 #include "CreateFile.h"
@@ -1495,6 +1496,17 @@ void GModel::createTopologyFromFaces(std::vector<discreteFace*> &discFaces)
     (*it)->parametrize(face2Vert, region2Vert);
   }
 
+  //fill edgeLoops of Faces or correct sign of l_edges
+  // for (std::vector<discreteFace*>::iterator itF = discFaces.begin();
+  //      itF != discFaces.end(); itF++){
+  //   //EMI, TODO
+  //   std::list<GEdgeLoop> edgeLoops = (*itF)->edgeLoops; 
+  //   edgeLoops.clear();
+  //   GEdgeLoop el((*itF)->edges());
+  //   edgeLoops.push_back(el);
+  // }
+
+
  //we need to recreate lines, triangles and tets
   //that contain those new MEdgeVertices
   for (std::map<GFace*, std::map<MVertex*, MVertex*, std::less<MVertex*> > >::iterator