30
31
32
33
34#include "implicit_f.inc"
35
36
37
38#include "com04_c.inc"
39
40
41
42 INTEGER IXS(NIXS,*), IXC(NIXC,*), IXTG(NIXTG,*),
43 . IADD(*), INVC(*)
44
45
46
47 INTEGER I, J, N
48
49 iadd(1) = 1
50 iadd(2) = 1
51
52 DO i=3,numnod+1
53 iadd(i)=0
54 ENDDO
55
56 DO j=2,5
57 DO i=1,numelc
58 n = ixc(j,i) + 2
59 iadd(n)=iadd(n)+1
60 ENDDO
61 ENDDO
62
63 DO j=2,4
64 DO i=1,numeltg
65 n = ixtg(j,i) + 2
66 iadd(n)=iadd(n)+1
67 ENDDO
68 ENDDO
69
70 DO j=2,9
71 DO i=1,numels
72 n = ixs(j,i) + 2
73 iadd(n)=iadd(n)+1
74 ENDDO
75 ENDDO
76
77 DO i=3,numnod+1
78 iadd(i)=iadd(i)+iadd(i-1)
79 ENDDO
80
81 DO i=1,numelc
82
83 n = ixc(2,i) + 1
84 invc(iadd(n)) = i
85 iadd(n) = iadd(n) + 1
86
87 n = ixc(3,i) + 1
88 invc(iadd(n)) = i
89 iadd(n) = iadd(n) + 1
90
91 n = ixc(4,i) + 1
92 invc(iadd(n)) = i
93 iadd(n) = iadd(n) + 1
94
95 n = ixc(5,i) + 1
96 invc(iadd(n)) = i
97 iadd(n) = iadd(n) + 1
98
99 ENDDO
100
101 DO i=1,numeltg
102 n = ixtg(2,i) + 1
103 invc(iadd(n)) = i+numelc
104 iadd(n) = iadd(n) + 1
105
106 n = ixtg(3,i) + 1
107 invc(iadd(n)) = i+numelc
108 iadd(n) = iadd(n) + 1
109
110 n = ixtg(4,i) + 1
111 invc(iadd(n)) = i+numelc
112 iadd(n) = iadd(n) + 1
113
114 ENDDO
115
116 DO i=1,numels
117
118 n = ixs(2,i) + 1
119 invc(iadd(n)) = i+numelc+numeltg
120 iadd(n) = iadd(n) + 1
121
122 n = ixs(3,i) + 1
123 invc(iadd(n)) = i+numelc+numeltg
124 iadd(n) = iadd(n) + 1
125
126 n = ixs(4,i) + 1
127 invc(iadd(n)) = i+numelc+numeltg
128 iadd(n) = iadd(n) + 1
129
130 n = ixs(5,i) + 1
131 invc(iadd(n)) = i+numelc+numeltg
132 iadd(n) = iadd(n) + 1
133
134 n = ixs(6,i) + 1
135 invc(iadd(n)) = i+numelc+numeltg
136 iadd(n) = iadd(n) + 1
137
138 n = ixs(7,i) + 1
139 invc(iadd(n)) = i+numelc+numeltg
140 iadd(n) = iadd(n) + 1
141
142 n = ixs(8,i) + 1
143 invc(iadd(n)) = i+numelc+numeltg
144 iadd(n) = iadd(n) + 1
145
146 n = ixs(9,i) + 1
147 invc(iadd(n)) = i+numelc+numeltg
148 iadd(n) = iadd(n) + 1
149 ENDDO
150
151
152
153
154 RETURN
1.15.0