Merge pull request #14855 from marcosvanella/master

FDS Source: Discard RC faces when bearing cartesian cell is solid.

Author: marcosvanella

Source/ccib.f90

@@ -7444,7 +7444,7 @@ SUBROUTINE INIT_CUTCELL_DATA(T,DT,FIRST_CALL)
    WHERE(FCVAR(0:IBP1,0:JBP1,0:KBP1,CC_FGSC,JAXIS)==CC_SOLID) V(0:IBP1,0:JBP1,0:KBP1) = 0._EB
    WHERE(FCVAR(0:IBP1,0:JBP1,0:KBP1,CC_FGSC,KAXIS)==CC_SOLID) W(0:IBP1,0:JBP1,0:KBP1) = 0._EB
    US = U; VS = V; WS = W
-   
+
    ! External mesh CFACEs initialize P1 BCs:
    DO ICF=1,N_EXTERNAL_CFACE_CELLS+N_INTWALL_CFACE_CELLS
       CFA  => CFACE(ICF)
@@ -21080,6 +21080,7 @@ SUBROUTINE GET_GASPHASE_REGRCFACES_DATA
    DO ICC=1,M%N_CUTCELL_MESH
       CC => CUT_CELL(ICC); IJK(IAXIS:KAXIS) = CC%IJK(IAXIS:KAXIS)
       DO JCC=1,CC%NCELL
+         IF(CC%UNKZ(JCC) < 1) CYCLE
          ! Loop faces and test:
          DO IFC=1,CC%CCELEM(1,JCC)
             IFACE = CC%CCELEM(IFC+1,JCC)
@@ -21372,6 +21373,7 @@ SUBROUTINE GET_GASPHASE_REGRCFACES_DATA
    DO ICC=1,M%N_CUTCELL_MESH
       CC => CUT_CELL(ICC); IJK(IAXIS:KAXIS) = CC%IJK(IAXIS:KAXIS)
       DO JCC=1,CC%NCELL
+         IF(CC%UNKZ(JCC) < 1) CYCLE
          ! Loop faces and test:
          IFC_LOOP : DO IFC=1,CC%CCELEM(1,JCC)
             IFACE = CC%CCELEM(IFC+1,JCC)