Fix writing quadrilateral faces, the index order in reference elements changed a while ago

dune/fufem/boundarywriter.hh

@@ -125,9 +125,9 @@ void writeBoundary(const BoundaryPatch<GridView>& surface, const std::string& fi
 
         if (n==4)
             fprintf(fp, "  %d %d %d\n",
-                    globalToLocal[indexSet.subIndex(inside, refElement.subEntity(it.indexInInside(), 1, 2, dim), dim)] + 1,
                     globalToLocal[indexSet.subIndex(inside, refElement.subEntity(it.indexInInside(), 1, 3, dim), dim)] + 1,
-                    globalToLocal[indexSet.subIndex(inside, refElement.subEntity(it.indexInInside(), 1, 0, dim), dim)] + 1);
+                    globalToLocal[indexSet.subIndex(inside, refElement.subEntity(it.indexInInside(), 1, 2, dim), dim)] + 1,
+                    globalToLocal[indexSet.subIndex(inside, refElement.subEntity(it.indexInInside(), 1, 1, dim), dim)] + 1);
 
     }