66
67
68
69#include "implicit_f.inc"
70#include "comlock.inc"
71
72
73
74#include "scr17_c.inc"
75
76
77
78 INTEGER, INTENT(IN) :: NEL
79 INTEGER ICP,IDEG(*)
80
82 . off(*),offg(*),vol(*),ksi,eta,zeta,wi,
83 . px1(*), px2(*), px3(*), px4(*),
84 . px5(*), px6(*), px7(*), px8(*),
85 . py1(*), py2(*), py3(*), py4(*),
86 . py5(*), py6(*), py7(*), py8(*),
87 . pz1(*), pz2(*), pz3(*), pz4(*),
88 . pz5(*), pz6(*), pz7(*), pz8(*),
89 . aj4(*),aj5(*),aj6(*),
90 . aj7(*),aj8(*),aj9(*),
91 . rx(*) ,ry(*) ,rz(*) ,
92 . sx(*) ,sy(*) ,sz(*) ,
93 . smax(*),deltax(*),
94 . pxc1(*), pxc2(*), pxc3(*), pxc4(*),
95 . pyc1(*), pyc2(*), pyc3(*), pyc4(*),
96 . pzc1(*), pzc2(*), pzc3(*), pzc4(*),
97 . bxy1(*),bxy2(*),bxy3(*),bxy4(*),
98 . bxy5(*),bxy6(*),bxy7(*),bxy8(*),
99 . byx1(*),byx2(*),byx3(*),byx4(*),
100 . byx5(*),byx6(*),byx7(*),byx8(*),
101 . bxz1(*),bxz2(*),bxz3(*),bxz4(*),
102 . bxz5(*),bxz6(*),bxz7(*),bxz8(*),
103 . bzx1(*),bzx2(*),bzx3(*),bzx4(*),
104 . bzx5(*),bzx6(*),bzx7(*),bzx8(*),
105 . byz1(*),byz2(*),byz3(*),byz4(*),
106 . byz5(*),byz6(*),byz7(*),byz8(*),
107 . bzy1(*),bzy2(*),bzy3(*),bzy4(*),
108 . bzy5(*),bzy6(*),bzy7(*),bzy8(*),
109 . bxx1(*),bxx2(*),bxx3(*),bxx4(*),
110 . bxx5(*),bxx6(*),bxx7(*),bxx8(*),
111 . byy1(*),byy2(*),byy3(*),byy4(*),
112 . byy5(*),byy6(*),byy7(*),byy8(*),
113 . bzz1(*),bzz2(*),bzz3(*),bzz4(*),
114 . bzz5(*),bzz6(*),bzz7(*),bzz8(*),nu(*)
115 DOUBLE PRECISION
116 . VOLP(*)
117
118
119
120 INTEGER I, J ,ICOR
121
122
124 . a1pr1,a1pr3,a1pr5,a1pr7,
125 . a4pr1,a4pr3,a4pr5,a4pr7,
126 . a7pr1,a7pr3,a7pr5,a7pr7,
127 . a2ps1,a2ps2,a2ps5,a2ps6,
128 . a5ps1,a5ps2,a5ps5,a5ps6,
129 . a8ps1,a8ps2,a8ps5,a8ps6,
130 . a3pt1,a3pt2,a3pt3,a3pt4,
131 . a6pt1,a6pt2,a6pt3,a6pt4,
132 . a9pt1,a9pt2,a9pt3,a9pt4,
133 . bxhi,byhi,bzhi,tmp,deti
134
135 DO i=1,nel
136 rx(i) = aj4(i)
137 ry(i) = aj5(i)
138 rz(i) = aj6(i)
139 sx(i) = aj7(i)
140 sy(i) = aj8(i)
141 sz(i) = aj9(i)
142 ENDDO
143
144 DO i=1,nel
145 off(i) = offg(i)
146 vol(i) = volp(i)
147 ENDDO
148 IF (idts6==0) THEN
149 DO i=1,nel
150 deltax(i)=
min(deltax(i),hundred28*vol(i)*smax(i)/wi)
151 ENDDO
152 ELSE
153 DO i=1,nel
154 IF (ideg(i)==0) deltax(i)=
min(deltax(i),hundred28*vol(i)*smax(i)/wi)
155 ENDDO
156 END IF
157
158
159
160
161 IF (icp==11) THEN
163 1 px1, px2, px3, px4,
164 2 px5, px6, px7, px8,
165 3 py1, py2, py3, py4,
166 4 py5, py6, py7, py8,
167 5 pz1, pz2, pz3, pz4,
168 6 pz5, pz6, pz7, pz8,
169 7 pxc1, pxc2, pxc3, pxc4,
170 8 pyc1, pyc2, pyc3, pyc4,
171 9 pzc1, pzc2, pzc3, pzc4,
172 a bxy1, bxy2, bxy3, bxy4,
173 b bxy5, bxy6, bxy7, bxy8,
174 c byx1, byx2, byx3, byx4,
175 d byx5, byx6, byx7, byx8,
176 e bxz1, bxz2, bxz3, bxz4,
177 f bxz5, bxz6, bxz7, bxz8,
178 g bzx1, bzx2, bzx3, bzx4,
179 h bzx5, bzx6, bzx7, bzx8,
180 i byz1, byz2, byz3, byz4,
181 j byz5, byz6, byz7, byz8,
182 k bzy1, bzy2, bzy3, bzy4,
183 l bzy5, bzy6, bzy7, bzy8,
184 m bxx1, bxx2, bxx3, bxx4,
185 n bxx5, bxx6, bxx7, bxx8,
186 o byy1, byy2, byy3, byy4,
187 p byy5, byy6, byy7, byy8,
188 q bzz1, bzz2, bzz3, bzz4,
189 r bzz5, bzz6, bzz7, bzz8,
190 s nu, nel)
191 ENDIF
192
193 RETURN
subroutine s8bicp11(px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, bxx1, bxx2, bxx3, bxx4, bxx5, bxx6, bxx7, bxx8, byy1, byy2, byy3, byy4, byy5, byy6, byy7, byy8, bzz1, bzz2, bzz3, bzz4, bzz5, bzz6, bzz7, bzz8, nu, nel)