i7cor3t.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i7cor3t (jlt, x, irect, nsv, cand_e, cand_n, igap, gap, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, stif, ix1, ix2, ix3, ix4, nsn, gap_s, gap_m, gapv, gapmax, gapmin, curv_max, ityp, nin, gap_s_l, gap_m_l, intth, drad, dgapload)

Function/Subroutine Documentation

i7cor3t()

subroutine i7cor3t	(	integer	jlt,
			x,
		integer, dimension(4,*)	irect,
		integer, dimension(*)	nsv,
		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,
			stif,
		integer, dimension(mvsiz)	ix1,
		integer, dimension(mvsiz)	ix2,
		integer, dimension(mvsiz)	ix3,
		integer, dimension(mvsiz)	ix4,
		integer	nsn,
			gap_s,
			gap_m,
			gapv,
			gapmax,
			gapmin,
			curv_max,
		integer	ityp,
		integer	nin,
			gap_s_l,
			gap_m_l,
		integer	intth,
		intent(in)	drad,
		intent(in)	dgapload )

Definition at line 31 of file i7cor3t.F.

C============================================================================
C  cette routine est appelee par : I10STO(/int10/i10sto.F)
C                                  I7STO(/int7/i7sto.F)
C----------------------------------------------------------------------------
C  cette routine appelle : -
C============================================================================
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-----------------------------------------------
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(*),
     .    JLT,IDT, NOINT, NSN, IGAP ,ITYP, NIN, INTTH
      INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ)
      my_real
     .   x(3,*), gapv(*), gap_s(*), gap_m(*),curv_max(*),
     .   gap, gapmax, gapmin
      my_real , INTENT(IN) :: dgapload ,drad
      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), stif(mvsiz),
     .   gap_s_l(*),gap_m_l(*)
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,IADD
C-----------------------------------------------
C
       IF(igap==0)THEN
         DO i=1,jlt
           gapv(i)=max(gap+dgapload,drad)
         END DO
       ELSEIF(igap == 3)THEN 
         iadd = 9 
         DO i=1,jlt
           j = cand_n(i)
           IF(j<=nsn) THEN 
             gapv(i)=gap_s(j)+gap_m(cand_e(i))
             gapv(i)=min(gap_s_l(j)+gap_m_l(cand_e(i)),gapv(i))
           ELSE
             ig = j-nsn
             gapv(i)=xrem(9,ig)+gap_m(cand_e(i))
             gapv(i)=min(xrem(10,ig)+gap_m_l(cand_e(i)),gapv(i))
           END IF
           gapv(i)=min(gapv(i),gapmax)
           gapv(i)=max(gapmin,gapv(i))
           gapv(i)=max(drad,gapv(i)+dgapload)
         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))
           gapv(i)=max(drad,gapv(i)+dgapload)
         END DO
       END IF
       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))
       ENDDO
 
      IF(ityp == 7)THEN
        DO i=1,jlt
          gapv(i) = gapv(i) + curv_max(cand_e(i))
        END DO
      ENDIF
C
      RETURN