38
40
41
42
43#include "implicit_f.inc"
44
45
46
47#include "mvsiz_p.inc"
48
49
50
51#include "vect01_c.inc"
52#include "com04_c.inc"
53
54
55
56 INTEGER NEL
57 INTEGER ,INTENT(IN) :: NINTEMP
58 INTEGER NC(MVSIZ,20),MXT(*),NGL(*),NGEO(*),IXS(NIXS,*),IXS20(12,*)
59
61 . x(3,*), v(3,*), xx(mvsiz,20), yy(mvsiz,20), zz(mvsiz,20),
62 . vx(mvsiz,20), vy(mvsiz,20), vz(mvsiz,20),
63
64 . sigg(nel,6),eintg
65
66
67
68 INTEGER I, IPERM1(20),IPERM2(20),N,N1,N2,NN,IUN
69
70 DATA iperm1/0,0,0,0,0,0,0,0,1,2,3,4,1,2,3,4,5,6,7,8/
71 DATA iperm2/0,0,0,0,0,0,0,0,2,3,4,1,5,6,7,8,6,7,8,5/
72
73
74
77
78 iun=1
79
80 DO i=lft,llt
81 ngeo(i)=ixs(10,i)
82 ngl(i) =ixs(11,i)
83 mxt
84 nc(i,1)=ixs(2,i)
85 nc(i,2)=ixs(3,i)
86 nc(i,3)=ixs(4,i)
87 nc(i,4)=ixs(5,i)
88 nc(i,5)=ixs(6,i)
89 nc(i,6)=ixs(7,i)
90 nc(i,7)=ixs(8,i)
91 nc(i,8)=ixs(9,i)
92 nc(i,9) =ixs20(1,i)
93 nc(i,10) =ixs20(2,i)
94 nc(i,11) =ixs20(3,i)
95 nc(i,12) =ixs20(4,i)
96 nc(i,13) =ixs20(5,i)
97 nc(i,14) =ixs20(6,i)
98 nc(i,15) =ixs20(7,i)
99 nc(i,16) =ixs20(8,i)
100 nc(i,17) =ixs20(9,i)
101 nc(i,18) =ixs20(10,i)
102 nc(i,19) =ixs20(11,i)
103 nc(i,20) =ixs20(12,i)
104 dtelem(i)=ep30
105 sti(i)=zero
106 eintg(i)=zero
107 rhog(i)=zero
108 qg(i)=zero
109 sigg(i,1)=zero
110 sigg(i,2)=zero
111 sigg(i,3)=zero
112 sigg(i,4)=zero
113 sigg(i,5)=zero
114 sigg(i,6)=zero
116
117 nc(i,1) =ixs(6,i)
118 nc(i,2) =ixs(7,i)
119 nc(i,3) =ixs(8,i)
120 nc(i,4) =ixs(9,i)
121 nc(i,5) =ixs(2,i)
122 nc(i,6) =ixs(3,i)
123 nc(i,7) =ixs(4,i)
124 nc(i,8) =ixs(5,i)
125 nc(i,9) =ixs20(9,i)
126 nc(i,10)=ixs20(10,i)
127 nc(i,11)=ixs20(11,i)
128 nc(i,12)=ixs20(12,i)
129 nc(i,17)=ixs20(1,i)
130 nc(i,18)=ixs20(2,i)
131 nc(i,19)=ixs20(3,i)
132 nc(i,20)=ixs20(4,i)
133 ixs(2,i) = nc(i,1)
134 ixs(3,i) = nc(i,2)
135 ixs(4,i) = nc(i,3)
136 ixs(5,i) = nc(i,4)
137 ixs(6,i) = nc(i,5
138 ixs(7,i) = nc(i,6)
139 ixs(8,i) = nc(i,7)
140 ixs(9,i) = nc(i,8)
141 ixs20(1,i) = nc(i,9)
142 ixs20(2,i) = nc(i,10)
143 ixs20(3,i) = nc(i,11)
144 ixs20(4,i) = nc(i,12)
145 ixs20(5,i) = nc(i,13)
146 ixs20(6,i) = nc(i,14)
147 ixs20(7,i) = nc(i,15)
148 ixs20(8,i) = nc(i,16)
149 ixs20(9,i) = nc(i,17)
150 ixs20(10,i) = nc(i,18)
151 ixs20(11,i) = nc(i,19)
152 ixs20(12,i) = nc(i,20)
153 ENDIF
154 ENDDO
155
156
157
158
159 DO n=1,20
160 DO i=lft,llt
161 nn =
max(iun,nc(i,n))
162 xx(i,n)=x(1,nn)
163 yy(i,n)=x(2,nn)
164 zz(i,n)=x(3,nn)
165 vx(i,n)=v(1,nn)
166 vy(i,n)=v(2,nn)
167 vz(i,n)=v(3,nn)
168 ENDDO
169 ENDDO
170
171 DO i=lft,llt
172 mass(i)=zero
173 ENDDO
174
175 DO n=9,20
176 n1=iperm1(n)
177 n2=iperm2(n)
178 DO i=lft,llt
179 IF(nc(i,n)==0)THEN
180 xx(i,n) = half*(xx(i,n1)+xx(i,n2))
181 yy(i,n) = half*(yy(i,n1)+yy(i,n2))
182 zz(i,n) = half*(zz(i,n1)+zz(i,n2))
183 vx(i,n) = half*(vx(i,n1)+vx(i,n2))
184 vy(i,n) = half*(vy(i,n1)+vy(i,n2))
185 vz(i,n) = half*(vz(i,n1)+vz(i,n2))
186 ENDIF
187 ENDDO
188 ENDDO
189
190 IF (jthe < 0 .or. nintemp > 0) THEN
191 IF(nintemp > 0 ) THEN
192 DO n =1,20
193 DO i=lft,llt
194 nn =
max(iun,nc(i,n))
195 IF(temp(nn)== zero) temp(nn) = temp0(i)
196 ENDDO
197 ENDDO
198 ELSE
199 DO n =1,20
200 DO i=lft,llt
201 nn =
max(iun,nc(i,n))
202 temp(nn) = temp0(i)
203 ENDDO
204 ENDDO
205 ENDIF
206 ENDIF
207
208 RETURN
function checkvolume_8n(x, ixs)
1.15.0