i23cor3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i23cor3 (jlt, nin, x, irect, nsn, nsv, cand_e, cand_n, stf, stfn, msr, ms, v, xi, yi, zi, ix1, ix2, ix3, ix4, nsvg, igsti, stif, kmin, kmax, igap, gap, gap_s, gapv, gapmax, gapmin, gap_m, vxi, vyi, vzi, msi, nodnx_sms, nsms, kinet, x1, y1, z1, x2, y2, z2, x3, y3, z3, x4, y4, z4, nx1, nx2, nx3, nx4, ny1, ny2, ny3, ny4, nz1, nz2, nz3, nz4, kini, index)

Function/Subroutine Documentation

i23cor3()

subroutine i23cor3	(	integer	jlt,
		integer	nin,
			x,
		integer, dimension(4,*)	irect,
		integer	nsn,
		integer, dimension(*)	nsv,
		integer, dimension(*)	cand_e,
		integer, dimension(*)	cand_n,
			stf,
			stfn,
		integer, dimension(*)	msr,
			ms,
			v,
			xi,
			yi,
			zi,
		integer, dimension(mvsiz)	ix1,
		integer, dimension(mvsiz)	ix2,
		integer, dimension(mvsiz)	ix3,
		integer, dimension(mvsiz)	ix4,
		integer, dimension(mvsiz)	nsvg,
		integer	igsti,
			stif,
			kmin,
			kmax,
		integer	igap,
			gap,
			gap_s,
			gapv,
			gapmax,
			gapmin,
			gap_m,
			vxi,
			vyi,
			vzi,
			msi,
		integer, dimension(*)	nodnx_sms,
		integer, dimension(mvsiz)	nsms,
		integer, dimension(*)	kinet,
			x1,
			y1,
			z1,
			x2,
			y2,
			z2,
			x3,
			y3,
			z3,
			x4,
			y4,
			z4,
			nx1,
			nx2,
			nx3,
			nx4,
			ny1,
			ny2,
			ny3,
			ny4,
			nz1,
			nz2,
			nz3,
			nz4,
		integer, dimension(*)	kini,
		integer, dimension(*)	index )

Definition at line 30 of file i23cor3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE tri7box
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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "sms_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IRECT(4,*), NSV(*), CAND_E(*), CAND_N(*),KINET(*),KINI(*),
     .        MSR(*), NODNX_SMS(*), INDEX(*),
     .        JLT, NOINT, IGAP , NSN, NIN, IGSTI
      INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
     .        NSVG(MVSIZ), NSMS(MVSIZ)
C     REAL
      my_real
     .   x(3,*), stf(*), stfn(*), gap_s(*), gap_m(*),
     .   ms(*), v(3,*),
     .   gap, kmin, kmax, gapmax, gapmin
C     REAL
      my_real
     .   x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
     .   y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
     .   z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
     .   xi(mvsiz), yi(mvsiz), zi(mvsiz), 
     .   nx1(mvsiz), nx2(mvsiz), nx3(mvsiz), nx4(mvsiz),
     .   ny1(mvsiz), ny2(mvsiz), ny3(mvsiz), ny4(mvsiz),
     .   nz1(mvsiz), nz2(mvsiz), nz3(mvsiz), nz4(mvsiz),
     .   pene(mvsiz),stif(mvsiz) ,gapv(mvsiz),
     .   vxi(mvsiz), vyi(mvsiz), vzi(mvsiz), msi(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,J ,IL, L, IG, ITMP, NN, NI
C-----------------------------------------------
      IF(igap==0)THEN
       DO i=1,jlt
         gapv(i)=gap
       ENDDO
      ELSE
         DO i=1,jlt
          IF(cand_n(i)<=nsn) THEN
            gapv(i)=gap_s(cand_n(i))+gap_m(cand_e(i))
          ELSE
            gapv(i)=gapfi(nin)%P(cand_n(i)-nsn)+gap_m(cand_e(i))
          ENDIF
           gapv(i)=min(gapv(i),gapmax)
           gapv(i)=max(gapmin,gapv(i))
         END DO
      ENDIF
C
      DO i=1,jlt
        ni = cand_n(i)
        IF(ni<=nsn)THEN
          ig = nsv(ni)
          nsvg(i) = ig
C Kini is not used in TYPE23
          kini(i) = kinet(ig)
          xi(i) = x(1,ig)
          yi(i) = x(2,ig)
          zi(i) = x(3,ig)
          vxi(i) = v(1,ig)
          vyi(i) = v(2,ig)
          vzi(i) = v(3,ig)
          msi(i)= ms(ig)
        ELSE
          nn = ni - nsn
          nsvg(i) = -nn
          kini(i) = kinfi(nin)%P(nn)
          xi(i) = xfi(nin)%P(1,nn)
          yi(i) = xfi(nin)%P(2,nn)
          zi(i) = xfi(nin)%P(3,nn)
          vxi(i)= vfi(nin)%P(1,nn)
          vyi(i)= vfi(nin)%P(2,nn)
          vzi(i)= vfi(nin)%P(3,nn)
          msi(i)= msfi(nin)%P(nn)
        END IF
      END DO
C
       DO i=1,jlt
        l=cand_e(i)
C
        ix1(i)=irect(1,l)
        ix2(i)=irect(2,l)
        ix3(i)=irect(3,l)
        ix4(i)=irect(4,l)
       END DO
C
       DO i=1,jlt
        l=cand_e(i)
C
        x1(i)=x(1,ix1(i))
        y1(i)=x(2,ix1(i))
        z1(i)=x(3,ix1(i))
C
        x2(i)=x(1,ix2(i))
        y2(i)=x(2,ix2(i))
        z2(i)=x(3,ix2(i))
C
        x3(i)=x(1,ix3(i))
        y3(i)=x(2,ix3(i))
        z3(i)=x(3,ix3(i))
C
        x4(i)=x(1,ix4(i))
        y4(i)=x(2,ix4(i))
        z4(i)=x(3,ix4(i))
C
       END DO
C
      DO i=1,jlt
       l  = cand_e(i)
       ni = cand_n(i)
       IF(ni<=nsn)THEN
         stif(i)=stf(l)*abs(stfn(ni))
       ELSE
         nn = ni - nsn
         stif(i)=stf(l)*abs(stifi(nin)%P(nn))
       END IF
      ENDDO
C
      IF(idtmins==2)THEN
       DO i=1,jlt
        IF(nsvg(i)>0)THEN
          nsms(i)=nodnx_sms(nsvg(i))
     .                       +nodnx_sms(ix1(i))+nodnx_sms(ix2(i))
     .                       +nodnx_sms(ix3(i))+nodnx_sms(ix4(i))
        ELSE
          nn=-nsvg(i)
          nsms(i)=nodnxfi(nin)%P(nn)
     .                       +nodnx_sms(ix1(i))+nodnx_sms(ix2(i))
     .                       +nodnx_sms(ix3(i))+nodnx_sms(ix4(i))
        END IF
       ENDDO
       IF(idtmins_int/=0)THEN
         DO i=1,jlt
          IF(nsms(i)==0)nsms(i)=-1
         ENDDO
       END IF
      ELSEIF(idtmins_int/=0)THEN
        DO i=1,jlt
         nsms(i)=-1
        ENDDO
      ENDIF
C
      RETURN