-
Notifications
You must be signed in to change notification settings - Fork 142
bug fix for particle with periodic condition #345
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
LGTM, the next reviewer can squash and merge
src/particle.f90
Outdated
| ! local data | ||
| real(mytype) :: iface | ||
| real(mytype),save :: bcord(6) | ||
| real(mytype),save :: bcord(6),lenpe(6) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
For consistency with the rest of the code this should be
real(mytype), dimension(6), save :: bcord, lenpe
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
done
| lenpe(1)=xlx | ||
| lenpe(2)=-xlx | ||
| lenpe(3)=yly | ||
| lenpe(4)=-yly | ||
| lenpe(5)=zlz | ||
| lenpe(6)=-zlz |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Why is the length being defined as 2xL in each direction?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
perhaps add a comment to explain this
The implemetation of particle periodic condition is wrong. particles would be still out of the domain.
the bug is fixed