diff --git a/example/core_tutorial/04_aosoa_advanced_unmanaged/advanced_aosoa_unmanaged.cpp b/example/core_tutorial/04_aosoa_advanced_unmanaged/advanced_aosoa_unmanaged.cpp
index fc6315b95..ed476ca07 100644
--- a/example/core_tutorial/04_aosoa_advanced_unmanaged/advanced_aosoa_unmanaged.cpp
+++ b/example/core_tutorial/04_aosoa_advanced_unmanaged/advanced_aosoa_unmanaged.cpp
@@ -13,6 +13,7 @@
 #include <Cabana_Sort.hpp>
 
 #include <cstdio>
+#include <iostream>
 
 //---------------------------------------------------------------------------//
 // Example of an unmanaged AoSoA.
diff --git a/example/grid_tutorial/11_semi_structured_solver_multi_variate/hypre_semi_structured_solver_multi_example.cpp b/example/grid_tutorial/11_semi_structured_solver_multi_variate/hypre_semi_structured_solver_multi_example.cpp
index bc34b28f8..2e7a6ec26 100644
--- a/example/grid_tutorial/11_semi_structured_solver_multi_variate/hypre_semi_structured_solver_multi_example.cpp
+++ b/example/grid_tutorial/11_semi_structured_solver_multi_variate/hypre_semi_structured_solver_multi_example.cpp
@@ -16,6 +16,7 @@
 #include <mpi.h>
 
 #include <array>
+#include <iostream>
 #include <vector>
 
 //---------------------------------------------------------------------------//
diff --git a/example/grid_tutorial/11_structured_solver/structured_solver_example.cpp b/example/grid_tutorial/11_structured_solver/structured_solver_example.cpp
index 7cb002b2e..b0cc1ca4e 100644
--- a/example/grid_tutorial/11_structured_solver/structured_solver_example.cpp
+++ b/example/grid_tutorial/11_structured_solver/structured_solver_example.cpp
@@ -16,6 +16,7 @@
 #include <mpi.h>
 
 #include <array>
+#include <iostream>
 #include <vector>
 
 //---------------------------------------------------------------------------//
diff --git a/example/grid_tutorial/11_structured_solver_hypre/hypre_structured_solver_example.cpp b/example/grid_tutorial/11_structured_solver_hypre/hypre_structured_solver_example.cpp
index ff830a9b2..8724084fd 100644
--- a/example/grid_tutorial/11_structured_solver_hypre/hypre_structured_solver_example.cpp
+++ b/example/grid_tutorial/11_structured_solver_hypre/hypre_structured_solver_example.cpp
@@ -16,6 +16,7 @@
 #include <mpi.h>
 
 #include <array>
+#include <iostream>
 #include <vector>
 
 //---------------------------------------------------------------------------//
diff --git a/grid/src/Cabana_Grid_ReferenceStructuredSolver.hpp b/grid/src/Cabana_Grid_ReferenceStructuredSolver.hpp
index 3bbee97fc..ba4fad696 100644
--- a/grid/src/Cabana_Grid_ReferenceStructuredSolver.hpp
+++ b/grid/src/Cabana_Grid_ReferenceStructuredSolver.hpp
@@ -29,6 +29,7 @@
 #include <Kokkos_Core.hpp>
 
 #include <array>
+#include <iostream>
 #include <memory>
 #include <numeric>
 #include <set>
diff --git a/grid/unit_test/tstLocalMesh2d.hpp b/grid/unit_test/tstLocalMesh2d.hpp
index a4a180fa6..83b747c3d 100644
--- a/grid/unit_test/tstLocalMesh2d.hpp
+++ b/grid/unit_test/tstLocalMesh2d.hpp
@@ -24,6 +24,7 @@
 
 #include <array>
 #include <cmath>
+#include <iostream>
 #include <vector>
 
 using namespace Cabana::Grid;