diff --git a/src/modules/mesh/moduleMesh.f90 b/src/modules/mesh/moduleMesh.f90 index 7c3c4c3..c23990a 100644 --- a/src/modules/mesh/moduleMesh.f90 +++ b/src/modules/mesh/moduleMesh.f90 @@ -1016,6 +1016,14 @@ MODULE moduleMesh W = partTemp%part%v - velocity normW = NORM2(W) + IF (normW > 1.D-10) THEN + !If relative velocity is too low, skip collision and move to next particle + partTemp => partTemp%next + + CYCLE + + END IF + lW = l * normW AW = coulombMatrix(k)%A_i/normW @@ -1082,6 +1090,14 @@ MODULE moduleMesh W = partTemp%part%v - velocity normW = NORM2(W) + IF (normW > 1.D-10) THEN + !If relative velocity is too low, skip collision and move to next particle + partTemp => partTemp%next + + CYCLE + + END IF + lW = l * normW AW = coulombMatrix(k)%A_j/normW diff --git a/src/modules/solver/moduleSolver.f90 b/src/modules/solver/moduleSolver.f90 index 12f9e9a..e539fed 100644 --- a/src/modules/solver/moduleSolver.f90 +++ b/src/modules/solver/moduleSolver.f90 @@ -70,7 +70,6 @@ MODULE moduleSolver CHARACTER(:), ALLOCATABLE:: pusherType REAL(8):: tau, tauSp - !TODO: Reorganize if Cart pushers are combined SELECT CASE(mesh%dimen) CASE(0) self%pushParticle => push0D