diff --git a/Geo/MElement.cpp b/Geo/MElement.cpp
index 5d0fdc7ef712876035d7d179fec8a1bab90a7ff8..66fd267c4eb731292970a529f7a86be003e314f5 100644
--- a/Geo/MElement.cpp
+++ b/Geo/MElement.cpp
@@ -971,6 +971,7 @@ MElement *MElementFactory::create(int type, std::vector<MVertex*> &v,
   case MSH_QUA_4:  return new MQuadrangle(v, num, part);
   case MSH_QUA_8:  return new MQuadrangle8(v, num, part);
   case MSH_QUA_9:  return new MQuadrangle9(v, num, part);
+  case MSH_QUA_12: return new MQuadrangleN(v, 3, num, part);
   case MSH_QUA_16: return new MQuadrangleN(v, 3, num, part);
   case MSH_QUA_25: return new MQuadrangleN(v, 4, num, part);
   case MSH_QUA_36: return new MQuadrangleN(v, 5, num, part);
diff --git a/tutorial/t6.geo b/tutorial/t6.geo
index b0eb0a19b8ea128ef6487bff05f481e411d36b28..07f850eef488c67ba4cfef9a25f62e2322bd8258 100644
--- a/tutorial/t6.geo
+++ b/tutorial/t6.geo
@@ -53,3 +53,5 @@ Recombine Surface{1};
 
 // Apply an elliptic smoother to the grid
 Mesh.Smoothing = 100;
+
+Physical Surface(1) = 1;