34
35
36
37#include "implicit_f.inc"
38
39
40
41#include "mvsiz_p.inc"
42
43
44
45#include "scr18_c.inc"
46
47
48
49 INTEGER NCTRL,G_PLA,G_EPSD,NEL
50
52 . volm(*),sigm(nel,6),eintm(*),rhom(*),
53 . qm(*), fx(mvsiz,*),fy(mvsiz,*),
54 . fz(mvsiz,*),btdba(3*nctrl,*),stig(mvsiz,*),
55 . mass(nctrl,mvsiz),mmunk(mvsiz),aface(6,*),
56 . detmin(mvsiz),eplasm(*),epsdg(*)
57
58
59
60 INTEGER I, J, ITEST
61
62 DO i=1,nel
63 sigm(i,1) = zero
64 sigm(i,2) = zero
65 sigm(i,3) = zero
66 sigm(i,4) = zero
67 sigm(i,5) = zero
68 sigm(i,6) = zero
69 rhom(i) = zero
70 eintm(i) = zero
71 qm(i) = zero
72 volm(i) = zero
73 DO j=1,nctrl
74 fx(i,j)=zero
75 fy(i,j)=zero
76 fz(i,j)=zero
77 stig(i,j)= zero
78 ENDDO
79 ENDDO
80
81 IF(idtmin(101)==1)THEN
82 DO i=1,nel
83 DO j=1,nctrl
84 mmunk(i)=ep10
85 mass(j,i)=zero
86 ENDDO
87 ENDDO
88 ENDIF
89
90 IF(idtmin(101)==2)THEN
91 DO i=1,nel
92 aface(1,i)=zero
93 aface(2,i)=zero
94 aface(3,i)=zero
95 aface(4,i)=zero
96 aface(5,i)=zero
97 aface(6,i)=zero
98 detmin(i)=ep10
99 ENDDO
100 ENDIF
101
102 IF (g_pla > 0) THEN
103 DO i=1,nel
104 eplasm(i)=zero
105 ENDDO
106 ENDIF
107 IF (g_epsd > 0) THEN
108 DO i=1,nel
109 epsdg(i)=zero
110 ENDDO
111 ENDIF
112
113 RETURN