29
30
31
32
33
34
35#include "implicit_f.inc"
36
37
38
39#include "com01_c.inc"
40#include "com08_c.inc"
41#include "com04_c.inc"
42#include "inter22.inc"
43
44
45
46 INTEGER NODFT, NODLT, NDDIM,ITAB(*),NALE(*)
47 my_real fzero(3,numnod) ,a(3,numnod) ,v(3,numnod) ,ar(3,numnod) ,vr(3,numnod)
48
49
50
51 INTEGER I,J,K,N
52
53
54
55
56
57 DO n=1,numnod
58 v(1,n)=v(1,n)+dt12*a(1,n)
59 v(2,n)=v(2,n)+dt12*a(2,n)
60 v(3,n)=v(3,n)+dt12*a(3,n)
61 a(1,n)=zero
62 a(2,n)=zero
63 a(3,n)=zero
64 ENDDO
65
66
67 IF(isigi == 2 .OR. isigi == 4) THEN
68
69 DO i=1,numnod
70 a(1,i)=fzero(1,i)
71 a(2,i)=fzero(2,i)
72 a(3,i)=fzero(3,i)
73 ENDDO
74
75 ENDIF
76
77 IF(iroddl /= 0) THEN
78
79 DO n=1,numnod
80 vr(1,n)=vr(1,n)+dt12*ar(1,n)
81 vr(2,n)=vr(2,n)+dt12*ar(2,n)
82 vr(3,n)=vr(3,n)+dt12*ar(3,n)
83 ar(1,n)=zero
84 ar(2,n)=zero
85 ar(3,n)=zero
86 ENDDO
87
88 ENDIF
89
90
91
92
93 IF(int22 > 0) THEN
94
95 DO n=1,numnod
96 IF(nale(n) /= 0)THEN
97 v(1,n)=zero
98 v(2,n)=zero
99 v(3,n)=zero
100 vr(1,n)=zero
101 vr(2,n)=zero
102 vr(3,n)=zero
103 ENDIF
104 ENDDO
105
106 ENDIF
107
108 RETURN