Corrections and some basic printing

src/modules/mesh/moduleMesh@boundaryParticle.f90

@@ -109,10 +109,9 @@ submodule(moduleMesh) boundaryParticle
 
         boundary%s_incident = s_incident
 
-        allocate(boundary%alpha(1:mesh%numEdges)) ! TODO: Change this so only the edges associated to the boundary are here
-        boundary%alpha = 0.d0
 
         boundary%update => quasiNeutrality_update
+        boundary%print  => quasiNeutrality_print
 
       end select
 
@@ -391,7 +390,7 @@ submodule(moduleMesh) boundaryParticle
           alpha = 1.d0 - density_incident/density_rest
 
           ! Apply correction with a factor of 0.1 to avoid fast changes
-          self%alpha(edge%n) = self%alpha(edge%n) + 0.1d0 * alpha
+          self%alpha(edge%n) = self%alpha(edge%n) + 1.0d-2 * alpha
 
           ! Limit alpha between 0 and 1
           self%alpha(edge%n) = min(self%alpha(edge%n), 1.d0)
@@ -405,6 +404,21 @@ submodule(moduleMesh) boundaryParticle
 
     end subroutine quasiNeutrality_update
 
+    subroutine quasiNeutrality_print(self)
+      implicit none
+
+      class(boundaryParticleGeneric), intent(inout):: self
+
+      print *, 'test'
+
+      select type(self)
+      type is(boundaryQuasiNeutrality)
+        print*, self%alpha
+
+      end select
+
+    end subroutine quasiNeutrality_print
+
     ! Generic boundary conditions for internal use
     module subroutine genericReflection(edge, part)
       use moduleCaseParam