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