s6proj3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s6proj3 (x1, x2, x3, x4, x5, x6, y1, y2, y3, y4, y5, y6, z1, z2, z3, z4, z5, z6, vx1, vx2, vx3, vx4, vx5, vx6, vy1, vy2, vy3, vy4, vy5, vy6, vz1, vz2, vz3, vz4, vz5, vz6, di, nel)

Function/Subroutine Documentation

s6proj3()

subroutine s6proj3	(		x1,
			x2,
			x3,
			x4,
			x5,
			x6,
			y1,
			y2,
			y3,
			y4,
			y5,
			y6,
			z1,
			z2,
			z3,
			z4,
			z5,
			z6,
			vx1,
			vx2,
			vx3,
			vx4,
			vx5,
			vx6,
			vy1,
			vy2,
			vy3,
			vy4,
			vy5,
			vy6,
			vz1,
			vz2,
			vz3,
			vz4,
			vz5,
			vz6,
			di,
		integer, intent(in)	nel )

Definition at line 29 of file s6proj3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
C     REAL
      my_real
     .   x1(*), x2(*), x3(*), x4(*), x5(*), x6(*), 
     .   y1(*), y2(*), y3(*), y4(*), y5(*), y6(*), 
     .   z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), 
     .  vx1(*), vx2(*), vx3(*), vx4(*), vx5(*), vx6(*), 
     .  vy1(*), vy2(*), vy3(*), vy4(*), vy5(*), vy6(*), 
     .  vz1(*), vz2(*), vz3(*), vz4(*), vz5(*), vz6(*), 
     .  di(mvsiz,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   xvy,xvz,yvx,yvz,zvx,zvy, 
     .   ar(3),alr(3)
C-----------------------------------------------
      DO i=1,nel
        xvy= x1(i)*vy1(i)+x2(i)*vy2(i)+x3(i)*vy3(i)
     .      +x4(i)*vy4(i)+x5(i)*vy5(i)+x6(i)*vy6(i)
        xvz= x1(i)*vz1(i)+x2(i)*vz2(i)+x3(i)*vz3(i)
     .      +x4(i)*vz4(i)+x5(i)*vz5(i)+x6(i)*vz6(i)
        yvx= y1(i)*vx1(i)+y2(i)*vx2(i)+y3(i)*vx3(i)
     .      +y4(i)*vx4(i)+y5(i)*vx5(i)+y6(i)*vx6(i)
        yvz= y1(i)*vz1(i)+y2(i)*vz2(i)+y3(i)*vz3(i)
     .      +y4(i)*vz4(i)+y5(i)*vz5(i)+y6(i)*vz6(i)
        zvy= z1(i)*vy1(i)+z2(i)*vy2(i)+z3(i)*vy3(i)
     .      +z4(i)*vy4(i)+z5(i)*vy5(i)+z6(i)*vy6(i)
        zvx= z1(i)*vx1(i)+z2(i)*vx2(i)+z3(i)*vx3(i)
     .      +z4(i)*vx4(i)+z5(i)*vx5(i)+z6(i)*vx6(i)
        ar(1)=-zvy+yvz
        ar(2)= zvx-xvz
        ar(3)=-yvx+xvy
        alr(1) =di(i,1)*ar(1)+di(i,4)*ar(2)+di(i,5)*ar(3)
        alr(2) =di(i,4)*ar(1)+di(i,2)*ar(2)+di(i,6)*ar(3)
        alr(3) =di(i,5)*ar(1)+di(i,6)*ar(2)+di(i,3)*ar(3)
C
        vx1(i) = vx1(i)-(z1(i)*alr(2)-y1(i)*alr(3))
        vy1(i) = vy1(i)-(x1(i)*alr(3)-z1(i)*alr(1))
        vz1(i) = vz1(i)-(y1(i)*alr(1)-x1(i)*alr(2))
        vx2(i) = vx2(i)-(z2(i)*alr(2)-y2(i)*alr(3))
        vy2(i) = vy2(i)-(x2(i)*alr(3)-z2(i)*alr(1))
        vz2(i) = vz2(i)-(y2(i)*alr(1)-x2(i)*alr(2))
        vx3(i) = vx3(i)-(z3(i)*alr(2)-y3(i)*alr(3))
        vy3(i) = vy3(i)-(x3(i)*alr(3)-z3(i)*alr(1))
        vz3(i) = vz3(i)-(y3(i)*alr(1)-x3(i)*alr(2))
        vx4(i) = vx4(i)-(z4(i)*alr(2)-y4(i)*alr(3))
        vy4(i) = vy4(i)-(x4(i)*alr(3)-z4(i)*alr(1))
        vz4(i) = vz4(i)-(y4(i)*alr(1)-x4(i)*alr(2))
        vx5(i) = vx5(i)-(z5(i)*alr(2)-y5(i)*alr(3))
        vy5(i) = vy5(i)-(x5(i)*alr(3)-z5(i)*alr(1))
        vz5(i) = vz5(i)-(y5(i)*alr(1)-x5(i)*alr(2))
        vx6(i) = vx6(i)-(z6(i)*alr(2)-y6(i)*alr(3))
        vy6(i) = vy6(i)-(x6(i)*alr(3)-z6(i)*alr(1))
        vz6(i) = vz6(i)-(y6(i)*alr(1)-x6(i)*alr(2))
C
      ENDDO
      RETURN