Fix injection

Half-Maxwellian distribution should inject particles correctly
This commit is contained in:
Jorge Gonzalez 2025-07-27 17:14:38 +02:00
commit 7f73b69dc2
2 changed files with 18 additions and 4 deletions

View file

@ -298,10 +298,7 @@ MODULE moduleInject
REAL(8):: v
v = 0.D0
DO WHILE (v <= 0.D0)
v = self%vTh*randomMaxwellian()
END DO
v = self%vTh*randomHalfMaxwellian()
END FUNCTION randomVelHalfMaxwellian