Merge pull request #14762 from marcosvanella/master

FDS Source : SET_CUTCELLS_3D, drop zero length body-plane intersection segments.

Author: marcosvanella

Source/geom.f90

@@ -37,7 +37,7 @@ MODULE COMPLEX_GEOMETRY
 ! Threshold cut-cell volume ratio used to define very small cut-cells, tied to NOADVANCE.
 REAL(EB), PARAMETER :: MIN_VOL_FACTOR   = 5.E-4_EB
 REAL(EB), PARAMETER :: ADIFF_INFO_FACTOR= 1.E-1_EB
-REAL(EB), PARAMETER :: SNAP_DIST_FACTOR = 1.E-5_EB
+REAL(EB), PARAMETER :: SNAP_DIST_FACTOR = 1.E-4_EB
 REAL(EB), PARAMETER :: MIN_LENGTH_FACTOR= 1.E-2_EB
 
 INTEGER,  SAVE ::      NGUARD = 5        ! Layers of guard-cells.
@@ -202,7 +202,7 @@ MODULE COMPLEX_GEOMETRY
 INTEGER,  ALLOCATABLE, DIMENSION(:,:):: CC_IS_CRS2
 REAL(EB), ALLOCATABLE, DIMENSION(:,:):: CC_SEG_TAN
 INTEGER :: X1NOC, X2NOC, X3NOC
-INTEGER, PARAMETER :: MAX_CELL_POLYLINES =100
+INTEGER, PARAMETER :: MAX_CELL_POLYLINES = 200
 
 REAL(EB):: VAL_TESTX_LOW,VAL_TESTX_HIGH,VAL_TESTY_LOW,VAL_TESTY_HIGH,VAL_TESTZ_LOW,VAL_TESTZ_HIGH
 
@@ -11363,6 +11363,10 @@ SUBROUTINE GET_BODINT_PLANE(X1AXIS,X1PLN,INDX1,PLNORMAL,X2AXIS,X3AXIS,&
 ISEG_NEW = 0
 IF(.NOT.TRI_ONPLANE_ONLY) THEN
    DO ISEG=1,BODINT_PLANE%NSEGS
+       SEG(NOD1:NOD2) = BODINT_PLANE%SEGS(NOD1:NOD2,ISEG)
+       XYZ1(IAXIS:KAXIS) = BODINT_PLANE%XYZ(IAXIS:KAXIS,SEG(NOD1))
+       XYZ2(IAXIS:KAXIS) = BODINT_PLANE%XYZ(IAXIS:KAXIS,SEG(NOD2))
+       IF( NORM2(XYZ2((/X2AXIS,X3AXIS/))-XYZ1((/X2AXIS,X3AXIS/))) < 0.1_EB*GEOMEPS) CYCLE
        IF ( (BODINT_PLANE%SEGTYPE(NOD1,ISEG) == CC_SOLID) .AND. &
             (BODINT_PLANE%SEGTYPE(NOD2,ISEG) == CC_SOLID) ) CYCLE