Skip to content
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.

Sign up for GitHub

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

FDS Source: also update BC%II/JJ/KK for particles wherever BC%IIG/JJG… #13686

Merged
merged 1 commit into from
Nov 4, 2024
Merged

FDS Source: also update BC%II/JJ/KK for particles wherever BC%IIG/JJG… #13686

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
14 changes: 8 additions & 6 deletions Source/part.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -529,9 +529,8 @@ SUBROUTINE INSERT_SPRAY_PARTICLES
ELSE
CALL GET_IJK(BC%X,BC%Y,BC%Z,NM,XI,YJ,ZK,II,JJ,KK)
IC = CELL_INDEX(II,JJ,KK)
BC%IIG = II
BC%JJG = JJ
BC%KKG = KK
BC%IIG = II; BC%JJG = JJ; BC%KKG = KK
BC%II = II; BC%JJ = JJ; BC%KK = KK
IF (.NOT.CELL(IC)%SOLID) EXIT CHOOSE_COORDS
ENDIF

Expand Down Expand Up @@ -865,6 +864,7 @@ SUBROUTINE PARTICLE_FACE_INSERT(WALL_INDEX,CFACE_INDEX)

! Update idicies in case offset puts location in a different cell
CALL GET_IJK(BC%X,BC%Y,BC%Z,NM,XI,YJ,ZK,BC%IIG,BC%JJG,BC%KKG)
BC%II = BC%IIG; BC%JJ = BC%JJG; BC%KK = BC%KKG

! Save the insertion time (TP) and scalar property (SP) for the particle

Expand Down Expand Up @@ -1472,9 +1472,8 @@ SUBROUTINE VOLUME_INIT_PARTICLE
INTEGER :: ND
TYPE (PROFILE_TYPE), POINTER :: PF

BC%IIG = II
BC%JJG = JJ
BC%KKG = KK
BC%IIG = II; BC%JJG = JJ; BC%KKG = KK
BC%II = II; BC%JJ = JJ; BC%KK = KK
LP%U = IN%U0
LP%V = IN%V0
LP%W = IN%W0
Expand Down Expand Up @@ -1886,6 +1885,7 @@ SUBROUTINE MOVE_PARTICLES(T,DT,NM)
IF (LP%PATH_PARTICLE) THEN
CALL MOVE_IN_GAS
CALL GET_IJK(BC%X,BC%Y,BC%Z,NM,XI,YJ,ZK,BC%IIG,BC%JJG,BC%KKG)
BC%II = BC%IIG; BC%JJ = BC%JJG; BC%KK = BC%KKG
ENDIF
EXIT TIME_STEP_LOOP
ENDIF
Expand Down Expand Up @@ -1939,6 +1939,7 @@ SUBROUTINE MOVE_PARTICLES(T,DT,NM)
! Determine the cell indices of the new particle location.

CALL GET_IJK(BC%X,BC%Y,BC%Z,NM,XI,YJ,ZK,BC%IIG,BC%JJG,BC%KKG)
BC%II = BC%IIG; BC%JJ = BC%JJG; BC%KK = BC%KKG

! If the particle is not near a boundary cell, cycle.

Expand Down Expand Up @@ -2520,6 +2521,7 @@ SUBROUTINE MOVE_PARTICLES(T,DT,NM)

! Store containing volume at new location
CALL GET_IJK(BC%X,BC%Y,BC%Z,NM,XI,YJ,ZK,BC%IIG,BC%JJG,BC%KKG)
BC%II = BC%IIG; BC%JJ = BC%JJG; BC%KK = BC%KKG
CALL GET_RVC(NM,BC%IIG,BC%JJG,BC%KKG,LP%RVC)

ENDDO TIME_STEP_LOOP
Expand Down