pdefo3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	pdefo3 (v, exx, exy, exz, al, nc1, nc2, nc3, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, vx1g, vx2g, vy1g, vy2g, vz1g, vz2g, nel)

Function/Subroutine Documentation

pdefo3()

subroutine pdefo3	(		v,
			exx,
			exy,
			exz,
			al,
		integer, dimension(mvsiz)	nc1,
		integer, dimension(mvsiz)	nc2,
		integer, dimension(mvsiz)	nc3,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z,
			vx1g,
			vx2g,
			vy1g,
			vy2g,
			vz1g,
			vz2g,
		integer, intent(in)	nel )

Definition at line 28 of file pdefo3.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
      INTEGER NC1(MVSIZ),NC2(MVSIZ),NC3(MVSIZ)
      my_real
     .   v(3,*),exx(mvsiz),exy(mvsiz),exz(mvsiz),al(mvsiz),
     .   e1x(mvsiz),e1y(mvsiz),e1z(mvsiz),e2x(mvsiz),e2y(mvsiz),
     .   e2z(mvsiz),e3x(mvsiz),e3y(mvsiz),e3z(mvsiz),vx1g(mvsiz),
     .   vx2g(mvsiz),vy1g(mvsiz),vy2g(mvsiz),vz1g(mvsiz),vz2g(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   vx1(mvsiz),vx2(mvsiz),vy1(mvsiz),vy2(mvsiz),vz1(mvsiz),vz2(mvsiz)
C-----------------------------------------------
      DO i=1,nel
        vx1g(i)=v(1,nc1(i))
        vy1g(i)=v(2,nc1(i))
        vz1g(i)=v(3,nc1(i))
        vx2g(i)=v(1,nc2(i))
        vy2g(i)=v(2,nc2(i))
        vz2g(i)=v(3,nc2(i))
      ENDDO
C
      DO i=1,nel
        vx1(i)=e1x(i)*vx1g(i)+e1y(i)*vy1g(i)+e1z(i)*vz1g(i)
        vy1(i)=e2x(i)*vx1g(i)+e2y(i)*vy1g(i)+e2z(i)*vz1g(i)
        vz1(i)=e3x(i)*vx1g(i)+e3y(i)*vy1g(i)+e3z(i)*vz1g(i)
        vx2(i)=e1x(i)*vx2g(i)+e1y(i)*vy2g(i)+e1z(i)*vz2g(i)
        vy2(i)=e2x(i)*vx2g(i)+e2y(i)*vy2g(i)+e2z(i)*vz2g(i)
        vz2(i)=e3x(i)*vx2g(i)+e3y(i)*vy2g(i)+e3z(i)*vz2g(i)
      ENDDO
C
      DO i=1,nel
        exx(i)=(vx2(i)-vx1(i))/al(i)
        exy(i)=(vy2(i)-vy1(i))/al(i)
        exz(i)=(vz2(i)-vz1(i))/al(i)
      ENDDO
C-----------------------------------------------
      RETURN