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