OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
scoor3.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "vect01_c.inc"
#include "scr05_c.inc"
#include "scr18_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine scoor3 (x, ixs, v, w, gama0, gama, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, vdx1, vdx2, vdx3, vdx4, vdx5, vdx6, vdx7, vdx8, vdy1, vdy2, vdy3, vdy4, vdy5, vdy6, vdy7, vdy8, vdz1, vdz2, vdz3, vdz4, vdz5, vdz6, vdz7, vdz8, vdx, vdy, vdz, vd2, vis, offg, off, sav, rho, rhoo, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, ngl, mxt, ngeo, vr, vxr, vyr, vzr, fr_wave, fr_wav, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, xdp, iparg, ng, nel)

Function/Subroutine Documentation

◆ scoor3()

subroutine scoor3 ( x,
integer, dimension(nixs,*) ixs,
v,
w,
gama0,
gama,
x1,
x2,
x3,
x4,
x5,
x6,
x7,
x8,
y1,
y2,
y3,
y4,
y5,
y6,
y7,
y8,
z1,
z2,
z3,
z4,
z5,
z6,
z7,
z8,
vx1,
vx2,
vx3,
vx4,
vx5,
vx6,
vx7,
vx8,
vy1,
vy2,
vy3,
vy4,
vy5,
vy6,
vy7,
vy8,
vz1,
vz2,
vz3,
vz4,
vz5,
vz6,
vz7,
vz8,
vdx1,
vdx2,
vdx3,
vdx4,
vdx5,
vdx6,
vdx7,
vdx8,
vdy1,
vdy2,
vdy3,
vdy4,
vdy5,
vdy6,
vdy7,
vdy8,
vdz1,
vdz2,
vdz3,
vdz4,
vdz5,
vdz6,
vdz7,
vdz8,
vdx,
vdy,
vdz,
vd2,
vis,
offg,
off,
double precision, dimension(nel,21) sav,
rho,
rhoo,
integer, dimension(*) nc1,
integer, dimension(*) nc2,
integer, dimension(*) nc3,
integer, dimension(*) nc4,
integer, dimension(*) nc5,
integer, dimension(*) nc6,
integer, dimension(*) nc7,
integer, dimension(*) nc8,
integer, dimension(*) ngl,
integer, dimension(*) mxt,
integer, dimension(*) ngeo,
vr,
vxr,
vyr,
vzr,
fr_wave,
fr_wav,
double precision, dimension(*) xd1,
double precision, dimension(*) xd2,
double precision, dimension(*) xd3,
double precision, dimension(*) xd4,
double precision, dimension(*) xd5,
double precision, dimension(*) xd6,
double precision, dimension(*) xd7,
double precision, dimension(*) xd8,
double precision, dimension(*) yd1,
double precision, dimension(*) yd2,
double precision, dimension(*) yd3,
double precision, dimension(*) yd4,
double precision, dimension(*) yd5,
double precision, dimension(*) yd6,
double precision, dimension(*) yd7,
double precision, dimension(*) yd8,
double precision, dimension(*) zd1,
double precision, dimension(*) zd2,
double precision, dimension(*) zd3,
double precision, dimension(*) zd4,
double precision, dimension(*) zd5,
double precision, dimension(*) zd6,
double precision, dimension(*) zd7,
double precision, dimension(*) zd8,
double precision, dimension(3,*) xdp,
integer, dimension(nparg,*) iparg,
integer ng,
integer nel )

Definition at line 35 of file scoor3.F.

52C-----------------------------------------------
53C I m p l i c i t T y p e s
54C-----------------------------------------------
55#include "implicit_f.inc"
56C-----------------------------------------------
57C G l o b a l P a r a m e t e r s
58C-----------------------------------------------
59#include "mvsiz_p.inc"
60C-----------------------------------------------
61C C o m m o n B l o c k s
62C-----------------------------------------------
63#include "vect01_c.inc"
64#include "scr05_c.inc"
65#include "scr18_c.inc"
66#include "param_c.inc"
67C-----------------------------------------------
68C D u m m y A r g u m e n t s
69C-----------------------------------------------
70 INTEGER IXS(NIXS,*),IPARG(NPARG,*),NG,NEL
71 my_real
72 . x(3,*),v(3,*),w(3,*), vis(*),vr(3,*),vxr(*),vyr(*),vzr(*),
73 . x1(*), x2(*), x3(*), x4(*), x5(*), x6(*), x7(*), x8(*),
74 . y1(*), y2(*), y3(*), y4(*), y5(*), y6(*), y7(*), y8(*),
75 . z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), z7(*), z8(*),
76 . vx1(*), vx2(*), vx3(*), vx4(*), vx5(*), vx6(*), vx7(*), vx8(*),
77 . vy1(*), vy2(*), vy3(*), vy4(*), vy5(*), vy6(*), vy7(*), vy8(*),
78 . vz1(*), vz2(*), vz3(*), vz4(*), vz5(*), vz6(*), vz7(*), vz8(*),
79 . vdx1(*),vdx2(*),vdx3(*),vdx4(*),vdx5(*),vdx6(*),vdx7(*),vdx8(*),
80 . vdy1(*),vdy2(*),vdy3(*),vdy4(*),vdy5(*),vdy6(*),vdy7(*),vdy8(*),
81 . vdz1(*),vdz2(*),vdz3(*),vdz4(*),vdz5(*),vdz6(*),vdz7(*),vdz8(*),
82 . vdx(*), vdy(*), vdz(*),vd2(*),offg(*),off(*),rho(*),
83 . rhoo(*),gama0(nel,6),gama(mvsiz,6),fr_wave(*),fr_wav(*)
84
85 INTEGER NC1(*), NC2(*), NC3(*), NC4(*),
86 . NC5(*), NC6(*), NC7(*), NC8(*), MXT(*), NGL(*),NGEO(*)
87
88 double precision
89 . xdp(3,*),sav(nel,21),
90 . xd1(*), xd2(*), xd3(*), xd4(*), xd5(*), xd6(*), xd7(*), xd8(*),
91 . yd1(*), yd2(*), yd3(*), yd4(*), yd5(*), yd6(*), yd7(*), yd8(*),
92 . zd1(*), zd2(*), zd3(*), zd4(*), zd5(*), zd6(*), zd7(*), zd8(*)
93C-----------------------------------------------
94C L o c a l V a r i a b l e s
95C-----------------------------------------------
96 INTEGER I, J,MXT_1
97 my_real, DIMENSION(MVSIZ) ::
98 . rx, ry, rz, sx, sy, sz, tx, ty, tz,
99 . e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z
100 my_real
101 . off_l
102C=======================================================================
103 mxt_1 = ixs(1,1)
104 DO i=1,nel
105 vis(i)=zero
106 ngeo(i)=ixs(10,i)
107 ngl(i)=ixs(11,i)
108 mxt(i)=mxt_1
109 nc1(i)=ixs(2,i)
110 nc2(i)=ixs(3,i)
111 nc3(i)=ixs(4,i)
112 nc4(i)=ixs(5,i)
113 nc5(i)=ixs(6,i)
114 nc6(i)=ixs(7,i)
115 nc7(i)=ixs(8,i)
116 nc8(i)=ixs(9,i)
117 rhoo(i)=rho(i)
118 ENDDO
119
120 off_l = zero
121 !----------------------------
122 ! NODAL COORDINATES |
123 !----------------------------
124 IF(iresp==1)THEN
125 DO i=1,nel
126 xd1(i)=xdp(1,nc1(i))
127 yd1(i)=xdp(2,nc1(i))
128 zd1(i)=xdp(3,nc1(i))
129 xd2(i)=xdp(1,nc2(i))
130 yd2(i)=xdp(2,nc2(i))
131 zd2(i)=xdp(3,nc2(i))
132 xd3(i)=xdp(1,nc3(i))
133 yd3(i)=xdp(2,nc3(i))
134 zd3(i)=xdp(3,nc3(i))
135 xd4(i)=xdp(1,nc4(i))
136 yd4(i)=xdp(2,nc4(i))
137 zd4(i)=xdp(3,nc4(i))
138 xd5(i)=xdp(1,nc5(i))
139 yd5(i)=xdp(2,nc5(i))
140 zd5(i)=xdp(3,nc5(i))
141 xd6(i)=xdp(1,nc6(i))
142 yd6(i)=xdp(2,nc6(i))
143 zd6(i)=xdp(3,nc6(i))
144 xd7(i)=xdp(1,nc7(i))
145 yd7(i)=xdp(2,nc7(i))
146 zd7(i)=xdp(3,nc7(i))
147 xd8(i)=xdp(1,nc8(i))
148 yd8(i)=xdp(2,nc8(i))
149 zd8(i)=xdp(3,nc8(i))
150 ENDDO
151 ELSE
152 DO i=1,nel
153 xd1(i)=x(1,nc1(i))
154 yd1(i)=x(2,nc1(i))
155 zd1(i)=x(3,nc1(i))
156 xd2(i)=x(1,nc2(i))
157 yd2(i)=x(2,nc2(i))
158 zd2(i)=x(3,nc2(i))
159 xd3(i)=x(1,nc3(i))
160 yd3(i)=x(2,nc3(i))
161 zd3(i)=x(3,nc3(i))
162 xd4(i)=x(1,nc4(i))
163 yd4(i)=x(2,nc4(i))
164 zd4(i)=x(3,nc4(i))
165 xd5(i)=x(1,nc5(i))
166 yd5(i)=x(2,nc5(i))
167 zd5(i)=x(3,nc5(i))
168 xd6(i)=x(1,nc6(i))
169 yd6(i)=x(2,nc6(i))
170 zd6(i)=x(3,nc6(i))
171 xd7(i)=x(1,nc7(i))
172 yd7(i)=x(2,nc7(i))
173 zd7(i)=x(3,nc7(i))
174 xd8(i)=x(1,nc8(i))
175 yd8(i)=x(2,nc8(i))
176 zd8(i)=x(3,nc8(i))
177 ENDDO
178 ENDIF
179
180 CALL srepiso3(
181 1 xd1, xd2, xd3, xd4,
182 2 xd5, xd6, xd7, xd8,
183 3 yd1, yd2, yd3, yd4,
184 4 yd5, yd6, yd7, yd8,
185 5 zd1, zd2, zd3, zd4,
186 6 zd5, zd6, zd7, zd8,
187 7 rx, ry, rz, sx,
188 8 sy, sz, tx, ty,
189 9 tz, nel)
190 CALL sreploc3(
191 1 rx, ry, rz, sx,
192 2 sy, sz, tx, ty,
193 3 tz, e1x, e2x, e3x,
194 4 e1y, e2y, e3y, e1z,
195 5 e2z, e3z, llt)
196 IF (isorth == 0) THEN
197 DO i=1,nel
198 gama(i,1) = one
199 gama(i,2) = zero
200 gama(i,3) = zero
201 gama(i,4) = zero
202 gama(i,5) = one
203 gama(i,6) = zero
204 ENDDO
205 ELSE
206 CALL sorthdir3(
207 1 rx, ry, rz, sx,
208 2 sy, sz, tx, ty,
209 3 tz, e1x, e2x, e3x,
210 4 e1y, e2y, e3y, e1z,
211 5 e2z, e3z, gama0, gama,
212 6 nel, irep)
213 ENDIF
214 IF(ismstr<=4.AND.jlag>0.OR.(ismstr==12.AND.idtmin(1)==3)) THEN
215C
216 DO i=1,nel
217 IF(abs(offg(i))>one)THEN
218 xd1(i)=sav(i,1)
219 yd1(i)=sav(i,2)
220 zd1(i)=sav(i,3)
221 xd2(i)=sav(i,4)
222 yd2(i)=sav(i,5)
223 zd2(i)=sav(i,6)
224 xd3(i)=sav(i,7)
225 yd3(i)=sav(i,8)
226 zd3(i)=sav(i,9)
227 xd4(i)=sav(i,10)
228 yd4(i)=sav(i,11)
229 zd4(i)=sav(i,12)
230 xd5(i)=sav(i,13)
231 yd5(i)=sav(i,14)
232 zd5(i)=sav(i,15)
233 xd6(i)=sav(i,16)
234 yd6(i)=sav(i,17)
235 zd6(i)=sav(i,18)
236 xd7(i)=sav(i,19)
237 yd7(i)=sav(i,20)
238 zd7(i)=sav(i,21)
239 xd8(i)=zero
240 yd8(i)=zero
241 zd8(i)=zero
242 off(i) = abs(offg(i))-one
243 off_l = min(off_l,offg(i))
244 ELSE
245 off(i) = abs(offg(i))
246 off_l = min(off_l,offg(i))
247 ENDIF
248 ENDDO
249
250 ELSE
251 DO i=1,nel
252 off(i) = abs(offg(i))
253 off_l = min(off_l,offg(i))
254 ENDDO
255 ENDIF
256
257 !copy and cast XD (DP) to X (SP) to assure coherence between XD and X
258 DO i=1,nel
259 x1(i)= xd1(i)
260 y1(i)= yd1(i)
261 z1(i)= zd1(i)
262 x2(i)= xd2(i)
263 y2(i)= yd2(i)
264 z2(i)= zd2(i)
265 x3(i)= xd3(i)
266 y3(i)= yd3(i)
267 z3(i)= zd3(i)
268 x4(i)= xd4(i)
269 y4(i)= yd4(i)
270 z4(i)= zd4(i)
271 x5(i)= xd5(i)
272 y5(i)= yd5(i)
273 z5(i)= zd5(i)
274 x6(i)= xd6(i)
275 y6(i)= yd6(i)
276 z6(i)= zd6(i)
277 x7(i)= xd7(i)
278 y7(i)= yd7(i)
279 z7(i)= zd7(i)
280 x8(i)= xd8(i)
281 y8(i)= yd8(i)
282 z8(i)= zd8(i)
283 ENDDO
284
285 DO i=1,nel
286 vx1(i)=v(1,nc1(i))
287 vy1(i)=v(2,nc1(i))
288 vz1(i)=v(3,nc1(i))
289 vx2(i)=v(1,nc2(i))
290 vy2(i)=v(2,nc2(i))
291 vz2(i)=v(3,nc2(i))
292 vx3(i)=v(1,nc3(i))
293 vy3(i)=v(2,nc3(i))
294 vz3(i)=v(3,nc3(i))
295 vx4(i)=v(1,nc4(i))
296 vy4(i)=v(2,nc4(i))
297 vz4(i)=v(3,nc4(i))
298 vx5(i)=v(1,nc5(i))
299 vy5(i)=v(2,nc5(i))
300 vz5(i)=v(3,nc5(i))
301 vx6(i)=v(1,nc6(i))
302 vy6(i)=v(2,nc6(i))
303 vz6(i)=v(3,nc6(i))
304 vx7(i)=v(1,nc7(i))
305 vy7(i)=v(2,nc7(i))
306 vz7(i)=v(3,nc7(i))
307 vx8(i)=v(1,nc8(i))
308 vy8(i)=v(2,nc8(i))
309 vz8(i)=v(3,nc8(i))
310 ENDDO
311
312 IF(isrot/=0)THEN
313 DO i=1,nel
314 vxr(i)=(vr(1,nc1(i))+vr(1,nc2(i))+vr(1,nc3(i))+vr(1,nc4(i)) + vr(1,nc5(i))+vr(1,nc6(i))+vr(1,nc7(i))+vr(1,nc8(i)))/eight
315 vyr(i)=(vr(2,nc1(i))+vr(2,nc2(i))+vr(2,nc3(i))+vr(2,nc4(i)) + vr(2,nc5(i))+vr(2,nc6(i))+vr(2,nc7(i))+vr(2,nc8(i)))/eight
316 vzr(i)=(vr(3,nc1(i))+vr(3,nc2(i))+vr(3,nc3(i))+vr(3,nc4(i)) + vr(3,nc5(i))+vr(3,nc6(i))+vr(3,nc7(i))+vr(3,nc8(i)))/eight
317 fr_wav(i)=max(
318 + fr_wave(nc1(i)),fr_wave(nc2(i)),fr_wave(nc3(i)),fr_wave(nc4(i)),
319 + fr_wave(nc5(i)),fr_wave(nc6(i)),fr_wave(nc7(i)),fr_wave(nc8(i)))
320 ENDDO
321 ENDIF
322
323 IF(off_l<zero)THEN
324 DO i=1,nel
325 IF(offg(i)<zero)THEN
326 vx1(i)=zero
327 vy1(i)=zero
328 vz1(i)=zero
329 vx2(i)=zero
330 vy2(i)=zero
331 vz2(i)=zero
332 vx3(i)=zero
333 vy3(i)=zero
334 vz3(i)=zero
335 vx4(i)=zero
336 vy4(i)=zero
337 vz4(i)=zero
338 vx5(i)=zero
339 vy5(i)=zero
340 vz5(i)=zero
341 vx6(i)=zero
342 vy6(i)=zero
343 vz6(i)=zero
344 vx7(i)=zero
345 vy7(i)=zero
346 vz7(i)=zero
347 vx8(i)=zero
348 vy8(i)=zero
349 vz8(i)=zero
350 ENDIF
351 ENDDO
352 ENDIF
353
354 IF (jlag/=0)THEN
355 DO i=1,nel
356 vd2(i)=zero
357 vdx1(i)=zero
358 vdy1(i)=zero
359 vdz1(i)=zero
360 vdx2(i)=zero
361 vdy2(i)=zero
362 vdz2(i)=zero
363 vdx3(i)=zero
364 vdy3(i)=zero
365 vdz3(i)=zero
366 vdx4(i)=zero
367 vdy4(i)=zero
368 vdz4(i)=zero
369 vdx5(i)=zero
370 vdy5(i)=zero
371 vdz5(i)=zero
372 vdx6(i)=zero
373 vdy6(i)=zero
374 vdz6(i)=zero
375 vdx7(i)=zero
376 vdy7(i)=zero
377 vdz7(i)=zero
378 vdx8(i)=zero
379 vdy8(i)=zero
380 vdz8(i)=zero
381 ENDDO
382 RETURN
383 ELSEIF(jale/=0)THEN
384 IF(iparg(64,ng)==1)THEN
385 DO i=1,nel
386 vd2(i) =zero
387 vdx1(i)=zero
388 vdy1(i)=zero
389 vdz1(i)=zero
390 vdx2(i)=zero
391 vdy2(i)=zero
392 vdz2(i)=zero
393 vdx3(i)=zero
394 vdy3(i)=zero
395 vdz3(i)=zero
396 vdx4(i)=zero
397 vdy4(i)=zero
398 vdz4(i)=zero
399 vdx5(i)=zero
400 vdy5(i)=zero
401 vdz5(i)=zero
402 vdx6(i)=zero
403 vdy6(i)=zero
404 vdz6(i)=zero
405 vdx7(i)=zero
406 vdy7(i)=zero
407 vdz7(i)=zero
408 vdx8(i)=zero
409 vdy8(i)=zero
410 vdz8(i)=zero
411 ENDDO
412 RETURN
413 ENDIF
414 DO i=1,nel
415 vdx1(i)=vx1(i)-w(1,nc1(i))
416 vdy1(i)=vy1(i)-w(2,nc1(i))
417 vdz1(i)=vz1(i)-w(3,nc1(i))
418 vdx2(i)=vx2(i)-w(1,nc2(i))
419 vdy2(i)=vy2(i)-w(2,nc2(i))
420 vdz2(i)=vz2(i)-w(3,nc2(i))
421 vdx3(i)=vx3(i)-w(1,nc3(i))
422 vdy3(i)=vy3(i)-w(2,nc3(i))
423 vdz3(i)=vz3(i)-w(3,nc3(i))
424 vdx4(i)=vx4(i)-w(1,nc4(i))
425 vdy4(i)=vy4(i)-w(2,nc4(i))
426 vdz4(i)=vz4(i)-w(3,nc4(i))
427 vdx5(i)=vx5(i)-w(1,nc5(i))
428 vdy5(i)=vy5(i)-w(2,nc5(i))
429 vdz5(i)=vz5(i)-w(3,nc5(i))
430 vdx6(i)=vx6(i)-w(1,nc6(i))
431 vdy6(i)=vy6(i)-w(2,nc6(i))
432 vdz6(i)=vz6(i)-w(3,nc6(i))
433 vdx7(i)=vx7(i)-w(1,nc7(i))
434 vdy7(i)=vy7(i)-w(2,nc7(i))
435 vdz7(i)=vz7(i)-w(3,nc7(i))
436 vdx8(i)=vx8(i)-w(1,nc8(i))
437 vdy8(i)=vy8(i)-w(2,nc8(i))
438 vdz8(i)=vz8(i)-w(3,nc8(i))
439 ENDDO
440 ELSEIF(jeul/=0)THEN
441 IF(iparg(64,ng)==1)THEN
442 DO i=1,nel
443 vd2(i) =zero
444 vdx1(i)=zero
445 vdy1(i)=zero
446 vdz1(i)=zero
447 vdx2(i)=zero
448 vdy2(i)=zero
449 vdz2(i)=zero
450 vdx3(i)=zero
451 vdy3(i)=zero
452 vdz3(i)=zero
453 vdx4(i)=zero
454 vdy4(i)=zero
455 vdz4(i)=zero
456 vdx5(i)=zero
457 vdy5(i)=zero
458 vdz5(i)=zero
459 vdx6(i)=zero
460 vdy6(i)=zero
461 vdz6(i)=zero
462 vdx7(i)=zero
463 vdy7(i)=zero
464 vdz7(i)=zero
465 vdx8(i)=zero
466 vdy8(i)=zero
467 vdz8(i)=zero
468 ENDDO
469 RETURN
470 ENDIF
471 DO i=1,nel
472 vdx1(i)=vx1(i)
473 vdy1(i)=vy1(i)
474 vdz1(i)=vz1(i)
475 vdx2(i)=vx2(i)
476 vdy2(i)=vy2(i)
477 vdz2(i)=vz2(i)
478 vdx3(i)=vx3(i)
479 vdy3(i)=vy3(i)
480 vdz3(i)=vz3(i)
481 vdx4(i)=vx4(i)
482 vdy4(i)=vy4(i)
483 vdz4(i)=vz4(i)
484 vdx5(i)=vx5(i)
485 vdy5(i)=vy5(i)
486 vdz5(i)=vz5(i)
487 vdx6(i)=vx6(i)
488 vdy6(i)=vy6(i)
489 vdz6(i)=vz6(i)
490 vdx7(i)=vx7(i)
491 vdy7(i)=vy7(i)
492 vdz7(i)=vz7(i)
493 vdx8(i)=vx8(i)
494 vdy8(i)=vy8(i)
495 vdz8(i)=vz8(i)
496 ENDDO
497 ENDIF
498
499 DO i=1,nel
500 vdx(i)=one_over_8*(vdx1(i)+vdx2(i)+vdx3(i)+vdx4(i)+vdx5(i)+vdx6(i)+vdx7(i)+vdx8(i))
501 vdy(i)=one_over_8*(vdy1(i)+vdy2(i)+vdy3(i)+vdy4(i)+vdy5(i)+vdy6(i)+vdy7(i)+vdy8(i))
502 vdz(i)=one_over_8*(vdz1(i)+vdz2(i)+vdz3(i)+vdz4(i)+vdz5(i)+vdz6(i)+vdz7(i)+vdz8(i))
503 vd2(i)=(vdx(i)**2+vdy(i)**2+vdz(i)**2)
504 ENDDO
505C-----------
506 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
sorthdir3
subroutine sorthdir3(rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, gama0, gama, nel, irep)
Definition sorthdir3.F:42
srepiso3
subroutine srepiso3(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, f1x, f1y, f1z, f2x, f2y, f2z)
Definition srepiso3.F:35
sreploc3
subroutine sreploc3(rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition sreploc3.F:32