diff --git a/Isosurfacing_3/examples/Isosurfacing_3/contouring_image.cpp b/Isosurfacing_3/examples/Isosurfacing_3/contouring_image.cpp
index 18d02a792b7..707c6061070 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/contouring_image.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/contouring_image.cpp
@@ -95,17 +95,23 @@ int main(int argc, char* argv[])
   }
 
   // convert image to a Cartesian grid
-  auto [grid, values] = IS::IO::read_Image_3<Kernel>(image);
+  Grid grid;
+  Values values { grid };
+  if(!IS::IO::read_Image_3(image, grid, values))
+  {
+    std::cerr << "Error: Cannot convert image to Cartesian grid" << std::endl;
+    return EXIT_FAILURE;
+  }
+
+  std::cout << "Bbox: " << grid.bbox() << std::endl;
+  std::cout << "Cell dimensions: " << grid.spacing()[0] << " " << grid.spacing()[1] << " " << grid.spacing()[2] << std::endl;
+  std::cout << "Cell #: " << grid.xdim() << ", " << grid.ydim() << ", " << grid.zdim() << std::endl;
 
   for (std::size_t i=0; i<grid.xdim(); ++i)
     for (std::size_t j=0; j<grid.ydim(); ++j)
       for (std::size_t k=0; k<grid.zdim(); ++k)
         values(i, j, k) = - values(i, j, k); // inside out
 
-  std::cout << "Bbox: " << grid.bbox() << std::endl;
-  std::cout << "Cell dimensions: " << grid.spacing()[0] << " " << grid.spacing()[1] << " " << grid.spacing()[2] << std::endl;
-  std::cout << "Cell #: " << grid.xdim() << ", " << grid.ydim() << ", " << grid.zdim() << std::endl;
-
   run_marching_cubes(grid, isovalue, values);
 
   run_dual_contouring(grid, isovalue, values);