i10corp3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i10corp3 (jlt, x, irect, cand_e, cand_n, stif, h1, h2, h3, h4, n1, n2, n3, nt1, nt2, nt3, ix1, ix2, ix3, ix4, nsvg, vxi, vyi, vzi, msi, dxi, dyi, dzi, nsn, nin, jlt_new, lrem)

Function/Subroutine Documentation

i10corp3()

subroutine i10corp3	(	integer	jlt,
			x,
		integer, dimension(4,*)	irect,
		integer, dimension(*)	cand_e,
		integer, dimension(*)	cand_n,
			stif,
			h1,
			h2,
			h3,
			h4,
			n1,
			n2,
			n3,
			nt1,
			nt2,
			nt3,
		integer, dimension(mvsiz)	ix1,
		integer, dimension(mvsiz)	ix2,
		integer, dimension(mvsiz)	ix3,
		integer, dimension(mvsiz)	ix4,
		integer, dimension(mvsiz)	nsvg,
			vxi,
			vyi,
			vzi,
			msi,
			dxi,
			dyi,
			dzi,
		integer	nsn,
		integer	nin,
		integer	jlt_new,
		integer	lrem )

Definition at line 31 of file i10corp3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE tri7box
      USE imp_intm
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 IRECT(4,*), CAND_E(*), CAND_N(*),
     .         JLT,JLT_NEW,NIN,NSN,LREM
      INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
     .        NSVG(MVSIZ)
C     REAL
C     REAL
      my_real
     .   vxi(mvsiz), vyi(mvsiz), vzi(mvsiz), msi(mvsiz),
     .   x(3,*), stif(*),n1(mvsiz), n2(mvsiz), n3(mvsiz),
     .   nt1(mvsiz), nt2(mvsiz), nt3(mvsiz),
     .   h1(mvsiz), h2(mvsiz), h3(mvsiz), h4(mvsiz),
     .   dxi(mvsiz),dyi(mvsiz),dzi(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,J  ,NI, L, NN, NS ,NE
      my_real
     .   t1,t2,t3,norm
C-----------------------------------------------
C
C
        DO i=1,jlt
         ni = cand_n(i)
         l  = cand_e(i)
         IF(ni>nsn)THEN
           nn = ni - nsn
           jlt_new = jlt_new + 1
           vxi(jlt_new)= vfi(nin)%P(1,nn)
           vyi(jlt_new)= vfi(nin)%P(2,nn)
           vzi(jlt_new)= vfi(nin)%P(3,nn)
           msi(jlt_new)= msfi(nin)%P(nn)
           ix1(jlt_new)  = irect(1,l)
           ix2(jlt_new)  = irect(2,l)
           ix3(jlt_new)  = irect(3,l)
           ix4(jlt_new)  = irect(4,l)
             ns=ind_int(nin)%P(nn)
             ne=shf_int(nin)+jlt_new+lrem
           nsvg(jlt_new) = ns
           h1(jlt_new)  = h_e(1,ne)
           h2(jlt_new)  = h_e(2,ne)
           h3(jlt_new)  = h_e(3,ne)
           h4(jlt_new)  = h_e(4,ne)
             stif(jlt_new)=stifs(ne)
             n1(jlt_new)=n_e(1,ne)
             n2(jlt_new)=n_e(2,ne)
             n3(jlt_new)=n_e(3,ne)
C----------displacement
             dxi(jlt_new)=dfi(1,ns)
             dyi(jlt_new)=dfi(2,ns)
             dzi(jlt_new)=dfi(3,ns)
           t1 = x(1,ix3(jlt_new)) - x(1,ix1(jlt_new))
           t2 = x(2,ix3(jlt_new)) - x(2,ix1(jlt_new))
           t3 = x(3,ix3(jlt_new)) - x(3,ix1(jlt_new))
           norm = one/sqrt(t1*t1+t2*t2+t3*t3)
           nt1(jlt_new) = t1*norm
           nt2(jlt_new) = t2*norm
           nt3(jlt_new) = t3*norm
C
         END IF
C
        END DO
C
      RETURN