i12loc3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i12loc3 (xs, xm, irect, lmsr, msr, nsv, iloc, nseg, lft, llt, nft)

Function/Subroutine Documentation

i12loc3()

subroutine i12loc3	(		xs,
			xm,
		integer, dimension(4,*)	irect,
		integer, dimension(*)	lmsr,
		integer, dimension(*)	msr,
		integer, dimension(*)	nsv,
		integer, dimension(*)	iloc,
		integer, dimension(*)	nseg,
		integer, intent(inout)	lft,
		integer, intent(inout)	llt,
		integer, intent(inout)	nft )

Definition at line 28 of file i12loc3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: LFT
      INTEGER, INTENT(INOUT) :: LLT
      INTEGER, INTENT(INOUT) :: NFT
      INTEGER IRECT(4,*), LMSR(*), MSR(*), NSV(*), ILOC(*), NSEG(*)
      my_real
     .   xs(3,*),xm(3,*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, IL,  JL, JLNEW, LL2, LL1, LL, LG, J, K, M, N
      my_real
     .   gms, cms,
     .   dms, ems, fms
C-----------------------------------------------
      DO i=lft,llt
        il=i+nft
        jl=iloc(il)
        jlnew=jl
        ll2=nseg(jl+1)-1
        ll1=nseg(jl)
        gms=1.e30
        DO ll=ll1,ll2
          lg=lmsr(ll)
          j=irect(1,lg)
          k=irect(2,lg)
          m=irect(3,lg)
          n=irect(4,lg)
          cms=(xs(1,il)-xm(1,j))**2+(xs(2,il)-xm(2,j))**2
     &                             +(xs(3,il)-xm(3,j))**2
          dms=(xs(1,il)-xm(1,k))**2+(xs(2,il)-xm(2,k))**2
     &                             +(xs(3,il)-xm(3,k))**2
          ems=(xs(1,il)-xm(1,m))**2+(xs(2,il)-xm(2,m))**2
     &                             +(xs(3,il)-xm(3,m))**2
          fms=(xs(1,il)-xm(1,n))**2+(xs(2,il)-xm(2,n))**2
     &                             +(xs(3,il)-xm(3,n))**2
          IF(cms<gms) THEN
            gms=cms
            jlnew=j
          ENDIF
          IF(dms<gms) THEN
            gms=dms
            jlnew=k
          ENDIF
          IF(ems<gms) THEN
            gms=ems
            jlnew=m
          ENDIF
          IF(fms<gms) THEN
            gms=fms
            jlnew=n
          ENDIF
        END DO !LL=LL1,LL2
        iloc(il)=jlnew      
      END DO !I=LFT,LLT
 
      RETURN