diff --git a/src/SimNDT/core/scenario.py b/src/SimNDT/core/scenario.py
index 070515a..c70a321 100644
--- a/src/SimNDT/core/scenario.py
+++ b/src/SimNDT/core/scenario.py
@@ -1,3 +1,9 @@
+#
+# Create the simulation model: Scenario
+#
+# scenario.I is the material label matrix
+# scenario.Iabs is the material label matrix with boundaries
+#
 __author__ = 'Miguel Molero'
 
 
@@ -64,7 +70,7 @@ class Scenario:
         self.Width = int(self.N/float(self.Pixel_mm))
         self.Height = int(self.N/float(self.Pixel_mm))
         self.resetBoundary()
-
+        Label = 255
         for boundary in boundaries:
             if boundary.BC == BC.AirLayer:
                 size = 1
@@ -131,6 +137,7 @@ class Scenario:
 
         if obj.Name == "ellipse":
             print("Adding ellipse to Scenario.")
+            print(obj)
             a = obj.a
             b = obj.b
             theta = obj.theta