Fix an issue with different time steps
Sometimes the division between number of steps and time step was not giving the right results. Nevertheless, this just indicates that the species have to be in separated arrays so that the assignment of particles in list, pushing and scattering can be dealt independently. Thus, this is the first step in creating separate arrays of particles per species.
This commit is contained in:
parent
541a6ff97a
commit
b73c8531e8
3 changed files with 7 additions and 3 deletions
|
|
@ -816,9 +816,8 @@ MODULE moduleMesh
|
|||
|
||||
IF (MOD(t, everyColl) == 0) THEN
|
||||
!Collisions need to be performed in this iteration
|
||||
!$OMP DO SCHEDULE(DYNAMIC) PRIVATE(part_i, part_j, partTemp_i, partTemp_j)
|
||||
!$OMP DO SCHEDULE(DYNAMIC) PRIVATE(part_i, part_j, partTemp_i, partTemp_j, cell)
|
||||
DO e=1, self%numCells
|
||||
|
||||
cell => self%cells(e)%obj
|
||||
|
||||
!TODO: Simplify this, to many sublevels
|
||||
|
|
|
|||
|
|
@ -123,7 +123,12 @@ MODULE moduleSolver
|
|||
END SELECT
|
||||
|
||||
self%pushSpecies = .FALSE.
|
||||
self%every = INT(tauSp/tau)
|
||||
self%every = FLOOR(tauSp/tau)
|
||||
!Correct value if not fulfilled
|
||||
IF (tau*REAL(self%every) < tauSp) THEN
|
||||
self%every = self%every + 1
|
||||
|
||||
END IF
|
||||
|
||||
END SUBROUTINE initPusher
|
||||
|
||||
|
|
|
|||
Loading…
Add table
Add a link
Reference in a new issue