OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
scoor3.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| scoor3 ../engine/source/elements/solid/solide/scoor3.F
25!||--- called by ------------------------------------------------------
26!|| s8cforc3 ../engine/source/elements/thickshell/solide8c/s8cforc3.F
27!|| s8eforc3 ../engine/source/elements/solid/solide8e/s8eforc3.F
28!|| s8zforc3 ../engine/source/elements/solid/solide8z/s8zforc3.F
29!|| sforc3 ../engine/source/elements/solid/solide/sforc3.F
30!||--- calls -----------------------------------------------------
31!|| sorthdir3 ../engine/source/elements/solid/solide/sorthdir3.F
32!|| srepiso3 ../engine/source/elements/solid/solide/srepiso3.F
33!|| sreploc3 ../engine/source/elements/solid/solide/sreploc3.F
34!||--- uses -----------------------------------------------------
35!|| element_mod ../common_source/modules/elements/element_mod.F90
36!||====================================================================
37 SUBROUTINE scoor3(X,IXS,V,W,GAMA0,GAMA,
38 . X1, X2, X3, X4, X5, X6, X7, X8,
39 . Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8,
40 . Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8,
41 . VX1, VX2, VX3, VX4, VX5, VX6, VX7, VX8,
42 . VY1, VY2, VY3, VY4, VY5, VY6, VY7, VY8,
43 . VZ1, VZ2, VZ3, VZ4, VZ5, VZ6, VZ7, VZ8,
44 . VDX1, VDX2, VDX3, VDX4, VDX5, VDX6, VDX7, VDX8,
45 . VDY1, VDY2, VDY3, VDY4, VDY5, VDY6, VDY7, VDY8,
46 . VDZ1, VDZ2, VDZ3, VDZ4, VDZ5, VDZ6, VDZ7, VDZ8,
47 . VDX,VDY,VDZ,VD2,VIS,OFFG,OFF,SAV,RHO,
48 . RHOO,NC1,NC2,NC3,NC4,NC5,NC6,NC7,NC8,NGL,MXT,NGEO,
49 . VR,VXR,VYR,VZR,FR_WAVE,FR_WAV,
50 . XD1, XD2 , XD3, XD4, XD5, XD6, XD7, XD8,
51 . YD1, YD2 , YD3, YD4, YD5, YD6, YD7, YD8,
52 . ZD1, ZD2 , ZD3, ZD4, ZD5, ZD6, ZD7, ZD8,
53 . XDP, IPARG, NG , NEL)
54 use element_mod , only : nixs
55C-----------------------------------------------
56C I m p l i c i t T y p e s
57C-----------------------------------------------
58#include "implicit_f.inc"
59C-----------------------------------------------
60C G l o b a l P a r a m e t e r s
61C-----------------------------------------------
62#include "mvsiz_p.inc"
63C-----------------------------------------------
64C C o m m o n B l o c k s
65C-----------------------------------------------
66#include "vect01_c.inc"
67#include "scr05_c.inc"
68#include "scr18_c.inc"
69#include "param_c.inc"
70C-----------------------------------------------
71C D u m m y A r g u m e n t s
72C-----------------------------------------------
73 INTEGER IXS(NIXS,*),IPARG(NPARG,*),NG,NEL
74 my_real
75 . X(3,*),V(3,*),W(3,*), VIS(*),VR(3,*),VXR(*),VYR(*),VZR(*),
76 . X1(*), X2(*), X3(*), X4(*), X5(*), X6(*), X7(*), X8(*),
77 . Y1(*), Y2(*), Y3(*), Y4(*), Y5(*), Y6(*), Y7(*), Y8(*),
78 . Z1(*), Z2(*), Z3(*), Z4(*), Z5(*), Z6(*), Z7(*), Z8(*),
79 . VX1(*), VX2(*), VX3(*), VX4(*), VX5(*), VX6(*), VX7(*), VX8(*),
80 . VY1(*), VY2(*), VY3(*), VY4(*), VY5(*), VY6(*), VY7(*), VY8(*),
81 . VZ1(*), VZ2(*), VZ3(*), VZ4(*), VZ5(*), VZ6(*), VZ7(*), VZ8(*),
82 . VDX1(*),VDX2(*),VDX3(*),VDX4(*),VDX5(*),VDX6(*),VDX7(*),VDX8(*),
83 . VDY1(*),VDY2(*),VDY3(*),VDY4(*),VDY5(*),VDY6(*),VDY7(*),VDY8(*),
84 . VDZ1(*),VDZ2(*),VDZ3(*),VDZ4(*),VDZ5(*),VDZ6(*),VDZ7(*),VDZ8(*),
85 . VDX(*), VDY(*), VDZ(*),VD2(*),OFFG(*),OFF(*),RHO(*),
86 . RHOO(*),GAMA0(NEL,6),GAMA(MVSIZ,6),FR_WAVE(*),FR_WAV(*)
87
88 INTEGER NC1(*), NC2(*), NC3(*), NC4(*),
89 . nc5(*), nc6(*), nc7(*), nc8(*), mxt(*), ngl(*),ngeo(*)
90
91 double precision
92 . xdp(3,*),sav(nel,21),
93 . xd1(*), xd2(*), xd3(*), xd4(*), xd5(*), xd6(*), xd7(*), xd8(*),
94 . yd1(*), yd2(*), yd3(*), yd4(*), yd5(*), yd6(*), yd7(*), yd8(*),
95 . zd1(*), zd2(*), zd3(*), zd4(*), zd5(*), zd6(*), zd7(*), zd8(*)
96C-----------------------------------------------
97C L o c a l V a r i a b l e s
98C-----------------------------------------------
99 INTEGER I,MXT_1
100 my_real, DIMENSION(MVSIZ) ::
101 . RX, RY, RZ, SX, SY, SZ, TX, TY, TZ,
102 . E1X, E1Y, E1Z, E2X, E2Y, E2Z, E3X, E3Y, E3Z
103 my_real
104 . off_l
105C=======================================================================
106 mxt_1 = ixs(1,1)
107 DO i=1,nel
108 vis(i)=zero
109 ngeo(i)=ixs(10,i)
110 ngl(i)=ixs(11,i)
111 mxt(i)=mxt_1
112 nc1(i)=ixs(2,i)
113 nc2(i)=ixs(3,i)
114 nc3(i)=ixs(4,i)
115 nc4(i)=ixs(5,i)
116 nc5(i)=ixs(6,i)
117 nc6(i)=ixs(7,i)
118 nc7(i)=ixs(8,i)
119 nc8(i)=ixs(9,i)
120 rhoo(i)=rho(i)
121 ENDDO
122
123 off_l = zero
124 !----------------------------
125 ! NODAL COORDINATES |
126 !----------------------------
127 IF(iresp==1)THEN
128 DO i=1,nel
129 xd1(i)=xdp(1,nc1(i))
130 yd1(i)=xdp(2,nc1(i))
131 zd1(i)=xdp(3,nc1(i))
132 xd2(i)=xdp(1,nc2(i))
133 yd2(i)=xdp(2,nc2(i))
134 zd2(i)=xdp(3,nc2(i))
135 xd3(i)=xdp(1,nc3(i))
136 yd3(i)=xdp(2,nc3(i))
137 zd3(i)=xdp(3,nc3(i))
138 xd4(i)=xdp(1,nc4(i))
139 yd4(i)=xdp(2,nc4(i))
140 zd4(i)=xdp(3,nc4(i))
141 xd5(i)=xdp(1,nc5(i))
142 yd5(i)=xdp(2,nc5(i))
143 zd5(i)=xdp(3,nc5(i))
144 xd6(i)=xdp(1,nc6(i))
145 yd6(i)=xdp(2,nc6(i))
146 zd6(i)=xdp(3,nc6(i))
147 xd7(i)=xdp(1,nc7(i))
148 yd7(i)=xdp(2,nc7(i))
149 zd7(i)=xdp(3,nc7(i))
150 xd8(i)=xdp(1,nc8(i))
151 yd8(i)=xdp(2,nc8(i))
152 zd8(i)=xdp(3,nc8(i))
153 ENDDO
154 ELSE
155 DO i=1,nel
156 xd1(i)=x(1,nc1(i))
157 yd1(i)=x(2,nc1(i))
158 zd1(i)=x(3,nc1(i))
159 xd2(i)=x(1,nc2(i))
160 yd2(i)=x(2,nc2(i))
161 zd2(i)=x(3,nc2(i))
162 xd3(i)=x(1,nc3(i))
163 yd3(i)=x(2,nc3(i))
164 zd3(i)=x(3,nc3(i))
165 xd4(i)=x(1,nc4(i))
166 yd4(i)=x(2,nc4(i))
167 zd4(i)=x(3,nc4(i))
168 xd5(i)=x(1,nc5(i))
169 yd5(i)=x(2,nc5(i))
170 zd5(i)=x(3,nc5(i))
171 xd6(i)=x(1,nc6(i))
172 yd6(i)=x(2,nc6(i))
173 zd6(i)=x(3,nc6(i))
174 xd7(i)=x(1,nc7(i))
175 yd7(i)=x(2,nc7(i))
176 zd7(i)=x(3,nc7(i))
177 xd8(i)=x(1,nc8(i))
178 yd8(i)=x(2,nc8(i))
179 zd8(i)=x(3,nc8(i))
180 ENDDO
181 ENDIF
182
183 CALL srepiso3(
184 1 xd1, xd2, xd3, xd4,
185 2 xd5, xd6, xd7, xd8,
186 3 yd1, yd2, yd3, yd4,
187 4 yd5, yd6, yd7, yd8,
188 5 zd1, zd2, zd3, zd4,
189 6 zd5, zd6, zd7, zd8,
190 7 rx, ry, rz, sx,
191 8 sy, sz, tx, ty,
192 9 tz, nel)
193 CALL sreploc3(
194 1 rx, ry, rz, sx,
195 2 sy, sz, tx, ty,
196 3 tz, e1x, e2x, e3x,
197 4 e1y, e2y, e3y, e1z,
198 5 e2z, e3z, llt)
199 IF (isorth == 0) THEN
200 DO i=1,nel
201 gama(i,1) = one
202 gama(i,2) = zero
203 gama(i,3) = zero
204 gama(i,4) = zero
205 gama(i,5) = one
206 gama(i,6) = zero
207 ENDDO
208 ELSE
209 CALL sorthdir3(
210 1 rx, ry, rz, sx,
211 2 sy, sz, tx, ty,
212 3 tz, e1x, e2x, e3x,
213 4 e1y, e2y, e3y, e1z,
214 5 e2z, e3z, gama0, gama,
215 6 nel, irep)
216 ENDIF
217 IF(ismstr<=4.AND.jlag>0.OR.(ismstr==12.AND.idtmin(1)==3)) THEN
218C
219 DO i=1,nel
220 IF(abs(offg(i))>one)THEN
221 xd1(i)=sav(i,1)
222 yd1(i)=sav(i,2)
223 zd1(i)=sav(i,3)
224 xd2(i)=sav(i,4)
225 yd2(i)=sav(i,5)
226 zd2(i)=sav(i,6)
227 xd3(i)=sav(i,7)
228 yd3(i)=sav(i,8)
229 zd3(i)=sav(i,9)
230 xd4(i)=sav(i,10)
231 yd4(i)=sav(i,11)
232 zd4(i)=sav(i,12)
233 xd5(i)=sav(i,13)
234 yd5(i)=sav(i,14)
235 zd5(i)=sav(i,15)
236 xd6(i)=sav(i,16)
237 yd6(i)=sav(i,17)
238 zd6(i)=sav(i,18)
239 xd7(i)=sav(i,19)
240 yd7(i)=sav(i,20)
241 zd7(i)=sav(i,21)
242 xd8(i)=zero
243 yd8(i)=zero
244 zd8(i)=zero
245 off(i) = abs(offg(i))-one
246 off_l = min(off_l,offg(i))
247 ELSE
248 off(i) = abs(offg(i))
249 off_l = min(off_l,offg(i))
250 ENDIF
251 ENDDO
252
253 ELSE
254 DO i=1,nel
255 off(i) = abs(offg(i))
256 off_l = min(off_l,offg(i))
257 ENDDO
258 ENDIF
259
260 !copy and cast XD (DP) to X (SP) to assure coherence between XD and X
261 DO i=1,nel
262 x1(i)= xd1(i)
263 y1(i)= yd1(i)
264 z1(i)= zd1(i)
265 x2(i)= xd2(i)
266 y2(i)= yd2(i)
267 z2(i)= zd2(i)
268 x3(i)= xd3(i)
269 y3(i)= yd3(i)
270 z3(i)= zd3(i)
271 x4(i)= xd4(i)
272 y4(i)= yd4(i)
273 z4(i)= zd4(i)
274 x5(i)= xd5(i)
275 y5(i)= yd5(i)
276 z5(i)= zd5(i)
277 x6(i)= xd6(i)
278 y6(i)= yd6(i)
279 z6(i)= zd6(i)
280 x7(i)= xd7(i)
281 y7(i)= yd7(i)
282 z7(i)= zd7(i)
283 x8(i)= xd8(i)
284 y8(i)= yd8(i)
285 z8(i)= zd8(i)
286 ENDDO
287
288 DO i=1,nel
289 vx1(i)=v(1,nc1(i))
290 vy1(i)=v(2,nc1(i))
291 vz1(i)=v(3,nc1(i))
292 vx2(i)=v(1,nc2(i))
293 vy2(i)=v(2,nc2(i))
294 vz2(i)=v(3,nc2(i))
295 vx3(i)=v(1,nc3(i))
296 vy3(i)=v(2,nc3(i))
297 vz3(i)=v(3,nc3(i))
298 vx4(i)=v(1,nc4(i))
299 vy4(i)=v(2,nc4(i))
300 vz4(i)=v(3,nc4(i))
301 vx5(i)=v(1,nc5(i))
302 vy5(i)=v(2,nc5(i))
303 vz5(i)=v(3,nc5(i))
304 vx6(i)=v(1,nc6(i))
305 vy6(i)=v(2,nc6(i))
306 vz6(i)=v(3,nc6(i))
307 vx7(i)=v(1,nc7(i))
308 vy7(i)=v(2,nc7(i))
309 vz7(i)=v(3,nc7(i))
310 vx8(i)=v(1,nc8(i))
311 vy8(i)=v(2,nc8(i))
312 vz8(i)=v(3,nc8(i))
313 ENDDO
314
315 IF(isrot/=0)THEN
316 DO i=1,nel
317 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
318 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
319 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
320 fr_wav(i)=max(
321 + fr_wave(nc1(i)),fr_wave(nc2(i)),fr_wave(nc3(i)),fr_wave(nc4(i)),
322 + fr_wave(nc5(i)),fr_wave(nc6(i)),fr_wave(nc7(i)),fr_wave(nc8(i)))
323 ENDDO
324 ENDIF
325
326 IF(off_l<zero)THEN
327 DO i=1,nel
328 IF(offg(i)<zero)THEN
329 vx1(i)=zero
330 vy1(i)=zero
331 vz1(i)=zero
332 vx2(i)=zero
333 vy2(i)=zero
334 vz2(i)=zero
335 vx3(i)=zero
336 vy3(i)=zero
337 vz3(i)=zero
338 vx4(i)=zero
339 vy4(i)=zero
340 vz4(i)=zero
341 vx5(i)=zero
342 vy5(i)=zero
343 vz5(i)=zero
344 vx6(i)=zero
345 vy6(i)=zero
346 vz6(i)=zero
347 vx7(i)=zero
348 vy7(i)=zero
349 vz7(i)=zero
350 vx8(i)=zero
351 vy8(i)=zero
352 vz8(i)=zero
353 ENDIF
354 ENDDO
355 ENDIF
356
357 IF (jlag/=0)THEN
358 DO i=1,nel
359 vd2(i)=zero
360 vdx1(i)=zero
361 vdy1(i)=zero
362 vdz1(i)=zero
363 vdx2(i)=zero
364 vdy2(i)=zero
365 vdz2(i)=zero
366 vdx3(i)=zero
367 vdy3(i)=zero
368 vdz3(i)=zero
369 vdx4(i)=zero
370 vdy4(i)=zero
371 vdz4(i)=zero
372 vdx5(i)=zero
373 vdy5(i)=zero
374 vdz5(i)=zero
375 vdx6(i)=zero
376 vdy6(i)=zero
377 vdz6(i)=zero
378 vdx7(i)=zero
379 vdy7(i)=zero
380 vdz7(i)=zero
381 vdx8(i)=zero
382 vdy8(i)=zero
383 vdz8(i)=zero
384 ENDDO
385 RETURN
386 ELSEIF(jale/=0)THEN
387 IF(iparg(64,ng)==1)THEN
388 DO i=1,nel
389 vd2(i) =zero
390 vdx1(i)=zero
391 vdy1(i)=zero
392 vdz1(i)=zero
393 vdx2(i)=zero
394 vdy2(i)=zero
395 vdz2(i)=zero
396 vdx3(i)=zero
397 vdy3(i)=zero
398 vdz3(i)=zero
399 vdx4(i)=zero
400 vdy4(i)=zero
401 vdz4(i)=zero
402 vdx5(i)=zero
403 vdy5(i)=zero
404 vdz5(i)=zero
405 vdx6(i)=zero
406 vdy6(i)=zero
407 vdz6(i)=zero
408 vdx7(i)=zero
409 vdy7(i)=zero
410 vdz7(i)=zero
411 vdx8(i)=zero
412 vdy8(i)=zero
413 vdz8(i)=zero
414 ENDDO
415 RETURN
416 ENDIF
417 DO i=1,nel
418 vdx1(i)=vx1(i)-w(1,nc1(i))
419 vdy1(i)=vy1(i)-w(2,nc1(i))
420 vdz1(i)=vz1(i)-w(3,nc1(i))
421 vdx2(i)=vx2(i)-w(1,nc2(i))
422 vdy2(i)=vy2(i)-w(2,nc2(i))
423 vdz2(i)=vz2(i)-w(3,nc2(i))
424 vdx3(i)=vx3(i)-w(1,nc3(i))
425 vdy3(i)=vy3(i)-w(2,nc3(i))
426 vdz3(i)=vz3(i)-w(3,nc3(i))
427 vdx4(i)=vx4(i)-w(1,nc4(i))
428 vdy4(i)=vy4(i)-w(2,nc4(i))
429 vdz4(i)=vz4(i)-w(3,nc4(i))
430 vdx5(i)=vx5(i)-w(1,nc5(i))
431 vdy5(i)=vy5(i)-w(2,nc5(i))
432 vdz5(i)=vz5(i)-w(3,nc5(i))
433 vdx6(i)=vx6(i)-w(1,nc6(i))
434 vdy6(i)=vy6(i)-w(2,nc6(i))
435 vdz6(i)=vz6(i)-w(3,nc6(i))
436 vdx7(i)=vx7(i)-w(1,nc7(i))
437 vdy7(i)=vy7(i)-w(2,nc7(i))
438 vdz7(i)=vz7(i)-w(3,nc7(i))
439 vdx8(i)=vx8(i)-w(1,nc8(i))
440 vdy8(i)=vy8(i)-w(2,nc8(i))
441 vdz8(i)=vz8(i)-w(3,nc8(i))
442 ENDDO
443 ELSEIF(jeul/=0)THEN
444 IF(iparg(64,ng)==1)THEN
445 DO i=1,nel
446 vd2(i) =zero
447 vdx1(i)=zero
448 vdy1(i)=zero
449 vdz1(i)=zero
450 vdx2(i)=zero
451 vdy2(i)=zero
452 vdz2(i)=zero
453 vdx3(i)=zero
454 vdy3(i)=zero
455 vdz3(i)=zero
456 vdx4(i)=zero
457 vdy4(i)=zero
458 vdz4(i)=zero
459 vdx5(i)=zero
460 vdy5(i)=zero
461 vdz5(i)=zero
462 vdx6(i)=zero
463 vdy6(i)=zero
464 vdz6(i)=zero
465 vdx7(i)=zero
466 vdy7(i)=zero
467 vdz7(i)=zero
468 vdx8(i)=zero
469 vdy8(i)=zero
470 vdz8(i)=zero
471 ENDDO
472 RETURN
473 ENDIF
474 DO i=1,nel
475 vdx1(i)=vx1(i)
476 vdy1(i)=vy1(i)
477 vdz1(i)=vz1(i)
478 vdx2(i)=vx2(i)
479 vdy2(i)=vy2(i)
480 vdz2(i)=vz2(i)
481 vdx3(i)=vx3(i)
482 vdy3(i)=vy3(i)
483 vdz3(i)=vz3(i)
484 vdx4(i)=vx4(i)
485 vdy4(i)=vy4(i)
486 vdz4(i)=vz4(i)
487 vdx5(i)=vx5(i)
488 vdy5(i)=vy5(i)
489 vdz5(i)=vz5(i)
490 vdx6(i)=vx6(i)
491 vdy6(i)=vy6(i)
492 vdz6(i)=vz6(i)
493 vdx7(i)=vx7(i)
494 vdy7(i)=vy7(i)
495 vdz7(i)=vz7(i)
496 vdx8(i)=vx8(i)
497 vdy8(i)=vy8(i)
498 vdz8(i)=vz8(i)
499 ENDDO
500 ENDIF
501
502 DO i=1,nel
503 vdx(i)=one_over_8*(vdx1(i)+vdx2(i)+vdx3(i)+vdx4(i)+vdx5(i)+vdx6(i)+vdx7(i)+vdx8(i))
504 vdy(i)=one_over_8*(vdy1(i)+vdy2(i)+vdy3(i)+vdy4(i)+vdy5(i)+vdy6(i)+vdy7(i)+vdy8(i))
505 vdz(i)=one_over_8*(vdz1(i)+vdz2(i)+vdz3(i)+vdz4(i)+vdz5(i)+vdz6(i)+vdz7(i)+vdz8(i))
506 vd2(i)=(vdx(i)**2+vdy(i)**2+vdz(i)**2)
507 ENDDO
508C-----------
509 RETURN
510 END
scoor3
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)
Definition scoor3.F:54
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:43
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