diff --git a/src/mesh.h b/src/mesh.h
--- a/src/mesh.h
+++ b/src/mesh.h
@@ -83,6 +83,7 @@ class Mesh {
     time = 0.;
     plugin = 0;
     boundary_polygon=0;
+    iterations = 0;
   };
   ~Mesh(void) {
     if (boundary_polygon) {
@@ -391,6 +392,10 @@ class Mesh {
   
   Node* findNextBoundaryNode(Node*);
 
+  int iterations;
+  int getIterations(){return this->iterations;}
+  void incrementIterations(){this->iterations++;}
+
  private:
 
   // Data members