From e825221b0b40ad57a378a00ae7e00eff909ec8ae Mon Sep 17 00:00:00 2001
From: marcosvanella <marcos_vanella@hotmail.com>
Date: Wed, 10 Jan 2024 09:46:41 -0500
Subject: [PATCH] FDS Source : Cycle particle-boundary interaction when
 cut-cell is outside of mesh (no CFACEs).

---
 Source/part.f90 | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/Source/part.f90 b/Source/part.f90
index fe574b50f9f..33842671823 100644
--- a/Source/part.f90
+++ b/Source/part.f90
@@ -1976,6 +1976,10 @@ SUBROUTINE MOVE_PARTICLES(T,DT,NM)
          ! Check if there are CFACEs in this cell:
          INDCF_POS : IF ( INDCF > 0 ) THEN  ! Current grid cell has CFACEs
 
+            ! Drop if cell with boundary cut-faces found is outside mesh, no CFACEs defined.
+            ICF=CELL_INDEX(CUT_FACE(INDCF)%IJK(1),CUT_FACE(INDCF)%IJK(2),CUT_FACE(INDCF)%IJK(3))
+            IF (CELL(ICF)%EXTERIOR) CYCLE PARTICLE_LOOP
+
             DIST2_MIN = 1.E6_EB
             DO IFACE=1,CUT_FACE(INDCF)%NFACE  ! Loop through CFACEs and find the one closest to the particle
                ICF = CUT_FACE(INDCF)%CFACE_INDEX(IFACE)