diff --git a/rod-eoc.cc b/rod-eoc.cc
index e7f0f0f808a20da6220e1495e6f4f1ef0b27e708..d99a6dcc36445d4f6ede36bbe561ed0a67a20ff1 100644
--- a/rod-eoc.cc
+++ b/rod-eoc.cc
@@ -17,6 +17,7 @@
 #include "src/rodassembler.hh"
 #include "src/riemanniantrsolver.hh"
 #include "src/rodrefine.hh"
+#include "src/rodwriter.hh"
 
 typedef Dune::OneDGrid GridType;
 
@@ -191,6 +192,10 @@ int main (int argc, char *argv[]) try
         for (int j=i; j<numLevels; j++)
             globalRodRefine(grid, solution);
 
+        std::stringstream numberAsAscii;
+        numberAsAscii << i;
+        writeRod(solution, "rodGrid_" + numberAsAscii.str());
+
         assert(referenceSolution.size() == solution.size());
 
         // Compute max-norm difference