iqel02.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	iqel02 (x, irect, lmsr, msr, nsv, iloc, irtl, nsn, nseg, crst, nor)

Function/Subroutine Documentation

iqel02()

subroutine iqel02	(	dimension(3,numnod), intent(in)	x,
		integer, dimension(4,*), intent(in)	irect,
		integer, dimension(*), intent(in)	lmsr,
		integer, dimension(*), intent(in)	msr,
		integer, dimension(*), intent(in)	nsv,
		integer, dimension(*), intent(in)	iloc,
		integer, dimension(*), intent(in)	irtl,
		integer, intent(in)	nsn,
		integer, dimension(*), intent(in)	nseg,
		dimension(2,*), intent(in)	crst,
		dimension(3,*), intent(inout)	nor )

Definition at line 28 of file iqel02.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER,INTENT(IN)    :: NSN
      INTEGER,INTENT(IN)    :: IRECT(4,*), LMSR(*), MSR(*), NSV(*), ILOC(*), IRTL(*),NSEG(*)
      my_real,INTENT(IN)    :: x(3,numnod), crst(2,*)
      my_real,INTENT(INOUT) :: nor(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER II, I, J, L, M1, M2, LK, NUM, NPT, JJ
      my_real n2, n3, ym1, ym2, zm1, zm2, ss, sck, xmg
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      DO ii=1,nsn
        i   = nsv(ii)
        j   = iloc(ii)
        l   = irtl(ii)
        m1  = msr(irect(1,l))
        m2  = msr(irect(2,l))
        ym1 = x(2,m1)
        ym2 = x(2,m2)
        zm1 = x(3,m1)
        zm2 = x(3,m2)
        !C-------------------------
        !C     NORMAL VECTOR
        !C-------------------------
        n2  = zm1-zm2
        n3  = ym2-ym1
        ss  = crst(1,ii)
        sck = abs(ss)-one
        IF(abs(sck) <= fiveem2) THEN
          !C----------------------------------
          !C     EXTREMITIES
          !C     WEIGHTED MEAN VALUES FROM NORMAL VECTORS
          !C----------------------------------
          lk=l
          num=nseg(j+1)-nseg(j)
          IF(num == 2) THEN
            npt=nseg(j)-1
            DO jj=1,num
              l=lmsr(npt+jj)
              IF(l /= lk) THEN
                m1  = msr(irect(1,l))
                m2  = msr(irect(2,l))
                ym1 = x(2,m1)
                ym2 = x(2,m2)
                zm1 = x(3,m1)
                zm2 = x(3,m2)
                n2  = n2+zm1-zm2
                n3  = n3+ym2-ym1
              ENDIF
            ENDDO !next JJ
          ENDIF
        ENDIF
        xmg       = sqrt(n2*n2+n3*n3)
        n2        = n2/xmg
        n3        = n3/xmg
        nor(1,ii) = zero
        nor(2,ii) = n2
        nor(3,ii) = n3
        ENDDO !next II
C-----------------------------------------------        
      RETURN