diff --git a/test/data/sciver/crab_cube.fits.gz b/test/data/sciver/crab_cube.fits.gz
new file mode 100644
index 0000000..db2ff62
--- /dev/null
+++ b/test/data/sciver/crab_cube.fits.gz
diff --git a/test/data/sciver/crab_cube.xml b/test/data/sciver/crab_cube.xml
new file mode 100644
index 0000000..979f1db
--- /dev/null
+++ b/test/data/sciver/crab_cube.xml
@@ -0,0 +1,20 @@
+<?xml version="1.0" standalone="no"?>
+<source_library title="source library">
+  <source name="Crab" type="DiffuseSource">
+    <spectrum type="PowerLaw">
+      <parameter name="Prefactor" scale="1.0" value="1.0" min="1e-07" max="1000.0" free="1"/>
+      <parameter name="Index"     scale="1.0" value="0.0" min="-5.0"  max="+5.0"   free="1"/>
+      <parameter name="Scale"     scale="1e6" value="1.0" min="0.01"  max="1000.0" free="0"/>
+    </spectrum>
+    <spatialModel type="MapCubeFunction" file="data/sciver/crab_cube.fits">
+      <parameter name="Normalization" scale="1" value="1" min="0.1" max="10" free="0"/>
+    </spatialModel>
+  </source>
+  <source name="Background" type="CTAIrfBackground" instrument="CTA">
+    <spectrum type="PowerLaw">
+      <parameter name="Prefactor" scale="1.0" value="1.0" min="1e-3" max="1e3"    free="1"/>
+      <parameter name="Index"     scale="1.0" value="0.0" min="-5.0" max="+5.0"   free="1"/>
+      <parameter name="Scale"     scale="1e6" value="1.0" min="0.01" max="1000.0" free="0"/>
+    </spectrum>
+  </source>
+</source_library>
diff --git a/test/science_verification.py b/test/science_verification.py
index ab6ae78..717c9c5 100755
--- a/test/science_verification.py
+++ b/test/science_verification.py
@@ -197,8 +197,9 @@ class sciver(gammalib.GPythonTestSuite):
         #self.append(self.spat_rshell, "Test radial shell model")
         #self.append(self.spat_edisk, "Test elliptical disk model")
         #self.append(self.spat_egauss, "Test elliptical Gaussian model")
-        #self.append(self.spat_map, "Test diffuse map model")
-        self.append(self.spat_const, "Test diffuse const model")
+        self.append(self.spat_const, "Test diffuse isotropic model")
+        self.append(self.spat_map, "Test diffuse map model")
+        self.append(self.spat_cube, "Test diffuse cube model")
 
         # Return
         return
@@ -435,6 +436,18 @@ class sciver(gammalib.GPythonTestSuite):
         self.test("Pull_Background_Index")
         return
 
+    # Test diffuse isotropic model
+    def spat_const(self):
+        """
+        Test diffuse isotropic model.
+        """
+        self.pull("data/sciver/crab_const")
+        self.test("Pull_Crab_Prefactor")
+        self.test("Pull_Crab_Index")
+        self.test("Pull_Background_Prefactor")
+        self.test("Pull_Background_Index")
+        return
+
     # Test diffuse map model
     def spat_map(self):
         """
@@ -447,12 +460,12 @@ class sciver(gammalib.GPythonTestSuite):
         self.test("Pull_Background_Index")
         return
 
-    # Test diffuse isotropic model
-    def spat_const(self):
+    # Test diffuse cube model
+    def spat_cube(self):
         """
-        Test diffuse isotropic model.
+        Test diffuse cube model.
         """
-        self.pull("data/sciver/crab_const")
+        self.pull("data/sciver/crab_cube")
         self.test("Pull_Crab_Prefactor")
         self.test("Pull_Crab_Index")
         self.test("Pull_Background_Prefactor")