i22datainit.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "subvolumes.inc"
#include "inter22.inc"
#include "task_c.inc"
#include "com01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i22datainit ()

Function/Subroutine Documentation

i22datainit()

subroutine i22datainit

Definition at line 32 of file i22datainit.F.

      USE i22tri_mod
      USE i22bufbric_mod
      USE multimat_param_mod , ONLY : m51_n0phas, m51_nvphas
C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C Interface Type22 (/INTER/TYPE22) is an FSI coupling method based on cut cell method. 
C   This experimental cut cell method is not completed, abandoned, and is not an official option.
C
C=======================================================================   
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
#include      "subvolumes.inc"
#include      "inter22.inc"
#include      "task_c.inc"
#include      "com01_c.inc"
C-----------------------------------------------
 
 
C=======================================================================  
C               D A T A   D E F I N I T I O N S
C=======================================================================  
 
C=======================================================================
C definitions from subvolumes.inc
C=======================================================================
      INTEGER I,NIN,J,K, IFLG
 
      DATA lsectyp/14/
 
      !----------------------------------------------------!
      !  combinations number                               !
      !----------------------------------------------------!
      DATA s_tetra /08/
      DATA s_penta /12/
      DATA s_poly3 /24/
      DATA s_hexae /03/
      DATA s_poly4 /04/
      DATA s_poly4a/06/      
      DATA s_poly4b/06/      
      DATA s_polyc /24/      
      DATA s_polyh(1:8) /08,12,24,03,04,06,06,24/ !generic form                  
      !----------------------------------------------------!
      s22_max = s_tetra  + s_penta  + s_poly3  +  s_hexae  +  s_poly4
     .        + s_poly4a + s_poly4b + s_polyc
 
      !----------------------------------------------------!
      !  MULTIPLICITIES                                    !
      !----------------------------------------------------!
      DATA m_tetra /01/
      DATA m_penta /01/
      DATA m_poly3 /01/
      DATA m_hexae /02/ !-> the same intersection points configuration leads to 2 different (opposite) graphes
      DATA m_poly4 /02/ !-> the same intersection points configuration leads to 2 different (opposite) graphes
      DATA m_poly4a/02/ !-> the same intersection points configuration leads to 2 different (opposite) graphes
      DATA m_poly4b/02/ !-> the same intersection points configuration leads to 2 different (opposite) graphes            
      DATA m_polyc /01/
      DATA m_polyh(1:8) /01,01,01,02,02,02,02,01/ !generic form                  
      
      !----------------------------------------------------!
      !  COMBINATIONS INDEXES  (first index from 1:52)     !
      !----------------------------------------------------!
      DATA d_tetra /01/
      DATA d_penta /09/
      DATA d_poly3 /21/
      DATA d_hexae /45/
      DATA d_poly4 /51/
      DATA d_poly4a/59/
      DATA d_poly4b/71/            
      DATA d_polyc /83/      
      DATA d_polyh(1:8) /01,09,21,45,51,59,71,83/ !generic form                  
 
      !----------------------------------------------------!
      !  main NODES NUMBER                               !
      !----------------------------------------------------!
      DATA n_tetra /01/
      DATA n_penta /02/
      DATA n_poly3 /03/
      DATA n_hexae /04/
      DATA n_poly4 /04/
      DATA n_poly4a/04/
      DATA n_poly4b/04/            
      DATA n_polyc /04/      
      DATA n_polyh(1:8) /01,02,03,04,04,04,04,04/ !generic form            
 
      !----------------------------------------------------!
      !  NUMBER OF FACES                                   !
      !----------------------------------------------------!
      DATA f_tetra /03/
      DATA f_penta /04/
      DATA f_poly3 /05/
      DATA f_hexae /05/
      DATA f_poly4 /06/
      DATA f_poly4a/06/
      DATA f_poly4b/06/            
      DATA f_polyc /05/      
      DATA f_polyh(1:8) /03,04,05,05,06,06,06,05/ !generic form      
 
      !----------------------------------------------------!
      !  NUMBER OF CORNER(EDGE)                            !
      !----------------------------------------------------!
      DATA c_tetra /03/
      DATA c_penta /04/
      DATA c_poly3 /05/
      DATA c_hexae /04/
      DATA c_poly4 /06/
      DATA c_poly4a/06/
      DATA c_poly4b/06/            
      DATA c_polyc /06/      
      DATA c_polyh(1:8) /03,04,05,04,06,06,06,06/ !generic form
 
      !----------------------------------------------------!
      !  SUBVOLUMES COMBINATIONS                           !
      !----------------------------------------------------!
      !TETRA
      DATA bcode(01) /b'111000000000'/    !3584 !1
      DATA bcode(02) /b'000111000000'/    !0448 !2
      DATA bcode(03) /b'000000111000'/    !0056 !3
      DATA bcode(04) /b'000000000111'/    !0007 !4
      DATA bcode(05) /b'100100100000'/    !2336 !5
      DATA bcode(06) /b'010010000100'/    !1156 !6  
      DATA bcode(07) /b'001000010010'/    !0530 !7
      DATA bcode(08) /b'000001001001'/    !0073 !8
      !----------------------------------------------------!
      !PENTA
      DATA bcode(09) /b'011100100000'/    !1824 !9
      DATA bcode(10) /b'101010000100'/    !2692 !10
      DATA bcode(11) /b'110000010010'/    !3090 !11
      DATA bcode(12) /b'100011100000'/    !2272 !12
      DATA bcode(13) /b'010101000100'/    !1348 !13
      DATA bcode(14) /b'000110001001'/    !0393 !14
      DATA bcode(15) /b'100100011000'/    !2328 !15
      DATA bcode(16) /b'001000101010'/    !0554 !16
      DATA bcode(17) /b'000001110001'/    !0113 !17
      DATA bcode(18) /b'010010000011'/    !1155 !18
      DATA bcode(19) /b'001000010101'/    !0533 !19
      DATA bcode(20) /b'000001001110'/    !0078 !20
      !----------------------------------------------------!
      !POLY3
      DATA bcode(21) /b'101101000100'/    !2884 !21
      DATA bcode(22) /b'001110100100'/    !0932 !22
      DATA bcode(23) /b'011011100000'/    !1760 !23
      DATA bcode(24) /b'110001100100'/    !3172 !24
      DATA bcode(25) /b'010101000011'/    !1347 !25
      DATA bcode(26) /b'010011001010'/    !1226 !26
      DATA bcode(27) /b'010100001101'/    !1293 !27
      DATA bcode(28) /b'000110001110'/    !0398 !28
      DATA bcode(29) /b'001001100011'/    !0611 !29
      DATA bcode(30) /b'000001110110'/    !0118 !30
      DATA bcode(31) /b'001001011100'/    !0604 !31
      DATA bcode(32) /b'001000101101'/    !0557 !32
      DATA bcode(33) /b'011100011000'/    !1816 !33
      DATA bcode(34) /b'101100001010'/    !2826 !34
      DATA bcode(35) /b'010100110010'/    !1330 !35
      DATA bcode(36) /b'110000101010'/    !3114 !36
      DATA bcode(37) /b'100010101001'/    !2217 !37
      DATA bcode(38) /b'000110110001'/    !0433 !38
      DATA bcode(39) /b'100011011000'/    !2264 !39
      DATA bcode(40) /b'100101010001'/    !2385 !40
      DATA bcode(41) /b'110000010101'/    !3093 !41
      DATA bcode(42) /b'100010010110'/    !2198 !42
      DATA bcode(43) /b'011010010001'/    !1681 !43
      DATA bcode(44) /b'101010000011'/    !2691 !44
      !----------------------------------------------------!
      !HEXAE
      DATA bcode(45) /b'100010010001'/    !2193 !45-46
      DATA bcode(46) /b'100010010001'/    !2193      
      DATA bcode(47) /b'010100001010'/    !1290 !47-48
      DATA bcode(48) /b'010100001010'/    !1290      
      DATA bcode(49) /b'001001100100'/    !0612 !49-50
      DATA bcode(50) /b'001001100100'/    !0612      
      !----------------------------------------------------!
      !POLY4
      DATA bcode(51) /b'011011011000'/    !1752 !51-52
      DATA bcode(52) /b'011011011000'/    !1752      
      DATA bcode(53) /b'101101000011'/    !2883 !53-54
      DATA bcode(54) /b'101101000011'/    !2883      
      DATA bcode(55) /b'110000101101'/    !3117 !55-56
      DATA bcode(56) /b'110000101101'/    !3117 
      DATA bcode(57) /b'000110110110'/    !0438 !57-58     
      DATA bcode(58) /b'000110110110'/    !0438
      !----------------------------------------------------!
      !POLY4A
      DATA bcode(59) /b'101011001010'/    !POLY4A 45 - 2762
      DATA bcode(60) /b'101011001010'/    !POLY4A 26 - 2762      
      DATA bcode(61) /b'010100110101'/    !POLY4A 46 - 1333 
      DATA bcode(62) /b'010100110101'/    !POLY4A 25 - 1333    
      DATA bcode(63) /b'011010101001'/    !POLY4A 36 - 1705
      DATA bcode(64) /b'011010101001'/    !POLY4A 15 - 1705
      DATA bcode(65) /b'100101010110'/    !POLY4A 35 - 2390
      DATA bcode(66) /b'100101010110'/    !POLY4A 16 - 2390
      DATA bcode(67) /b'110001100011'/    !POLY4A 34 - 3171
      DATA bcode(68) /b'110001100011'/    !POLY4A 12 - 3171
      DATA bcode(69) /b'001110011100'/    !POLY4A 23 - 0924
      DATA bcode(70) /b'001110011100'/    !POLY4A 14 - 0924
      !POLY4B
      DATA bcode(71) /b'011101010001'/    !POLY4B 45 - 1873
      DATA bcode(72) /b'011101010001'/    !POLY4B 26 - 1873
      DATA bcode(73) /b'100010101110'/    !POLY4B 46 - 2222
      DATA bcode(74) /b'100010101110'/    !POLY4B 25 - 2222
      DATA bcode(75) /b'110001011100'/    !POLY4B 36 - 3148
      DATA bcode(76) /b'110001011100'/    !POLY4B 15 - 3148
      DATA bcode(77) /b'001110100011'/    !POLY4B 35 - 0931 
      DATA bcode(78) /b'001110100011'/    !POLY4B 16 - 0931
      DATA bcode(79) /b'101100001101'/    !POLY4B 34 - 2829
      DATA bCODE(80) /B'101100001101'/    !POLY4B 12 - 2829
      DATA bCODE(81) /B'010011110010'/    !POLY4B 23 - 1266
      DATA bCODE(82) /B'010011110010'/    !POLY4B 14 - 1266
      !----------------------------------------------------!
      !POLYC
      DATA bCODE(83)  /B'100110010001'/    !POLYC R+.4
      DATA bCODE(84)  /B'100010110001'/    !POLYC R+.7   
      DATA bCODE(85)  /B'100010011001'/    !POLYC R+.9
      DATA bCODE(86)  /B'100011010001'/    !POLYC R+.6   
                        
      DATA bCODE(87)  /B'110010010001'/    !POLYC R-.2
      DATA bCODE(88)  /B'101010010001'/    !POLYC R-.3 
      DATA bCODE(89)  /B'100010010011'/    !POLYC R-.11
      DATA bCODE(90)  /B'100010010101'/    !POLYC R-.10
                        
      DATA bCODE(91)  /B'110100001010'/    !POLYC S+.1
      DATA bCODE(92)  /B'011100001010'/    !POLYC S+.3
      DATA bCODE(93)  /B'010100011010'/    !POLYC S+.8
      DATA bCODE(94)  /B'010100101010'/    !POLYC S+.7
                      
      DATA bCODE(95)  /B'010110001010'/    !POLYC S-.5
      DATA bCODE(96)  /B'010100001110'/    !POLYC S-.10
      DATA bCODE(97)  /B'010100001011'/    !POLYC S-.12
      DATA bCODE(98)  /B'010101001010'/    !POLYC S-.6
                        
      DATA bCODE(99)  /B'001011100100'/    !POLYC T+.5
      DATA bCODE(100) /B'001101100100'/    !POLYC T+.4
      DATA bCODE(101) /B'101001100100'/    !POLYC T+.1
      DATA bCODE(102) /B'011001100100'/    !POLYC T+.2
                        
      DATA bCODE(103) /B'001001100101'/    !POLYC T-.12
      DATA bCODE(104) /B'001001101100'/    !POLYC T-.9
      DATA bCODE(105) /B'001001110100'/    !POLYC T-.8
      DATA bCODE(106) /B'001001100110'/    !POLYC T-.11
      !----------------------------------------------------!
 
      !----------------------------------------------------!
      !  RETRIEVING POLYHEDRA FAMILY (1-5) FROM A 1-58 CODE!
      !----------------------------------------------------!
       GetPolyhedraType(01:08)  = 1 !Type1=TETRA
       GetPolyhedraType(09:20)  = 2 !Type2=PENTA       
       GetPolyhedraType(21:44)  = 3 !Type2=POLY3
       GetPolyhedraType(45:50)  = 4 !Type2=HEXAE
       GetPolyhedraType(51:58)  = 5 !Type2=POLY4 
       GetPolyhedraType(59:70)  = 5 !Type2=POLY4A 
       GetPolyhedraType(71:82)  = 5 !Type2=POLY4B              
       GetPolyhedraType(83:106) = 6 !Type2=POLYC
 
       GetNumTria(1)           = 1 !Type1=TETRA
       GetNumTria(2)           = 2 !Type2=PENTA       
       GetNumTria(3)           = 3 !Type2=POLY3
       GetNumTria(4)           = 2 !Type2=HEXAE
       GetNumTria(5)           = 4 !Type2=POLY4
       GetNumTria(6)           = 4 !Type2=POLY4A
       GetNumTria(7)           = 4 !Type2=POLY4B             
       GetNumTria(8)           = 4 !Type2=POLYC       
       
       Gtria(:, 1, 1) = (/1,2,3/)  !TETRA (1 triangle )
 
       Gtria(:, 1, 2) = (/1,2,3/)  !PENTA (2 triangles)
       Gtria(:, 2, 2) = (/1,3,4/)       
 
       Gtria(:, 1, 3) = (/2,3,4/)  !POLY3 (3 triangles)
       Gtria(:, 2, 3) = (/1,2,4/)                      
       Gtria(:, 3, 3) = (/1,4,5/)                             
 
       Gtria(:, 1, 4) = (/1,2,3/)  !HEXAE (2 triangles)
       Gtria(:, 2, 4) = (/1,3,4/) 
       
       Gtria(:, 1, 5) = (/1,6,5/)  !POLY4 (4 triangles)
       Gtria(:, 2, 5) = (/1,5,4/)                      
       Gtria(:, 3, 5) = (/1,3,4/)                             
       Gtria(:, 4, 5) = (/1,2,3/)
 
       Gtria(:, 1, 6) = (/2,3,4/)  !POLY4A (4 triangles)
       Gtria(:, 2, 6) = (/2,4,5/)                      
       Gtria(:, 3, 6) = (/2,5,1/)                             
       Gtria(:, 4, 6) = (/1,5,6/)
       
       Gtria(:, 1, 7) = (/2,3,4/)  !POLY4B (4 triangles)
       Gtria(:, 2, 7) = (/2,4,5/)                      
       Gtria(:, 3, 7) = (/2,5,1/)                             
       Gtria(:, 4, 7) = (/1,5,6/)
              
       Gtria(:, 1, 8) = (/1,5,4/)  !POLYC (4 triangles)
       Gtria(:, 2, 8) = (/1,5,6/)                      
       Gtria(:, 3, 8) = (/1,6,2/)                             
       Gtria(:, 4, 8) = (/2,3,6/)                                 
       
      !----------------------------------------------------!
      !  COMBINATIONS NAMES                                !
      !----------------------------------------------------!
      !TETRA
      StrCODE( 1) ='tetra 1       '
      StrCODE( 2) ='tetra 3       '
      StrCODE( 3) ='tetra 6       '
      StrCODE( 4) ='tetra 8       '
      StrCODE( 5) ='tetra 2       '
      StrCODE( 6) ='tetra 4       '
      StrCODE( 7) ='tetra 5       '
      StrCODE( 8) ='tetra 7       '
      !----------------------------------------------------!
      !PENTA
      StrCODE( 9) ='penta 1-2     '
      StrCODE(10) ='penta 1-4     '
      StrCODE(11) ='penta 1-5     '
      StrCODE(12) ='penta 2-3     '
      StrCODE(13) ='penta 3-4     '
      StrCODE(14) ='penta 3-7     '
      StrCODE(15) ='penta 2-6     '
      StrCODE(16) ='penta 5-6     '
      StrCODE(17) ='penta 6-7     '
      StrCODE(18) ='penta 4-8     '
      StrCODE(19) ='penta 5-8     '
      StrCODE(20) ='penta 7-8     '
      !----------------------------------------------------!
      !POLY3
      StrCODE(21) ='poly3 1-4-3   '
      StrCODE(22) ='poly3 2-1-4   '
      StrCODE(23) ='poly3 3-2-1   '
      StrCODE(24) ='poly3 4-3-2   '
      StrCODE(25) ='poly3 3-4-8   '
      StrCODE(26) ='poly3 4-8-7   '
      StrCODE(27) ='poly3 7-3-4   '
      StrCODE(28) ='poly3 8-7-3   '
      StrCODE(29) ='poly3 5-6-7   '
      StrCODE(30) ='poly3 6-7-8   '
      StrCODE(31) ='poly3 7-8-5   '
      StrCODE(32) ='poly3 8-5-6   '
      StrCODE(33) ='poly3 1-2-6   '
      StrCODE(34) ='poly3 2-6-5   '
      StrCODE(35) ='poly3 5-1-2   '
      StrCODE(36) ='poly3 6-5-1   '
      StrCODE(37) ='poly3 2-3-7   '
      StrCODE(38) ='poly3 3-7-6   '
      StrCODE(39) ='poly3 6-2-3   '
      StrCODE(40) ='poly3 7-6-2   '
      StrCODE(41) ='poly3 1-5-8   '
      StrCODE(42) ='poly3 4-1-5   '
      StrCODE(43) ='poly3 5-8-4   '
      StrCODE(44) ='poly3 8-4-1   '
      !----------------------------------------------------!
      !HEXAE
      StrCODE(45) ='hexae r+      '
      StrCODE(46) ='hexae r-      '      
      StrCODE(47) ='hexae s+      '
      StrCODE(48) ='hexae s-      '      
      StrCODE(49) ='hexae t+      '
      StrCODE(50) ='hexae t-      '      
      !----------------------------------------------------!
      !POLY4
      StrCODE(51) ='poly4 1-2-6-3 ' 
      StrCODE(52) ='poly4 4-8-7-5 '      
      StrCODE(53) ='poly4 3-4-8-1 ' 
      StrCODE(54) ='poly4 5-6-7-5 ' 
      StrCODE(55) ='poly4 4-3-2-7 '           
      StrCODE(56) ='poly4 2-3-7-4 '      
      StrCODE(57) ='poly4 4-1-2-5 ' 
      StrCODE(58) ='poly4 8-7-3-6 ' 
      !----------------------------------------------------!
      !POLY4A
      StrCODE(59) ='poly4a 4-5    ' 
      StrCODE(60) ='poly4a 2-6    '      
      StrCODE(61) ='poly4a 4-6    ' 
      StrCODE(62) ='poly4a 2-5    ' 
      StrCODE(63) ='poly4a 3-6    '           
      StrCODE(64) ='poly4a 1-5    '      
      StrCODE(65) ='poly4a 3-5    ' 
      StrCODE(66) ='poly4a 1-6    ' 
      StrCODE(67) ='poly4a 3-4    '           
      StrCODE(68) ='poly4a 1-2    '      
      StrCODE(69) ='poly4a 2-3    ' 
      StrCODE(70) ='poly4a 1-4    '   
      !----------------------------------------------------!
      !POLY4B
      StrCODE(71) ='poly4b 4-5    ' 
      StrCODE(72) ='poly4b 2-6    '      
      StrCODE(73) ='poly4b 4-6    ' 
      StrCODE(74) ='poly4b 2-5    ' 
      StrCODE(75) ='poly4b 3-6    '           
      StrCODE(76) ='poly4b 1-5    '      
      StrCODE(77) ='poly4b 3-5    ' 
      StrCODE(78) ='poly4b 1-6    ' 
      StrCODE(79) ='poly4b 3-4    '           
      StrCODE(80) ='poly4b 1-2    '      
      StrCODE(81) ='poly4b 2-3    ' 
      StrCODE(82) ='poly4b 1-4    '              
      !----------------------------------------------------!
      !POLYC
      StrCODE(83)  ='polyc r+.04   '
      StrCODE(84)  ='polyc r+.07   '
      StrCODE(85)  ='polyc r+.09   '
      StrCODE(86)  ='polyc r+.06   '
      StrCODE(87)  ='polyc r-.02   ' 
      StrCODE(88)  ='polyc r-.03   '       
      StrCODE(89)  ='polyc r-.11   '
      StrCODE(90)  ='polyc r-.10   ' 
      StrCODE(91)  ='polyc s+.01   '
      StrCODE(92)  ='polyc s+.03   '
      StrCODE(93)  ='polyc s+.08   '
      StrCODE(94)  ='polyc s+.07   ' 
      StrCODE(95)  ='polyc s-.05   '
      StrCODE(96)  ='polyc s-.10   '
      StrCODE(97)  ='polyc s-.12   ' 
      StrCODE(98)  ='polyc s-.06   '  
      StrCODE(99)  ='polyc t+.05   '
      StrCODE(100) ='polyc t+.04   '
      StrCODE(101) ='polyc t+.01   '
      StrCODE(102) ='polyc t+.02   ' 
      StrCODE(103) ='polyc t-.12   '
      StrCODE(104) ='polyc t-.09   '
      StrCODE(105) ='polyc t-.08   ' 
      StrCODE(106) ='polyc t-.11   '                   
      
      !----------------------------------------------------!
      !  SUBVOLUME 'g'RAPH : 'f'ACES                       !
      !----------------------------------------------------!
      !TETRA
      Gface(:, 1)=(/1,6,4,0,0,0/)
      Gface(:, 2)=(/2,1,5,0,0,0/)
      Gface(:, 3)=(/4,3,5,0,0,0/)
      Gface(:, 4)=(/2,3,6,0,0,0/)
      Gface(:, 5)=(/5,1,4,0,0,0/)
      Gface(:, 6)=(/6,1,2,0,0,0/)
      Gface(:, 7)=(/6,3,4,0,0,0/)
      Gface(:, 8)=(/5,3,2,0,0,0/)
      !----------------------------------------------------!
      !PENTA
      Gface(:, 9)=(/6,5,4,1,0,0/)
      Gface(:,10)=(/4,2,1,6,0,0/)
      Gface(:,11)=(/1,3,6,4,0,0/)
      Gface(:,12)=(/4,2,5,1,0,0/)
      Gface(:,13)=(/5,6,2,1,0,0/)
      Gface(:,14)=(/1,3,5,2,0,0/)
      Gface(:,15)=(/1,3,4,5,0,0/)
      Gface(:,16)=(/6,5,3,4,0,0/)
      Gface(:,17)=(/4,2,3,5,0,0/)
      Gface(:,18)=(/1,3,2,6,0,0/)
      Gface(:,19)=(/4,2,6,3,0,0/)
      Gface(:,20)=(/5,6,3,2,0,0/)
      !----------------------------------------------------!
      !POLY3
      Gface(:,21)=(/1,5,2,6,4,0/)
      Gface(:,22)=(/1,2,6,4,5,0/)
      Gface(:,23)=(/1,6,4,5,2,0/)
      Gface(:,24)=(/1,4,5,2,6,0/)
      Gface(:,25)=(/2,3,6,1,5,0/)
      Gface(:,26)=(/2,5,3,6,1,0/)
      Gface(:,27)=(/2,6,1,5,3,0/)
      Gface(:,28)=(/2,1,5,3,6,0/)
      Gface(:,29)=(/3,2,5,4,6,0/)
      Gface(:,30)=(/3,6,2,5,4,0/)
      Gface(:,31)=(/3,4,6,2,5,0/)
      Gface(:,32)=(/3,5,4,6,2,0/)
      Gface(:,33)=(/4,3,5,1,6,0/)
      Gface(:,34)=(/4,6,3,5,1,0/)
      Gface(:,35)=(/4,5,1,6,3,0/)
      Gface(:,36)=(/4,1,6,3,5,0/)
      Gface(:,37)=(/5,3,2,1,4,0/)
      Gface(:,38)=(/5,4,3,2,1,0/)
      Gface(:,39)=(/5,2,1,4,3,0/)
      Gface(:,40)=(/5,1,4,3,2,0/)
      Gface(:,41)=(/6,2,3,4,1,0/)
      Gface(:,42)=(/6,3,4,1,2,0/)
      Gface(:,43)=(/6,1,2,3,4,0/)
      Gface(:,44)=(/6,4,1,2,3,0/)
      !----------------------------------------------------!
      !HEXA
      Gface(:,45)=(/5,3,2,1,4,0/)
      Gface(:,46)=(/6,3,4,1,2,0/)
      Gface(:,47)=(/4,3,5,1,6,0/)
      Gface(:,48)=(/2,6,1,5,3,0/)
      Gface(:,49)=(/3,4,6,2,5,0/)
      Gface(:,50)=(/1,4,5,2,6,0/)
      !----------------------------------------------------!
      !POLY4
      Gface(:,51)=(/6,5,1,4,2,3/)
      Gface(:,52)=(/1,3,6,2,4,5/)     
      Gface(:,53)=(/5,6,1,2,4,3/)
      Gface(:,54)=(/6,5,4,3,1,2/)
      Gface(:,55)=(/6,5,2,1,3,4/)
      Gface(:,56)=(/4,2,1,5,6,3/)     
      Gface(:,57)=(/2,4,1,6,5,3/)
      Gface(:,58)=(/6,5,3,2,4,1/)
      !----------------------------------------------------!
      !POLY4A
      Gface(:,59)=(/6,3,5,2,1,4/)
      Gface(:,60)=(/4,1,2,5,3,6/)     
      Gface(:,61)=(/5,1,6,2,3,4/)
      Gface(:,62)=(/4,3,2,6,1,5/)
      Gface(:,63)=(/5,4,6,1,2,3/)
      Gface(:,64)=(/6,4,5,3,2,1/)     
      Gface(:,65)=(/1,4,3,6,2,5/)
      Gface(:,66)=(/3,4,1,5,2,6/)
      Gface(:,67)=(/1,6,3,2,5,4/)
      Gface(:,68)=(/4,5,2,3,5,1/)     
      Gface(:,69)=(/4,5,2,1,5,3/)
      Gface(:,70)=(/2,6,4,3,5,1/)      
      !----------------------------------------------------!
      !POLY4B
      Gface(:,71)=(/6,1,5,2,3,4/)
      Gface(:,72)=(/5,1,6,4,3,2/)     
      Gface(:,73)=(/2,1,4,5,3,6/)
      Gface(:,74)=(/4,1,2,6,3,5/)
      Gface(:,75)=(/1,4,3,5,2,6/)
      Gface(:,76)=(/3,4,1,6,2,5/)     
      Gface(:,77)=(/6,4,5,1,2,3/)
      Gface(:,78)=(/5,4,6,3,2,1/)
      Gface(:,79)=(/1,5,3,2,6,4/)
      Gface(:,80)=(/4,6,2,3,5,1/)     
      Gface(:,81)=(/1,6,3,4,5,2/)
      Gface(:,82)=(/3,6,1,2,5,4/)      
      !----------------------------------------------------!
      !POLYC
      Gface(:,83) =(/5,3,2,4,1,0/)     
      Gface(:,84) =(/5,2,1,3,4,0/)     
      Gface(:,85) =(/5,1,4,2,3,0/)     
      Gface(:,86) =(/5,4,3,1,2,0/)                       
      Gface(:,87) =(/6,3,2,4,1,0/)     
      Gface(:,88) =(/6,2,1,3,4,0/)     
      Gface(:,89) =(/6,1,4,2,3,0/)     
      Gface(:,90) =(/6,4,3,1,2,0/)      
      Gface(:,91) =(/4,3,5,6,1,0/)     
      Gface(:,92) =(/4,5,1,3,6,0/)     
      Gface(:,93) =(/4,1,6,5,3,0/)     
      Gface(:,94) =(/4,6,3,1,5,0/) 
      Gface(:,95) =(/2,3,5,6,1,0/)     
      Gface(:,96) =(/2,5,1,3,6,0/)     
      Gface(:,97) =(/2,1,6,5,3,0/)     
      Gface(:,98) =(/2,6,3,1,5,0/) 
      Gface(:,99) =(/1,4,5,6,2,0/)     
      Gface(:,100)=(/1,5,2,4,6,0/)      
      Gface(:,101)=(/1,2,6,5,4,0/)      
      Gface(:,102)=(/1,6,4,2,5,0/)  
      Gface(:,103)=(/3,4,5,6,2,0/)      
      Gface(:,104)=(/3,6,4,2,5,0/)      
      Gface(:,105)=(/3,2,6,5,4,0/)      
      Gface(:,106)=(/3,5,2,4,6,0/)                   
      !--------------------------- -------------------------!
      !  SUBVOLUME 'g'RAPH : 'n'NUMBERR of 'pt'POINT        !  !enables to have number of point on local face of a given polyhedron. Check graph index1-6 corresponds to Fa,Fb,Fc,Fd,Fe,Ff
      !----------------------------------------------------!  
      DO I=01,08 ;GNPt(:,I)  = (/3,3,3,0,0,0/);ENDDO  
      DO I=09,20 ;GNPt(:,I)  = (/3,3,4,4,0,0/);ENDDO  
      DO I=21,44 ;GNPt(:,I)  = (/5,3,4,4,3,0/);ENDDO  
      DO I=45,50 ;GNPt(:,I)  = (/4,4,4,4,4,0/);ENDDO  
      DO I=51,58 ;GNPt(:,I)  = (/3,5,5,5,3,3/);ENDDO  
      DO I=59,70 ;GNPt(:,I)  = (/3,4,5,3,4,5/);ENDDO  
      DO I=71,82 ;GNPt(:,I)  = (/3,4,5,3,4,5/);ENDDO                
      DO I=83,106;GNPt(:,I)  = (/6,4,4,4,3,3/);ENDDO  
      
      DO I=01,08;GNnod(:,I)  = (/1,1,1,0,0,0/);ENDDO  
      DO I=09,20;GNnod(:,I)  = (/1,1,2,2,0,0/);ENDDO  
      DO I=21,44;GNnod(:,I)  = (/3,1,2,2,1,0/);ENDDO  
      DO I=45,50;GNnod(:,I)  = (/4,2,2,2,2,0/);ENDDO  
      DO I=51,58;GNnod(:,I)  = (/1,3,3,3,1,1/);ENDDO  
      DO I=59,70;GNnod(:,I)  = (/1,2,3,1,2,3/);ENDDO  
      DO I=71,82;GNnod(:,I)  = (/1,2,3,1,2,3/);ENDDO  
      DO I=83,106;GNnod(:,I) = (/4,2,2,2,2,0/);ENDDO                                  
      !----------------------------------------------------!
      !  SUBVOLUME 'g'RAPH : 'c'ORNER                      !  !EDGE id => intersection node id (CORNER).
      !----------------------------------------------------!
      !TETRA
      Gcorner(:, 1)=(/1,2,3,1,0,0,0/)      * -1
      Gcorner(:, 2)=(/6,5,4,6,0,0,0/)      * -1
      Gcorner(:, 3)=(/7,8,9,7,0,0,0/)      * -1
      Gcorner(:, 4)=(/10,12,11,10,0,0,0/)  * -1
      Gcorner(:, 5)=(/7,4,1,7,0,0,0/)
      Gcorner(:, 6)=(/10,2,5,10,0,0,0/)
      Gcorner(:, 7)=(/3,11,8,3,0,0,0/)
      Gcorner(:, 8)=(/6,9,12,6,0,0,0/)
      !----------------------------------------------------!
      !PENTA
      Gcorner(:, 9)=(/-2,-3,+7,+4,-2,0,0/)
      Gcorner(:,10)=(/-3,-1,+5,+10,-3,0,0/)
      Gcorner(:,11)=(/-1,-2,+11,+8,-1,0,0/)
      Gcorner(:,12)=(/+1,+7,-6,-5,+1,0,0/)
      Gcorner(:,13)=(/-4,-6,+10,+2,-4,0,0/)
      Gcorner(:,14)=(/-5,-4,+9,+12,-5,0,0/)
      Gcorner(:,15)=(/+4,+1,-8,-9,+4,0,0/)
      Gcorner(:,16)=(/+3,+11,-9,-7,+3,0,0/)
      Gcorner(:,17)=(/-7,-8,+12,+6,-7,0,0/)
      Gcorner(:,18)=(/+2,+5,-12,-11,+2,0,0/)
      Gcorner(:,19)=(/+8,+3,-10,-12,+8,0,0/)
      Gcorner(:,20)=(/+6,+9,-11,-10,+6,0,0/)
      !----------------------------------------------------!
      !POLY3
      Gcorner(:,21)=(/-3,-1,-4,-6,+10,-3,0/)
      Gcorner(:,22)=(/+7,+4,+5,+10,-3,+7,0/)
      Gcorner(:,23)=(/-6,-5,-2,-3,+7,-6,0/)
      Gcorner(:,24)=(/+10,+2,+1,+7,-6,+10,0/)
      Gcorner(:,25)=(/-4,-6,-12,-11,+2,-4,0/)
      Gcorner(:,26)=(/+2,+5,+6,+9,-11,+2,0/)
      Gcorner(:,27)=(/+9,+12,+10,+2,-4,+9,0/)
      Gcorner(:,28)=(/-11,-10,-5,-4,+9,-11,0/)
      Gcorner(:,29)=(/+3,+11,+12,+6,-7,+3,0/)
      Gcorner(:,30)=(/-7,-8,-11,-10,+6,-7,0/)
      Gcorner(:,31)=(/+6,+9,+8,+3,-10,+6,0/)
      Gcorner(:,32)=(/-10,-12,-9,-7,+3,-10,0/)
      Gcorner(:,33)=(/-2,-3,-8,-9,+4,-2,0/)
      Gcorner(:,34)=(/+4,+1,+3,+11,-9,+4,0/)
      Gcorner(:,35)=(/+11,8,+7,+4,-2,+11,0/)
      Gcorner(:,36)=(/-9,-7,-1,-2,+11,-9,0/)
      Gcorner(:,37)=(/+1,+7,+9,+12,-5,+1,0/)
      Gcorner(:,38)=(/-5,-4,-7,-8,+12,-5,0/)
      Gcorner(:,39)=(/-8,-9,-6,-5,+1,-8,0/)
      Gcorner(:,40)=(/+12,+6,+4,+1,-8,+12,0/)
      Gcorner(:,41)=(/-1,-2,-10,-12,+8,-1,0/)
      Gcorner(:,42)=(/+5,+10,+11,+8,-1,+5,0/)
      Gcorner(:,43)=(/+8,+3,+2,+5,-12,+8,0/)
      Gcorner(:,44)=(/-12,-11,-3,-1,+5,-12,0/)
      !----------------------------------------------------!
      !HEXAE
      Gcorner(:,45)=(/-8 ,+12,-5,+1,-8 ,0 ,0/)
      Gcorner(:,46)=(/-12,+8 ,-1,+5,-12,0 ,0/)
      Gcorner(:,47)=(/+11,-9,+4,-2,+11,0,0/)
      Gcorner(:,48)=(/-11,+2,-4,+9,-11,0,0/)
      Gcorner(:,49)=(/-7,+3,-10,+6,-7,0,0/)
      Gcorner(:,50)=(/-3,+7,-6,+10,-3,0,0/)
      !----------------------------------------------------!
      !POLY4
      Gcorner(:,51)=(/-02,-03,-08,-09,-06,-05,-02/)
      Gcorner(:,52)=(/+02,+05,+06,+09,+08,+03,+02/)
      Gcorner(:,53)=(/-04,-06,-12,-11,-03,-01,-04/)
      Gcorner(:,54)=(/+03,+11,+12,+06,+04,+01,+03/)
      Gcorner(:,55)=(/+10,+02,+01,+07,+09,+12,+10/)
      Gcorner(:,56)=(/+01,+07,+09,+12,+10,+02,+01/) 
      Gcorner(:,57)=(/+05,+10,+11,+08,+07,+04,+05/)
      Gcorner(:,58)=(/-11,-10,-05,-04,-07,-08,-11/)
      !----------------------------------------------------!
      !POLY4A
      Gcorner(:,59)=(/+11,-09,-06,-05,+01,+03,0/)    
      Gcorner(:,60)=(/-01,+05,+06,+09,-11,-03,0/)      
      Gcorner(:,61)=(/+04,-02,-10,-12,+08,+07,0/)    
      Gcorner(:,62)=(/-08,+12,+10,+02,-04,-07,0/)                    
      Gcorner(:,63)=(/-07,+03,+02,+05,-12,-09,0/)    
      Gcorner(:,64)=(/-03,+07,+09,+12,-05,-02,0/)      
      Gcorner(:,65)=(/+01,-08,-11,-10,+06,+04,0/)    
      Gcorner(:,66)=(/+08,-01,-04,-06,+10,+11,0/)                    
      Gcorner(:,67)=(/-02,+11,+12,+06,-07,-01,0/)    
      Gcorner(:,68)=(/+07,-06,-12,-11,+02,+01,0/)    
      Gcorner(:,69)=(/+03,-10,-05,-04,+09,+08,0/)    
      Gcorner(:,70)=(/+10,-03,-08,-09,+04,+05,0/)    
      !----------------------------------------------------!
      !POLY4B
      Gcorner(:,71)=(/-02,+04,+06,+12,-08,-03,0/)          
      Gcorner(:,72)=(/-04,+02,+03,+08,-12,-06,0/)    
      Gcorner(:,73)=(/+05,-01,-07,-09,+11,+10,0/)      
      Gcorner(:,74)=(/+01,-05,-10,-11,+09,+07,0/)    
      Gcorner(:,75)=(/-01,+08,+09,+06,-10,-02,0/)                    
      Gcorner(:,76)=(/-08,+01,+02,+10,-06,-09,0/)    
      Gcorner(:,77)=(/+03,-07,-04,-05,+12,+11,0/)      
      Gcorner(:,78)=(/+07,-03,-11,-12,+05,+04,0/)    
      Gcorner(:,79)=(/+04,-09,-12,-10,+03,+01,0/)                    
      Gcorner(:,80)=(/-03,+10,+12,+09,-04,-01,0/)    
      Gcorner(:,81)=(/+02,-11,-08,-07,+06,+05,0/)    
      Gcorner(:,82)=(/+11,-02,-05,-06,+07,+08,0/)    
      !----------------------------------------------------!
      !POLYC
      Gcorner(:,83) =(/-08,+12,-05,+01,+04,-04,0/)    !POLYC R+.4
      Gcorner(:,84) =(/+12,-05,+01,-08,-07,+07,0/)    !POLYC R+.7    
      Gcorner(:,85) =(/-05,+01,-08,+12,+09,-09,0/)    !POLYC R+.9
      Gcorner(:,86) =(/+01,-08,+12,-05,-06,+06,0/)    !POLYC R+.6                  
      Gcorner(:,87) =(/+08,-12,+05,-01,-02,+02,0/)    !POLYC R-.2
      Gcorner(:,88) =(/-12,+05,-01,+08,+03,-03,0/)    !POLYC R-.3    
      Gcorner(:,89) =(/+05,-01,+08,-12,-11,+11,0/)    !POLYC R-.11
      Gcorner(:,90) =(/-01,+08,-12,+05,+10,-10,0/)    !POLYC R-.10                 
      Gcorner(:,91) =(/+11,-09,+04,-02,-01,+01,0/)    !POLYC S+.1
      Gcorner(:,92) =(/-09,+04,-02,+11,+03,-03,0/)    !POLYC S+.3
      Gcorner(:,93) =(/+04,-02,+11,-09,-08,+08,0/)    !POLYC S+.8
      Gcorner(:,94) =(/-02,+11,-09,+04,+07,-07,0/)    !POLYC S+.7
      Gcorner(:,95) =(/-11,+09,-04,+02,+05,-05,0/)    !POLYC S-.5
      Gcorner(:,96) =(/+09,-04,+02,-11,-10,+10,0/)    !POLYC S-.10
      Gcorner(:,97) =(/-04,+02,-11,+09,+12,-12,0/)    !POLYC S-.12
      Gcorner(:,98) =(/+02,-11,+09,-04,-04,+06,0/)    !POLYC S-.6
      Gcorner(:,99) =(/-03,+07,-06,+10,+05,-05,0/)    !POLYC T+.5
      Gcorner(:,100)=(/+07,-06,+10,-03,-04,+04,0/)    !POLYC T+.4
      Gcorner(:,101)=(/-06,+10,-03,+07,+01,-01,0/)    !POLYC T+.1
      Gcorner(:,102)=(/+10,-03,+07,-06,-02,+02,0/)    !POLYC T+.2
      Gcorner(:,103)=(/+03,-07,+06,-10,-12,+12,0/)    !POLYC T-.12
      Gcorner(:,104)=(/-07,+06,-10,+03,+09,-09,0/)    !POLYC T-.9
      Gcorner(:,105)=(/+06,-10,+03,-07,-08,+08,0/)    !POLYC T-.8
      Gcorner(:,106)=(/-10,+03,-07,+06,+11,-11,0/)    !POLYC T-.11     
      
      !----------------------------------------------------!
      !  SUBVOLUME 'g'RAPH : 'main 'N'odes               !
      !----------------------------------------------------!
      !TETRA
      gnode(:, 1)=(/1,0,0,0/) ; bgnode( 1)=b'10000000'
      gnode(:, 2)=(/3,0,0,0/) ; bgnode( 2)=b'00100000'
      gnode(:, 3)=(/6,0,0,0/) ; bgnode( 3)=b'00000100'
      gnode(:, 4)=(/8,0,0,0/) ; bgnode( 4)=b'00000001'
      gnode(:, 5)=(/2,0,0,0/) ; bgnode( 5)=b'01000000'
      gnode(:, 6)=(/4,0,0,0/) ; bgnode( 6)=b'00010000'
      gnode(:, 7)=(/5,0,0,0/) ; bgnode( 7)=b'00001000'
      gnode(:, 8)=(/7,0,0,0/) ; bgnode( 8)=b'00000010'
      !----------------------------------------------------!
      !PENTA
      gnode(:, 9)=(/1,2,0,0/) ; bgnode( 9)=b'11000000'
      gnode(:,10)=(/1,4,0,0/) ; bgnode(10)=b'10010000'
      gnode(:,11)=(/1,5,0,0/) ; bgnode(11)=b'10001000'
      gnode(:,12)=(/2,3,0,0/) ; bgnode(12)=b'01100000'
      gnode(:,13)=(/3,4,0,0/) ; bgnode(13)=b'00110000'
      gnode(:,14)=(/3,7,0,0/) ; bgnode(14)=b'00100010'
      gnode(:,15)=(/2,6,0,0/) ; bgnode(15)=b'01000100'
      gnode(:,16)=(/5,6,0,0/) ; bgnode(16)=b'00001100'
      gnode(:,17)=(/6,7,0,0/) ; bgnode(17)=b'00000110'
      gnode(:,18)=(/4,8,0,0/) ; bgnode(18)=b'00010001'
      gnode(:,19)=(/5,8,0,0/) ; bgnode(19)=b'00001001'
      gnode(:,20)=(/7,8,0,0/) ; bgnode(20)=b'00000011'
      !----------------------------------------------------!
      !POLY3
      gnode(:,21) =(/1,4,3,0/) ; bgnode(21)=b'10110000'
      gnode(:,22) =(/2,1,4,0/) ; bgnode(22)=b'11010000'
      gnode(:,23) =(/3,2,1,0/) ; bgnode(23)=b'11100000'
      gnode(:,24) =(/4,3,2,0/) ; bgnode(24)=b'01110000'
      gnode(:,25) =(/3,4,8,0/) ; bgnode(25)=b'00110001'
      gnode(:,26) =(/4,8,7,0/) ; bgnode(26)=b'00010011'
      gnode(:,27) =(/7,3,4,0/) ; bgnode(27)=b'00110010'
      gnode(:,28) =(/8,7,3,0/) ; bgnode(28)=b'00100011'      
      gnode(:,29) =(/5,6,7,0/) ; bgnode(29)=b'00001110'
      gnode(:,30) =(/6,7,8,0/) ; bgnode(30)=b'00000111'
      gnode(:,31) =(/7,8,5,0/) ; bgnode(31)=b'00001011'
      gnode(:,32) =(/8,5,6,0/) ; bgnode(32)=b'00001101'      
      gnode(:,33) =(/1,2,6,0/) ; bgnode(33)=b'11000100'
      gnode(:,34) =(/2,6,5,0/) ; bgnode(34)=b'01001100'
      gnode(:,35) =(/5,1,2,0/) ; bgnode(35)=b'11001000'
      gnode(:,36) =(/6,5,1,0/) ; bgnode(36)=b'10001100'
      gnode(:,37) =(/2,3,7,0/) ; bgnode(37)=b'01100010'
      gnode(:,38) =(/3,7,6,0/) ; bgnode(38)=b'00100110'
      gnode(:,39) =(/6,2,3,0/) ; bgnode(39)=b'01100100'
      gnode(:,40) =(/7,6,2,0/) ; bgnode(40)=b'01000110'
      gnode(:,41) =(/1,5,8,0/) ; bgnode(41)=b'10001001'
      gnode(:,42) =(/4,1,5,0/) ; bgnode(42)=b'10011000'
      gnode(:,43) =(/5,8,4,0/) ; bgnode(43)=b'00011001'
      gnode(:,44) =(/8,4,1,0/) ; bgnode(44)=b'10010001'
      !----------------------------------------------------!
      !HEXAE
      gnode(:,45)=(/6,7,3,2/) ; bgnode(45)=b'01100110'
      gnode(:,46)=(/8,5,1,4/) ; bgnode(46)=b'10011001'
      gnode(:,47)=(/5,6,2,1/) ; bgnode(47)=b'11001100'
      gnode(:,48)=(/8,4,3,7/) ; bgnode(48)=b'00110011'
      gnode(:,49)=(/6,5,8,7/) ; bgnode(49)=b'00001111'
      gnode(:,50)=(/1,2,3,4/) ; bgnode(50)=b'11110000' 
      !----------------------------------------------------!
      !POLY4
      gnode(:,51)=(/1,2,6,3/) ; bgnode(51)=b'11100100' 
      gnode(:,52)=(/4,8,7,5/) ; bgnode(52)=b'00011011' 
      gnode(:,53)=(/3,4,8,1/) ; bgnode(53)=b'10110001' 
      gnode(:,54)=(/5,6,7,2/) ; bgnode(54)=b'01001110' 
      gnode(:,55)=(/4,3,2,7/) ; bgnode(55)=b'01110010' 
      gnode(:,56)=(/2,3,7,4/) ; bgnode(56)=b'01110010' 
      gnode(:,57)=(/4,1,5,2/) ; bgnode(57)=b'11011000' 
      gnode(:,58)=(/8,7,3,6/) ; bgnode(58)=b'00100111' 
      !----------------------------------------------------!
      !POLY4A
      gnode(:,59)=(/5,6,2,3/) ; bgnode(59)=b'01101100' 
      gnode(:,60)=(/1,4,8,7/) ; bgnode(60)=b'10010011'  
      gnode(:,61)=(/2,1,5,8/) ; bgnode(61)=b'11001001' 
      gnode(:,62)=(/6,7,3,4/) ; bgnode(62)=b'00110110'                                                                          
      gnode(:,63)=(/6,5,8,4/) ; bgnode(63)=b'00011101' 
      gnode(:,64)=(/1,2,3,7/) ; bgnode(64)=b'11100010' 
      gnode(:,65)=(/2,6,7,8/) ; bgnode(65)=b'01000111' 
      gnode(:,66)=(/5,1,4,3/) ; bgnode(66)=b'10111000'                                                      
      gnode(:,67)=(/1,5,6,7/) ; bgnode(67)=b'10001110' 
      gnode(:,68)=(/2,3,4,8/) ; bgnode(68)=b'01110001' 
      gnode(:,69)=(/5,8,7,3/) ; bgnode(69)=b'00101011' 
      gnode(:,70)=(/4,1,2,6/) ; bgnode(70)=b'11010100'                                                     
      !----------------------------------------------------!
      !POLY4B
      gnode(:,71)=(/1,2,6,7/) ; bgnode(71)=b'11000110'   
      gnode(:,72)=(/3,4,8,5/) ; bgnode(72)=b'00111001'   
      gnode(:,73)=(/4,1,5,6/) ; bgnode(73)=b'10011100'   
      gnode(:,74)=(/2,3,7,8/) ; bgnode(74)=b'01100011'                                               
      gnode(:,75)=(/1,5,8,7/) ; bgnode(75)=b'10001011'   
      gnode(:,76)=(/6,2,3,4/) ; bgnode(76)=b'01110100'   
      gnode(:,77)=(/5,6,7,4/) ; bgnode(77)=b'00011110'   
      gnode(:,78)=(/2,1,4,8/) ; bgnode(78)=b'11010001'                           
      gnode(:,79)=(/2,6,5,8/) ; bgnode(79)=b'01001101'   
      gnode(:,80)=(/1,4,3,7/) ; bgnode(80)=b'10110010'   
      gnode(:,81)=(/4,8,7,6/) ; bgnode(81)=b'00010111'   
      gnode(:,82)=(/5,1,2,3/) ; bgnode(82)=b'11101000'   
      !----------------------------------------------------!
      !POLYC
      gnode(:,83) =(/6,7,3,2/) ; bgnode(83)  =b'01100110'          !polyc r+.4
      gnode(:,84) =(/7,3,2,6/) ; bgnode(84)  =b'01100110'          !POLYC R+.7    
      gnode(:,85) =(/3,2,6,7/) ; bgnode(85)  =b'01100110'          !POLYC R+.9
      gnode(:,86) =(/2,6,7,3/) ; bgnode(86)  =b'01100110'          !POLYC R+.6                                                                            
      gnode(:,87) =(/5,8,4,1/) ; bgnode(87)  =b'10011001'          !POLYC R-.2
      gnode(:,88) =(/8,4,1,5/) ; bgnode(88)  =b'10011001'          !POLYC R-.3 
      gnode(:,89) =(/4,1,5,8/) ; bgnode(89)  =b'10011001'          !POLYC R-.11
      gnode(:,90) =(/1,5,8,4/) ; bgnode(90)  =b'10011001'          !POLYC R-.10                                                       
      gnode(:,91) =(/5,6,2,1/) ; bgnode(91)  =b'11001100'          !POLYC S+.1
      gnode(:,92) =(/6,2,1,5/) ; bgnode(92)  =b'11001100'          !POLYC S+.3
      gnode(:,93) =(/2,1,5,6/) ; bgnode(93)  =b'11001100'          !POLYC S+.8
      gnode(:,94) =(/1,5,6,2/) ; bgnode(94)  =b'11001100'          !POLYC S+.7                                                       
      gnode(:,95) =(/8,7,3,4/) ; bgnode(95)  =b'00110011'          !POLYC S-.5
      gnode(:,96) =(/7,3,4,8/) ; bgnode(96)  =b'00110011'          !POLYC S-.10
      gnode(:,97) =(/3,4,8,7/) ; bgnode(97)  =b'00110011'          !POLYC S-.12
      gnode(:,98) =(/4,8,7,3/) ; bgnode(98)  =b'00110011'          !POLYC S-.6                                              
      gnode(:,99) =(/1,2,3,4/) ; bgnode(99)  =b'11110000'          !POLYC T+.5
      gnode(:,100)=(/2,3,4,1/) ; bgnode(100) =b'11110000'          !POLYC T+.4
      gnode(:,101)=(/3,4,1,2/) ; bgnode(101) =b'11110000'          !POLYC T+.1
      gnode(:,102)=(/4,1,2,3/) ; bgnode(102) =b'11110000'          !POLYC T+.2                          
      gnode(:,103)=(/5,6,7,8/) ; bgnode(103) =b'00001111'          !POLYC T-.12
      gnode(:,104)=(/8,5,6,7/) ; bgnode(104) =b'00001111'          !POLYC T-.9
      gnode(:,105)=(/7,8,5,6/) ; bgnode(105) =b'00001111'          !POLYC T-.8
      gnode(:,106)=(/6,7,8,5/) ; bgnode(106) =b'00001111'          !POLYC T-.11
 
      !----------------------------------------------------!
      !  NODE TAGS FOR EACH OF 58 COMBINATIONS             !
      !----------------------------------------------------!
                            !  1        10        20        30        40        50
                            !  |        |         |         |         |         |  
                            ! '1234567890123456789012345678901234567890123456789012345678'/ 
c      DATA b22CombiNode(01) /B'0000000000000000000000000000000000000000000000000000000000'/  !1,9,10,11,21,22,23,24,25,26,32,33,44,46,47,49,52,54,56,58
c      DATA b22CombiNode(02) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiNode(03) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiNode(04) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiNode(05) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiNode(06) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiNode(07) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiNode(08) /B'0000000000000000000000000000000000000000000000000000000000'/                                    
 
c      DATA b22CombiEdge(01) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(02) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(03) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(04) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(05) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(06) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(07) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(08) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(09) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(10) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(11) /B'0000000000000000000000000000000000000000000000000000000000'/
c      DATA b22CombiEdge(12) /B'0000000000000000000000000000000000000000000000000000000000'/                    
 
C=======================================================================
C definitions
C=======================================================================
 
      !diagonal vectors
      int22_buf%iFACE(1:2,1,1)=(/1,3/); int22_buf%iFACE(1:2,2,1)=(/2,4/)   
      int22_buf%iFACE(1:2,1,2)=(/3,8/); int22_buf%iFACE(1:2,2,2)=(/4,7/) 
      int22_buf%iFACE(1:2,1,3)=(/5,7/); int22_buf%iFACE(1:2,2,3)=(/6,8/) 
      int22_buf%iFACE(1:2,1,4)=(/1,6/); int22_buf%iFACE(1:2,2,4)=(/2,5/) 
      int22_buf%iFACE(1:2,1,5)=(/2,7/); int22_buf%iFACE(1:2,2,5)=(/3,6/) 
      int22_buf%iFACE(1:2,1,6)=(/4,5/); int22_buf%iFACE(1:2,2,6)=(/1,8/)
      
      !face_id -> (node_id, node_id, node_id, node_id)
      int22_buf%nodFACE(1,1:4) = (/1,2,3,4/) 
      int22_buf%nodFACE(2,1:4) = (/4,3,7,8/) 
      int22_buf%nodFACE(3,1:4) = (/8,7,6,5/) 
      int22_buf%nodFACE(4,1:4) = (/5,6,2,1/) 
      int22_buf%nodFACE(5,1:4) = (/3,2,6,7/) 
      int22_buf%nodFACE(6,1:4) = (/5,1,4,8/)
      
      !same in binary
      int22_buf%bNodFACE(01) = b'000011110' !last digit is position 0 and is unused
      int22_buf%bNodFACE(02) = b'110011000'      
      int22_buf%bNodFACE(03) = b'111100000'
      int22_buf%bNodFACE(04) = b'001100110'
      int22_buf%bNodFACE(05) = b'011001100'
      int22_buf%bNodFACE(06) = b'100110010'
      
      !(node_id, face_id) -> Is_On_face
      int22_buf%IsNodeOnFace(1,1:6) = (/.true. ,.false.,.false.,.true .,.false.,.true. /)
      int22_buf%IsNodeOnFace(2,1:6) = (/.true. ,.false.,.false.,.true. ,.true. ,.false./)
      int22_buf%IsNodeOnFace(3,1:6) = (/.true. ,.true. ,.false.,.false.,.true. ,.false./)
      int22_buf%IsNodeOnFace(4,1:6) = (/.true. ,.true. ,.false.,.false.,.false.,.true. /)
      int22_buf%IsNodeOnFace(5,1:6) = (/.false.,.false.,.true. ,.true. ,.false.,.true. /)
      int22_buf%IsNodeOnFace(6,1:6) = (/.false.,.false.,.true. ,.true. ,.true. ,.false./)
      int22_buf%IsNodeOnFace(7,1:6) = (/.false.,.true. ,.true. ,.false.,.true. ,.false./)
      int22_buf%IsNodeOnFace(8,1:6) = (/.false.,.true. ,.true. ,.false.,.false.,.true. /)                                    
 
      !edge number -> (node_id, node_id)            
      int22_buf%iEDGE(1:2, 1)=(/1,2/) ; int22_buf%iEDGE(1:2,  4)=(/3,2/);
      int22_buf%iEDGE(1:2, 2)=(/1,4/) ; int22_buf%iEDGE(1:2,  5)=(/3,4/);
      int22_buf%iEDGE(1:2, 3)=(/1,5/) ; int22_buf%iEDGE(1:2,  6)=(/3,7/);
      int22_buf%iEDGE(1:2, 7)=(/6,2/) ; int22_buf%iEDGE(1:2, 10)=(/8,4/);
      int22_buf%iEDGE(1:2, 8)=(/6,5/) ; int22_buf%iEDGE(1:2, 11)=(/8,5/);
      int22_buf%iEDGE(1:2, 9)=(/6,7/) ; int22_buf%iEDGE(1:2, 12)=(/8,7/);
      
!      !diagonals Identification (if needed)
!      iEDGE(1:2, 13)=(/1,3/) ; iEDGE(1:2, 16)=(/2,4/);
!      iEDGE(1:2, 14)=(/1,6/) ; iEDGE(1:2, 17)=(/2,5/);
!      iEDGE(1:2, 15)=(/1,8/) ; iEDGE(1:2, 18)=(/2,7/);
!      !
!      iEDGE(1:2, 19)=(/3,6/) ; iEDGE(1:2, 21)=(/4,7/)
!      iEDGE(1:2, 20)=(/3,8/) ; iEDGE(1:2, 22)=(/4,5/);
!      !
!      iEDGE(1:2, 23)=(/5,7/) ; iEDGE(1:2, 24)=(/6,8/); 
     
      ! node_id -> 3 adjacent edges
      int22_buf%i22WhichEdge(1, 1:3)=(/-01,-02,-03/)
      int22_buf%i22WhichEdge(2, 1:3)=(/ 01, 04, 07/)
      int22_buf%i22WhichEdge(3, 1:3)=(/-04,-05,-06/)
      int22_buf%i22WhichEdge(4, 1:3)=(/ 02, 05, 10/)
      int22_buf%i22WhichEdge(5, 1:3)=(/ 03, 08, 11/)
      int22_buf%i22WhichEdge(6, 1:3)=(/-07,-08,-09/)
      int22_buf%i22WhichEdge(7, 1:3)=(/ 06, 09, 12/)
      int22_buf%i22WhichEdge(8, 1:3)=(/-10,-11,-12/) 
      
      ! node_id -> 3 adjacent faces
      int22_buf%iFacesFromNode(1,1:3) = (/ 1,4,6 /)                                         
      int22_buf%iFacesFromNode(2,1:3) = (/ 1,4,5 /)  
      int22_buf%iFacesFromNode(3,1:3) = (/ 1,2,5 /)  
      int22_buf%iFacesFromNode(4,1:3) = (/ 1,2,6 /)  
      int22_buf%iFacesFromNode(5,1:3) = (/ 3,4,6 /)  
      int22_buf%iFacesFromNode(6,1:3) = (/ 3,4,5 /)  
      int22_buf%iFacesFromNode(7,1:3) = (/ 2,3,5 /)  
      int22_buf%iFacesFromNode(8,1:3) = (/ 2,3,6 /)                                            
      
      !(face_id,node_id) -> edge_id
      int22_buf%iGetEdge(1,1)=+03; int22_buf%iGetEdge(1,2)=-07; int22_buf%iGetEdge(1,3)=+06; int22_buf%iGetEdge(1,4)=-10
      int22_buf%iGetEdge(2,4)=-02; int22_buf%iGetEdge(2,3)=+04; int22_buf%iGetEdge(2,7)=-09; int22_buf%iGetEdge(2,8)=+11
      int22_buf%iGetEdge(3,8)=+10; int22_buf%iGetEdge(3,7)=-06; int22_buf%iGetEdge(3,6)=+07; int22_buf%iGetEdge(3,5)=-03
      int22_buf%iGetEdge(4,5)=-11; int22_buf%iGetEdge(4,6)=+09; int22_buf%iGetEdge(4,2)=-04; int22_buf%iGetEdge(4,1)=+02
      int22_buf%iGetEdge(5,3)=+05; int22_buf%iGetEdge(5,2)=-01; int22_buf%iGetEdge(5,6)=+08; int22_buf%iGetEdge(5,7)=-12
      int22_buf%iGetEdge(6,5)=-08; int22_buf%iGetEdge(6,1)=+01; int22_buf%iGetEdge(6,4)=-05; int22_buf%iGetEdge(6,8)=+12                             
 
      !(face_id) -> {edge_ids}
      int22_buf%iNormalEdgesList(1,1:4)= (/ +03,+06,-07,-10 /) 
      int22_buf%iNormalEdgesList(2,1:4)= (/ -02,+04,-09,+11 /)
      int22_buf%iNormalEdgesList(3,1:4)= (/ -03,-06,+07,+10 /)
      int22_buf%iNormalEdgesList(4,1:4)= (/ +02,-04,+09,-11 /)
      int22_buf%iNormalEdgesList(5,1:4)= (/ -01,+05,+08,-12 /)
      int22_buf%iNormalEdgesList(6,1:4)= (/ +01,-05,-08,+12 /)                              
 
 
      !(face_id,node_id) -> opposite node
      int22_buf%IGetOppositeNode(1,1)=05; int22_buf%IGetOppositeNode(1,2)=06; 
      int22_buf%IGetOppositeNode(2,4)=01; int22_buf%IGetOppositeNode(2,3)=02; 
      int22_buf%IGetOppositeNode(3,8)=04; int22_buf%IGetOppositeNode(3,7)=03; 
      int22_buf%IGetOppositeNode(4,5)=08; int22_buf%IGetOppositeNode(4,6)=07; 
      int22_buf%IGetOppositeNode(5,3)=04; int22_buf%IGetOppositeNode(5,2)=01; 
      int22_buf%IGetOppositeNode(6,5)=06; int22_buf%IGetOppositeNode(6,1)=02;
       
      int22_buf%IGetOppositeNode(1,3)=07; int22_buf%IGetOppositeNode(1,4)=08;
      int22_buf%IGetOppositeNode(2,7)=06; int22_buf%IGetOppositeNode(2,8)=05;
      int22_buf%IGetOppositeNode(3,6)=02; int22_buf%IGetOppositeNode(3,5)=01;
      int22_buf%IGetOppositeNode(4,2)=03; int22_buf%IGetOppositeNode(4,1)=04;
      int22_buf%IGetOppositeNode(5,6)=05; int22_buf%IGetOppositeNode(5,7)=08;
      int22_buf%IGetOppositeNode(6,4)=03; int22_buf%IGetOppositeNode(6,8)=07; 
      
      !(face_id,node_id) -> (tangent node 1, tangent node 2)
      int22_buf%IGetTransvNodes(:,:,:)   = 0
      !     
      int22_buf%IGetTransvNodes(1,1,1:2) = (/4,2/) ; 
      int22_buf%IGetTransvNodes(1,2,1:2) = (/1,3/) ; 
      int22_buf%IGetTransvNodes(1,3,1:2) = (/2,4/) ; 
      int22_buf%IGetTransvNodes(1,4,1:2) = (/3,1/) ; 
      
      int22_buf%IGetTransvNodes(2,4,1:2) = (/3,8/) ; 
      int22_buf%IGetTransvNodes(2,3,1:2) = (/7,4/) ; 
      int22_buf%IGetTransvNodes(2,7,1:2) = (/8,3/) ; 
      int22_buf%IGetTransvNodes(2,8,1:2) = (/4,7/) ; 
 
      int22_buf%IGetTransvNodes(3,8,1:2) = (/5,7/) ;
      int22_buf%IGetTransvNodes(3,7,1:2) = (/8,6/) ;
      int22_buf%IGetTransvNodes(3,6,1:2) = (/7,5/) ;
      int22_buf%IGetTransvNodes(3,5,1:2) = (/6,8/) ;      
 
      int22_buf%IGetTransvNodes(4,5,1:2) = (/1,6/) ; 
      int22_buf%IGetTransvNodes(4,6,1:2) = (/5,2/) ; 
      int22_buf%IGetTransvNodes(4,2,1:2) = (/6,1/) ; 
      int22_buf%IGetTransvNodes(4,1,1:2) = (/2,5/) ; 
      
      int22_buf%IGetTransvNodes(5,3,1:2) = (/2,7/) ; 
      int22_buf%IGetTransvNodes(5,2,1:2) = (/6,3/) ; 
      int22_buf%IGetTransvNodes(5,6,1:2) = (/7,2/) ; 
      int22_buf%IGetTransvNodes(5,7,1:2) = (/3,6/) ; 
 
      int22_buf%IGetTransvNodes(6,1,1:2) = (/4,5/) ;
      int22_buf%IGetTransvNodes(6,4,1:2) = (/8,1/) ;
      int22_buf%IGetTransvNodes(6,8,1:2) = (/5,4/) ;
      int22_buf%IGetTransvNodes(6,5,1:2) = (/1,8/) ;            
 
      ! face_id -> (edge_id, edge_id, edge_id, edge_id)
      int22_buf%iGetEdgesFromFace(1,1:4) = (/01,02,05,04/)
      int22_buf%iGetEdgesFromFace(2,1:4) = (/05,10,12,06/)
      int22_buf%iGetEdgesFromFace(3,1:4) = (/12,11,08,09/)
      int22_buf%iGetEdgesFromFace(4,1:4) = (/01,03,08,07/)
      int22_buf%iGetEdgesFromFace(5,1:4) = (/04,06,09,07/)
      int22_buf%iGetEdgesFromFace(6,1:4) = (/02,03,11,10/) 
      
      ! edge_id -> edge_id, edge_id     
      int22_buf%iGetAdjFace( 01, 1:2 ) = (/ 01,04 /)
      int22_buf%iGetAdjFace( 02, 1:2 ) = (/ 01,06 /)
      int22_buf%iGetAdjFace( 03, 1:2 ) = (/ 04,06 /)
      int22_buf%iGetAdjFace( 04, 1:2 ) = (/ 01,05 /)
      int22_buf%iGetAdjFace( 05, 1:2 ) = (/ 01,02 /)
      int22_buf%iGetAdjFace( 06, 1:2 ) = (/ 02,05 /)
      int22_buf%iGetAdjFace( 07, 1:2 ) = (/ 04,05 /)
      int22_buf%iGetAdjFace( 08, 1:2 ) = (/ 03,04 /)
      int22_buf%iGetAdjFace( 09, 1:2 ) = (/ 03,05 /)
      int22_buf%iGetAdjFace( 10, 1:2 ) = (/ 02,06 /)
      int22_buf%iGetAdjFace( 11, 1:2 ) = (/ 03,06 /)
      int22_buf%iGetAdjFace( 12, 1:2 ) = (/ 02,03 /)
      
      
      !FACE_ID,EDGE_ID- > NODE_ID
      int22_buf%iGetNodeFromEdgeAndFace(:,:) =  0
      
      int22_buf%iGetNodeFromEdgeAndFace(1,03) = 1 ; int22_buf%iGetNodeFromEdgeAndFace(2,02) = 4
      int22_buf%iGetNodeFromEdgeAndFace(1,06) = 3 ; int22_buf%iGetNodeFromEdgeAndFace(2,04) = 3
      int22_buf%iGetNodeFromEdgeAndFace(1,07) = 6 ; int22_buf%iGetNodeFromEdgeAndFace(2,09) = 7
      int22_buf%iGetNodeFromEdgeAndFace(1,10) = 8 ; int22_buf%iGetNodeFromEdgeAndFace(2,11) = 8           
 
      int22_buf%iGetNodeFromEdgeAndFace(3,03) = 5 ; int22_buf%iGetNodeFromEdgeAndFace(4,02) = 1
      int22_buf%iGetNodeFromEdgeAndFace(3,06) = 7 ; int22_buf%iGetNodeFromEdgeAndFace(4,04) = 2
      int22_buf%iGetNodeFromEdgeAndFace(3,07) = 6 ; int22_buf%iGetNodeFromEdgeAndFace(4,09) = 6
      int22_buf%iGetNodeFromEdgeAndFace(3,10) = 8 ; int22_buf%iGetNodeFromEdgeAndFace(4,11) = 5
      
      int22_buf%iGetNodeFromEdgeAndFace(5,01) = 2 ; int22_buf%iGetNodeFromEdgeAndFace(6,01) = 1
      int22_buf%iGetNodeFromEdgeAndFace(5,05) = 3 ; int22_buf%iGetNodeFromEdgeAndFace(6,05) = 4 
      int22_buf%iGetNodeFromEdgeAndFace(5,08) = 6 ; int22_buf%iGetNodeFromEdgeAndFace(6,08) = 5 
      int22_buf%iGetNodeFromEdgeAndFace(5,12) = 7 ; int22_buf%iGetNodeFromEdgeAndFace(6,12) = 8 
      
 
 
      !iRightEdge(Face,nodeID)
      !next edge is clockwise edge when face is looked from cell centroid.
      ! "-": Means centripetal
      ! "+" : means centrifugal
      int22_buf%iLeftEdge(1,1)=+02; int22_buf%iLeftEdge(1,2)=-01; int22_buf%iLeftEdge(1,3)=+04; int22_buf%iLeftEdge(1,4)=-05;  !1234  <-  nodFACE(1,1:4) = (/1,2,3,4/) 
      int22_buf%iLeftEdge(2,1)=-10; int22_buf%iLeftEdge(2,2)=+05; int22_buf%iLeftEdge(2,3)=-06; int22_buf%iLeftEdge(2,4)=+12;  !4378  <-  nodFACE(2,1:4) = (/4,3,7,8/) 
      int22_buf%iLeftEdge(3,1)=+11; int22_buf%iLeftEdge(3,2)=-12; int22_buf%iLeftEdge(3,3)=+09; int22_buf%iLeftEdge(3,4)=-08;  !8765  <-  nodFACE(3,1:4) = (/8,7,6,5/) 
      int22_buf%iLeftEdge(4,1)=-03; int22_buf%iLeftEdge(4,2)=+08; int22_buf%iLeftEdge(4,3)=-07; int22_buf%iLeftEdge(4,4)=+01;  !5621  <-  nodFACE(4,1:4) = (/5,6,2,1/) 
      int22_buf%iLeftEdge(5,1)=+06; int22_buf%iLeftEdge(5,2)=-04; int22_buf%iLeftEdge(5,3)=+07; int22_buf%iLeftEdge(5,4)=-09;  !3267  <-  nodFACE(5,1:4) = (/3,2,6,7/) 
      int22_buf%iLeftEdge(6,1)=-11; int22_buf%iLeftEdge(6,2)=+03; int22_buf%iLeftEdge(6,3)=-02; int22_buf%iLeftEdge(6,4)=+10;  !5148  <-  nodFACE(6,1:4) = (/5,1,4,8/)             
      
      IF(int22>0)THEN
        DO nin=1,int22
          DO j=1,12*i22len
            edge_list(nin,j)%NBCUT = 0
          ENDDO
        ENDDO
      ENDIF
            
C=======================================================================
      idt_int22                                                       = 0
      dt22_min                                                        = ep30
      dx22_min                                                        = ep30
      dx22min_l(0:nthread-1)                                          = ep30
      v22_max                                                         = zero
      v22max_l(0:nthread-1)                                           = zero
      ncross22                                                        = five
      critmerge22                                                     = half                                                                         
      critdemerge22                                                   = half                                                                         
      critdvol22                                                      = fourth                                                                          
      nadj_f                                                          = 5                                                                              
      nbold(1)                                                        = nb                                                                             
      iiad22(:,:)                                                     = 0                                                                              
      nbcut_max                                                       = 8 + 6  !cut plane in cut cell, maximum reached with 8 tetra (all 8 nodes isolated)  +6 closure surface
      i22_degenerated                                                 = 0                                                                              
      uvarl(0:nthread-1,1:i22len,1:m51_n0phas+trimat*m51_nvphas+i22law37) = zero
      supercellvol_l(0:nthread-1,0:trimat+i22law37,1:i22len)          = zero
      eint_l(0:nthread-1,1:i22len)                                    = zero
      mom_l(1:3,0:nthread-1,1:i22len)                                 = zero
      rho_l(0:nthread-1,1:i22len)                                     = zero
      sig_l(0:nthread-1,1:6,1:i22len)                                 = zero
      vold_l(0:nthread-1,0:trimat+i22law37,1:i22len)                  = zero      
      unlinked_cells_l(0:nthread-1,1:2,1:i22len)                      = 0         
      n_unlinked_l(0:nthread-1)                                       = 0         
      DO nin=1,int22
        old_secndlist(nin,1:i22len)%Num                               = 0         
        old_secndlist(nin,1:i22len)%NumSecndNodes                     = 0         
        DO j=1,24
          old_secndlist(nin,1:i22len)%FM(j)                           = 0         
          old_secndlist(nin,1:i22len)%FV(j)                           = 0         
          old_secndlist(nin,1:i22len)%IV(j)                           = 0         
          old_secndlist(nin,1:i22len)%IBV(j)                          = 0         
          old_secndlist(nin,1:i22len)%ICELLv(j)                       = 0         
          old_secndlist(nin,1:i22len)%VOL(j)                          = zero      
          old_secndlist(nin,1:i22len)%NumNOD_Cell(j)                  = 0
          DO k=1,8
          old_secndlist(nin,1:i22len)%ListNodID(j,k)                  = 0
          ENDDO
        enddo!next J
      enddo!next NIN
C=======================================================================
 
      !--------------------!
      !         ALL        !
      !   ON(1) / OFF(0)   !      
      !--------------------!      
      iflg = 0 !IOUTP22
 
      IF(int22 > 0)THEN
        ibug22               =  0
        !-----------------------!                                                                             
        ibug22_ccbuflist     =  0
        ibug22_antidiff      = -0 
        ibug22_upwind        = -0
        ibug22_tri           =  0
        ibug22_trivox        =  0                                                                                   
        ibug22_intersect     = -0
        ibug22_subvol        =  12781
        ibug22_aera          =  0
        ibug22_ident         =  12781
        ibug22_convec        =  0
        ibug22_spe_convec    =  0
        ibug22_eos           =  0
        ibug22_dvol          = -0
        ibug22_sigeps51      =  0
        ibug22_ssp           =  0
        ibug22_flux          = -0
        ibug22_flux22        = -0
        ibug22_volfrac       =  0
        ibug22_conti         =  0
        ibug22_fint          =  12781
        ibug22_sinit         =  0
        ibug22_nodalp        =  0
        ibug22_amouil        =  0
        ibug22_fcont         = -0
        ibug22_fcontass      = -0
        ibug22_swet          = -0
        ibug22_destroy       = -0
        ibug22_merge         = -0
        ibug22_link_switch   = -0
        ibug22_tracking      = -0
        ibug22_prediction    = -0
        ibug22_undirectlink  = -0
        ibug22_truss         =  0
        ibug22_orphannodes   =  0
        ibug22_dtmin         =  0
        ibug22_outp_ccbuffer =  0    !enables to output bricks in cut cell buffer (HyperMesh tcl script)            
        ibug22_outp_intpoint =  0    !enables to output intersection points (HyperMesh tcl script)                  
        ibug22_nvar          =  2    !-1: all, 1:mass, 2:energy                                                     
        ibug22_itrimat       =  1    !-1: all, 0:global, 1:submat1, 2:submat2, 3:submat3, 4:submat4   
        ibug22_cc_length     =  0              
      ENDIF
      
      
 
      !---------------!            
      IF(int22 > 0)THEN
        ibug22_ccbuflist     = ibug22_ccbuflist       * iflg
        ibug22_antidiff      = ibug22_antidiff        * iflg    
        ibug22_upwind        = ibug22_upwind          * iflg    
        ibug22_tri           = ibug22_tri             * iflg    
        ibug22_trivox        = ibug22_trivox          * iflg    
        ibug22_intersect     = ibug22_intersect       * iflg    
        ibug22_subvol        = ibug22_subvol          * iflg    
        ibug22_aera          = ibug22_aera            * iflg    
        ibug22_ident         = ibug22_ident           * iflg    
        ibug22_convec        = ibug22_convec          * iflg    
        ibug22_spe_convec    = ibug22_spe_convec      * iflg    
        ibug22_eos           = ibug22_eos             * iflg    
        ibug22_dvol          = ibug22_dvol            * iflg    
        ibug22_sigeps51      = ibug22_sigeps51        * iflg    
        ibug22_ssp           = ibug22_ssp             * iflg    
        ibug22_flux          = ibug22_flux            * iflg    
        ibug22_flux22        = ibug22_flux22          * iflg    
        ibug22_volfrac       = ibug22_volfrac         * iflg    
        ibug22_conti         = ibug22_conti           * iflg    
        ibug22_fint          = ibug22_fint            * iflg    
        ibug22_sinit         = ibug22_sinit           * iflg    
        ibug22_nodalp        = ibug22_nodalp          * iflg    
        ibug22_amouil        = ibug22_amouil          * iflg    
        ibug22_fcont         = ibug22_fcont           * iflg    
        ibug22_fcontass      = ibug22_fcontass        * iflg    
        ibug22_swet          = ibug22_swet            * iflg    
        ibug22_destroy       = ibug22_destroy         * iflg    
        ibug22_merge         = ibug22_merge           * iflg    
        ibug22_link_switch   = ibug22_link_switch     * iflg    
        ibug22_tracking      = ibug22_tracking        * iflg            
        ibug22_prediction    = ibug22_prediction      * iflg            
        ibug22_undirectlink  = ibug22_undirectlink    * iflg    
        ibug22_truss         = ibug22_truss           * iflg    
        ibug22_dtmin         = ibug22_dtmin           * iflg    
        ibug22_orphannodes   = ibug22_orphannodes     * iflg         
        ibug22_outp_ccbuffer = ibug22_outp_ccbuffer   * iflg    
        ibug22_outp_intpoint = ibug22_outp_intpoint   * iflg    
        ibug22_nvar          = ibug22_nvar            * iflg    
        ibug22_itrimat       = ibug22_itrimat         * iflg   
        ibug22_cc_length     = ibug22_cc_length       * iflg
      ENDIF      
      
      IF(int22>0 .AND. ibug22 == 1 .AND. iflg/=0)THEN
        print *, "=================="        
        print *, "= Output Flags : ="
        print *, "=================="        
        IF(ibug22              /=0)print *,"IBUG22              "
        IF(ibug22_tri          /=0)print *,"IBUG22_tri          "
        IF(ibug22_trivox       /=0)print *,"IBUG22_trivox       "
        IF(ibug22_intersect    /=0)print *,"IBUG22_intersect    "
        IF(ibug22_subvol       /=0)print *,"IBUG22_subvol       "
        IF(ibug22_aera         /=0)print *,"IBUG22_aera         "
        IF(ibug22_ident        /=0)print *,"IBUG22_ident        "
        IF(ibug22_convec       /=0)print *,"IBUG22_convec       "
        IF(ibug22_eos          /=0)print *,"IBUG22_eos          "
        IF(ibug22_dvol         /=0)print *,"IBUG22_dvol         "
        IF(ibug22_sigeps51     /=0)print *,"IBUG22_sigeps51     "
        IF(ibug22_ssp          /=0)print *,"IBUG22_ssp          "
        IF(ibug22_flux         /=0)print *,"IBUG22_flux         "
        IF(ibug22_flux22       /=0)print *,"IBUG22_flux22       "
        IF(ibug22_volfrac      /=0)print *,"IBUG22_volfrac      "
        IF(ibug22_conti        /=0)print *,"IBUG22_conti        "
        IF(ibug22_fint         /=0)print *,"IBUG22_fint         "
        IF(ibug22_sinit        /=0)print *,"IBUG22_sinit        "
        IF(ibug22_nodalp       /=0)print *,"IBUG22_nodalp       "
        IF(ibug22_amouil       /=0)print *,"IBUG22_amouil       "
        IF(ibug22_spe_convec   /=0)print *,"IBUG22_spe_convec   "
        IF(ibug22_fcont        /=0)print *,"IBUG22_fcont        " 
        IF(ibug22_swet         /=0)print *,"IBUG22_Swet         " 
        IF(ibug22_destroy      /=0)print *,"IBUG22_destroy      "               
        IF(ibug22_merge        /=0)print *,"IBUG22_merge        " 
        IF(ibug22_undirectlink /=0)print *,"IBUG22_UndirectLink "         
        IF(ibug22_truss        /=0)print *,"IBUG22_Truss        " 
        IF(ibug22_dtmin        /=0)print *,"IBUG22_DtMin        "                 
        IF(ibug22_outp_ccbuffer/=0)print *,"IBUG22_OUTP_CCbuffer" 
        IF(ibug22_outp_intpoint/=0)print *,"IBUG22_OUTP_IntPoint"                 
        print *, "=================="                                                                                                                     
      ENDIF
 
!      IF(NTHREAD>1)THEN
!        print *, "**error inter22 : Multithreading not yet validated"
!        stop
!      ENDIF
 
      RETURN