Init for edges changed. Now, input

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

@@ -140,7 +140,7 @@ MODULE moduleMesh2DCyl
 
     !EDGE FUNCTIONS
     !Init edge element
-    SUBROUTINE initEdge2DCyl(self, n, p, bt, physicalSurface)
+    SUBROUTINE initEdge2DCyl(self, n, p, bt)
       USE moduleSpecies
       USE moduleErrors
       USE moduleConstParam, ONLY: PI
@@ -149,8 +149,7 @@ MODULE moduleMesh2DCyl
       CLASS(meshEdge2DCyl), INTENT(out):: self
       INTEGER, INTENT(in):: n
       INTEGER, INTENT(in):: p(:)
-      INTEGER, INTENT(in):: bt
-      INTEGER, INTENT(in):: physicalSurface
+      INTEGER, INTENT(in):: bt(1:nsPecies)
       REAL(8), DIMENSION(1:3):: r1, r2
       INTEGER:: s
 
@@ -178,17 +177,13 @@ MODULE moduleMesh2DCyl
                          0.D0 /)
       self%normal = self%normal/NORM2(self%normal)
 
-      !Boundary index
-      self%boundary => boundaries(bt)
-      ALLOCATE(self%fboundary(1:nSpecies))
+      !Boundary particle linking
+      allocate(self%boundariesParticle(1:nSpecies))
       !Assign functions to boundary
-      DO s = 1, nSpecies
-        CALL pointBoundaryFunction(self, s)
+      do s = 1, nSpecies
+        self%boundariesParticle(s)%obj => boundariesParticle(bt(s))%obj
 
-      END DO
-
-      !Physical surface
-      self%physicalSurface = physicalSurface
+      end do
 
     END SUBROUTINE initEdge2DCyl