diff --git a/src/modules/mesh/moduleMesh@boundaryEM.f90 b/src/modules/mesh/moduleMesh@boundaryEM.f90
index 26363b5..f20804a 100644
--- a/src/modules/mesh/moduleMesh@boundaryEM.f90
+++ b/src/modules/mesh/moduleMesh@boundaryEM.f90
@@ -224,9 +224,6 @@ submodule(moduleMesh) boundaryEM
 
             mom_center = edge%gatherF(edge%centerXi(), edge%nNodes, mom_nodes)
 
-            ! Only account for charge exiting the surface
-            mom_center = max(mom_center, 0.0d0)
-            
             edgeDensityCurrent = edgeDensityCurrent + qSpecies(s) * mom_center
 
           end do