facepernode.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	facepernode (iponoelfa, inoelfa, lakonfa, ipkonfa, konfa, nsurfs, inoelsize)

Function/Subroutine Documentation

facepernode()

subroutine facepernode	(	integer, dimension(*)	iponoelfa,
		integer, dimension(3,*)	inoelfa,
		character*8, dimension(*)	lakonfa,
		integer, dimension(*)	ipkonfa,
		integer, dimension(*)	konfa,
		integer	nsurfs,
		integer	inoelsize
	)

!
      implicit none
!
      character*8 lakonfa(*)
!
      integer iponoelfa(*),inoelfa(3,*),ipkonfa(*),konfa(*),i,j,nsurfs,
     &  inoelfree,nope,indexe,node,inoelsize
!
!     determining the external faces belonging to the nodes on
!     the surface of the structure
!
      inoelfree=1
      do i=1,nsurfs
         if(ipkonfa(i).lt.0) cycle
         if(lakonfa(i)(2:2).eq.'8') then
            nope=8
         elseif(lakonfa(i)(2:2).eq.'4') then
            nope=4
         elseif(lakonfa(i)(2:2).eq.'6') then
            nope=6
         elseif(lakonfa(i)(2:2).eq.'3') then
            nope=3
         endif
         indexe=ipkonfa(i)
         do j=1,nope
            node=konfa(indexe+j)
            inoelfa(1,inoelfree)=i
            inoelfa(2,inoelfree)=j
            inoelfa(3,inoelfree)=iponoelfa(node)
            iponoelfa(node)=inoelfree
            inoelfree=inoelfree+1
         enddo
      enddo
!
!     size of field inoelfa
!
      inoelsize=inoelfree-1
!
      return