Calculates a reflection parameter (alpha) per edge based on the ratio of densities between the incident species and the rest

src/modules/init/moduleInput.f90

@@ -826,6 +826,7 @@ MODULE moduleInput
       real(8):: eThreshold !Energy threshold
       integer:: speciesID, electronSecondaryID
       character(:), allocatable:: speciesName, crossSection, electronSecondary
+      integer:: s_incident
 
       ! Read models of particles
       object = 'boundaries.particles'
@@ -900,7 +901,12 @@ MODULE moduleInput
             END IF
 
           case('quasiNeutrality')
-            call initQuasiNeutrality(self)
+            call config%get(object // '.incident', speciesName, found)
+            if (.not. found) call criticalError("Incident species name not found for quasiNeutrality boundary model", 'readBoundary')
+
+            s_incident = speciesName2Index(speciesName)
+
+            call initQuasiNeutrality(self, s_incident)
 
           CASE DEFAULT
             CALL criticalError('Boundary type ' // bType // ' undefined', 'readBoundary')