29 . G1X,G1Y,G1Z,G2X,G2Y,G2Z,G3X,G3Y,G3Z)
33#include "implicit_f.inc"
40 . mxx(*), myy(*), mzz(*),g1x(*),g1y(*),g1z(*),
41 . g2x(*),g2y(*),g2z(*),g3x(*),g3y(*),g3z(*)
48 . sx,sy,sz,fxx,fxy,fxz,fyx,fyy,fyz,fzx,fzy,fzz
51 fxx = mxx(i)*g1x(i)*g1x(i)+myy(i)*g1y(i)*g1y(i)+mzz(i)*g1z(i)*g1z(i)
52 fyy = mxx(i)*g2x(i)*g2x(i)+myy(i)*g2y(i)*g2y(i)+mzz(i)*g2z(i)*g2z(i)
53 fzz = mxx(i)*g3x(i)*g3x(i)+myy(i)*g3y(i)*g3y(i)+mzz(i)*g3z(i)*g3z(i)
subroutine sordeftb(jft, jlt, mxx, myy, mzz, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z)