Proper output file for particle boundaries. Missing input toggle

src/modules/solver/moduleSolver.f90

@@ -521,10 +521,9 @@ MODULE moduleSolver
       USE moduleCompTime
       USE moduleProbe
       USE moduleCaseParam, ONLY: timeStep
+      use moduleMeshInoutCommon, only: boundariesParticle_print
       IMPLICIT NONE
 
-      integer:: b ! TEMPORARY
-
       CALL outputProbes()
 
       counterOutput = counterOutput + 1
@@ -537,6 +536,10 @@ MODULE moduleSolver
         CALL mesh%printOutput()
         IF (ASSOCIATED(meshForMCC)) CALL meshForMCC%printColl()
         CALL mesh%printEM()
+
+        ! Output of boundaries. TODO: Add an if
+        call boundariesParticle_print()
+
         WRITE(*, "(5X,A21,I10,A1,I10)") "t/tFinal: ", timeStep, "/", tFinal
         WRITE(*, "(5X,A21,I10)")        "Particles: ", nPartOld
         IF (timeStep == 0) THEN
@@ -553,17 +556,6 @@ MODULE moduleSolver
         END IF
         WRITE(*,*)
 
-        ! TEMPORARY
-        ! Output of boundaries. TODO: Move this to mesh module
-        do b = 1, nBoundariesParticle
-          if (associated(boundariesParticle(b)%obj%print)) then
-            print *, b
-            call boundariesParticle(b)%obj%print()
-
-          end if
-
-        end do
-
       END IF
 
       counterCPUTime = counterCPUTime + 1