@@ -4795,7 +4795,7 @@ SUBROUTINE REASSIGN_WALL_CELLS(T,NM)
47954795USE COMP_FUNCTIONS, ONLY : CURRENT_TIME
47964796INTEGER , INTENT (IN ) :: NM
47974797REAL (EB), INTENT (IN ) :: T
4798- INTEGER :: N, I1,I2,J1,J2,K1,K2,I,J,K,IW,ICG,IC,OBST_INDEX,NOM,IIO,JJO,KKO
4798+ INTEGER :: I1,I2,J1,J2,K1,K2,I,J,K,IW,ICG,IC,OBST_INDEX,NOM,IIO,JJO,KKO
47994799REAL (EB) :: TNOW
48004800LOGICAL :: CREATE,REMOVE
48014801TYPE (OBSTRUCTION_TYPE), POINTER :: OB
@@ -4822,17 +4822,18 @@ SUBROUTINE REASSIGN_WALL_CELLS(T,NM)
48224822 IC = CELL_INDEX(BC% II ,BC% JJ ,BC% KK )
48234823 ICG = CELL_INDEX(BC% IIG,BC% JJG,BC% KKG)
48244824 IF (CELL(ICG)% SOLID) CYCLE
4825- OBST_INDEX= 0
4825+ OBST_INDEX = 0
4826+ REMOVE = .FALSE.
4827+ CREATE = .FALSE.
48264828 CALL GET_BOUNDARY_TYPE
48274829 ENDIF
48284830ENDDO
48294831
48304832! Loop over all obstructions in the current mesh and initialize newly exposed or covered wall cell faces
48314833
4832- OBSTRUCTION_LOOP: DO N = 1 ,N_OBST
4834+ OBSTRUCTION_LOOP: DO OBST_INDEX = 1 ,N_OBST
48334835
4834- OB = > OBSTRUCTION(N)
4835- OBST_INDEX= N
4836+ OB = > OBSTRUCTION(OBST_INDEX)
48364837
48374838IF (.NOT. OB% SCHEDULED_FOR_REMOVAL .AND. .NOT. OB% SCHEDULED_FOR_CREATION) CYCLE OBSTRUCTION_LOOP
48384839
@@ -4964,7 +4965,8 @@ SUBROUTINE GET_BOUNDARY_TYPE
49644965WC = > MESHES(NM)% WALL(IW)
49654966BC = > MESHES(NM)% BOUNDARY_COORD(WC% BC_INDEX)
49664967
4967- IF (IW<= N_EXTERNAL_WALL_CELLS .AND. REMOVE) THEN
4968+ IF (IW<= N_EXTERNAL_WALL_CELLS) THEN
4969+ ! IF (IW<=N_EXTERNAL_WALL_CELLS .AND. REMOVE) THEN
49684970 EWC = > MESHES(NM)% EXTERNAL_WALL(IW)
49694971 WC% BOUNDARY_TYPE = SOLID_BOUNDARY
49704972 CELL(IC)% SOLID = .TRUE.
@@ -5023,12 +5025,6 @@ SUBROUTINE GET_BOUNDARY_TYPE
50235025 CALL INIT_WALL_CELL(NM,BC% II,BC% JJ,BC% KK,WC% OBST_INDEX,IW,BC% IOR,WC% SURF_INDEX,IERR,T)
50245026 WC = > MESHES(NM)% WALL(IW)
50255027 IF (IW<= N_EXTERNAL_WALL_CELLS) EWC% PRESSURE_BC_TYPE = PRESSURE_BC_TYPE
5026- ! This code is under construction
5027- ! SF => SURFACE(WC%SURF_INDEX)
5028- ! IF (SF%VARIABLE_THICKNESS .OR. SF%HT_DIM>1) THEN
5029- ! CALL FIND_WALL_BACK_INDEX(NM,IW)
5030- ! CALL REALLOCATE_ONE_D_ARRAYS(NM,WALL_CELL=IW)
5031- ! ENDIF
50325028ENDIF
50335029
50345030! Special cases 1: BURNed_AWAY obstruction exposes a surface that also burns, in which case the surface is to ignite immediately.
0 commit comments