Merge pull request #15437 from mcgratta/master

FDS Source: Fix logic for particle position

Author: mcgratta

Source/part.f90

@@ -1013,9 +1013,9 @@ SUBROUTINE INSERT_VOLUMETRIC_PARTICLES
 ! If the INIT volume is outside the current mesh, return
 
 IF (IN%SHAPE/='RING' .AND. IN%SHAPE/='LINE') THEN
-   IF ((IN_X1-XF)>50._EB*TWO_EPSILON_EB .OR. (IN_X2-XS)<-50._EB*TWO_EPSILON_EB .OR. &
-       (IN_Y1-YF)>50._EB*TWO_EPSILON_EB .OR. (IN_Y2-YS)<-50._EB*TWO_EPSILON_EB .OR. &
-       (IN_Z1-ZF)>50._EB*TWO_EPSILON_EB .OR. (IN_Z2-ZS)<-50._EB*TWO_EPSILON_EB) RETURN
+   IF ((IN_X1-XF)>-50._EB*TWO_EPSILON_EB .OR. (IN_X2-XS)<50._EB*TWO_EPSILON_EB .OR. &
+       (IN_Y1-YF)>-50._EB*TWO_EPSILON_EB .OR. (IN_Y2-YS)<50._EB*TWO_EPSILON_EB .OR. &
+       (IN_Z1-ZF)>-50._EB*TWO_EPSILON_EB .OR. (IN_Z2-ZS)<50._EB*TWO_EPSILON_EB) RETURN
 ELSEIF (IN%SHAPE=='RING') THEN
    IF (RING_MESH_INTERSECTION_ARC(NM,IN%X0,IN%Y0,IN%RADIUS)<TWO_EPSILON_EB) RETURN
 ELSE

Source/read.f90

@@ -12886,12 +12886,12 @@ SUBROUTINE READ_INIT
       MESH_LOOP: DO NM=1,NMESHES
          M=>MESHES(NM)
          IF (XYZ(1)>=M%XS .AND. XYZ(1)<=M%XF .AND. XYZ(2)>=M%YS .AND.  XYZ(2)<=M%YF .AND. XYZ(3)>=M%ZS .AND. XYZ(3)<=M%ZF) THEN
-            IF (ABS(XYZ(1)-M%XS)<TWO_EPSILON_EB) XYZ(1) = XYZ(1) + 0.01_EB*M%DXI
-            IF (ABS(XYZ(1)-M%XF)<TWO_EPSILON_EB) XYZ(1) = XYZ(1) - 0.01_EB*M%DXI
-            IF (ABS(XYZ(2)-M%YS)<TWO_EPSILON_EB) XYZ(2) = XYZ(2) + 0.01_EB*M%DETA
-            IF (ABS(XYZ(2)-M%YF)<TWO_EPSILON_EB) XYZ(2) = XYZ(2) - 0.01_EB*M%DETA
-            IF (ABS(XYZ(3)-M%ZS)<TWO_EPSILON_EB) XYZ(3) = XYZ(3) + 0.01_EB*M%DZETA
-            IF (ABS(XYZ(3)-M%ZF)<TWO_EPSILON_EB) XYZ(3) = XYZ(3) - 0.01_EB*M%DZETA
+            IF (ABS(XYZ(1)-M%XS)<50._EB*TWO_EPSILON_EB) XYZ(1) = XYZ(1) + 0.01_EB*M%DXI
+            IF (ABS(XYZ(1)-M%XF)<50._EB*TWO_EPSILON_EB) XYZ(1) = XYZ(1) - 0.01_EB*M%DXI
+            IF (ABS(XYZ(2)-M%YS)<50._EB*TWO_EPSILON_EB) XYZ(2) = XYZ(2) + 0.01_EB*M%DETA
+            IF (ABS(XYZ(2)-M%YF)<50._EB*TWO_EPSILON_EB) XYZ(2) = XYZ(2) - 0.01_EB*M%DETA
+            IF (ABS(XYZ(3)-M%ZS)<50._EB*TWO_EPSILON_EB) XYZ(3) = XYZ(3) + 0.01_EB*M%DZETA
+            IF (ABS(XYZ(3)-M%ZF)<50._EB*TWO_EPSILON_EB) XYZ(3) = XYZ(3) - 0.01_EB*M%DZETA
             EXIT MESH_LOOP
          ENDIF
       ENDDO MESH_LOOP