networkneighbor.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	networkneighbor (nelem, node, nelemnei, nodenei, ibranch, iponoel, inoel, ipkon, kon)

Function/Subroutine Documentation

networkneighbor()

subroutine networkneighbor	(	integer	nelem,
		integer	node,
		integer	nelemnei,
		integer	nodenei,
		integer	ibranch,
		integer, dimension(*)	iponoel,
		integer, dimension(2,*)	inoel,
		integer, dimension(*)	ipkon,
		integer, dimension(*)	kon
	)

!
!     looks for the neighboring element and neighboring end node of
!     node "node" of element "nelem". If the neighboring end node
!     belongs to more than 2 elements ibranch=1, else ibranch=0
!     
      implicit none
!
      integer nelem,node,nelemnei,nodenei,ibranch,index,indexe,
     &  iponoel(*),inoel(2,*),ipkon(*),kon(*),iel
!
      nelemnei=0
      ibranch=0
!
      index=iponoel(node)
      if(index.eq.0) then
         write(*,*) '*ERROR in networkneighbor:node',node
         write(*,*) '       does not belong to network element',nelem
         call exit(201)
      endif
!
      do
         iel=inoel(1,index)
!
         if(iel.eq.nelem) then
            index=inoel(2,index)
            if(index.eq.0) exit
            cycle
         endif
!
!        neighboring element; check whether a neighboring element
!        was already found
!
         if(nelemnei.ne.0) then
            ibranch=1
            exit
         endif
!
         nelemnei=iel
         indexe=ipkon(iel)
         if(kon(indexe+1).eq.node) then
            nodenei=kon(indexe+3)
         else
            nodenei=kon(indexe+1)
         endif
         index=inoel(2,index)
         if(index.eq.0) exit
      enddo
!
      return