addizdofcload.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	addizdofcload (nodeforc, ndirforc, nactdof, mi, izdof, nzdof, iforc, iznode, nznode, nk, imdnode, nmdnode, xforc)

Function/Subroutine Documentation

addizdofcload()

subroutine addizdofcload	(	integer, dimension(2,*)	nodeforc,
		integer, dimension(*)	ndirforc,
		integer, dimension(0:mi(2),*)	nactdof,
		integer, dimension(*)	mi,
		integer, dimension(*)	izdof,
		integer	nzdof,
		integer	iforc,
		integer, dimension(*)	iznode,
		integer	nznode,
		integer	nk,
		integer, dimension(*)	imdnode,
		integer	nmdnode,
		real*8, dimension(*)	xforc
	)

!
!     adds the dof in which a point force was applied to iznode, izdof
!     and to ** imdnode if user-defined load **
!     (needed in dyna.c and steadystate.c)
!
      implicit none
!
      integer nodeforc(2,*),ndirforc(*),iforc,node,j,jdof,mi(*),nk,
     &  nactdof(0:mi(2),*),izdof(*),nzdof,iznode(*),nznode,nodebasis,
     &  imdnode(*),nmdnode
!
      real*8 xforc(*)
!
      node=nodeforc(1,iforc)
!
!     adding the nodes in the basis sector to iznode
!
      nodebasis=mod(node,nk)
      call addimd(iznode,nznode,nodebasis)
!
!     user-defined load
!
      if((xforc(iforc).lt.1.2357111318d0).and.
     &     (xforc(iforc).gt.1.2357111316d0)) then
         call addimd(imdnode,nmdnode,node)
      endif
!
      j=ndirforc(iforc)
!
!     C-convention!
!
      jdof=nactdof(j,node)-1
!
      if(jdof.gt.0) call addimd(izdof,nzdof,jdof)
!
      return