i23cor3t.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i23cor3t (jlt, x, irect, cand_e, cand_n, igap, gap, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, ix1, ix2, ix3, ix4, nsn, gap_s, gapv, gapmax, gapmin, curv_max, nsv, msr, gap_m)

Function/Subroutine Documentation

i23cor3t()

subroutine i23cor3t	(	integer	jlt,
			x,
		integer, dimension(4,*)	irect,
		integer, dimension(*)	cand_e,
		integer, dimension(*)	cand_n,
		integer	igap,
			gap,
			x1,
			x2,
			x3,
			x4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			xi,
			yi,
			zi,
		integer, dimension(mvsiz)	ix1,
		integer, dimension(mvsiz)	ix2,
		integer, dimension(mvsiz)	ix3,
		integer, dimension(mvsiz)	ix4,
		integer	nsn,
			gap_s,
			gapv,
			gapmax,
			gapmin,
			curv_max,
		integer, dimension(*)	nsv,
		integer, dimension(*)	msr,
			gap_m )

Definition at line 30 of file i23cor3t.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IRECT(4,*), CAND_E(*), CAND_N(*), NSV(*), MSR(*),
     .    JLT,IDT, NOINT, NSN, IGAP ,ITYP
      INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ)
C     REAL
      my_real
     .   x(3,*), gapv(*), gap_s(*),curv_max(*), gap_m(*),
     .   gap, 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)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,J  ,IL, L, NN, IG,JFT
C-----------------------------------------------
      IF(igap==0)THEN
        DO i=1,jlt
          gapv(i)=gap
        END DO
      ELSE
        DO i=1,jlt
          j = cand_n(i)
          IF(j<=nsn) THEN
            gapv(i)=gap_s(j)+gap_m(cand_e(i))
          ELSE
            ig = j-nsn
             gapv(i)=xrem(9,ig)+gap_m(cand_e(i))
          END IF
          gapv(i)=min(gapv(i),gapmax)
          gapv(i)=max(gapmin,gapv(i))
        END DO
      END IF
C
      DO i=1,jlt
        gapv(i)=sqrt(three)*gapv(i)
      END DO
C
      DO i=1,jlt
       j = cand_n(i)
       IF(j<=nsn) THEN
         ig = nsv(j)
         xi(i) = x(1,ig)
         yi(i) = x(2,ig)
         zi(i) = x(3,ig)
       ELSE
         ig = j-nsn
         xi(i) = xrem(1,ig)
         yi(i) = xrem(2,ig)
         zi(i) = xrem(3,ig)
       ENDIF
C
       l  = cand_e(i)
C
       ix1(i)=irect(1,l)
       x1(i)=x(1,ix1(i))
       y1(i)=x(2,ix1(i))
       z1(i)=x(3,ix1(i))
C
       ix2(i)=irect(2,l)
       x2(i)=x(1,ix2(i))
       y2(i)=x(2,ix2(i))
       z2(i)=x(3,ix2(i))
C
       ix3(i)=irect(3,l)
       x3(i)=x(1,ix3(i))
       y3(i)=x(2,ix3(i))
       z3(i)=x(3,ix3(i))
C
       ix4(i)=irect(4,l)
       x4(i)=x(1,ix4(i))
       y4(i)=x(2,ix4(i))
       z4(i)=x(3,ix4(i))
C
      END DO
C
      RETURN