30 1 PXC1, PXC2, PXC3, PXC4,
31 2 PYC1, PYC2, PYC3, PYC4,
32 3 PZC1, PZC2, PZC3, PZC4,
44#include "implicit_f.inc"
51 INTEGER,
INTENT(IN) :: NEL
54 . VX1(*), VX2(*), VX3(*), VX4(*), VX5(*), VX6(*), VX7(*), VX8(*),
55 . VY1(*), VY2(*), VY3(*), VY4(*), VY5(*), VY6(*), VY7(*), VY8(*),
56 . VZ1(*), VZ2(*), VZ3(*), VZ4(*), VZ5(*), VZ6(*), VZ7(*), VZ8(*),
57 . PXC1(*), PXC2(*), PXC3(*), PXC4(*),
58 . PYC1(*), PYC2(*), PYC3(*), PYC4(*),
59 . PZC1(*), PZC2(*), PZC3(*), PZC4(*), DVCA(*)
66 dvca(i)=pxc1(i)*(vx1(i)-vx7(i))+pxc2(i)*(vx2(i)-vx8(i))
67 . +pxc3(i)*(vx3(i)-vx5(i))+pxc4(i)*(vx4(i)-vx6(i))
68 . +pyc1(i)*(vy1(i)-vy7(i))+pyc2(i)*(vy2(i)-vy8(i))
69 . +pyc3(i)*(vy3(i)-vy5(i))+pyc4(i)*(vy4(i)-vy6(i))
70 . +pzc1(i)*(vz1(i)-vz7(i))+pzc2(i)*(vz2(i)-vz8(i))
71 . +pzc3(i)*(vz3(i)-vz5(i))+pzc4(i)*(vz4(i)-vz6(i))
subroutine s8edefoc3(pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, dvca, nel)