inint0.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	inint0 (x, irect, nseg, nod2seg, nsv, msr, iloc, nmn, nsn, nrt, sirect, s_n2seg)

Function/Subroutine Documentation

inint0()

subroutine inint0	(	dimension(3,*), intent(in)	x,
		integer, dimension(4,sirect/4), intent(in)	irect,
		integer, dimension(*), intent(in)	nseg,
		integer, dimension(s_n2seg), intent(inout)	nod2seg,
		integer, dimension(*), intent(in)	nsv,
		integer, dimension(*), intent(in)	msr,
		integer, dimension(*), intent(inout)	iloc,
		integer, intent(in)	nmn,
		integer, intent(in)	nsn,
		integer, intent(in)	nrt,
		integer, intent(in)	sirect,
		integer, intent(in)	s_n2seg )

Parameters

[in]	sirect	array size IRECT(SIRECT)
[in]	s_n2seg	array size NOD2SEG(S_N2SEG)

Definition at line 31 of file inint0.F.

C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
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(IN) :: SIRECT !< array size IRECT(SIRECT)
      INTEGER,INTENT(IN) :: S_N2SEG !< array size NOD2SEG(S_N2SEG)
      INTEGER, INTENT(IN) :: NMN, NSN, NRT
      INTEGER, INTENT(IN) :: IRECT(4,SIRECT/4)
      INTEGER, INTENT(IN) :: NSEG(*), NSV(*), MSR(*)
      INTEGER, INTENT(INOUT) :: NOD2SEG(S_N2SEG)
      INTEGER, INTENT(INOUT) :: ILOC(*)
      my_real, INTENT(IN) :: x(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER II, I, JJ, J, KK, L1, L2, L
      my_real lower_dist, dist
      LOGICAL lFOUND
C-----------------------------------------------
C   P r e - C o n d i t i o n
C-----------------------------------------------
      IF(nrt == 0) RETURN
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      kk = 0
C  
C   Compute nearest secnd node from main node (stored in ILOC)
C      inter type1 : main<->secnd ; nearest main node (lag) from ale node (ale)
      DO ii=1,nmn
         i=msr(ii)
         lower_dist = ep20
         DO jj=1,nsn
           j=nsv(jj)
           dist = (x(1,i)-x(1,j))**2 + (x(2,i)-x(2,j))**2 + (x(3,i)-x(3,j))**2
           IF(dist <= lower_dist)THEN
             lower_dist = dist
             iloc(ii)= jj ! nearest secnd node more main node ii
           ENDIF
         ENDDO
      ENDDO
C  
C   Compute NOD2SEG : List of Connected Segments to secnd nodes
C       inter1 : main <-> secnd ; connected main segment to main nodes (lag)
      DO i=1,nrt
         DO j=1,4
           IF(irect(j,i) /= 0) THEN
             lfound=.false.
             DO kk=1,nsn
               IF(nsv(kk) == irect(j,i)) THEN
                 lfound = .true.
                 EXIT
               ENDIF
             ENDDO
 
             IF(lfound)THEN
               ii=0
               l1=nseg(kk)
               l2=nseg(kk+1)-1
               DO l=l1,l2
                 IF(nod2seg(l) /= 0) ii=ii+1
               ENDDO
               nod2seg(l1+ii)=i
             ENDIF
           ENDIF
         ENDDO
      ENDDO
C
      RETURN