diff --git a/test/mixeddimension/embedded/1p2c_richards2c/richardstestproblem.hh b/test/mixeddimension/embedded/1p2c_richards2c/richardstestproblem.hh
index 4350a2bd431408d57bbec81eca6030f7103d399c..97737e10a87c65021b26f933841942a45ad56374 100644
--- a/test/mixeddimension/embedded/1p2c_richards2c/richardstestproblem.hh
+++ b/test/mixeddimension/embedded/1p2c_richards2c/richardstestproblem.hh
@@ -81,7 +81,7 @@ SET_BOOL_PROP(RichardsTestProblem, VtkAddVelocity, true);
 // Set the grid parameter group
 SET_STRING_PROP(RichardsTestProblem, GridParameterGroup, "SoilGrid");
 
-// Use mass fractions
+// Use mole fractions
 SET_BOOL_PROP(RichardsTestProblem, UseMoles, true);
 }
 
@@ -235,6 +235,7 @@ public:
         sourceValues[conti0EqIdx] = 2* M_PI *rootRadius * Kr *(lowDimVolVars.pressure(wPhaseIdx) - bulkVolVars.pressure(wPhaseIdx))
                                     *bulkVolVars.molarDensity(wPhaseIdx);
 
+        //! advective transport over root wall
         // compute correct upwind concentration
         if (sourceValues[conti0EqIdx] > 0)
             sourceValues[transportEqIdx] = sourceValues[conti0EqIdx]
@@ -243,6 +244,7 @@ public:
             sourceValues[transportEqIdx] = sourceValues[conti0EqIdx]
                                             * bulkVolVars.moleFraction(wPhaseIdx, transportCompIdx);
 
+        //! diffusive transport over root wall
         sourceValues[transportEqIdx] += 2* M_PI *rootRadius * 1.0e-8
                                         *(lowDimVolVars.moleFraction(wPhaseIdx, transportCompIdx) - bulkVolVars.moleFraction(wPhaseIdx, transportCompIdx))
                                         *0.5*(bulkVolVars.molarDensity(wPhaseIdx) + lowDimVolVars.molarDensity(wPhaseIdx));
diff --git a/test/mixeddimension/embedded/1p2c_richards2c/rootsystemtestproblem.hh b/test/mixeddimension/embedded/1p2c_richards2c/rootsystemtestproblem.hh
index 26a95ff07d0785d304216c5f17ccc5f4923c8871..b9343bbca670e486fb1b205418e4edb2baa4cc7c 100644
--- a/test/mixeddimension/embedded/1p2c_richards2c/rootsystemtestproblem.hh
+++ b/test/mixeddimension/embedded/1p2c_richards2c/rootsystemtestproblem.hh
@@ -77,7 +77,7 @@ SET_BOOL_PROP(RootsystemTestProblem, ProblemEnableGravity, true);
 // Enable velocity output
 SET_BOOL_PROP(RootsystemTestProblem, VtkAddVelocity, true);
 
-// Use mass fractions
+// Use mole fractions
 SET_BOOL_PROP(RootsystemTestProblem, UseMoles, true);
 }
 
@@ -315,6 +315,7 @@ public:
         sourceValues[conti0EqIdx] = 2* M_PI *rootRadius * Kr *(bulkVolVars.pressure(wPhaseIdx) - lowDimVolVars.pressure(wPhaseIdx))
                                     *bulkVolVars.molarDensity(wPhaseIdx);
 
+        //! advective transport over root wall
         // compute correct upwind concentration
         if (sourceValues[conti0EqIdx] < 0)
             sourceValues[transportEqIdx] = sourceValues[conti0EqIdx]
@@ -323,6 +324,7 @@ public:
             sourceValues[transportEqIdx] = sourceValues[conti0EqIdx]
                                             * bulkVolVars.moleFraction(wPhaseIdx, transportCompIdx);
 
+        //! diffusive transport over root wall
         sourceValues[transportEqIdx] += 2* M_PI *rootRadius * 1.0e-8
                                         *(bulkVolVars.moleFraction(wPhaseIdx, transportCompIdx) - lowDimVolVars.moleFraction(wPhaseIdx, transportCompIdx))
                                         *0.5*(bulkVolVars.molarDensity(wPhaseIdx) + lowDimVolVars.molarDensity(wPhaseIdx));