diff --git a/src/modules/mesh/moduleMesh@boundaryEM.f90 b/src/modules/mesh/moduleMesh@boundaryEM.f90
index 5d4fa58..26363b5 100644
--- a/src/modules/mesh/moduleMesh@boundaryEM.f90
+++ b/src/modules/mesh/moduleMesh@boundaryEM.f90
@@ -207,6 +207,7 @@ submodule(moduleMesh) boundaryEM
         do e = 1, size(self%edges)
           edge => self%edges(e)%obj
 
+          chargeTau = 0.0d0
           edgeDensityCurrent = 0.0d0
 
           nodes = edge%getNodes(edge%nNodes)
@@ -257,7 +258,7 @@ submodule(moduleMesh) boundaryEM
       select type(self)
       type is(boundaryEMFloating)
         write(fileID, '(A,",",A)') 'Total charge (C)', 'Potential (V)'
-        write(fileID, '(*('//fmtColReal//'))') self%charge*qe*v_ref*n_ref*L_ref**2*ti_ref, self%potential * Volt_ref
+        write(fileID, '(*('//fmtColReal//'))') self%charge*qe*v_ref*n_ref*L_ref**2*ti_ref, self%potential*Volt_ref
 
       end select