diff --git a/test/science_verification.py b/test/science_verification.py
index 98ca005..ab6ae78 100755
--- a/test/science_verification.py
+++ b/test/science_verification.py
@@ -191,13 +191,14 @@ class sciver(gammalib.GPythonTestSuite):
         #self.append(self.spec_nodes, "Test nodes model")
 
         # Append spatial tests
-        #self.append(self.spec_ptsrc, "Test point source model")
-        #self.append(self.spec_rdisk, "Test radial disk model")
-        #self.append(self.spec_rgauss, "Test radial Gaussian model")
-        #self.append(self.spec_rshell, "Test radial shell model")
-        #self.append(self.spec_edisk, "Test elliptical disk model")
-        #self.append(self.spec_egauss, "Test elliptical Gaussian model")
-        self.append(self.spec_map, "Test diffuse map model")
+        #self.append(self.spat_ptsrc, "Test point source model")
+        #self.append(self.spat_rdisk, "Test radial disk model")
+        #self.append(self.spat_rgauss, "Test radial Gaussian model")
+        #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")
 
         # Return
         return
@@ -341,7 +342,7 @@ class sciver(gammalib.GPythonTestSuite):
         return
 
     # Test point source model
-    def spec_ptsrc(self):
+    def spat_ptsrc(self):
         """
         Test point source model.
         """
@@ -355,7 +356,7 @@ class sciver(gammalib.GPythonTestSuite):
         return
 
     # Test radial disk model
-    def spec_rdisk(self):
+    def spat_rdisk(self):
         """
         Test radial disk model.
         """
@@ -370,7 +371,7 @@ class sciver(gammalib.GPythonTestSuite):
         return
 
     # Test radial Gaussian model
-    def spec_rgauss(self):
+    def spat_rgauss(self):
         """
         Test radial Gaussian model.
         """
@@ -385,7 +386,7 @@ class sciver(gammalib.GPythonTestSuite):
         return
 
     # Test radial shell model
-    def spec_rshell(self):
+    def spat_rshell(self):
         """
         Test radial shell model.
         """
@@ -401,7 +402,7 @@ class sciver(gammalib.GPythonTestSuite):
         return
 
     # Test elliptical disk model
-    def spec_edisk(self):
+    def spat_edisk(self):
         """
         Test elliptical disk model.
         """
@@ -418,7 +419,7 @@ class sciver(gammalib.GPythonTestSuite):
         return
 
     # Test elliptical Gaussian model
-    def spec_egauss(self):
+    def spat_egauss(self):
         """
         Test elliptical Gaussian model.
         """
@@ -435,7 +436,7 @@ class sciver(gammalib.GPythonTestSuite):
         return
 
     # Test diffuse map model
-    def spec_map(self):
+    def spat_map(self):
         """
         Test diffuse map model.
         """
@@ -446,6 +447,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
+
 
 # ======================== #
 # Main routine entry point #