diff --git a/src/modules/mesh/moduleMesh@boundaryEM.f90 b/src/modules/mesh/moduleMesh@boundaryEM.f90
index e86ee5c..43750d3 100644
--- a/src/modules/mesh/moduleMesh@boundaryEM.f90
+++ b/src/modules/mesh/moduleMesh@boundaryEM.f90
@@ -262,7 +262,7 @@ submodule(moduleMesh) boundaryEM
           allocate(mom_nodes(1:edge%nNodes))
           do s = 1, nSpecies
             mom_center = 0.0d0
-            do n = 1, self%nNodes
+            do n = 1, edge%nNodes
               node => mesh%nodes(nodes(n))%obj
 
                 ! Minus sign to get the flux exiting the surface (opposite to the normal) 
@@ -405,7 +405,7 @@ submodule(moduleMesh) boundaryEM
             allocate(mom_nodes(1:edge%nNodes))
 
             do s = 1, nSpecies
-              do n = 1, self%nNodes
+              do n = 1, edge%nNodes
                 node => mesh%nodes(nodes(n))%obj
 
                 ! Minus sign to get the flux exiting the surface (opposite to the normal)