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