Forgot to link edge to boundary procedures. 1D case working right now

src/modules/mesh/1DCart/moduleMesh1DCart.f90

@@ -130,6 +130,13 @@ MODULE moduleMesh1DCart
       call meshEdgePointer_add(physicalSurfaces(ps_index)%edges, self%n)
       call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n1%n)
 
+      ! Link edge to particle boundaries
+      allocate(self%boundariesParticle(1:nSpecies))
+      do s = 1, nSpecies
+        self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj
+
+      end do
+
     end subroutine initEdge1DCart
 
     !Get nodes from edge

src/modules/mesh/1DRad/moduleMesh1DRad.f90

@@ -130,6 +130,13 @@ MODULE moduleMesh1DRad
       call meshEdgePointer_add(physicalSurfaces(ps_index)%edges, self%n)
       call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n1%n)
 
+      ! Link edge to particle boundaries
+      allocate(self%boundariesParticle(1:nSpecies))
+      do s = 1, nSpecies
+        self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj
+
+      end do
+
     end subroutine initEdge1DRad
 
     !Get nodes from edge

src/modules/mesh/2DCart/moduleMesh2DCart.f90

@@ -170,6 +170,13 @@ MODULE moduleMesh2DCart
       call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n1%n)
       call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n2%n)
 
+      ! Link edge to particle boundaries
+      allocate(self%boundariesParticle(1:nSpecies))
+      do s = 1, nSpecies
+        self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj
+
+      end do
+
     END SUBROUTINE initEdge2DCart
 
     !Get nodes from edge

src/modules/mesh/2DCyl/moduleMesh2DCyl.f90

@@ -179,6 +179,13 @@ MODULE moduleMesh2DCyl
       call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n1%n)
       call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n2%n)
 
+      ! Link edge to particle boundaries
+      allocate(self%boundariesParticle(1:nSpecies))
+      do s = 1, nSpecies
+        self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj
+
+      end do
+
     END SUBROUTINE initEdge2DCyl
 
     !Get nodes from edge

src/modules/mesh/3DCart/moduleMesh3DCart.f90

@@ -153,6 +153,13 @@ MODULE moduleMesh3DCart
       call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n2%n)
       call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n3%n)
 
+      ! Link edge to particle boundaries
+      allocate(self%boundariesParticle(1:nSpecies))
+      do s = 1, nSpecies
+        self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj
+
+      end do
+
     END SUBROUTINE initEdge3DCartTria
 
     !Get nodes from surface