i3dis3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i3dis3 (igimp, nty, dist, x1, x2, x3, x4, xi, y1, y2, y3, y4, yi, z1, z2, z3, z4, zi, xface, n1, n2, n3, ssc, ttc, alp, ans, xp, yp, zp, h1, h2, h3, h4, lft, llt)

Function/Subroutine Documentation

i3dis3()

subroutine i3dis3	(	integer	igimp,
		integer	nty,
			dist,
		intent(in)	x1,
		intent(in)	x2,
		intent(in)	x3,
		intent(in)	x4,
		intent(in)	xi,
		intent(in)	y1,
		intent(in)	y2,
		intent(in)	y3,
		intent(in)	y4,
		intent(in)	yi,
		intent(in)	z1,
		intent(in)	z2,
		intent(in)	z3,
		intent(in)	z4,
		intent(in)	zi,
		intent(inout)	xface,
		intent(in)	n1,
		intent(in)	n2,
		intent(in)	n3,
		intent(inout)	ssc,
		intent(inout)	ttc,
		intent(in)	alp,
		intent(inout)	ans,
		intent(inout)	xp,
		intent(inout)	yp,
		intent(inout)	zp,
		intent(inout)	h1,
		intent(inout)	h2,
		intent(inout)	h3,
		intent(inout)	h4,
		integer, intent(inout)	lft,
		integer, intent(inout)	llt )

Definition at line 32 of file i3dis3.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 IGIMP,NTY
      my_real dist(*)
      my_real, DIMENSION(MVSIZ), INTENT(IN) ::  x1,x2,x3,x4,xi
      my_real, DIMENSION(MVSIZ), INTENT(IN) ::  y1,y2,y3,y4,yi
      my_real, DIMENSION(MVSIZ), INTENT(IN) ::  z1,z2,z3,z4,zi
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) ::  xface,ans
      my_real, DIMENSION(MVSIZ), INTENT(IN) ::  n1,n2,n3,alp
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: xp,yp,zp,ssc,ttc
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: h1,h2,h3,h4
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C-----------------------------------------------
      DO i=lft,llt
        IF (xface(i) == zero) cycle
        IF (abs(ssc(i))>one+alp(i) .OR. abs(ttc(i))>one+alp(i)) THEN
          xface(i)=zero
        ELSE
          IF(abs(ssc(i)) > one) ssc(i)=ssc(i)/abs(ssc(i))
          IF(abs(ttc(i)) > one) ttc(i)=ttc(i)/abs(ttc(i))
        ENDIF
      END DO
C
      igimp=0
      DO i=lft,llt
        igimp = igimp + abs(xface(i))
      END DO
      IF (igimp == 0) RETURN
C
      DO i=lft,llt
        h1(i) = fourth*(one-ttc(i))*(one-ssc(i))
        h2(i) = fourth*(one-ttc(i))*(one+ssc(i))
        h3(i) = fourth*(one+ttc(i))*(one+ssc(i))
        h4(i) = fourth*(one+ttc(i))*(one-ssc(i))
      END DO
C
      DO i=lft,llt
        xp(i)=h1(i)*x1(i)+h2(i)*x2(i)+h3(i)*x3(i)+h4(i)*x4(i)
        yp(i)=h1(i)*y1(i)+h2(i)*y2(i)+h3(i)*y3(i)+h4(i)*y4(i)
        zp(i)=h1(i)*z1(i)+h2(i)*z2(i)+h3(i)*z3(i)+h4(i)*z4(i)
      END DO
C
      DO i=lft,llt
        ans(i)= n1(i)*(xi(i)-xp(i))
     .         +n2(i)*(yi(i)-yp(i))
     .         +n3(i)*(zi(i)-zp(i))
      END DO
      IF (nty == 8) THEN
        DO i=lft,llt
          dist(i) = ans(i)
        ENDDO
      ENDIF
C-----------
      RETURN