diff --git a/src/modules/mesh/moduleMesh@boundaryEM.f90 b/src/modules/mesh/moduleMesh@boundaryEM.f90
index 094d671..e86ee5c 100644
--- a/src/modules/mesh/moduleMesh@boundaryEM.f90
+++ b/src/modules/mesh/moduleMesh@boundaryEM.f90
@@ -265,7 +265,8 @@ submodule(moduleMesh) boundaryEM
             do n = 1, self%nNodes
               node => mesh%nodes(nodes(n))%obj
 
-              ! Minus sign to get the flux exiting the edge
+                ! Minus sign to get the flux exiting the surface (opposite to the normal) 
+                ! Momentum is reescaled as value at the node has no data about the node volume
               mom_nodes(n) = - dot_product(node%output(s)%mom, edge%normal)/(node%v*Vol_ref*n_ref)
 
             end do
@@ -407,7 +408,7 @@ submodule(moduleMesh) boundaryEM
               do n = 1, self%nNodes
                 node => mesh%nodes(nodes(n))%obj
 
-                ! Minus sign as we look at the values exiting the surface (opposite to the normal) 
+                ! Minus sign to get the flux exiting the surface (opposite to the normal) 
                 ! Momentum is reescaled as value at the node has no data about the node volume
                 mom_nodes(n) = - dot_product(node%output(s)%mom, edge%normal)/(node%v*Vol_ref*n_ref)