diff --git a/ProcessLib/BoundaryCondition/BoundaryCondition.cpp b/ProcessLib/BoundaryCondition/BoundaryCondition.cpp
index ece058652cab9c7b95924506e5534f02731bc08b..0727ea6fe93cdb9eaef3cce8cdc25a26000756e3 100644
--- a/ProcessLib/BoundaryCondition/BoundaryCondition.cpp
+++ b/ProcessLib/BoundaryCondition/BoundaryCondition.cpp
@@ -102,6 +102,9 @@ BoundaryConditionBuilder::createDirichletBoundaryCondition(
         });
     ids.erase(ids_new_end_iterator, std::end(ids));
 
+    DBUG("Found %d nodes for Dirichlet BCs for the variable %d and componet %d",
+         ids.size(), variable_id, config.component_id);
+
     return ProcessLib::createDirichletBoundaryCondition(
         config.config, std::move(ids), dof_table, mesh.getID(), variable_id,
         config.component_id, parameters);
diff --git a/ProcessLib/BoundaryCondition/GenericNaturalBoundaryCondition-impl.h b/ProcessLib/BoundaryCondition/GenericNaturalBoundaryCondition-impl.h
index 25e8988090a56fc94d5194228ff59f66ab5508a3..a6681a76cbfe4817d6f535f3cc7d1255374d0d31 100644
--- a/ProcessLib/BoundaryCondition/GenericNaturalBoundaryCondition-impl.h
+++ b/ProcessLib/BoundaryCondition/GenericNaturalBoundaryCondition-impl.h
@@ -39,6 +39,8 @@ GenericNaturalBoundaryCondition<BoundaryConditionData,
            static_cast<int>(dof_table_bulk.getNumberOfComponents()));
 
     std::vector<MeshLib::Node*> nodes = MeshLib::getUniqueNodes(_elements);
+    DBUG("Found %d nodes for Natural BCs for the varialbe %d and component %d",
+         nodes.size(), variable_id, component_id);
 
     auto const& mesh_subsets =
         dof_table_bulk.getMeshSubsets(variable_id, component_id);
diff --git a/ProcessLib/LIE/BoundaryCondition/GenericNaturalBoundaryCondition-impl.h b/ProcessLib/LIE/BoundaryCondition/GenericNaturalBoundaryCondition-impl.h
index 7c8815f61a8e3349c4d0fa93b9eccdc45c802532..aa9e9cd0417c9e4dbd7246070797bde610a6351e 100644
--- a/ProcessLib/LIE/BoundaryCondition/GenericNaturalBoundaryCondition-impl.h
+++ b/ProcessLib/LIE/BoundaryCondition/GenericNaturalBoundaryCondition-impl.h
@@ -47,6 +47,8 @@ GenericNaturalBoundaryCondition<BoundaryConditionData,
            static_cast<int>(dof_table_bulk.getNumberOfComponents()));
 
     std::vector<MeshLib::Node*> nodes = MeshLib::getUniqueNodes(_elements);
+    DBUG("Found %d nodes for Natural BCs for the varialbe %d and component %d",
+         nodes.size(), variable_id, component_id);
 
     auto const& mesh_subsets =
         dof_table_bulk.getMeshSubsets(variable_id, component_id);