uboun.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	uboun (boun, kstep, kinc, time, node, idof, coords, vold, mi, iponoel, inoel, ipobody, xbody, ibody)

Function/Subroutine Documentation

uboun()

subroutine uboun	(	real*8, intent(out)	boun,
		integer, intent(in)	kstep,
		integer, intent(in)	kinc,
		real*8, dimension(2), intent(in)	time,
		integer, intent(in)	node,
		integer, intent(in)	idof,
		real*8, dimension(3), intent(in)	coords,
		real*8, dimension(0:mi(2),*), intent(in)	vold,
		integer, dimension(*), intent(in)	mi,
		integer, dimension(*), intent(in)	iponoel,
		integer, dimension(2,*), intent(in)	inoel,
		integer, dimension(2,*), intent(in)	ipobody,
		real*8, dimension(7,*), intent(in)	xbody,
		integer, dimension(3,*), intent(in)	ibody
	)

!
!     user subroutine uboun
!
!
!     INPUT:
!
!     kstep              step number
!     kinc               increment number
!     time(1)            current step time
!     time(2)            current total time
!     node               node number
!     idof               degree of freedom
!     coords  (1..3)     global coordinates of the node
!     vold(0..4,1..nk)   solution field in all nodes
!                        (not available for CFD-calculations)
!                        0: temperature
!                        1: displacement in global x-direction
!                           (or mass flow rate for fluid nodes)
!                        2: displacement in global y-direction
!                        3: displacement in global z-direction
!                        4: not used
!     mi(1)              max # of integration points per element (max
!                        over all elements)
!     mi(2)              max degree of freedomm per node (max over all
!                        nodes) in fields like v(0:mi(2))...
!     iponoel(i)         the network elements to which node i belongs
!                        are stored in inoel(1,iponoel(i)),
!                        inoel(1,inoel(2,iponoel(i)))...... until
!                        inoel(2,inoel(2,inoel(2......)=0
!     inoel(1..2,*)      field containing the network elements
!     ipobody(1,i)       points to an entry in fields ibody and xbody 
!                        containing the body load applied to element i, 
!                        if any, else 0
!     ipobody(2,i)       index referring to the line in field ipobody
!                        containing a pointer to the next body load
!                        applied to element i, else 0
!     ibody(1,i)         code identifying the kind of body load i:
!                        1=centrifugal, 2=gravity, 3=generalized gravity
!     ibody(2,i)         amplitude number for load i
!     ibody(3,i)         load case number for load i
!     xbody(1,i)         size of body load i
!     xbody(2..4,i)      for centrifugal loading: point on the axis,
!                        for gravity loading with known gravity vector:
!                          normalized gravity vector
!     xbody(5..7,i)      for centrifugal loading: normalized vector on the
!                          rotation axis
!
!     OUTPUT:
!
!     boun               boundary value for degree of freedom idof
!                        in node "node"
!           
      implicit none
!
      integer kstep,kinc,node,idof,mi(*),iponoel(*),inoel(2,*),
     &  ipobody(2,*),ibody(3,*)
! 
      real*8 boun,time(2),coords(3),vold(0:mi(2),*),xbody(7,*)
!
      intent(in) kstep,kinc,time,node,idof,coords,vold,mi,
     &                 iponoel,inoel,ipobody,xbody,ibody
!
      intent(out) boun
!
      boun=300.d0/8.d0*(4.d0-coords(2)*coords(2))
!
      return