OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
scumu3.F File Reference
#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "scr18_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine scumu3 (offg, e, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, stifn, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nvc, ar, fr_wave, fr_wav, mx1, my1, mz1, mx2, my2, mz2, mx3, my3, mz3, mx4, my4, mz4, mx5, my5, mz5, mx6, my6, mz6, mx7, my7, mz7, mx8, my8, mz8, them, fthe, condn, conde, nel, jthe, isrot, ipartsph, nodadt_therm)

Function/Subroutine Documentation

◆ scumu3()

subroutine scumu3 ( offg,
e,
integer, dimension(*) nc1,
integer, dimension(*) nc2,
integer, dimension(*) nc3,
integer, dimension(*) nc4,
integer, dimension(*) nc5,
integer, dimension(*) nc6,
integer, dimension(*) nc7,
integer, dimension(*) nc8,
stifn,
sti,
f11,
f21,
f31,
f12,
f22,
f32,
f13,
f23,
f33,
f14,
f24,
f34,
f15,
f25,
f35,
f16,
f26,
f36,
f17,
f27,
f37,
f18,
f28,
f38,
integer nvc,
ar,
fr_wave,
fr_wav,
mx1,
my1,
mz1,
mx2,
my2,
mz2,
mx3,
my3,
mz3,
mx4,
my4,
mz4,
mx5,
my5,
mz5,
mx6,
my6,
mz6,
mx7,
my7,
mz7,
mx8,
my8,
mz8,
them,
fthe,
condn,
conde,
integer, intent(in) nel,
integer, intent(in) jthe,
integer, intent(in) isrot,
integer, intent(in) ipartsph,
integer, intent(in) nodadt_therm )

Definition at line 35 of file scumu3.F.

55C-----------------------------------------------
56C I m p l i c i t T y p e s
57C-----------------------------------------------
58#include "implicit_f.inc"
59#include "comlock.inc"
60C-----------------------------------------------
61C G l o b a l P a r a m e t e r s
62C-----------------------------------------------
63#include "mvsiz_p.inc"
64C-----------------------------------------------
65C C o m m o n B l o c k s
66C-----------------------------------------------
67#include "com01_c.inc"
68C-----------------------------------------------
69C D u m m y A r g u m e n t s
70C-----------------------------------------------
71 INTEGER, INTENT(IN) :: NEL
72 INTEGER, INTENT(IN) :: JTHE
73 INTEGER, INTENT(IN) :: ISROT
74 INTEGER, INTENT(IN) :: IPARTSPH
75 INTEGER, INTENT(IN) :: NODADT_THERM
76 INTEGER NC1(*), NC2(*), NC3(*), NC4(*), NC5(*), NC6(*), NC7(*),
77 . NC8(*), NVC
78C REAL
79 my_real
80 . offg(*),e(3,*),stifn(*),sti(*),
81 . f11(*),f21(*),f31(*),f12(*),f22(*),f32(*),
82 . f13(*),f23(*),f33(*),f14(*),f24(*),f34(*),
83 . f15(*),f25(*),f35(*),f16(*),f26(*),f36(*),
84 . f17(*),f27(*),f37(*),f18(*),f28(*),f38(*),
85 . ar(3,*),fr_wave(*),fr_wav(*),
86 . mx1(*),my1(*),mz1(*),mx2(*),my2(*),mz2(*),
87 . mx3(*),my3(*),mz3(*),mx4(*),my4(*),mz4(*),
88 . mx5(*),my5(*),mz5(*),mx6(*),my6(*),mz6(*),
89 . mx7(*),my7(*),mz7(*),mx8(*),my8(*),mz8(*),
90 . them(mvsiz,8), fthe(*),condn(*),conde(*)
91C-----------------------------------------------
92C C o m m o n B l o c k s
93C-----------------------------------------------
94#include "scr18_c.inc"
95C-----------------------------------------------
96C L o c a l V a r i a b l e s
97C-----------------------------------------------
98 INTEGER I, J
99 my_real
100 . off_l
101 INTEGER NVC1,NVC2,NVC3,NVC4,NVC5,NVC6,NVC7,NVC8
102C-----------------------------------------------
103 off_l = zero
104 DO i=1,nel
105 off_l = min(off_l,offg(i))
106 ENDDO
107 IF(ipartsph==0)THEN
108 IF(off_l<zero)THEN
109 DO i=1,nel
110 IF(offg(i)<zero)THEN
111 f11(i)=zero
112 f21(i)=zero
113 f31(i)=zero
114 f12(i)=zero
115 f22(i)=zero
116 f32(i)=zero
117 f13(i)=zero
118 f23(i)=zero
119 f33(i)=zero
120 f14(i)=zero
121 f24(i)=zero
122 f34(i)=zero
123 f15(i)=zero
124 f25(i)=zero
125 f35(i)=zero
126 f16(i)=zero
127 f26(i)=zero
128 f36(i)=zero
129 f17(i)=zero
130 f27(i)=zero
131 f37(i)=zero
132 f18(i)=zero
133 f28(i)=zero
134 f38(i)=zero
135 sti(i)=zero
136 ENDIF
137 ENDDO
138 ENDIF
139 ELSE
140 IF(off_l<=zero)THEN
141 DO i=1,nel
142 IF(offg(i)<=zero)THEN
143C
144C case /DT/BRICK/DEL => OFFG==0 & F,STI/=0
145 f11(i)=zero
146 f21(i)=zero
147 f31(i)=zero
148 f12(i)=zero
149 f22(i)=zero
150 f32(i)=zero
151 f13(i)=zero
152 f23(i)=zero
153 f33(i)=zero
154 f14(i)=zero
155 f24(i)=zero
156 f34(i)=zero
157 f15(i)=zero
158 f25(i)=zero
159 f35(i)=zero
160 f16(i)=zero
161 f26(i)=zero
162 f36(i)=zero
163 f17(i)=zero
164 f27(i)=zero
165 f37(i)=zero
166 f18(i)=zero
167 f28(i)=zero
168 f38(i)=zero
169 sti(i)=zero
170 ENDIF
171 ENDDO
172 ENDIF
173 ENDIF
174 IF(jthe < 0 ) THEN
175 IF(off_l<=zero)THEN
176 DO j=1,8
177 DO i=1,nel
178 IF(offg(i)<=zero)THEN
179 them(i,j)=zero
180 ENDIF
181 ENDDO
182 ENDDO
183 ENDIF
184 IF(nodadt_therm == 1) THEN
185 IF(off_l<zero)THEN
186 DO i=1,nel
187 IF(offg(i)<zero)THEN
188 conde(i)=zero
189 ENDIF
190 ENDDO
191 ENDIF
192 ENDIF
193 ENDIF
194 nvc1 = nvc / 128
195 nvc2 =(nvc-nvc1*128) / 64
196 nvc3 =(nvc-nvc1*128-nvc2*64) / 32
197 nvc4 =(nvc-nvc1*128-nvc2*64-nvc3*32)/16
198 nvc5 =(nvc-nvc1*128-nvc2*64-nvc3*32-nvc4*16)/8
199 nvc6 =(nvc-nvc1*128-nvc2*64-nvc3*32-nvc4*16-nvc5*8)/4
200 nvc7 =(nvc-nvc1*128-nvc2*64-nvc3*32-nvc4*16-nvc5*8-nvc6*4)/2
201 nvc8 = nvc-nvc1*128-nvc2*64-nvc3*32-nvc4*16-nvc5*8-nvc6*4-nvc7*2
202C
203C because 2*Mnodal = 2*Melement/8
204 DO i=1,nel
205 sti(i)=fourth*sti(i)
206 END DO
207 IF(nodadt_therm == 1 ) THEN
208 DO i=1,nel
209 conde(i)=one_over_8*conde(i)
210 END DO
211 ENDIF
212C
213C
214 IF(jthe >= 0) THEN
215 IF(nvc1==0)THEN
216 DO i=1,nel
217 e(1,nc1(i))=e(1,nc1(i))+f11(i)
218 e(2,nc1(i))=e(2,nc1(i))+f21(i)
219 e(3,nc1(i))=e(3,nc1(i))+f31(i)
220 stifn(nc1(i))=stifn(nc1(i))+sti(i)
221 ENDDO
222 ELSE
223 DO i=1,nel
224 e(1,nc1(i))=e(1,nc1(i))+f11(i)
225 e(2,nc1(i))=e(2,nc1(i))+f21(i)
226 e(3,nc1(i))=e(3,nc1(i))+f31(i)
227 stifn(nc1(i))=stifn(nc1(i))+sti(i)
228 ENDDO
229 ENDIF
230 IF(nvc2==0)THEN
231 DO i=1,nel
232 e(1,nc2(i))=e(1,nc2(i))+f12(i)
233 e(2,nc2(i))=e(2,nc2(i))+f22(i)
234 e(3,nc2(i))=e(3,nc2(i))+f32(i)
235 stifn(nc2(i))=stifn(nc2(i))+sti(i)
236 ENDDO
237 ELSE
238 DO i=1,nel
239 e(1,nc2(i))=e(1,nc2(i))+f12(i)
240 e(2,nc2(i))=e(2,nc2(i))+f22(i)
241 e(3,nc2(i))=e(3,nc2(i))+f32(i)
242 stifn(nc2(i))=stifn(nc2(i))+sti(i)
243 ENDDO
244 ENDIF
245 IF(nvc3==0)THEN
246 DO i=1,nel
247 e(1,nc3(i))=e(1,nc3(i))+f13(i)
248 e(2,nc3(i))=e(2,nc3(i))+f23(i)
249 e(3,nc3(i))=e(3,nc3(i))+f33(i)
250 stifn(nc3(i))=stifn(nc3(i))+sti(i)
251 ENDDO
252 ELSE
253 DO i=1,nel
254 e(1,nc3(i))=e(1,nc3(i))+f13(i)
255 e(2,nc3(i))=e(2,nc3(i))+f23(i)
256 e(3,nc3(i))=e(3,nc3(i))+f33(i)
257 stifn(nc3(i))=stifn(nc3(i))+sti(i)
258 ENDDO
259 ENDIF
260 IF(nvc4==0)THEN
261 DO i=1,nel
262 e(1,nc4(i))=e(1,nc4(i))+f14(i)
263 e(2,nc4(i))=e(2,nc4(i))+f24(i)
264 e(3,nc4(i))=e(3,nc4(i))+f34(i)
265 stifn(nc4(i))=stifn(nc4(i))+sti(i)
266 ENDDO
267 ELSE
268 DO i=1,nel
269 e(1,nc4(i))=e(1,nc4(i))+f14(i)
270 e(2,nc4(i))=e(2,nc4(i))+f24(i)
271 e(3,nc4(i))=e(3,nc4(i))+f34(i)
272 stifn(nc4(i))=stifn(nc4(i))+sti(i)
273 ENDDO
274 ENDIF
275 IF(nvc5==0)THEN
276 DO i=1,nel
277 e(1,nc5(i))=e(1,nc5(i))+f15(i)
278 e(2,nc5(i))=e(2,nc5(i))+f25(i)
279 e(3,nc5(i))=e(3,nc5(i))+f35(i)
280 stifn(nc5(i))=stifn(nc5(i))+sti(i)
281 ENDDO
282 ELSE
283 DO i=1,nel
284 e(1,nc5(i))=e(1,nc5(i))+f15(i)
285 e(2,nc5(i))=e(2,nc5(i))+f25(i)
286 e(3,nc5(i))=e(3,nc5(i))+f35(i)
287 stifn(nc5(i))=stifn(nc5(i))+sti(i)
288 ENDDO
289 ENDIF
290 IF(nvc6==0)THEN
291 DO i=1,nel
292 e(1,nc6(i))=e(1,nc6(i))+f16(i)
293 e(2,nc6(i))=e(2,nc6(i))+f26(i)
294 e(3,nc6(i))=e(3,nc6(i))+f36(i)
295 stifn(nc6(i))=stifn(nc6(i))+sti(i)
296 ENDDO
297 ELSE
298 DO i=1,nel
299 e(1,nc6(i))=e(1,nc6(i))+f16(i)
300 e(2,nc6(i))=e(2,nc6(i))+f26(i)
301 e(3,nc6(i))=e(3,nc6(i))+f36(i)
302 stifn(nc6(i))=stifn(nc6(i))+sti(i)
303 ENDDO
304 ENDIF
305 IF(nvc7==0)THEN
306 DO i=1,nel
307 e(1,nc7(i))=e(1,nc7(i))+f17(i)
308 e(2,nc7(i))=e(2,nc7(i))+f27(i)
309 e(3,nc7(i))=e(3,nc7(i))+f37(i)
310 stifn(nc7(i))=stifn(nc7(i))+sti(i)
311 ENDDO
312 ELSE
313 DO i=1,nel
314 e(1,nc7(i))=e(1,nc7(i))+f17(i)
315 e(2,nc7(i))=e(2,nc7(i))+f27(i)
316 e(3,nc7(i))=e(3,nc7(i))+f37(i)
317 stifn(nc7(i))=stifn(nc7(i))+sti(i)
318 ENDDO
319 ENDIF
320 IF(nvc8==0)THEN
321 DO i=1,nel
322 e(1,nc8(i))=e(1,nc8(i))+f18(i)
323 e(2,nc8(i))=e(2,nc8(i))+f28(i)
324 e(3,nc8(i))=e(3,nc8(i))+f38(i)
325 stifn(nc8(i))=stifn(nc8(i))+sti(i)
326 ENDDO
327 ELSE
328 DO i=1,nel
329 e(1,nc8(i))=e(1,nc8(i))+f18(i)
330 e(2,nc8(i))=e(2,nc8(i))+f28(i)
331 e(3,nc8(i))=e(3,nc8(i))+f38(i)
332 stifn(nc8(i))=stifn(nc8(i))+sti(i)
333 ENDDO
334 ENDIF
335CCC
336 ELSE
337 IF(nodadt_therm == 1 ) THEN
338 IF(nvc1==0)THEN
339 DO i=1,nel
340 e(1,nc1(i))=e(1,nc1(i))+f11(i)
341 e(2,nc1(i))=e(2,nc1(i))+f21(i)
342 e(3,nc1(i))=e(3,nc1(i))+f31(i)
343 stifn(nc1(i))=stifn(nc1(i))+sti(i)
344 fthe(nc1(i)) = fthe(nc1(i)) + them(i,1)
345 condn(nc1(i))= condn(nc1(i))+ conde(i)
346 ENDDO
347 ELSE
348 DO i=1,nel
349 e(1,nc1(i))=e(1,nc1(i))+f11(i)
350 e(2,nc1(i))=e(2,nc1(i))+f21(i)
351 e(3,nc1(i))=e(3,nc1(i))+f31(i)
352 stifn(nc1(i))=stifn(nc1(i))+sti(i)
353 fthe(nc1(i)) = fthe(nc1(i)) + them(i,1)
354 condn(nc1(i))= condn(nc1(i))+ conde(i)
355 ENDDO
356 ENDIF
357 IF(nvc2==0)THEN
358 DO i=1,nel
359 e(1,nc2(i))=e(1,nc2(i))+f12(i)
360 e(2,nc2(i))=e(2,nc2(i))+f22(i)
361 e(3,nc2(i))=e(3,nc2(i))+f32(i)
362 stifn(nc2(i))=stifn(nc2(i))+sti(i)
363 fthe(nc2(i)) = fthe(nc2(i)) + them(i,2)
364 condn(nc2(i))= condn(nc2(i))+ conde(i)
365 ENDDO
366 ELSE
367 DO i=1,nel
368 e(1,nc2(i))=e(1,nc2(i))+f12(i)
369 e(2,nc2(i))=e(2,nc2(i))+f22(i)
370 e(3,nc2(i))=e(3,nc2(i))+f32(i)
371 stifn(nc2(i))=stifn(nc2(i))+sti(i)
372 fthe(nc2(i)) = fthe(nc2(i)) + them(i,2)
373 condn(nc2(i))= condn(nc2(i))+ conde(i)
374 ENDDO
375 ENDIF
376 IF(nvc3==0)THEN
377 DO i=1,nel
378 e(1,nc3(i))=e(1,nc3(i))+f13(i)
379 e(2,nc3(i))=e(2,nc3(i))+f23(i)
380 e(3,nc3(i))=e(3,nc3(i))+f33(i)
381 stifn(nc3(i))=stifn(nc3(i))+sti(i)
382 fthe(nc3(i)) = fthe(nc3(i)) + them(i,3)
383 condn(nc3(i))= condn(nc3(i))+ conde(i)
384 ENDDO
385 ELSE
386 DO i=1,nel
387 e(1,nc3(i))=e(1,nc3(i))+f13(i)
388 e(2,nc3(i))=e(2,nc3(i))+f23(i)
389 e(3,nc3(i))=e(3,nc3(i))+f33(i)
390 stifn(nc3(i))=stifn(nc3(i))+sti(i)
391 fthe(nc3(i)) = fthe(nc3(i)) + them(i,3)
392 condn(nc3(i))= condn(nc3(i))+ conde(i)
393 ENDDO
394 ENDIF
395 IF(nvc4==0)THEN
396 DO i=1,nel
397 e(1,nc4(i))=e(1,nc4(i))+f14(i)
398 e(2,nc4(i))=e(2,nc4(i))+f24(i)
399 e(3,nc4(i))=e(3,nc4(i))+f34(i)
400 stifn(nc4(i))=stifn(nc4(i))+sti(i)
401 fthe(nc4(i)) = fthe(nc4(i)) + them(i,4)
402 condn(nc4(i))= condn(nc4(i))+ conde(i)
403 ENDDO
404 ELSE
405 DO i=1,nel
406 e(1,nc4(i))=e(1,nc4(i))+f14(i)
407 e(2,nc4(i))=e(2,nc4(i))+f24(i)
408 e(3,nc4(i))=e(3,nc4(i))+f34(i)
409 stifn(nc4(i))=stifn(nc4(i))+sti(i)
410 fthe(nc4(i)) = fthe(nc4(i)) + them(i,4)
411 condn(nc4(i))= condn(nc4(i))+ conde(i)
412 ENDDO
413 ENDIF
414 IF(nvc5==0)THEN
415 DO i=1,nel
416 e(1,nc5(i))=e(1,nc5(i))+f15(i)
417 e(2,nc5(i))=e(2,nc5(i))+f25(i)
418 e(3,nc5(i))=e(3,nc5(i))+f35(i)
419 stifn(nc5(i))=stifn(nc5(i))+sti(i)
420 fthe(nc5(i)) = fthe(nc5(i)) + them(i,5)
421 condn(nc5(i))= condn(nc5(i))+ conde(i)
422 ENDDO
423 ELSE
424 DO i=1,nel
425 e(1,nc5(i))=e(1,nc5(i))+f15(i)
426 e(2,nc5(i))=e(2,nc5(i))+f25(i)
427 e(3,nc5(i))=e(3,nc5(i))+f35(i)
428 stifn(nc5(i))=stifn(nc5(i))+sti(i)
429 fthe(nc5(i)) = fthe(nc5(i)) + them(i,5)
430 condn(nc5(i))= condn(nc5(i))+ conde(i)
431 ENDDO
432 ENDIF
433 IF(nvc6==0)THEN
434 DO i=1,nel
435 e(1,nc6(i))=e(1,nc6(i))+f16(i)
436 e(2,nc6(i))=e(2,nc6(i))+f26(i)
437 e(3,nc6(i))=e(3,nc6(i))+f36(i)
438 stifn(nc6(i))=stifn(nc6(i))+sti(i)
439 fthe(nc6(i)) = fthe(nc6(i)) + them(i,6)
440 condn(nc6(i))= condn(nc6(i))+ conde(i)
441 ENDDO
442 ELSE
443 DO i=1,nel
444 e(1,nc6(i))=e(1,nc6(i))+f16(i)
445 e(2,nc6(i))=e(2,nc6(i))+f26(i)
446 e(3,nc6(i))=e(3,nc6(i))+f36(i)
447 stifn(nc6(i))=stifn(nc6(i))+sti(i)
448 fthe(nc6(i)) = fthe(nc6(i)) + them(i,6)
449 condn(nc6(i))= condn(nc6(i))+ conde(i)
450 ENDDO
451 ENDIF
452 IF(nvc7==0)THEN
453 DO i=1,nel
454 e(1,nc7(i))=e(1,nc7(i))+f17(i)
455 e(2,nc7(i))=e(2,nc7(i))+f27(i)
456 e(3,nc7(i))=e(3,nc7(i))+f37(i)
457 stifn(nc7(i))=stifn(nc7(i))+sti(i)
458 fthe(nc7(i)) = fthe(nc7(i)) + them(i,7)
459 condn(nc7(i))= condn(nc7(i))+ conde(i)
460 ENDDO
461 ELSE
462 DO i=1,nel
463 e(1,nc7(i))=e(1,nc7(i))+f17(i)
464 e(2,nc7(i))=e(2,nc7(i))+f27(i)
465 e(3,nc7(i))=e(3,nc7(i))+f37(i)
466 stifn(nc7(i))=stifn(nc7(i))+sti(i)
467 fthe(nc7(i)) = fthe(nc7(i)) + them(i,7)
468 condn(nc7(i))= condn(nc7(i))+ conde(i)
469 ENDDO
470 ENDIF
471 IF(nvc8==0)THEN
472 DO i=1,nel
473 e(1,nc8(i))=e(1,nc8(i))+f18(i)
474 e(2,nc8(i))=e(2,nc8(i))+f28(i)
475 e(3,nc8(i))=e(3,nc8(i))+f38(i)
476 stifn(nc8(i))=stifn(nc8(i))+sti(i)
477 fthe(nc8(i)) = fthe(nc8(i)) + them(i,8)
478 condn(nc8(i))= condn(nc8(i))+ conde(i)
479 ENDDO
480 ELSE
481 DO i=1,nel
482 e(1,nc8(i))=e(1,nc8(i))+f18(i)
483 e(2,nc8(i))=e(2,nc8(i))+f28(i)
484 e(3,nc8(i))=e(3,nc8(i))+f38(i)
485 stifn(nc8(i))=stifn(nc8(i))+sti(i)
486 fthe(nc8(i)) = fthe(nc8(i)) + them(i,8)
487 condn(nc8(i))= condn(nc8(i))+ conde(i)
488 ENDDO
489 ENDIF
490C
491 ELSE ! NODADT_THERM
492C
493 IF(nvc1==0)THEN
494 DO i=1,nel
495 e(1,nc1(i))=e(1,nc1(i))+f11(i)
496 e(2,nc1(i))=e(2,nc1(i))+f21(i)
497 e(3,nc1(i))=e(3,nc1(i))+f31(i)
498 stifn(nc1(i))=stifn(nc1(i))+sti(i)
499 fthe(nc1(i)) = fthe(nc1(i)) + them(i,1)
500 ENDDO
501 ELSE
502 DO i=1,nel
503 e(1,nc1(i))=e(1,nc1(i))+f11(i)
504 e(2,nc1(i))=e(2,nc1(i))+f21(i)
505 e(3,nc1(i))=e(3,nc1(i))+f31(i)
506 stifn(nc1(i))=stifn(nc1(i))+sti(i)
507 fthe(nc1(i)) = fthe(nc1(i)) + them(i,1)
508 ENDDO
509 ENDIF
510 IF(nvc2==0)THEN
511 DO i=1,nel
512 e(1,nc2(i))=e(1,nc2(i))+f12(i)
513 e(2,nc2(i))=e(2,nc2(i))+f22(i)
514 e(3,nc2(i))=e(3,nc2(i))+f32(i)
515 stifn(nc2(i))=stifn(nc2(i))+sti(i)
516 fthe(nc2(i)) = fthe(nc2(i)) + them(i,2)
517 ENDDO
518 ELSE
519 DO i=1,nel
520 e(1,nc2(i))=e(1,nc2(i))+f12(i)
521 e(2,nc2(i))=e(2,nc2(i))+f22(i)
522 e(3,nc2(i))=e(3,nc2(i))+f32(i)
523 stifn(nc2(i))=stifn(nc2(i))+sti(i)
524 fthe(nc2(i)) = fthe(nc2(i)) + them(i,2)
525 ENDDO
526 ENDIF
527 IF(nvc3==0)THEN
528 DO i=1,nel
529 e(1,nc3(i))=e(1,nc3(i))+f13(i)
530 e(2,nc3(i))=e(2,nc3(i))+f23(i)
531 e(3,nc3(i))=e(3,nc3(i))+f33(i)
532 stifn(nc3(i))=stifn(nc3(i))+sti(i)
533 fthe(nc3(i)) = fthe(nc3(i)) + them(i,3)
534 ENDDO
535 ELSE
536 DO i=1,nel
537 e(1,nc3(i))=e(1,nc3(i))+f13(i)
538 e(2,nc3(i))=e(2,nc3(i))+f23(i)
539 e(3,nc3(i))=e(3,nc3(i))+f33(i)
540 stifn(nc3(i))=stifn(nc3(i))+sti(i)
541 fthe(nc3(i)) = fthe(nc3(i)) + them(i,3)
542 ENDDO
543 ENDIF
544 IF(nvc4==0)THEN
545 DO i=1,nel
546 e(1,nc4(i))=e(1,nc4(i))+f14(i)
547 e(2,nc4(i))=e(2,nc4(i))+f24(i)
548 e(3,nc4(i))=e(3,nc4(i))+f34(i)
549 stifn(nc4(i))=stifn(nc4(i))+sti(i)
550 fthe(nc4(i)) = fthe(nc4(i)) + them(i,4)
551 ENDDO
552 ELSE
553 DO i=1,nel
554 e(1,nc4(i))=e(1,nc4(i))+f14(i)
555 e(2,nc4(i))=e(2,nc4(i))+f24(i)
556 e(3,nc4(i))=e(3,nc4(i))+f34(i)
557 stifn(nc4(i))=stifn(nc4(i))+sti(i)
558 fthe(nc4(i)) = fthe(nc4(i)) + them(i,4)
559 ENDDO
560 ENDIF
561 IF(nvc5==0)THEN
562 DO i=1,nel
563 e(1,nc5(i))=e(1,nc5(i))+f15(i)
564 e(2,nc5(i))=e(2,nc5(i))+f25(i)
565 e(3,nc5(i))=e(3,nc5(i))+f35(i)
566 stifn(nc5(i))=stifn(nc5(i))+sti(i)
567 fthe(nc5(i)) = fthe(nc5(i)) + them(i,5)
568 ENDDO
569 ELSE
570 DO i=1,nel
571 e(1,nc5(i))=e(1,nc5(i))+f15(i)
572 e(2,nc5(i))=e(2,nc5(i))+f25(i)
573 e(3,nc5(i))=e(3,nc5(i))+f35(i)
574 stifn(nc5(i))=stifn(nc5(i))+sti(i)
575 fthe(nc5(i)) = fthe(nc5(i)) + them(i,5)
576 ENDDO
577 ENDIF
578 IF(nvc6==0)THEN
579 DO i=1,nel
580 e(1,nc6(i))=e(1,nc6(i))+f16(i)
581 e(2,nc6(i))=e(2,nc6(i))+f26(i)
582 e(3,nc6(i))=e(3,nc6(i))+f36(i)
583 stifn(nc6(i))=stifn(nc6(i))+sti(i)
584 fthe(nc6(i)) = fthe(nc6(i)) + them(i,6)
585 ENDDO
586 ELSE
587 DO i=1,nel
588 e(1,nc6(i))=e(1,nc6(i))+f16(i)
589 e(2,nc6(i))=e(2,nc6(i))+f26(i)
590 e(3,nc6(i))=e(3,nc6(i))+f36(i)
591 stifn(nc6(i))=stifn(nc6(i))+sti(i)
592 fthe(nc6(i)) = fthe(nc6(i)) + them(i,6)
593 ENDDO
594 ENDIF
595 IF(nvc7==0)THEN
596 DO i=1,nel
597 e(1,nc7(i))=e(1,nc7(i))+f17(i)
598 e(2,nc7(i))=e(2,nc7(i))+f27(i)
599 e(3,nc7(i))=e(3,nc7(i))+f37(i)
600 stifn(nc7(i))=stifn(nc7(i))+sti(i)
601 fthe(nc7(i)) = fthe(nc7(i)) + them(i,7)
602 ENDDO
603 ELSE
604 DO i=1,nel
605 e(1,nc7(i))=e(1,nc7(i))+f17(i)
606 e(2,nc7(i))=e(2,nc7(i))+f27(i)
607 e(3,nc7(i))=e(3,nc7(i))+f37(i)
608 stifn(nc7(i))=stifn(nc7(i))+sti(i)
609 fthe(nc7(i)) = fthe(nc7(i)) + them(i,7)
610 ENDDO
611 ENDIF
612 IF(nvc8==0)THEN
613 DO i=1,nel
614 e(1,nc8(i))=e(1,nc8(i))+f18(i)
615 e(2,nc8(i))=e(2,nc8(i))+f28(i)
616 e(3,nc8(i))=e(3,nc8(i))+f38(i)
617 stifn(nc8(i))=stifn(nc8(i))+sti(i)
618 fthe(nc8(i)) = fthe(nc8(i)) + them(i,8)
619 ENDDO
620 ELSE
621 DO i=1,nel
622 e(1,nc8(i))=e(1,nc8(i))+f18(i)
623 e(2,nc8(i))=e(2,nc8(i))+f28(i)
624 e(3,nc8(i))=e(3,nc8(i))+f38(i)
625 stifn(nc8(i))=stifn(nc8(i))+sti(i)
626 fthe(nc8(i)) = fthe(nc8(i)) + them(i,8)
627 ENDDO
628 ENDIF
629C
630 ENDIF
631 ENDIF
632C
633 IF(isrot/=0)THEN
634 IF(off_l<zero)THEN
635 DO i=1,nel
636 IF(offg(i)<zero)THEN
637 mx1(i)=0.
638 my1(i)=0.
639 mz1(i)=0.
640 mx2(i)=0.
641 my2(i)=0.
642 mz2(i)=0.
643 mx3(i)=0.
644 my3(i)=0.
645 mz3(i)=0.
646 mx4(i)=0.
647 my4(i)=0.
648 mz4(i)=0.
649 mx5(i)=0.
650 my5(i)=0.
651 mz5(i)=0.
652 mx6(i)=0.
653 my6(i)=0.
654 mz6(i)=0.
655 mx7(i)=0.
656 my7(i)=0.
657 mz7(i)=0.
658 mx8(i)=0.
659 my8(i)=0.
660 mz8(i)=0.
661 ENDIF
662 ENDDO
663 ENDIF
664#include "lockon.inc"
665 DO i=1,nel
666 ar(1,nc1(i))=ar(1,nc1(i))+mx1(i)
667 ar(2,nc1(i))=ar(2,nc1(i))+my1(i)
668 ar(3,nc1(i))=ar(3,nc1(i))+mz1(i)
669C
670 ar(1,nc2(i))=ar(1,nc2(i))+mx2(i)
671 ar(2,nc2(i))=ar(2,nc2(i))+my2(i)
672 ar(3,nc2(i))=ar(3,nc2(i))+mz2(i)
673C
674 ar(1,nc3(i))=ar(1,nc3(i))+mx3(i)
675 ar(2,nc3(i))=ar(2,nc3(i))+my3(i)
676 ar(3,nc3(i))=ar(3,nc3(i))+mz3(i)
677C
678 ar(1,nc4(i))=ar(1,nc4(i))+mx4(i)
679 ar(2,nc4(i))=ar(2,nc4(i))+my4(i)
680 ar(3,nc4(i))=ar(3,nc4(i))+mz4(i)
681C
682 ar(1,nc5(i))=ar(1,nc5(i))+mx5(i)
683 ar(2,nc5(i))=ar(2,nc5(i))+my5(i)
684 ar(3,nc5(i))=ar(3,nc5(i))+mz5(i)
685C
686 ar(1,nc6(i))=ar(1,nc6(i))+mx6(i)
687 ar(2,nc6(i))=ar(2,nc6(i))+my6(i)
688 ar(3,nc6(i))=ar(3,nc6(i))+mz6(i)
689C
690 ar(1,nc7(i))=ar(1,nc7(i))+mx7(i)
691 ar(2,nc7(i))=ar(2,nc7(i))+my7(i)
692 ar(3,nc7(i))=ar(3,nc7(i))+mz7(i)
693C
694 ar(1,nc8(i))=ar(1,nc8(i))+mx8(i)
695 ar(2,nc8(i))=ar(2,nc8(i))+my8(i)
696 ar(3,nc8(i))=ar(3,nc8(i))+mz8(i)
697 ENDDO
698#include "lockoff.inc"
699C
700C--------------------------------------------
701C Front wave
702C--------------------------------------------
703 IF (ifrwv/=0)THEN
704#include "lockon.inc"
705 DO i=1,nel
706 IF(fr_wave(nc1(i))==0.0)fr_wave(nc1(i))=-fr_wav(i)
707 IF(fr_wave(nc2(i))==0.0)fr_wave(nc2(i))=-fr_wav(i)
708 IF(fr_wave(nc3(i))==0.0)fr_wave(nc3(i))=-fr_wav(i)
709 IF(fr_wave(nc4(i))==0.0)fr_wave(nc4(i))=-fr_wav(i)
710 IF(fr_wave(nc5(i))==0.0)fr_wave(nc5(i))=-fr_wav(i)
711 IF(fr_wave(nc6(i))==0.0)fr_wave(nc6(i))=-fr_wav(i)
712 IF(fr_wave(nc7(i))==0.0)fr_wave(nc7(i))=-fr_wav(i)
713 IF(fr_wave(nc8(i))==0.0)fr_wave(nc8(i))=-fr_wav(i)
714c FR_WAVE(NC1(I)) = MAX(FR_WAVE(NC1(I)),FR_WAV(I))
715c FR_WAVE(NC2(I)) = MAX(FR_WAVE(NC2(I)),FR_WAV(I))
716c FR_WAVE(NC3(I)) = MAX(FR_WAVE(NC3(I)),FR_WAV(I))
717c FR_WAVE(NC4(I)) = MAX(FR_WAVE(NC4(I)),FR_WAV(I))
718c FR_WAVE(NC5(I)) = MAX(FR_WAVE(NC5(I)),FR_WAV(I))
719c FR_WAVE(NC6(I)) = MAX(FR_WAVE(NC6(I)),FR_WAV(I))
720c FR_WAVE(NC7(I)) = MAX(FR_WAVE(NC7(I)),FR_WAV(I))
721c FR_WAVE(NC8(I)) = MAX(FR_WAVE(NC8(I)),FR_WAV(I))
722 ENDDO
723#include "lockoff.inc"
724 ENDIF
725 ENDIF
726 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20