c3evec3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	c3evec3 (jft, jlt, area, x1, x2, x3, y1, y2, y3, z1, z2, z3, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x31, y31, z31, x2l, x3l, y3l)

Function/Subroutine Documentation

c3evec3()

subroutine c3evec3	(	integer	jft,
		integer	jlt,
			area,
			x1,
			x2,
			x3,
			y1,
			y2,
			y3,
			z1,
			z2,
			z3,
			e1x,
			e2x,
			e3x,
			e1y,
			e2y,
			e3y,
			e1z,
			e2z,
			e3z,
			x31,
			y31,
			z31,
			x2l,
			x3l,
			y3l )

Definition at line 34 of file c3evec3.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT
      my_real x1(mvsiz), x2(mvsiz), x3(mvsiz),
     .        y1(mvsiz), y2(mvsiz), y3(mvsiz),
     .        z1(mvsiz), z2(mvsiz), z3(mvsiz),
     .        e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
     .        e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
     .        e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),
     .        x31(mvsiz), y31(mvsiz), z31(mvsiz),
     .        x32(mvsiz), y32(mvsiz), z32(mvsiz),
     .        x21(mvsiz), y21(mvsiz), z21(mvsiz), area(mvsiz),
     .        x2l(mvsiz), x3l(mvsiz), y3l(mvsiz)
C-----------------------------------------------
      INTEGER I
      my_real sum
C-----------------------------------------------
      DO i=jft,jlt
        x21(i)=x2(i)-x1(i)
        y21(i)=y2(i)-y1(i)
        z21(i)=z2(i)-z1(i)
        x31(i)=x3(i)-x1(i)
        y31(i)=y3(i)-y1(i)
        z31(i)=z3(i)-z1(i)
        x32(i)=x3(i)-x2(i)
        y32(i)=y3(i)-y2(i)
        z32(i)=z3(i)-z2(i)
      ENDDO
C
      DO i=jft,jlt
        e1x(i)= x21(i)
        e1y(i)= y21(i)
        e1z(i)= z21(i)
        x2l(i) = sqrt(e1x(i)*e1x(i)+e1y(i)*e1y(i)+e1z(i)*e1z(i))
        e1x(i)=e1x(i)/x2l(i)
        e1y(i)=e1y(i)/x2l(i)
        e1z(i)=e1z(i)/x2l(i)
      ENDDO
C
      DO i=jft,jlt
        e3x(i)=y31(i)*z32(i)-z31(i)*y32(i)
        e3y(i)=z31(i)*x32(i)-x31(i)*z32(i)
        e3z(i)=x31(i)*y32(i)-y31(i)*x32(i)
        sum = sqrt(e3x(i)*e3x(i)+e3y(i)*e3y(i)+e3z(i)*e3z(i))
        e3x(i)=e3x(i)/sum
        e3y(i)=e3y(i)/sum
        e3z(i)=e3z(i)/sum
        area(i) = half * sum
      ENDDO
C
      DO i=jft,jlt
        e2x(i)=e3y(i)*e1z(i)-e3z(i)*e1y(i)
        e2y(i)=e3z(i)*e1x(i)-e3x(i)*e1z(i)
        e2z(i)=e3x(i)*e1y(i)-e3y(i)*e1x(i)
        sum = sqrt(e2x(i)*e2x(i)+e2y(i)*e2y(i)+e2z(i)*e2z(i))
        e2x(i)=e2x(i)/sum
        e2y(i)=e2y(i)/sum
        e2z(i)=e2z(i)/sum
        y3l(i)=e2x(i)*x31(i)+e2y(i)*y31(i)+e2z(i)*z31(i)
        x3l(i)=e1x(i)*x31(i)+e1y(i)*y31(i)+e1z(i)*z31(i)
      ENDDO
C-----------
      RETURN