From cca57d5ab049f6d505ce506d99f42eb69a5f96dd Mon Sep 17 00:00:00 2001 From: JGonzalez Date: Fri, 13 Jan 2023 15:32:25 +0100 Subject: [PATCH] Incorrect tau in electrostatic pushers During the improve performant step, an error in the electrostatic pushers was introduced, resulting in these using the minimum time step and not the species time step when calculating the acceleration. --- src/modules/solver/pusher/modulePusher.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/modules/solver/pusher/modulePusher.f90 b/src/modules/solver/pusher/modulePusher.f90 index bc02912..f75c733 100644 --- a/src/modules/solver/pusher/modulePusher.f90 +++ b/src/modules/solver/pusher/modulePusher.f90 @@ -24,7 +24,7 @@ MODULE modulePusher !Get the electric field at particle position qmEFT = mesh%cells(part%cell)%obj%gatherElectricField(part%Xi) - qmEFt = qmEFt*part%species%qm*tauMin + qmEFt = qmEFt*part%species%qm*tauIn !Update velocity part%v = part%v + qmEFt @@ -127,7 +127,7 @@ MODULE modulePusher part_temp = part !Get electric field at particle position qmEFT = mesh%cells(part_temp%cell)%obj%gatherElectricField(part_temp%Xi) - qmEFt = qmEFt*part_temp%species%qm*tauMin + qmEFt = qmEFt*part_temp%species%qm*tauIn !z part_temp%v(1) = part%v(1) + qmEFt(1) part_temp%r(1) = part%r(1) + part_temp%v(1)*tauIn @@ -203,7 +203,7 @@ MODULE modulePusher part_temp = part !Get electric field at particle position qmEFT = mesh%cells(part_temp%cell)%obj%gatherElectricField(part_temp%Xi) - qmEFt = qmEFt*part_temp%species%qm*tauMin + qmEFt = qmEFt*part_temp%species%qm*tauIn !r,theta v_p_oh_star(1) = part%v(1) + qmEFt(1) x_new = part%r(1) + v_p_oh_star(1)*tauIn