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