Correcting some issues with the floating boundary condition

src/modules/init/moduleInput.f90

@@ -1344,6 +1344,13 @@ MODULE moduleInput
 
             ! Specific assignments per type
             select type(bound)
+            type is(boundaryEMFloating)
+              if (associated(physicalSurfaces(ps)%EM, bound)) then
+                bound%edges = [bound%edges, physicalSurfaces(ps)%edges]
+                bound%nEdges = size(bound%edges)
+
+              end if
+
             type is(boundaryEMFreeCurrent)
               if (associated(physicalSurfaces(ps)%EM, bound)) then
                 bound%edges = [bound%edges, physicalSurfaces(ps)%edges]

src/modules/mesh/moduleMesh.f90

@@ -1089,6 +1089,7 @@ MODULE moduleMesh
     real(8):: invC ! Inverse of the capacitance of the surface
     real(8):: charge ! Charge accumulated on surface
     type(meshEdgePointer), allocatable:: edges(:) !Array with edges
+    integer:: nEdges
     contains
       ! boundaryEMGeneric DEFERRED PROCEDURES
       procedure, pass:: apply => applyFloating

src/modules/mesh/moduleMesh@boundaryEM.f90

@@ -252,7 +252,7 @@ submodule(moduleMesh) boundaryEM
 
       select type(self)
       type is(boundaryEMFloating)
-        do e = 1, size(self%edges)
+        do e = 1, self%nEdges
           edge => self%edges(e)%obj
 
           chargeTau = 0.0d0