diff --git a/ProcessLib/ProcessVariable.h b/ProcessLib/ProcessVariable.h
index 5ce2566e92d05d632e54accfeabd6e7a6d8633b0..3b2e05b06ee211bc31ffa5c82ea82c4b278bad1e 100644
--- a/ProcessLib/ProcessVariable.h
+++ b/ProcessLib/ProcessVariable.h
@@ -67,10 +67,14 @@ public:
         const int variable_id,
         const int component_id)
     {
+        // Find all boundary conditions matching the component id. There can be
+        // more than one such boundary condition.
         for (auto& bc_config : _dirichlet_bc_configs)
         {
             if (bc_config.second != component_id)
                 continue;
+            // Create/initialize the boundary condition with matching component
+            // id and output it through the OutputIterator.
             DirichletBc<GlobalIndexType> bc;
             bc_config.first->initialize(searcher, dof_table, variable_id,
                                         component_id, bc);
@@ -86,10 +90,15 @@ public:
                           const int variable_id,
                           const int component_id)
     {
+        // Find all boundary conditions matching the component id. There can be
+        // more than one such boundary condition.
         for (auto& bc_config : _neumann_bc_configs)
         {
             if (bc_config.second != component_id)
                 continue;
+
+            // Create/initialize the boundary condition with matching component
+            // id and output it through the OutputIterator.
             bc_config.first->initialize(searcher);
             bcs++ = std::unique_ptr<NeumannBc<GlobalSetup>>{
                 new NeumannBc<GlobalSetup>(*bc_config.first,