From a35923aa38881fa33ef72971095782fb63fb1f72 Mon Sep 17 00:00:00 2001
From: mcgratta <mcgratta@gmail.com>
Date: Fri, 17 Jan 2025 10:16:54 -0500
Subject: [PATCH] FDS Source: Issue #14072. Fix bug in BURN_AWAY

---
 Source/init.f90 | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/Source/init.f90 b/Source/init.f90
index fdc5dd78bbf..5a8633df791 100644
--- a/Source/init.f90
+++ b/Source/init.f90
@@ -5037,8 +5037,8 @@ SUBROUTINE GET_BOUNDARY_TYPE
 
 SF => SURFACE(WC%SURF_INDEX)
 
-IF (REMOVE .AND. ( (SF%THERMAL_BC_INDEX==THERMALLY_THICK.AND.(SF%VARIABLE_THICKNESS.OR.SF%HT_DIM>1)) &
-                  .OR. SF%PYROLYSIS_MODEL==PYROLYSIS_SPECIFIED ) ) THEN
+USE_OLD_WALL_CELL: IF (REMOVE .AND. ( (SF%THERMAL_BC_INDEX==THERMALLY_THICK.AND.(SF%VARIABLE_THICKNESS.OR.SF%HT_DIM>1)) &
+                      .OR. SF%PYROLYSIS_MODEL==PYROLYSIS_SPECIFIED ) ) THEN
 
    BC => MESHES(NM)%BOUNDARY_COORD(WC%BC_INDEX)
    IIG = BC%IIG
@@ -5059,6 +5059,7 @@ SUBROUTINE GET_BOUNDARY_TYPE
       EWC => MESHES(NM)%EXTERNAL_WALL(CELL(ICG)%WALL_INDEX(IOR))
       NOM = EWC%NOM
       IF (NOM>0) THEN
+         IF (.NOT.ALLOCATED(MESHES(NOM)%WALL)) EXIT USE_OLD_WALL_CELL
          IIO = EWC%IIO_MIN
          JJO = EWC%JJO_MIN
          KKO = EWC%KKO_MIN
@@ -5073,6 +5074,7 @@ SUBROUTINE GET_BOUNDARY_TYPE
    ENDIF
 
    SWAP: IF (IW_OLD>0) THEN
+      IF (MESHES(NOM)%CELL(ICG_OLD)%SOLID) EXIT SWAP
       WC_OLD => MESHES(NOM)%WALL(IW_OLD)
       IF (WC_OLD%OD_INDEX==0) EXIT SWAP
       IF (SF%PYROLYSIS_MODEL==PYROLYSIS_SPECIFIED) THEN
@@ -5096,7 +5098,7 @@ SUBROUTINE GET_BOUNDARY_TYPE
          ENDIF
       ENDIF
    ENDIF SWAP
-ENDIF
+ENDIF USE_OLD_WALL_CELL
 
 END SUBROUTINE GET_BOUNDARY_TYPE