From 38fa37c9957cc152439bbe9ab58f40b72167ef04 Mon Sep 17 00:00:00 2001 From: JGonzalez Date: Mon, 17 Jul 2023 16:29:52 +0200 Subject: [PATCH] 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 | 4 ++++ src/modules/mesh/moduleMeshBoundary.f90 | 5 ++++- 2 files changed, 8 insertions(+), 1 deletion(-) create mode 100644 data/see/constant.dat diff --git a/data/see/constant.dat b/data/see/constant.dat new file mode 100644 index 0000000..060c8ad --- /dev/null +++ b/data/see/constant.dat @@ -0,0 +1,4 @@ +#Relative energy (eV) yield () +0.000 1.000E-01 +1.000 1.000E-01 + diff --git a/src/modules/mesh/moduleMeshBoundary.f90 b/src/modules/mesh/moduleMeshBoundary.f90 index 34dc8a7..ffd1af9 100644 --- a/src/modules/mesh/moduleMeshBoundary.f90 +++ b/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')