Merge pull request #14073 from mcgratta/master

FDS Source: Issue #14072. Fix bug in BURN_AWAY

Author: mcgratta

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