46
47
48
49#include "implicit_f.inc"
50
51
52
53#include "mvsiz_p.inc"
54
55
56
57
58
59
60 INTEGER, INTENT(IN) :: JCVT
61 INTEGER, INTENT(IN) :: NEL
62 INTEGER , IKSUP
64 . py1(*), py2(*), py3(*), py4(*),
65 . pz1(*), pz2(*), pz3(*), pz4(*),
66 . pyc1(*), pyc2(*), pzc1(*), pzc2(*), ay(*), r22(*), r23(*),
67 . k11(2,2,*), k12(2,2,*), k13(2,2,*), k14(2,2,*), k22(2,2,*),
68 . k23(2,2,*), k24(2,2,*), k33(2,2,*), k34(2,2,*), k44(2,2,*),
69 . k11u(2,2,*),k12u(2,2,*),k13u(2,2,*),k14u(2,2,*),k22u(2,2,*),
70 . k23u(2,2,*),k24u(2,2,*),k33u(2,2,*),k34u(2,2,*),k44u(2,2,*),
71 . k11l(2,2,*),k12l(2,2,*),k13l(2,2,*),k14l(2,2,*),k22l(2,2,*),
72 . k23l(2,2,*),k24l(2,2,*),k33l(2,2,*),k34l(2,2,*),k44l(2,2,*),
73 . dd(3,3,*), gg(*), dg(3,3,*), g33(3,3,*), off(*)
74
75
76
77 INTEGER I, J, IS
79 . pyc3(mvsiz), pyc4(mvsiz), pzc3(mvsiz), pzc4(mvsiz)
80
81
82
83 DO i=1,nel
84 pyc3(i) = -pyc1(i)
85 pyc4(i) = -pyc2(i)
86 pzc3(i) = -pzc1(i)
87 pzc4(i) = -pzc2(i)
88 ENDDO
89
90 IF (iksup==0) THEN
91 is = 1
93 1 py1, pz1, py1, pz1,
94 2 pyc1, pzc1, pyc1, pzc1,
95 3 ay, r22, r23, dd,
96 4 gg, k11, k11u, k11l,
97 5 is, nel, jcvt)
99 1 py2, pz2, py2, pz2,
100 2 pyc2, pzc2, pyc2, pzc2,
101 3 ay, r22, r23, dd,
102 4 gg, k22, k22u, k22l,
103 5 is, nel, jcvt)
105 1 py3, pz3, py3, pz3,
106 2 pyc3, pzc3, pyc3, pzc3,
107 3 ay, r22, r23, dd,
108 4 gg, k33, k33u, k33l,
109 5 is, nel, jcvt)
111 1 py4, pz4, py4, pz4,
112 2 pyc4, pzc4, pyc4, pzc4,
113 3 ay, r22, r23, dd,
114 4 gg, k44, k44u, k44l,
115 5 is, nel, jcvt)
116 is = 0
118 1 py1, pz1, py2, pz2,
119 2 pyc1, pzc1, pyc2, pzc2,
120 3 ay, r22, r23, dd,
121 4 gg, k12, k12u, k12l,
122 5 is, nel, jcvt)
124 1 py1, pz1, py3, pz3,
125 2 pyc1, pzc1, pyc3, pzc3,
126 3 ay, r22, r23, dd,
127 4 gg, k13, k13u, k13l,
128 5 is, nel, jcvt)
130 1 py1, pz1, py4, pz4,
131 2 pyc1, pzc1, pyc4, pzc4,
132 3 ay, r22, r23, dd,
133 4 gg, k14, k14u, k14l,
134 5 is, nel, jcvt)
136 1 py2, pz2, py3, pz3,
137 2 pyc2, pzc2, pyc3, pzc3,
138 3 ay, r22, r23, dd,
139 4 gg, k23, k23u, k23l,
140 5 is, nel, jcvt)
142 1 py2, pz2, py4, pz4,
143 2 pyc2, pzc2, pyc4, pzc4,
144 3 ay, r22, r23, dd,
145 4 gg, k24, k24u, k24l,
146 5 is, nel, jcvt)
148 1 py3, pz3, py4, pz4,
149 2 pyc3, pzc3, pyc4, pzc4,
150 3 ay, r22, r23, dd,
151 4 gg, k34, k34u, k34l,
152 5 is, nel, jcvt)
153 ELSE
154 is = 1
156 1 py1, pz1, py1, pz1,
157 2 pyc1, pzc1, pyc1, pzc1,
158 3 ay, r22, r23, dd,
159 4 dg, g33, k11, is,
160 5 nel, jcvt)
162 1 py2, pz2, py2, pz2,
163 2 pyc2, pzc2, pyc2, pzc2,
164 3 ay, r22, r23, dd,
165 4 dg, g33, k22, is,
166 5 nel, jcvt)
168 1 py3, pz3, py3, pz3,
169 2 pyc3, pzc3, pyc3, pzc3,
170 3 ay, r22, r23, dd,
171 4 dg, g33, k33, is,
172 5 nel, jcvt)
174 1 py4, pz4, py4, pz4,
175 2 pyc4, pzc4, pyc4, pzc4,
176 3 ay, r22, r23, dd,
177 4 dg, g33, k44, is,
178 5 nel, jcvt)
179 is = 0
181 1 py1, pz1, py2, pz2,
182 2 pyc1, pzc1, pyc2, pzc2,
183 3 ay, r22, r23, dd,
184 4 dg, g33, k12, is,
185 5 nel, jcvt)
187 1 py1, pz1, py3, pz3,
188 2 pyc1, pzc1, pyc3, pzc3,
189 3 ay, r22, r23, dd,
190 4 dg, g33, k13, is,
191 5 nel, jcvt)
193 1 py1, pz1, py4, pz4,
194 2 pyc1, pzc1, pyc4, pzc4,
195 3 ay, r22, r23, dd,
196 4 dg, g33, k14, is,
197 5 nel, jcvt)
199 1 py2, pz2, py3, pz3,
200 2 pyc2, pzc2, pyc3, pzc3,
201 3 ay, r22, r23, dd,
202 4 dg, g33, k23, is,
203 5 nel, jcvt)
205 1 py2, pz2, py4, pz4,
206 2 pyc2, pzc2, pyc4, pzc4,
207 3 ay, r22, r23, dd,
208 4 dg, g33, k24, is,
209 5 nel, jcvt)
211 1 py3, pz3, py4, pz4,
212 2 pyc3, pzc3, pyc4, pzc4,
213 3 ay, r22, r23, dd,
214 4 dg, g33, k34, is,
215 5 nel, jcvt)
216 ENDIF
217
218 RETURN
subroutine q4kelij2(pyi, pzi, pyj, pzj, pyci, pzci, pycj, pzcj, ay, r22, r23, dd, dg, g33, kij, is, nel, jcvt)
subroutine q4kelijs2(pyi, pzi, pyj, pzj, pyci, pzci, pycj, pzcj, ay, r22, r23, dd, gg, kij, kiju, kijl, is, nel, jcvt)