pevec3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	pevec3 (rloc, r, al, nc1, nc2, nc3, x1, x2, x3, y1, y2, y3, z1, z2, z3, rx1g, rx2g, ry1g, ry2g, rz1g, rz2g, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, nel)

Function/Subroutine Documentation

pevec3()

subroutine pevec3	(		rloc,
			r,
			al,
		integer, dimension(mvsiz)	nc1,
		integer, dimension(mvsiz)	nc2,
		integer, dimension(mvsiz)	nc3,
			x1,
			x2,
			x3,
			y1,
			y2,
			y3,
			z1,
			z2,
			z3,
			rx1g,
			rx2g,
			ry1g,
			ry2g,
			rz1g,
			rz2g,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z,
		integer, intent(in)	nel )

Definition at line 28 of file pevec3.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-----------------------------------------------
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ)
      my_real
     .   rloc(3,*),r(3,*),al(*),x1(mvsiz),x2(mvsiz),x3(mvsiz),
     .   y1(mvsiz),y2(mvsiz),y3(mvsiz),z1(mvsiz),z2(mvsiz),
     .   z3(mvsiz),rx1g(mvsiz),rx2g(mvsiz),ry1g(mvsiz),ry2g(mvsiz),
     . rz1g(mvsiz),rz2g(mvsiz),rx1(mvsiz),rx2(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
      my_real
     .   sum(mvsiz), 
     .   theta, sum2(mvsiz), sum3(mvsiz), cost(mvsiz),
     . sint(mvsiz)
C--------------------------------------------
      DO i=1,nel
        rx1g(i)=r(1,nc1(i))
        ry1g(i)=r(2,nc1(i))
        rz1g(i)=r(3,nc1(i))
        rx2g(i)=r(1,nc2(i))
        ry2g(i)=r(2,nc2(i))
        rz2g(i)=r(3,nc2(i))
      ENDDO
C
      DO i=1,nel
        e2x(i)=rloc(1,i)
        e2y(i)=rloc(2,i)
        e2z(i)=rloc(3,i)
      ENDDO
C
      DO i=1,nel
        e1x(i)=x2(i)-x1(i)
        e1y(i)=y2(i)-y1(i)
        e1z(i)=z2(i)-z1(i)
      ENDDO
C
      DO i=1,nel
        al(i)=sqrt(e1x(i)**2+e1y(i)**2+e1z(i)**2)
      ENDDO
C
      DO i=1,nel
        e1x(i)=e1x(i)/al(i)
        e1y(i)=e1y(i)/al(i)
        e1z(i)=e1z(i)/al(i)
      ENDDO
C
      DO i=1,nel
        e3x(i)=e1y(i)*e2z(i)-e1z(i)*e2y(i)
        e3y(i)=e1z(i)*e2x(i)-e1x(i)*e2z(i)
        e3z(i)=e1x(i)*e2y(i)-e1y(i)*e2x(i)
      ENDDO
C
      DO i=1,nel
        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)
      ENDDO
C--------------------------------------------
C     TORSION MOYENNE EN COORDONNEES GLOBALES
C--------------------------------------------
      DO i=1,nel
        rx1(i)=e1x(i)*rx1g(i)+e1y(i)*ry1g(i)+e1z(i)*rz1g(i)
        rx2(i)=e1x(i)*rx2g(i)+e1y(i)*ry2g(i)+e1z(i)*rz2g(i)
        theta=(rx1(i)+rx2(i))/two*dt1
        sum2(i)=sqrt(e2x(i)**2+e2y(i)**2+e2z(i)**2)
        sum3(i)=sqrt(e3x(i)**2+e3y(i)**2+e3z(i)**2)
        cost(i)=cos(theta)/sum2(i)
        sint(i) =sin(theta)/sum3(i)
      ENDDO
C
C ... it is modified.
C
      DO i=1,nel
        e2x(i)=e2x(i)*cost(i)+e3x(i)*sint(i)
        e2y(i)=e2y(i)*cost(i)+e3y(i)*sint(i)
        e2z(i)=e2z(i)*cost(i)+e3z(i)*sint(i)
      ENDDO
C
      DO i=1,nel
        sum(i)=sqrt(e2x(i)**2+e2y(i)**2+e2z(i)**2)
      ENDDO
C
      DO i=1,nel
        e2x(i)=e2x(i)/sum(i)
        e2y(i)=e2y(i)/sum(i)
        e2z(i)=e2z(i)/sum(i)
      ENDDO
C
      DO i=1,nel
        e3x(i)=e1y(i)*e2z(i)-e1z(i)*e2y(i)
        e3y(i)=e1z(i)*e2x(i)-e1x(i)*e2z(i)
        e3z(i)=e1x(i)*e2y(i)-e1y(i)*e2x(i)
      ENDDO
C
      DO i=1,nel
        sum(i)=sqrt(e3x(i)**2+e3y(i)**2+e3z(i)**2)
        e3x(i)=e3x(i)/sum(i)
        e3y(i)=e3y(i)/sum(i)
        e3z(i)=e3z(i)/sum(i)
      ENDDO
C
      DO i=1,nel
        rloc(1,i)=e2x(i)
        rloc(2,i)=e2y(i)
        rloc(3,i)=e2z(i)
      ENDDO
C--------------------------------------------
      RETURN