i11corp3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i11corp3 (jlt, x, irects, irectm, cand_e, cand_n, stif, hs1, hs2, hm1, hm2, nx, ny, nz, n1, n2, m1, m2, vxs1, vys1, vzs1, vxs2, vys2, vzs2, vxm1, vym1, vzm1, vxm2, vym2, vzm2, dxs1, dys1, dzs1, dxs2, dys2, dzs2, dxm1, dym1, dzm1, dxm2, dym2, dzm2, ms1, ms2, mm1, mm2, ms, d, v, nrts, nin, jlt_new, lrem)
subroutine	i11corkp3 (jlt, xxs1, xxs2, xys1, xys2, xzs1, xzs2, dxs1, dys1, dzs1, dxs2, dys2, dzs2, d, cand_n, nrts, n1, n2, nin, iupd)

Function/Subroutine Documentation

i11corkp3()

subroutine i11corkp3	(	integer	jlt,
			xxs1,
			xxs2,
			xys1,
			xys2,
			xzs1,
			xzs2,
			dxs1,
			dys1,
			dzs1,
			dxs2,
			dys2,
			dzs2,
			d,
		integer, dimension(*)	cand_n,
		integer	nrts,
		integer, dimension(mvsiz)	n1,
		integer, dimension(mvsiz)	n2,
		integer	nin,
		integer	iupd )

Definition at line 145 of file i11corp3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      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 CAND_N(*),JLT,NIN,NRTS,IUPD
      INTEGER N1(MVSIZ), N2(MVSIZ)
C     REAL
      my_real
     .   dxs1(mvsiz),dys1(mvsiz),dzs1(mvsiz),
     .   dxs2(mvsiz),dys2(mvsiz),dzs2(mvsiz),
     .   xxs1(mvsiz), xxs2(mvsiz), xys1(mvsiz), xys2(mvsiz),
     .   xzs1(mvsiz), xzs2(mvsiz),d(3,*)
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,NN1,NN2
C-----------------------------------------------
C
        DO i=1,jlt
         ni = cand_n(i)
         IF(ni>nrts)THEN
           nn = ni - nrts
             ns=ind_int(nin)%P(nn)
C----------displacement
             dxs1(i)=dfi(1,ns)
             dys1(i)=dfi(2,ns)
             dzs1(i)=dfi(3,ns)
             dxs2(i)=dfi(1,ns+1)
             dys2(i)=dfi(2,ns+1)
             dzs2(i)=dfi(3,ns+1)
           ELSE
            nn1 = n1(i)
            nn2 = n2(i)
            dxs1(i) = d(1,nn1)
            dys1(i) = d(2,nn1)
            dzs1(i) = d(3,nn1)
            dxs2(i) = d(1,nn2)
            dys2(i) = d(2,nn2)
            dzs2(i) = d(3,nn2)
         END IF
C----------update
         IF(iupd>0.AND.ni>nrts)THEN
          xxs1(i) = xxs1(i) + dxs1(i)
          xys1(i) = xys1(i) + dys1(i)
          xzs1(i) = xzs1(i) + dzs1(i)
          xxs2(i) = xxs2(i) + dxs2(i)
          xys2(i) = xys2(i) + dys2(i)
          xzs2(i) = xzs2(i) + dzs2(i)
         END IF
C
        END DO
C
      RETURN

i11corp3()

subroutine i11corp3	(	integer	jlt,
			x,
		integer, dimension(2,*)	irects,
		integer, dimension(2,*)	irectm,
		integer, dimension(*)	cand_e,
		integer, dimension(*)	cand_n,
			stif,
			hs1,
			hs2,
			hm1,
			hm2,
			nx,
			ny,
			nz,
		integer, dimension(mvsiz)	n1,
		integer, dimension(mvsiz)	n2,
		integer, dimension(mvsiz)	m1,
		integer, dimension(mvsiz)	m2,
			vxs1,
			vys1,
			vzs1,
			vxs2,
			vys2,
			vzs2,
			vxm1,
			vym1,
			vzm1,
			vxm2,
			vym2,
			vzm2,
			dxs1,
			dys1,
			dzs1,
			dxs2,
			dys2,
			dzs2,
			dxm1,
			dym1,
			dzm1,
			dxm2,
			dym2,
			dzm2,
			ms1,
			ms2,
			mm1,
			mm2,
			ms,
			d,
			v,
		integer	nrts,
		integer	nin,
		integer	jlt_new,
		integer	lrem )

Definition at line 31 of file i11corp3.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 IRECTS(2,*),IRECTM(2,*), CAND_E(*), CAND_N(*),
     .         JLT,JLT_NEW,NIN,NRTS,LREM
      INTEGER N1(MVSIZ), N2(MVSIZ), M1(MVSIZ), M2(MVSIZ)
C     REAL
      my_real
     .   x(3,*),d(3,*),  v(3,*), ms(*),
     .   hs1(mvsiz), hs2(mvsiz), hm1(mvsiz), hm2(mvsiz),
     .   nx(mvsiz), ny(mvsiz), nz(mvsiz), stif(mvsiz),
     .   ms1(mvsiz),ms2(mvsiz),mm1(mvsiz),mm2(mvsiz),
     .   vxs1(mvsiz),vys1(mvsiz),vzs1(mvsiz),vxs2(mvsiz),vys2(mvsiz),
     .   vzs2(mvsiz),vxm1(mvsiz),vym1(mvsiz),vzm1(mvsiz),vxm2(mvsiz),
     .   vym2(mvsiz),vzm2(mvsiz),dxs1(mvsiz),dys1(mvsiz),dzs1(mvsiz),
     .   dxs2(mvsiz),dys2(mvsiz),dzs2(mvsiz),dxm1(mvsiz),dym1(mvsiz),
     .   dzm1(mvsiz),dxm2(mvsiz),dym2(mvsiz),dzm2(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,NN1,NN2
C-----------------------------------------------
C
C
        DO i=1,jlt
         ni = cand_n(i)
         l  = cand_e(i)
         IF(ni>nrts)THEN
           nn = ni - nrts
           jlt_new = jlt_new + 1
           nn1=2*(nn-1)+1
           nn2=2*nn
           m1(jlt_new)=irectm(1,l)
           m2(jlt_new)=irectm(2,l)
           vxs1(jlt_new) = vfi(nin)%P(1,nn1)
           vys1(jlt_new) = vfi(nin)%P(2,nn1)
           vzs1(jlt_new) = vfi(nin)%P(3,nn1)
           vxs2(jlt_new) = vfi(nin)%P(1,nn2)
           vys2(jlt_new) = vfi(nin)%P(2,nn2)
           vzs2(jlt_new) = vfi(nin)%P(3,nn2)
           vxm1(jlt_new) = v(1,m1(jlt_new))
           vym1(jlt_new) = v(2,m1(jlt_new))
           vzm1(jlt_new) = v(3,m1(jlt_new))
           vxm2(jlt_new) = v(1,m2(jlt_new))
           vym2(jlt_new) = v(2,m2(jlt_new))
           vzm2(jlt_new) = v(3,m2(jlt_new))
           dxm1(jlt_new) = d(1,m1(jlt_new))
           dym1(jlt_new) = d(2,m1(jlt_new))
           dzm1(jlt_new) = d(3,m1(jlt_new))
           dxm2(jlt_new) = d(1,m2(jlt_new))
           dym2(jlt_new) = d(2,m2(jlt_new))
           dzm2(jlt_new) = d(3,m2(jlt_new))
           ms1(jlt_new) = msfi(nin)%P(nn1)
           ms2(jlt_new) = msfi(nin)%P(nn2)
           mm1(jlt_new) = ms(m1(jlt_new))
           mm2(jlt_new) = ms(m2(jlt_new))
C
             ns=ind_int(nin)%P(nn)
             ne=shf_int(nin)+jlt_new+lrem
           hs1(jlt_new)  = h_e(1,ne)
           hs2(jlt_new)  = h_e(2,ne)
           hm1(jlt_new)  = h_e(3,ne)
           hm2(jlt_new)  = h_e(4,ne)
             stif(jlt_new)=stifs(ne)
             nx(jlt_new)=n_e(1,ne)
             ny(jlt_new)=n_e(2,ne)
             nz(jlt_new)=n_e(3,ne)
C----------displacement
             dxs1(jlt_new)=dfi(1,ns)
             dys1(jlt_new)=dfi(2,ns)
             dzs1(jlt_new)=dfi(3,ns)
             dxs2(jlt_new)=dfi(1,ns+1)
             dys2(jlt_new)=dfi(2,ns+1)
             dzs2(jlt_new)=dfi(3,ns+1)
           n1(jlt_new)=ns
           n2(jlt_new)=ns+1
C
         END IF
C
        END DO
C
      RETURN