s8edefw3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s8edefw3 (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, wxx, wyy, wzz, nel)

Function/Subroutine Documentation

s8edefw3()

subroutine s8edefw3	(		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,
			wxx,
			wyy,
			wzz,
		integer, intent(in)	nel )

Definition at line 30 of file s8edefw3.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-----------------------------------------------
#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
C     REAL
      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(*), 
     .   wxx(*), wyy(*), wzz(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I, J
C     REAL
      my_real
     .   dxy, dyx, dxz,
     .   dzy, dyz, dzx,
     .   vx17, vy17, vz17,
     .   vx28, vy28, vz28,
     .   vx35, vy35, vz35,
     .   vx46, vy46, vz46, dt1d2, dt1d
C-----------------------------------------------
      dt1d2=half*dt1
       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)
         dxy=py1(i)*vx17+py2(i)*vx28+
     .    py3(i)*vx35+py4(i)*vx46
         dxz=pz1(i)*vx17+pz2(i)*vx28+
     .    pz3(i)*vx35+pz4(i)*vx46
         dyx=px1(i)*vy17+px2(i)*vy28+
     .    px3(i)*vy35+px4(i)*vy46
         dyz=pz1(i)*vy17+pz2(i)*vy28+
     .    pz3(i)*vy35+pz4(i)*vy46
         dzx=px1(i)*vz17+px2(i)*vz28+
     .    px3(i)*vz35+px4(i)*vz46
         dzy=py1(i)*vz17+py2(i)*vz28+
     .    py3(i)*vz35+py4(i)*vz46
C
        wzz(i)=dt1d2*(dyx-dxy)
        wyy(i)=dt1d2*(dxz-dzx)
        wxx(i)=dt1d2*(dzy-dyz)
       ENDDO
C
C
      RETURN
C