33
34
35
36
37
38
39#include "implicit_f.inc"
40
41
42
43#include "com04_c.inc"
44#include "param_c.inc"
45
46
47
48 INTEGER PROC, LLRBE2_L, LEN_IA,NBDDRBE2,NUMNOD_L,ITABRBE2M(*),
49 . IRBE2(NRBE2L,*), LRBE2(*), NODLOCAL(*)
50
51
52
53 INTEGER NLOCAL
55
56
57
58 INTEGER I, J, K, KL, NSNL, NSN, ILP, NLL, NN, M, N,
59 . IBVTMP(NRBE2L,NRBE2), LBVTMP(LLRBE2_L),IAD,II
60
61
62
63 kl=0
64 nsnl = 0
65 ii = 0
66 DO i = 1, nrbe2
67 k = irbe2(1,i)
68 m = irbe2(3,i)
69 nsn = irbe2(5,i)
70 iad = nsnl
71 IF (nsn > 0) ii = ii +1
72 IF (nsn ==0) cycle
73 DO j = 1, nsn
74 n = lrbe2(k+j)
75 IF (
nlocal(n,proc+1)==1)
THEN
76 nsnl = nsnl + 1
77 lbvtmp(nsnl) = nodlocal(n)
78 ENDIF
79 ENDDO
80 IF (
nlocal(m,proc+1)==1)
THEN
81 kl = kl + 1
82 DO j = 1, nrbe2l
83 ibvtmp(j,kl) = irbe2(j,i)
84 ENDDO
85 ibvtmp(1,kl) = iad
86 ibvtmp(5,kl) = nsnl-iad
87 ibvtmp(3,kl) = nodlocal(m)
88
89 ibvtmp(10,kl) = ii
90 ibvtmp(12,kl) = nsn
91
92 ENDIF
93 ENDDO
94
97 len_ia = len_ia + nrbe2l*kl + llrbe2_l
98
99 IF(nbddrbe2>0) THEN
100 itabrbe2m(1:numnod_l) = 0
101 ilp=0
102 DO k = 1, kl
103 n = ibvtmp(3,k)
104 IF(itabrbe2m(n)==0) THEN
105 ilp = ilp + 1
106
107 itabrbe2m(n) = n
108 ENDIF
109 END DO
110 ENDIF
111
112 RETURN
void write_i_c(int *w, int *len)