34
35
36
37#include "implicit_f.inc"
38
39
40
41#include "mvsiz_p.inc"
42
43
44
45 INTEGER, INTENT(IN) :: NEL
46
48 . vx(mvsiz,10),vy(mvsiz,10),vz(mvsiz,10),
49 . px(mvsiz,10),py(mvsiz,10),pz(mvsiz,10),
50 . dxx(*), dxy(*), dxz(*),
51 . dyx(*), dyy(*), dyz(*),
52 . dzx(*), dzy(*), dzz(*)
53
54
55
56
57
58
59 INTEGER I,N
60
61
62
63
64
65 DO i=1,nel
66 dxx(i)=px(i,1)*vx(i,1)+px(i,2)*vx(i,2)+px(i,3)*vx(i,3)+px(i,4)*vx(i,4)+px(i,5)*vx(i,5)
67 . +px(i,6)*vx(i,6)+px(i,7)*vx(i,7)+px(i,8)*vx(i,8)+px(i,9)*vx(i,9)+px(i,10)*vx(i,10)
68 dyy(i)=py(i,1)*vy(i,1)+py(i,2)*vy(i,2)+py(i,3)*vy(i,3)+py(i,4)*vy(i,4)+py(i,5)*vy(i,5)
69 . +py(i,6)*vy(i,6)+py(i,7)*vy(i,7)+py(i,8)*vy(i,8)+py(i,9)*vy(i,9)+py(i,10)*vy(i,10)
70 dzz(i)=pz(i,1)*vz(i,1)+pz(i,2)*vz(i,2)+pz(i,3)*vz(i,3)+pz(i,4)*vz(i,4)+pz(i,5)*vz(i,5)
71 . +pz(i,6)*vz(i,6)+pz(i,7)*vz(i,7)+pz(i,8)*vz(i,8)+pz(i,9)*vz(i,9)+pz(i,10)*vz(i,10)
72 dxy(i)=py(i,1)*vx(i,1)+py(i,2)*vx(i,2)+py(i,3)*vx(i,3)+py(i,4)*vx(i,4)+py(i,5)*vx(i,5)
73 . +py(i,6)*vx(i,6)+py(i,7)*vx(i,7)+py(i,8)*vx(i,8)+py(i,9)*vx(i,9)+py(i,10)*vx(i
74 dxz(i)=pz(i,1)*vx(i,1)+pz(i,2)*vx(i,2)+pz(i,3)*vx(i,3)+pz(i,4)*vx(i,4)+pz(i,5)*vx(i,5)
75 . +pz(i,6)*vx(i,6)+pz(i,7)*vx(i,7)+pz(i,8)*vx(i,8)+pz(i,9)*vx(i,9)+pz(i,10)*vx(i,10)
76 dyx(i)=px(i,1)*vy(i,1)+px(i,2)*vy(i,2)+px(i,3)*vy(i,3)+px(i,4)*vy(i,4)+px(i,5)*vy(i,5)
77 . +px(i,6)*vy(i,6)+px(i,7)*vy(i,7)+px(i,8)*vy(i,8)+px(i,9)*vy(i,9)+px
78 dyz(i)=pz(i,1)*vy(i,1)+pz(i,2)*vy(i,2)+pz(i,3)*vy(i,3)+pz(i,4)*vy(i,4)+pz(i,5)*vy(i,5)
79 . +pz(i,6)*vy(i,6)+pz(i,7)*vy(i,7)+pz(i,8)*vy(i,8)+pz(i,9)*vy(i,9)+pz(i,10)*vy(i,10)
80 dzx(i)=px(i,1)*vz(i,1)+px(i,2)*vz(i,2)+px(i,3)*vz(i,3)+px(i,4)*vz(i,4)+px(i,5)*vz(i,5)
81 . +px(i,6)*vz(i,6)+px(i,7)*vz(i,7)+px(i,8)*vz(i,8)+px(i,9)*vz(i,9)+px(i,10)*vz(i,10)
82 dzy(i)=py(i,1)*vz(i,1)+py(i,2)*vz(i,2)+py(i,3)*vz(i,3)+py(i,4)*vz(i,4)+py(i,5)*vz(i,5)
83 . +py(i,6)*vz(i,6)+py(i,7)*vz(i,7)+py(i,8)*vz(i,8)+py(i,9)*vz(i,9)+py(i,10)*vz(i,10)
84 ENDDO
85
86 RETURN
87