From bbf6d2faa3ad9ee967a75c65ff342bd88ab99233 Mon Sep 17 00:00:00 2001
From: eap <eap@opencascade.com>
Date: Mon, 17 Aug 2020 14:00:07 +0300
Subject: [PATCH] #19895 [CEA] Overlapping volume cells with viscous layer on
 shared faces

---
 src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx b/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
index 64825c0..1f11791 100644
--- a/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
+++ b/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
@@ -263,9 +263,15 @@ static int checkTmpFace(const SMDS_MeshNode* node1,
         bool isDomainInPyramid = ( i3base != i3 );
         return isDomainInPyramid ? HOLE_ID : vol->getshapeId();
       }
-      else
+      else // nbNodes == 6
       {
-        return vol->getshapeId(); // triangle is a prism top
+        int i1 = vol->GetNodeIndex( node1 );
+        if (( i1 == 5 && i2 == 4 && i3 == 3 ) ||
+            ( i1 == 4 && i2 == 3 && i3 == 5 ) ||
+            ( i1 == 3 && i2 == 5 && i3 == 4 ))
+          return HOLE_ID;
+        else
+          return vol->getshapeId(); // triangle is a prism top
       }
     }
   }
-- 
2.39.2