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