diff --git a/src/modules/mesh/moduleMesh@boundaryEM.f90 b/src/modules/mesh/moduleMesh@boundaryEM.f90
index 3ca875f..0be427a 100644
--- a/src/modules/mesh/moduleMesh@boundaryEM.f90
+++ b/src/modules/mesh/moduleMesh@boundaryEM.f90
@@ -413,20 +413,23 @@ submodule(moduleMesh) boundaryEM
                 ! Minus sign as we look at the values exiting the surface (opposite to the normal) 
                 ! Density and momentum are reescaled as 'den' has no data about the node volume
                 mom_nodes(n) = - dot_product(node%output(s)%mom, edge%normal)/(node%v*Vol_ref*n_ref)
-                den_nodes(n) =   node%output(s)%den/(node%v*Vol_ref*n_ref)
+                ! den_nodes(n) =   node%output(s)%den/(node%v*Vol_ref*n_ref)
 
               end do
 
-              mom_center    = mom_center    + edge%gatherF(edge%centerXi(), edge%nNodes, mom_nodes)
-              den_center    = den_center    + edge%gatherF(edge%centerXi(), edge%nNodes, den_nodes)
-              charge_center = charge_center + qSpecies(s)*edge%gatherF(edge%centerXi(), edge%nNodes, den_nodes)
+              ! mom_center    = mom_center    + edge%gatherF(edge%centerXi(), edge%nNodes, mom_nodes)
+              ! den_center    = den_center    + edge%gatherF(edge%centerXi(), edge%nNodes, den_nodes)
+              ! charge_center = charge_center + qSpecies(s)*edge%gatherF(edge%centerXi(), edge%nNodes, den_nodes)
+
+              mom_center    = mom_center    + qSpecies(s)*edge%gatherF(edge%centerXi(), edge%nNodes, mom_nodes)
 
             end do
 
-            if (den_center > 1.0e-12) then
-              self%deltaElectricField(e) = self%deltaElectricField(e) + mom_center/den_center * charge_center * tauMin
+            ! if (den_center > 1.0e-12) then
+              ! self%deltaElectricField(e) = self%deltaElectricField(e) + mom_center/den_center * charge_center * tauMin
 
-            end if
+            ! end if
+            self%deltaElectricField(e) = self%deltaElectricField(e) + mom_center * tauMin
 
             deallocate(nodes, mom_nodes, den_nodes)