cdkevec3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cdkevec3 (jft, jlt, area, x1, x2, x3, y1, y2, y3, z1, z2, z3, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)

Function/Subroutine Documentation

cdkevec3()

subroutine cdkevec3	(	integer	jft,
		integer	jlt,
			area,
			x1,
			x2,
			x3,
			y1,
			y2,
			y3,
			z1,
			z2,
			z3,
			e1x,
			e2x,
			e3x,
			e1y,
			e2y,
			e3y,
			e1z,
			e2z,
			e3z )

Definition at line 34 of file cdkevec3.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 area(mvsiz),
     .     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)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,I1
      my_real, DIMENSION(MVSIZ) ::  sum, rx, ry, rz, sx, sy, sz
C=======================================================================
      DO i=jft,jlt
        rx(i) = x2(i) - x1(i)
        ry(i) = y2(i) - y1(i)
        rz(i) = z2(i) - z1(i)
        sx(i) = x3(i) - x1(i)
        sy(i) = y3(i) - y1(i)
        sz(i) = z3(i) - z1(i)
      ENDDO
C----------------------------
C     LOCAL SYSTEM
C----------------------------
      i1 = 0
      CALL clskew3(jft,jlt,i1,
     .     rx, ry, rz, sx, sy, sz, 
     .     e1x,e2x,e3x,e1y,e2y,e3y,e1z,e2z,e3z,sum)
C
      DO i=jft,jlt
        area(i) = half*sum(i)
      ENDDO
c-----------
      RETURN