Skip to content

Commit c54f0a1

Browse files
authored
Merge pull request #15270 from mcgratta/master
FDS Source: Fix bug related to particles near boundaries
2 parents 96ca245 + 692cba8 commit c54f0a1

File tree

1 file changed

+3
-3
lines changed

1 file changed

+3
-3
lines changed

Source/part.f90

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1013,9 +1013,9 @@ SUBROUTINE INSERT_VOLUMETRIC_PARTICLES
10131013
! If the INIT volume is outside the current mesh, return
10141014

10151015
IF (IN%SHAPE/='RING' .AND. IN%SHAPE/='LINE') THEN
1016-
IF ((IN_X1-XF)>-50._EB*TWO_EPSILON_EB .OR. (IN_X2-XS)<50._EB*TWO_EPSILON_EB .OR. &
1017-
(IN_Y1-YF)>-50._EB*TWO_EPSILON_EB .OR. (IN_Y2-YS)<50._EB*TWO_EPSILON_EB .OR. &
1018-
(IN_Z1-ZF)>-50._EB*TWO_EPSILON_EB .OR. (IN_Z2-ZS)<50._EB*TWO_EPSILON_EB) RETURN
1016+
IF ((IN_X1-XF)>50._EB*TWO_EPSILON_EB .OR. (IN_X2-XS)<-50._EB*TWO_EPSILON_EB .OR. &
1017+
(IN_Y1-YF)>50._EB*TWO_EPSILON_EB .OR. (IN_Y2-YS)<-50._EB*TWO_EPSILON_EB .OR. &
1018+
(IN_Z1-ZF)>50._EB*TWO_EPSILON_EB .OR. (IN_Z2-ZS)<-50._EB*TWO_EPSILON_EB) RETURN
10191019
ELSEIF (IN%SHAPE=='RING') THEN
10201020
IF (RING_MESH_INTERSECTION_ARC(NM,IN%X0,IN%Y0,IN%RADIUS)<TWO_EPSILON_EB) RETURN
10211021
ELSE

0 commit comments

Comments
 (0)