sdefot3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sdefot3 (px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, nel)

Function/Subroutine Documentation

sdefot3()

subroutine sdefot3	(		px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			vx1,
			vx2,
			vx3,
			vx4,
			vx5,
			vx6,
			vx7,
			vx8,
			vy1,
			vy2,
			vy3,
			vy4,
			vy5,
			vy6,
			vy7,
			vy8,
			vz1,
			vz2,
			vz3,
			vz4,
			vz5,
			vz6,
			vz7,
			vz8,
			dxx,
			dxy,
			dxz,
			dyx,
			dyy,
			dyz,
			dzx,
			dzy,
			dzz,
		integer, intent(in)	nel )

Definition at line 30 of file sdefot3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.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
      my_real
     .   vx1(*), vx2(*), vx3(*), vx4(*), vx5(*), vx6(*), vx7(*), vx8(*),
     .   vy1(*), vy2(*), vy3(*), vy4(*), vy5(*), vy6(*), vy7(*), vy8(*),
     .   vz1(*), vz2(*), vz3(*), vz4(*), vz5(*), vz6(*), vz7(*), vz8(*),
     .   px1(*), px2(*), px3(*), px4(*),
     .   py1(*), py2(*), py3(*), py4(*),
     .   pz1(*), pz2(*), pz3(*), pz4(*), 
     .   dxx(*), dxy(*), dxz(*),
     .   dyx(*), dyy(*), dyz(*),
     .   dzx(*), dzy(*), dzz(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I, J
      my_real
     .   vx17, vy17, vz17,
     .   vx28, vy28, vz28,
     .   vx35, vy35, vz35,
     .   vx46, vy46, vz46
C-----------------------------------------------
      DO i=1,nel
        vx17=vx1(i)-vx7(i)
        vx28=vx2(i)-vx8(i)
        vx35=vx3(i)-vx5(i)
        vx46=vx4(i)-vx6(i)
        vy17=vy1(i)-vy7(i)
        vy28=vy2(i)-vy8(i)
        vy35=vy3(i)-vy5(i)
        vy46=vy4(i)-vy6(i)
        vz17=vz1(i)-vz7(i)
        vz28=vz2(i)-vz8(i)
        vz35=vz3(i)-vz5(i)
        vz46=vz4(i)-vz6(i)
!
        dxx(i)=px1(i)*vx17+px2(i)*vx28+
     .         px3(i)*vx35+px4(i)*vx46
        dyy(i)=py1(i)*vy17+py2(i)*vy28+
     .         py3(i)*vy35+py4(i)*vy46
        dzz(i)=pz1(i)*vz17+pz2(i)*vz28+
     .         pz3(i)*vz35+pz4(i)*vz46
        dxy(i)=py1(i)*vx17+py2(i)*vx28+
     .         py3(i)*vx35+py4(i)*vx46
        dxz(i)=pz1(i)*vx17+pz2(i)*vx28+
     .         pz3(i)*vx35+pz4(i)*vx46
        dyx(i)=px1(i)*vy17+px2(i)*vy28+
     .         px3(i)*vy35+px4(i)*vy46
        dyz(i)=pz1(i)*vy17+pz2(i)*vy28+
     .         pz3(i)*vy35+pz4(i)*vy46
        dzx(i)=px1(i)*vz17+px2(i)*vz28+
     .         px3(i)*vz35+px4(i)*vz46
        dzy(i)=py1(i)*vz17+py2(i)*vz28+
     .         py3(i)*vz35+py4(i)*vz46
      ENDDO
C-----------
      RETURN