diff --git a/inputFiles/wavePropagation/benchmarks/acouselas3D_flu_komatitsch.xml b/inputFiles/wavePropagation/benchmarks/acouselas3D_flu_komatitsch.xml
index 44cf59b73d9..fe41cbb0260 100644
--- a/inputFiles/wavePropagation/benchmarks/acouselas3D_flu_komatitsch.xml
+++ b/inputFiles/wavePropagation/benchmarks/acouselas3D_flu_komatitsch.xml
@@ -31,7 +31,7 @@
     <FiniteElements>
       <FiniteElementSpace
         name="FE1"
-        order="3"
+        order="1"
         formulation="SEM"/>
     </FiniteElements>
   </NumericalMethods>
diff --git a/src/coreComponents/physicsSolvers/wavePropagation/AcousticElasticWaveEquationSEM.cpp b/src/coreComponents/physicsSolvers/wavePropagation/AcousticElasticWaveEquationSEM.cpp
index 0537694f73e..576af99a37c 100644
--- a/src/coreComponents/physicsSolvers/wavePropagation/AcousticElasticWaveEquationSEM.cpp
+++ b/src/coreComponents/physicsSolvers/wavePropagation/AcousticElasticWaveEquationSEM.cpp
@@ -39,12 +39,23 @@ void AcousticElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups
   auto acousNodesSet = acousticSolver()->getSolverNodesSet();
   auto elasNodesSet = elasticSolver()->getSolverNodesSet();
 
+  if (false)  // debug
+  {
+    std::cout << "acoustic nodes" << std::endl;
+    for (auto val : acousNodesSet) std::cout << val << ' ';
+    std::cout << std::endl;
+    std::cout << "elastic nodes" << std::endl;
+    for (auto val : elasNodesSet) std::cout << val << ' ';
+    std::cout << std::endl;
+  }
+
   for (auto val : acousNodesSet)
   {
     if (elasNodesSet.contains(val)) m_interfaceNodesSet.insert(val);
   }
 
-  // std::cout << "\t[AcousticElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] size=" << m_solverTargetNodesSet.size() << std::endl;
+  std::cout << "\t[AcousticElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] "
+  << "m_interfaceNodesSet.size()=" << m_interfaceNodesSet.size() << std::endl;
   GEOS_THROW_IF( m_interfaceNodesSet.size() == 0, "Failed to compute interface: check xml input (solver order)", std::runtime_error );
 }
 
diff --git a/src/coreComponents/physicsSolvers/wavePropagation/AcousticWaveEquationSEM.cpp b/src/coreComponents/physicsSolvers/wavePropagation/AcousticWaveEquationSEM.cpp
index 17ba3ece579..0c60da01cde 100644
--- a/src/coreComponents/physicsSolvers/wavePropagation/AcousticWaveEquationSEM.cpp
+++ b/src/coreComponents/physicsSolvers/wavePropagation/AcousticWaveEquationSEM.cpp
@@ -320,7 +320,7 @@ void AcousticWaveEquationSEM::initializePostInitialConditionsPreSubGroups()
       } );
     } );
   } );
-  // std::cout << "\t[AcousticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] size=" << m_solverTargetNodesSet.size() << std::endl;
+  std::cout << "\t[AcousticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] m_solverTargetNodesSet.size()=" << m_solverTargetNodesSet.size() << std::endl;
 
 }
 
diff --git a/src/coreComponents/physicsSolvers/wavePropagation/ElasticWaveEquationSEM.cpp b/src/coreComponents/physicsSolvers/wavePropagation/ElasticWaveEquationSEM.cpp
index 4f5eac17f10..4e730e625fa 100644
--- a/src/coreComponents/physicsSolvers/wavePropagation/ElasticWaveEquationSEM.cpp
+++ b/src/coreComponents/physicsSolvers/wavePropagation/ElasticWaveEquationSEM.cpp
@@ -599,7 +599,7 @@ void ElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups()
     } );
   } );
   
-  // std::cout << "\t[ElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] size=" << m_solverTargetNodesSet.size() << std::endl;
+  std::cout << "\t[ElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] m_solverTargetNodesSet.size()=" << m_solverTargetNodesSet.size() << std::endl;
 }