cneveci.F File Reference

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

Go to the source code of this file.

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

Function/Subroutine Documentation

cneveci()

subroutine cneveci	(	integer	jft,
		integer	jlt,
		intent(out)	area,
		intent(in)	x1,
		intent(in)	x2,
		intent(in)	x3,
		intent(in)	x4,
		intent(in)	y1,
		intent(in)	y2,
		intent(in)	y3,
		intent(in)	y4,
		intent(in)	z1,
		intent(in)	z2,
		intent(in)	z3,
		intent(in)	z4,
		intent(out)	e1x,
		intent(out)	e2x,
		intent(out)	e3x,
		intent(out)	e1y,
		intent(out)	e2y,
		intent(out)	e3y,
		intent(out)	e1z,
		intent(out)	e2z,
		intent(out)	e3z )

Definition at line 32 of file cneveci.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, DIMENSION(MVSIZ), INTENT(IN)  :: 
     .                           x1,x2,x3,x4,y1,y2,y3,y4,z1,z2,z3,z4
      my_real, DIMENSION(MVSIZ), INTENT(OUT) :: area,
     .                   e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,I1
      my_real, DIMENSION(MVSIZ)  :: x21,y21,z21,x31,y31,z31,x42,y42,z42,
     .                              rx, ry, rz, sx, sy, sz,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)
        x42(i)=x4(i)-x2(i)
        y42(i)=y4(i)-y2(i)
        z42(i)=z4(i)-z2(i)
      ENDDO
C
      DO i=jft,jlt
        rx(i)=x31(i)-x42(i)
        ry(i)=y31(i)-y42(i)
        rz(i)=z31(i)-z42(i)
        sx(i)=x31(i)+x42(i)
        sy(i)=y31(i)+y42(i)
        sz(i)=z31(i)+z42(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) = fourth*sum(i)
      ENDDO
c-----------
      RETURN