Merge pull request #13365 from mcgratta/master

FDS Source: HT3D boundary conditions

Source/wall.f90

@@ -2767,7 +2767,6 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
    ONE_D%TMP(0)     = ONE_D%TMP(0)     + DELTA_TMP(0)
    DELTA_TMP(NWP+1) = MAX(TMPMIN,ONE_D%TMP(NWP)*RFACB2+QDXKB) - ONE_D%TMP(NWP+1)
    ONE_D%TMP(NWP+1) = ONE_D%TMP(NWP+1) + DELTA_TMP(NWP+1)
-
    ONE_D%DELTA_TMP(0:NWP+1) = ONE_D%DELTA_TMP(0:NWP+1) + DELTA_TMP(0:NWP+1)
 
    B1%Q_CON_F = B1%Q_CON_F + HTCF*DT_BC_SUB*(B1%TMP_G-0.5_EB*B1%TMP_F)
@@ -3645,7 +3644,6 @@ SUBROUTINE HT3D_TEMPERATURE_EXCHANGE(NM)
          IF (ABS(BC2%IOR)/=M%HT_3D_SWEEP_DIRECTION) CYCLE WEIGHT_LOOP
          ONE_D%TMP(I) = (ONE_D%RHO_C_S(I)*ONE_D%TMP(I) + &
                          THR_D%NODE(I)%ALTERNATE_WALL_WEIGHT(II)*ONE_D2%RHO_C_S(I_NODE)*ONE_D2%DELTA_TMP(I_NODE))/ONE_D%RHO_C_S(I)
-
       ENDDO WEIGHT_LOOP
       IF (THR_D%NODE(I)%MESH_NUMBER==NM) M%TMP(THR_D%NODE(I)%I,THR_D%NODE(I)%J,THR_D%NODE(I)%K) = ONE_D%TMP(I)
    ENDDO NODE_LOOP
@@ -3667,6 +3665,8 @@ SUBROUTINE HT3D_TEMPERATURE_EXCHANGE(NM)
    THR_D => M%BOUNDARY_THR_D(TW%TD_INDEX)
    IF (.NOT.ALLOCATED(THR_D%NODE)) CYCLE THIN_WALL_LOOP
 
+   TMP_1   = ONE_D%TMP(1)
+   TMP_NWP = ONE_D%TMP(NWP)
    NODE_LOOP_2: DO I=1,NWP  ! Nodes of the chosen wall cell.
       IF (THR_D%NODE(I)%ALTERNATE_WALL_COUNT==0) CYCLE NODE_LOOP_2
       IF (.NOT.ANY(ABS(THR_D%NODE(I)%ALTERNATE_WALL_IOR(:))==M%HT_3D_SWEEP_DIRECTION)) CYCLE NODE_LOOP_2
@@ -3687,6 +3687,8 @@ SUBROUTINE HT3D_TEMPERATURE_EXCHANGE(NM)
          ONE_D%TMP(I) = ONE_D%TMP(I) + THR_D%NODE(I)%ALTERNATE_WALL_WEIGHT(II)*ONE_D2%DELTA_TMP(I_NODE)
       ENDDO WEIGHT_LOOP_2
    ENDDO NODE_LOOP_2
+   ONE_D%TMP(0)     = ONE_D%TMP(0)     + ONE_D%TMP(1)   - TMP_1
+   ONE_D%TMP(NWP+1) = ONE_D%TMP(NWP+1) + ONE_D%TMP(NWP) - TMP_NWP
 
 ENDDO THIN_WALL_LOOP
 !$OMP END DO