36
37
38
39#include "implicit_f.inc"
40
41
42
43
44
45
46 INTEGER G_PLA,G_EPSD,NEL
47
49 . fx1(*), fx2(*), fx3(*), fx4(*), fx5(*), fx6(*),
50 . fy1(*), fy2(*), fy3(*), fy4(*), fy5(*), fy6(*),
51 . fz1(*), fz2(*), fz3(*), fz4(*), fz5(*), fz6(*),
52 . eintm(*),rhom(*),qm(*),epsdm(*),eplasm(*),sigm(nel,6)
53
54
55
56 INTEGER I
57
58 DO i=1,nel
59 sigm(i,1) = zero
60 sigm(i,2) = zero
61 sigm(i,3) = zero
62 sigm(i,4) = zero
63 sigm(i,5) = zero
64 sigm(i,6) = zero
65 rhom(i) = zero
66 eintm(i) = zero
67 qm(i) = zero
68 ENDDO
69 IF (g_epsd > 0) THEN
70 DO i=1,nel
71 epsdm(i) = zero
72 ENDDO
73 ENDIF
74 IF (g_pla > 0) THEN
75 DO i=1,nel
76 eplasm(i) = zero
77 ENDDO
78 ENDIF
79 DO i=1,nel
80 fx1(i)=zero
81 fy1(i)=zero
82 fz1(i)=zero
83 fx2(i)=zero
84 fy2(i)=zero
85 fz2(i)=zero
86 fx3(i)=zero
87 fy3(i)=zero
88 fz3(i)=zero
89 fx4(i)=zero
90 fy4(i)=zero
91 fz4(i)=zero
92 fx5(i)=zero
93 fy5(i)=zero
94 fz5(i)=zero
95 fx6(i)=zero
96 fy6(i)=zero
97 fz6(i)=zero
98 ENDDO
99
100 RETURN