Dummy velocity for testing

I have set the velocity of the particle to be a dummy Maxwellian
distribution to start testing things.

Also I added a constant yield table for testing.

Next step is to do the user input and run some tests.

data/see/constant.dat

@@ -0,0 +1,4 @@
+#Relative energy (eV) yield ()
+0.000 1.000E-01
+1.000 1.000E-01
+

src/modules/mesh/moduleMeshBoundary.f90

@@ -193,7 +193,7 @@ MODULE moduleMeshBoundary
 
       END SELECT
 
-      !Removes ionizing electron regardless the number of pair created
+      !Removes ionizing electron regardless the number of pairs created
       part%n_in = .FALSE.
 
     END SUBROUTINE ionization
@@ -214,6 +214,7 @@ MODULE moduleMeshBoundary
 
     SUBROUTINE secondaryEmission(edge, part)
       USE moduleSpecies
+      USE moduleRandom
       IMPLICIT NONE
 
       CLASS(meshEdge), INTENT(inout):: edge
@@ -270,6 +271,7 @@ MODULE moduleMeshBoundary
           secondaryElectron%n_in = .TRUE.
 
           !Assign velocity
+          secondaryElectron%v = 2.D0 * edge%normal + 1.D-1 * (/ randomMaxwellian(), randomMaxwellian(), randomMaxwellian() /)
 
           !Add particle to list
           CALL partSurfaces%setLock()
@@ -314,6 +316,7 @@ MODULE moduleMeshBoundary
 
       TYPE IS(boundarySEE)
         edge%fBoundary(s)%apply => secondaryEmission
+
       CLASS DEFAULT
         CALL criticalError("Boundary type not defined in this geometry", 'pointBoundaryFunction')