i6for3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i6for3 (irect, msr, nsv, irtl, stf, stfn, igimp, es, em, ansmx, fmx, fmy, fmz, peni, n1, n2, n3, xface, ans, h1, h2, h3, h4, thk, fni, fxi, fyi, fzi, fx1, fx2, fx3, fx4, fy1, fy2, fy3, fy4, fz1, fz2, fz3, fz4, lft, llt, nft)

Function/Subroutine Documentation

i6for3()

subroutine i6for3	(	integer, dimension(4,*)	irect,
		integer, dimension(*)	msr,
		integer, dimension(*)	nsv,
		integer, dimension(*)	irtl,
			stf,
			stfn,
		integer	igimp,
			es,
			em,
			ansmx,
			fmx,
			fmy,
			fmz,
			peni,
		intent(in)	n1,
		intent(in)	n2,
		intent(in)	n3,
		intent(inout)	xface,
		intent(inout)	ans,
		intent(in)	h1,
		intent(in)	h2,
		intent(in)	h3,
		intent(in)	h4,
		intent(in)	thk,
		intent(inout)	fni,
		intent(inout)	fxi,
		intent(inout)	fyi,
		intent(inout)	fzi,
		intent(inout)	fx1,
		intent(inout)	fx2,
		intent(inout)	fx3,
		intent(inout)	fx4,
		intent(inout)	fy1,
		intent(inout)	fy2,
		intent(inout)	fy3,
		intent(inout)	fy4,
		intent(inout)	fz1,
		intent(inout)	fz2,
		intent(inout)	fz3,
		intent(inout)	fz4,
		integer, intent(inout)	lft,
		integer, intent(inout)	llt,
		integer, intent(inout)	nft )

Definition at line 28 of file i6for3.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(INOUT) :: LFT
      INTEGER, INTENT(INOUT) :: LLT
      INTEGER, INTENT(INOUT) :: NFT
      INTEGER IGIMP
C     REAL
      my_real
     .   ansmx, fmx, fmy, fmz, peni
      INTEGER IRECT(4,*), MSR(*), NSV(*), IRTL(*)
C     REAL
      my_real
     .   stf(*), stfn(*), es(*), em(*)
      my_real, DIMENSION(MVSIZ), INTENT(IN) :: n1,n2,n3
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: ans,xface
      my_real, DIMENSION(MVSIZ), INTENT(IN) :: h1,h2,h3,h4,thk
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: fni
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: fxi,fyi,fzi
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: fx1,fx2,fx3,fx4
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: fy1,fy2,fy3,fy4
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: fz1,fz2,fz3,fz4
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, IL, L, J3, J2, J1, I3,I2, I1
C     REAL
      my_real
     .   penet
C-----------------------------------------------
      DO i=lft,llt
        penet = (ans(i)-thk(i))*xface(i)
        ansmx = min(ansmx,penet)
        ans(i)= min(zero,penet)
        IF (ans(i) == zero) xface(i)=zero
        ans(i) = xface(i)*ans(i)
      ENDDO
C
      igimp = 0
      DO i=lft,llt
        igimp = igimp+abs(xface(i))
      ENDDO
      IF (igimp == 0) RETURN
C
      DO i=lft,llt
        fni(i)=ans(i)
        fxi(i)=n1(i)*fni(i)
        fyi(i)=n2(i)*fni(i)
        fzi(i)=n3(i)*fni(i)
      ENDDO
C
      DO i=lft,llt
        fx1(i)=fxi(i)*h1(i)
        fy1(i)=fyi(i)*h1(i)
        fz1(i)=fzi(i)*h1(i)
C
        fx2(i)=fxi(i)*h2(i)
        fy2(i)=fyi(i)*h2(i)
        fz2(i)=fzi(i)*h2(i)
C
        fx3(i)=fxi(i)*h3(i)
        fy3(i)=fyi(i)*h3(i)
        fz3(i)=fzi(i)*h3(i)
C
        fx4(i)=fxi(i)*h4(i)
        fy4(i)=fyi(i)*h4(i)
        fz4(i)=fzi(i)*h4(i)
      ENDDO
C
      DO i=lft,llt
        il=i+nft
        l =irtl(il)
        j3=3*irect(1,l)
        j2=j3-1
        j1=j2-1
        em(j1)=em(j1)+fx1(i)
        em(j2)=em(j2)+fy1(i)
        em(j3)=em(j3)+fz1(i)
C
        j3=3*irect(2,l)
        j2=j3-1
        j1=j2-1
        em(j1)=em(j1)+fx2(i)
        em(j2)=em(j2)+fy2(i)
        em(j3)=em(j3)+fz2(i)
C
        j3=3*irect(3,l)
        j2=j3-1
        j1=j2-1
        em(j1)=em(j1)+fx3(i)
        em(j2)=em(j2)+fy3(i)
        em(j3)=em(j3)+fz3(i)
C
        j3=3*irect(4,l)
        j2=j3-1
        j1=j2-1
        em(j1)=em(j1)+fx4(i)
        em(j2)=em(j2)+fy4(i)
        em(j3)=em(j3)+fz4(i)
      ENDDO
C
      DO i=lft,llt
        i3=3*(i+nft)
        i2=i3-1
        i1=i2-1
        es(i1)=es(i1)-fxi(i)
        es(i2)=es(i2)-fyi(i)
        es(i3)=es(i3)-fzi(i)
        fmx = fmx + fxi(i)
        fmy = fmy + fyi(i)
        fmz = fmz + fzi(i)
      ENDDO
C-----------
      RETURN