Merge pull request #15053 from rmcdermo/master

FDS Source: fix 2nd ghost cell value of PBAR in wall

Author: rmcdermo

Source/wall.f90

@@ -273,7 +273,7 @@ SUBROUTINE ASSIGN_GHOST_VALUE(WALL_INDEX,BC,B1)
 USE COMPLEX_GEOMETRY, ONLY : CC_CGSC, CC_SOLID
 INTEGER, INTENT(IN) :: WALL_INDEX
 REAL(EB) :: ARO,ZZ_GET(1:N_TRACKED_SPECIES),RHO_OTHER,RHO_ZZ_OTHER(1:N_TOTAL_SCALARS),&
-            RSUM_TMP,RHO_OTHER_2,RHO_ZZ_OTHER_2(1:N_TOTAL_SCALARS)
+            RSUM_TMP,RHO_OTHER_2,RHO_ZZ_OTHER_2(1:N_TOTAL_SCALARS),PBAR_P_2
 INTEGER :: IIO,JJO,KKO,II2,JJ2,KK2,ICG,ICO
 LOGICAL :: CC_SOLID_FLAG,SECOND_ORDER_INTERPOLATED_BOUNDARY
 REAL(EB), POINTER, DIMENSION(:,:,:,:) :: OM_ZZP
@@ -365,7 +365,8 @@ SUBROUTINE ASSIGN_GHOST_VALUE(WALL_INDEX,BC,B1)
    ZZP(BC%II2,BC%JJ2,BC%KK2,1:N_TOTAL_SCALARS) = MAX(0._EB,MIN(1._EB,RHO_ZZ_OTHER_2(1:N_TOTAL_SCALARS)/RHO_OTHER_2))
    ZZ_GET(1:N_TRACKED_SPECIES) = ZZP(BC%II2,BC%JJ2,BC%KK2,1:N_TRACKED_SPECIES)
    CALL GET_SPECIFIC_GAS_CONSTANT(ZZ_GET,RSUM_TMP)
-   TMP(BC%II2,BC%JJ2,BC%KK2) = PBAR_P(BC%KK2,B1%PRESSURE_ZONE)/(RSUM_TMP*RHOP(BC%II2,BC%JJ2,BC%KK2))
+   PBAR_P_2 = 2._EB*PBAR_P(BC%KK,B1%PRESSURE_ZONE) - PBAR_P(BC%KKG,B1%PRESSURE_ZONE)
+   TMP(BC%II2,BC%JJ2,BC%KK2) = PBAR_P_2/(RSUM_TMP*RHOP(BC%II2,BC%JJ2,BC%KK2))
 ELSE
    RHOP(BC%II2,BC%JJ2,BC%KK2) = RHOP(BC%II,BC%JJ,BC%KK)
    ZZP(BC%II2,BC%JJ2,BC%KK2,1:N_TOTAL_SCALARS) = ZZP(BC%II,BC%JJ,BC%KK,1:N_TOTAL_SCALARS)