s20defo3.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	s20defo3 (npe, px, py, pz, vx, vy, vz, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, d4, d5, d6, wxx, wyy, wzz, rho, rhoo, volnp, voln, nel)

Function/Subroutine Documentation

s20defo3()

subroutine s20defo3	(	integer	npe,
			px,
			py,
			pz,
			vx,
			vy,
			vz,
			dxx,
			dxy,
			dxz,
			dyx,
			dyy,
			dyz,
			dzx,
			dzy,
			dzz,
			d4,
			d5,
			d6,
			wxx,
			wyy,
			wzz,
			rho,
			rhoo,
			volnp,
			voln,
		integer, intent(in)	nel )

Definition at line 29 of file s20defo3.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 NPE
      my_real
     .  vx(mvsiz,*),vy(mvsiz,*),vz(mvsiz,*),
     .  px(mvsiz,*),py(mvsiz,*),pz(mvsiz,*),
     .   dxx(*), dxy(*), dxz(*),
     .   dyx(*), dyy(*), dyz(*),
     .   dzx(*), dzy(*), dzz(*), d4(*), d5(*), d6(*),
     .   wxx(*), wyy(*), wzz(*),rho(*),
     .  rhoo(*),volnp(*),voln(*)
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,N
C     REAL
      my_real
     .   dt1d2
C-----------------------------------------------
C                                                                     12
      DO i=1,nel
        rhoo(i)=rho(i)
        voln(i)=volnp(i)
      ENDDO
      n = 1
      DO i=1,nel
        dxx(i)=px(i,n)*vx(i,n)
        dyy(i)=py(i,n)*vy(i,n)
        dzz(i)=pz(i,n)*vz(i,n)
        dxy(i)=py(i,n)*vx(i,n)
        dxz(i)=pz(i,n)*vx(i,n)
        dyx(i)=px(i,n)*vy(i,n)
        dyz(i)=pz(i,n)*vy(i,n)
        dzx(i)=px(i,n)*vz(i,n)
        dzy(i)=py(i,n)*vz(i,n)
      ENDDO
!
      DO n=2,npe
       DO i=1,nel
        dxx(i)=dxx(i)+px(i,n)*vx(i,n)
        dyy(i)=dyy(i)+py(i,n)*vy(i,n)
        dzz(i)=dzz(i)+pz(i,n)*vz(i,n)
        dxy(i)=dxy(i)+py(i,n)*vx(i,n)
        dxz(i)=dxz(i)+pz(i,n)*vx(i,n)
        dyx(i)=dyx(i)+px(i,n)*vy(i,n)
        dyz(i)=dyz(i)+pz(i,n)*vy(i,n)
        dzx(i)=dzx(i)+px(i,n)*vz(i,n)
        dzy(i)=dzy(i)+py(i,n)*vz(i,n)
       ENDDO
      ENDDO
C
      dt1d2=half*dt1
C
      DO i=1,nel
       d4(i)   =dxy(i)+dyx(i)
       d5(i)   =dyz(i)+dzy(i)
       d6(i)   =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
C
      RETURN
C