From 0b16a73a5a9606c5aef3e22e8c9f44249f0318e3 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Tue, 4 Aug 2009 16:43:08 +0100
Subject: [PATCH] better error checking if restart

---
 .../polyTopoChange/polyTopoChange/hexRef8.C          | 12 ++++++++++++
 1 file changed, 12 insertions(+)

diff --git a/src/dynamicMesh/polyTopoChange/polyTopoChange/hexRef8.C b/src/dynamicMesh/polyTopoChange/polyTopoChange/hexRef8.C
index af4a37bead..926f2f4913 100644
--- a/src/dynamicMesh/polyTopoChange/polyTopoChange/hexRef8.C
+++ b/src/dynamicMesh/polyTopoChange/polyTopoChange/hexRef8.C
@@ -350,6 +350,18 @@ void Foam::hexRef8::modFace
 // Bit complex way to determine the unrefined edge length.
 Foam::scalar Foam::hexRef8::getLevel0EdgeLength() const
 {
+    if (cellLevel_.size() != mesh_.nCells())
+    {
+        FatalErrorIn
+        (
+            "hexRef8::getLevel0EdgeLength() const"
+        )   << "Number of cells in mesh:" << mesh_.nCells()
+            << " does not equal size of cellLevel:" << cellLevel_.size()
+            << endl
+            << "This might be because of a restart with inconsistent cellLevel."
+            << abort(FatalError);
+    }
+
     // Determine minimum edge length per refinement level
     // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~