i23cor3t.F File Reference

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

Go to the source code of this file.

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

Function/Subroutine Documentation

i23cor3t()

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

Definition at line 29 of file i23cor3t.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 IGAP, IRECT(4,*), NSV(*), CAND_E(*), CAND_N(*), MSR(*)
C     REAL
      my_real
     .        gap, gapmin, gapmax, 
     .        x(3,*), gapv(*), gap_s(*), gap_m(*)
      INTEGER, DIMENSION(MVSIZ), INTENT(INOUT) :: IX1,IX2,IX3,IX4,NSVG
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: x1,x2,x3,x4
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: y1,y2,y3,y4
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: z1,z2,z3,z4
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: xi,yi,zi
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect07_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, IL, L, NN, IG
C-----------------------------------------------
      IF(igap==0)THEN
        DO i=lft,llt
          gapv(i)=gap
        END DO
      ELSE
        DO i=lft,llt
          gapv(i)=gap_s(cand_n(i))+gap_m(cand_e(i))
          gapv(i)=min(gapv(i),gapmax)
          gapv(i)=max(gapmin,gapv(i))
        END DO
      END IF
C
      DO i=lft,llt
        gapv(i)=sqrt(three)*gapv(i)
      END DO
C
      DO i=lft,llt
       ig = nsv(cand_n(i))
       nsvg(i) = ig
       xi(i) = x(1,ig)
       yi(i) = x(2,ig)
       zi(i) = x(3,ig)
      ENDDO
C
      DO i=lft,llt
C
        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)
C
      ENDDO
C
      DO i=lft,llt
C
        nn=ix1(i)
        x1(i)=x(1,nn)
        y1(i)=x(2,nn)
        z1(i)=x(3,nn)
C
        nn=ix2(i)
        x2(i)=x(1,nn)
        y2(i)=x(2,nn)
        z2(i)=x(3,nn)
C
        nn=ix3(i)
        x3(i)=x(1,nn)
        y3(i)=x(2,nn)
        z3(i)=x(3,nn)
C
        nn=ix4(i)
        x4(i)=x(1,nn)
        y4(i)=x(2,nn)
        z4(i)=x(3,nn)
C
      ENDDO
C
      RETURN