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 IXQ(NIXQ,NUMELQ), N4_VOIS(NUMELQ+NQVOIS,4), ITABM1(*),NV46
52 my_real flux(4,*),flux_vois(numelq+nqvois,nv46)
53 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
54
55
56
57 INTEGER I,II,K,KK,L,IV,NN,IPERM(2,4),TAG(0:NUMNOD), NI(2), IAD2
58
59
60
61 INTEGER SYSFUS2
62
63 DATA iperm/1,2,
64 . 2,3,
65 . 3,4,
66 . 4,1/
67
68
69
70 DO i=lft,llt
71 ii = i+nft
72 DO l=1,4
73 tag(ixq(1+l,ii)) = 0
74 ENDDO
75 ENDDO
76
77 DO i=lft,llt
78 ii = i+nft
79 iad2 = ale_connect%ee_connect%iad_connect(ii)
80 DO k=1,4
81 iv = ale_connect%ee_connect%connected(iad2 + k - 1)
82 IF(iv > numels)THEN
83 DO kk=1,4
84 IF(n4_vois(iv,1) /= 0.AND.
85 . flux_vois(iv,kk) /= -ep20)THEN
86 nn =0
87 DO l=1,2
88 ni(l) =
sysfus2(n4_vois(iv,iperm(l,kk)),itabm1,numnod )
89 tag(ni(l)) = 1
90 ENDDO
91 DO l=1,2
92 nn=nn+tag(ixq(1+iperm(l,k),ii))
93 ENDDO
94 DO l=1,2
95 tag(ni(l)) = 0
96 ENDDO
97 IF(nn == 2)THEN
98 flux(k,i)=-flux_vois(iv,kk)
99 EXIT
100 ENDIF
101 ENDIF
102 ENDDO
103 ENDIF
104 ENDDO
105 ENDDO
106
107
108 RETURN
integer function sysfus2(iu, itabm1, numnod)