35
36
37
39 use element_mod , only : nixs
40
41
42
43#include "implicit_f.inc"
44
45
46
47#include "spmd_c.inc"
48#include "com04_c.inc"
49#include "vect01_c.inc"
50
51
52
53 INTEGER IXS(NIXS,NUMELS), N4_VOIS(NUMELS+NSVOIS,8),ITABM1(*),NV46
54 my_real flux(6,*), flux_vois(numels+nsvois,nv46)
55 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
56
57
58
59 INTEGER I,II,K,KK,L,IV,NN,IPERM(4,6),TAG(0:NUMNOD), NI(4), IAD2
60
61
62
63 INTEGER SYSFUS2
64
65 DATA iperm/3,1,2,4,
66 . 7,4,3,8,
67 . 6,8,7,5,
68 . 2,5,6,1,
69 . 7,2,6,3,
70 . 8,1,4,5/
71
72
73
74 DO i=lft,llt
75 ii = i+nft
76 DO l=1,8
77 tag(ixs(1+l,ii)) = 0
78 ENDDO
79 ENDDO
80
81 DO i=lft,llt
82 ii = i+nft
83 iad2 = ale_connect%ee_connect%iad_connect(ii)
84 DO k=1,6
85 iv = ale_connect%ee_connect%connected(iad2 + k - 1)
86 IF(iv > numels)THEN
87 DO kk=1,6
88 IF(n4_vois(iv,1) /= 0.AND.
89 . flux_vois(iv,kk) /= -ep20)THEN
90 nn =0
91 DO l=1,4
92 ni(l) =
sysfus2(n4_vois(iv,iperm(l,kk)),itabm1,numnod)
93 tag(ni(l)) = 1
94 ENDDO
95 DO l=1,4
96 nn=nn+tag(ixs(1+iperm(l,k),ii))
97 ENDDO
98 DO l=1,4
99 tag(ni(l)) = 0
100 ENDDO
101 IF(nn==4)THEN
102 flux(k,i)=-flux_vois(iv,kk)
103 EXIT
104 ENDIF
105 ENDIF
106 ENDDO
107 ENDIF
108 enddo
109 enddo
110
111 RETURN
integer function sysfus2(iu, itabm1, numnod)