s8defo3.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	s8defo3 (vx, vy, vz, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, ipt, d1, d2, d3, d4, d5, d6, wxx, wyy, wzz, nel, jcvt)

Function/Subroutine Documentation

s8defo3()

subroutine s8defo3	(		vx,
			vy,
			vz,
			px1,
			px2,
			px3,
			px4,
			px5,
			px6,
			px7,
			px8,
			py1,
			py2,
			py3,
			py4,
			py5,
			py6,
			py7,
			py8,
			pz1,
			pz2,
			pz3,
			pz4,
			pz5,
			pz6,
			pz7,
			pz8,
		integer	ipt,
			d1,
			d2,
			d3,
			d4,
			d5,
			d6,
			wxx,
			wyy,
			wzz,
		integer, intent(in)	nel,
		integer, intent(in)	jcvt )

Definition at line 28 of file s8defo3.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, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: JCVT
C     REAL
      my_real
     .    vx(mvsiz,8), vy(mvsiz,8), vz(mvsiz,8),
     .   px1(mvsiz,8),px2(mvsiz,8),px3(mvsiz,8),px4(mvsiz,8),      
     .   px5(mvsiz,8),px6(mvsiz,8),px7(mvsiz,8),px8(mvsiz,8),      
     .   py1(mvsiz,8),py2(mvsiz,8),py3(mvsiz,8),py4(mvsiz,8),      
     .   py5(mvsiz,8),py6(mvsiz,8),py7(mvsiz,8),py8(mvsiz,8),      
     .   pz1(mvsiz,8),pz2(mvsiz,8),pz3(mvsiz,8),pz4(mvsiz,8),      
     .   pz5(mvsiz,8),pz6(mvsiz,8),pz7(mvsiz,8),pz8(mvsiz,8),
     .    d1(mvsiz,8), d2(mvsiz,8), d3(mvsiz,8), d4(mvsiz,8),
     .    d5(mvsiz,8), d6(mvsiz,8),
     .   wxx(mvsiz),wyy(mvsiz),wzz(mvsiz)  
      INTEGER IPT    
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com08_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   dxx(mvsiz), dyy(mvsiz), dzz(mvsiz),
     .   dxy(mvsiz), dyx(mvsiz), dyz(mvsiz),
     .   dzy(mvsiz), dzx(mvsiz), dxz(mvsiz), dt1d2
C-----------------------------------------------
C
      DO 30  i=1,nel
      dxx(i) =px1(i,ipt)*vx(i,1)+px2(i,ipt)*vx(i,2)
     +       +px3(i,ipt)*vx(i,3)+px4(i,ipt)*vx(i,4)
     +       +px5(i,ipt)*vx(i,5)+px6(i,ipt)*vx(i,6)
     +       +px7(i,ipt)*vx(i,7)+px8(i,ipt)*vx(i,8)
      dyy(i) =py1(i,ipt)*vy(i,1)+py2(i,ipt)*vy(i,2)
     +       +py3(i,ipt)*vy(i,3)+py4(i,ipt)*vy(i,4)
     +       +py5(i,ipt)*vy(i,5)+py6(i,ipt)*vy(i,6)
     +       +py7(i,ipt)*vy(i,7)+py8(i,ipt)*vy(i,8)
      dzz(i) =pz1(i,ipt)*vz(i,1)+pz2(i,ipt)*vz(i,2)
     +       +pz3(i,ipt)*vz(i,3)+pz4(i,ipt)*vz(i,4)
     +       +pz5(i,ipt)*vz(i,5)+pz6(i,ipt)*vz(i,6)
     +       +pz7(i,ipt)*vz(i,7)+pz8(i,ipt)*vz(i,8)
      dxy(i) =py1(i,ipt)*vx(i,1)+py2(i,ipt)*vx(i,2)
     +       +py3(i,ipt)*vx(i,3)+py4(i,ipt)*vx(i,4)
     +       +py5(i,ipt)*vx(i,5)+py6(i,ipt)*vx(i,6)
     +       +py7(i,ipt)*vx(i,7)+py8(i,ipt)*vx(i,8)
      dxz(i) =pz1(i,ipt)*vx(i,1)+pz2(i,ipt)*vx(i,2)
     +       +pz3(i,ipt)*vx(i,3)+pz4(i,ipt)*vx(i,4)
     +       +pz5(i,ipt)*vx(i,5)+pz6(i,ipt)*vx(i,6)
     +       +pz7(i,ipt)*vx(i,7)+pz8(i,ipt)*vx(i,8)
      dyx(i) =px1(i,ipt)*vy(i,1)+px2(i,ipt)*vy(i,2)
     +       +px3(i,ipt)*vy(i,3)+px4(i,ipt)*vy(i,4)
     +       +px5(i,ipt)*vy(i,5)+px6(i,ipt)*vy(i,6)
     +       +px7(i,ipt)*vy(i,7)+px8(i,ipt)*vy(i,8)
      dyz(i) =pz1(i,ipt)*vy(i,1)+pz2(i,ipt)*vy(i,2)
     +       +pz3(i,ipt)*vy(i,3)+pz4(i,ipt)*vy(i,4)
     +       +pz5(i,ipt)*vy(i,5)+pz6(i,ipt)*vy(i,6)
     +       +pz7(i,ipt)*vy(i,7)+pz8(i,ipt)*vy(i,8)
      dzx(i) =px1(i,ipt)*vz(i,1)+px2(i,ipt)*vz(i,2)
     +       +px3(i,ipt)*vz(i,3)+px4(i,ipt)*vz(i,4)
     +       +px5(i,ipt)*vz(i,5)+px6(i,ipt)*vz(i,6)
     +       +px7(i,ipt)*vz(i,7)+px8(i,ipt)*vz(i,8)
      dzy(i) =py1(i,ipt)*vz(i,1)+py2(i,ipt)*vz(i,2)
     +       +py3(i,ipt)*vz(i,3)+py4(i,ipt)*vz(i,4)
     +       +py5(i,ipt)*vz(i,5)+py6(i,ipt)*vz(i,6)
     +       +py7(i,ipt)*vz(i,7)+py8(i,ipt)*vz(i,8)
   30 CONTINUE
C
      dt1d2=half*dt1
C
      IF (jcvt/=0) THEN
       DO i=1,nel
         d4(i,ipt)= dxy(i)+dyx(i)
     .           -dt1*(dxx(i)*dxy(i)+dyx(i)*dyy(i)+dzx(i)*dzy(i))
         d5(i,ipt)= dyz(i)+dzy(i)
     .           -dt1*(dyy(i)*dyz(i)+dzy(i)*dzz(i)+dxy(i)*dxz(i))
         d6(i,ipt)= dxz(i)+dzx(i)
     .           -dt1*(dzz(i)*dzx(i)+dxz(i)*dxx(i)+dyz(i)*dyx(i))
         d1(i,ipt)= dxx(i)  
     .           -dt1d2*(dxx(i)*dxx(i)+dyx(i)*dyx(i)+dzx(i)*dzx(i))
         d2(i,ipt)= dyy(i)
     .           -dt1d2*(dyy(i)*dyy(i)+dzy(i)*dzy(i)+dxy(i)*dxy(i))
         d3(i,ipt)= dzz(i)
     .           -dt1d2*(dzz(i)*dzz(i)+dxz(i)*dxz(i)+dyz(i)*dyz(i))
        wxx(i)=zero
        wyy(i)=zero
        wzz(i)=zero
       ENDDO
      ELSE
       DO i=1,nel
        d1(i,ipt)=dxx(i)
        d2(i,ipt)=dyy(i)
        d3(i,ipt)=dzz(i)
        d4(i,ipt)=dxy(i)+dyx(i)
        d5(i,ipt)=dyz(i)+dzy(i)
        d6(i,ipt)=dxz(i)+dzx(i)
        wzz(i)=dt1d2*(dyx(i)-dxy(i))
        wyy(i)=dt1d2*(dxz(i)-dzx(i))
        wxx(i)=dt1d2*(dzy(i)-dyz(i))
       ENDDO
      ENDIF
      RETURN
C