thickness.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	thickness (dgdx, nobject, nodedesiboun, ndesiboun, objectset, xo, yo, zo, x, y, z, nx, ny, nz, co, ifree, ndesia, ndesib, iobject, ndesi, dgdxglob, nk)

Function/Subroutine Documentation

thickness()

subroutine thickness	(	real*8, dimension(ndesi,nobject)	dgdx,
		integer	nobject,
		integer, dimension(*)	nodedesiboun,
		integer	ndesiboun,
		character*81, dimension(4,*)	objectset,
		real*8, dimension(*)	xo,
		real*8, dimension(*)	yo,
		real*8, dimension(*)	zo,
		real*8, dimension(*)	x,
		real*8, dimension(*)	y,
		real*8, dimension(*)	z,
		integer, dimension(*)	nx,
		integer, dimension(*)	ny,
		integer, dimension(*)	nz,
		real*8, dimension(3,*)	co,
		integer	ifree,
		integer	ndesia,
		integer	ndesib,
		integer	iobject,
		integer	ndesi,
		real*8, dimension(2,nk,nobject)	dgdxglob,
		integer	nk
	)

!
!     calcualtion of the actual wall thickness      
!
      implicit none
!
      character*81 objectset(4,*)

      integer nobject,nodedesiboun(*),nnodesinside,i,ndesi,
     &        ndesiboun,j,k,neighbor(10),nx(*),ny(*),nz(*),
     &        istat,ndesia,ndesib,ifree,nnodes,irefnode,
     &        iactnode,iobject,nk
!
      real*8 dgdx(ndesi,nobject),xo(*),yo(*),zo(*),x(*),
     &       y(*),z(*),refdist,co(3,*),xdesi,ydesi,zdesi,      
     &       scale,actdist,dd,xrefnode,yrefnode,zrefnode,
     &       numericdist,dgdxglob(2,nk,nobject)
!   
      read(objectset(1,iobject)(61:80),'(f20.0)',iostat=istat) refdist
!      
      nnodes=1
      numericdist=0.02
!
!     check if upper or lower constraint is defined
!      
      if(objectset(1,iobject)(19:20).eq.'LE') then
         refdist=refdist*(1-numericdist)
      else
         refdist=refdist*(1+numericdist)
      endif
!
!     find minimum distance 
!
      do j=ndesia,ndesib
!         
         iactnode=nodedesiboun(j)    
         xdesi=co(1,iactnode)
         ydesi=co(2,iactnode)
         zdesi=co(3,iactnode)
!
         call near3d(xo,yo,zo,x,y,z,nx,ny,nz,xdesi,ydesi,zdesi,
     &        ifree,neighbor,nnodes)
!  
!        Calculate distance and check if node is active
!
         irefnode=neighbor(1)
         xrefnode=xo(irefnode)
         yrefnode=yo(irefnode)
         zrefnode=zo(irefnode)
         dd=(xrefnode-xdesi)**2+(yrefnode-ydesi)**2
     &      +(zrefnode-zdesi)**2
         actdist=dsqrt(dd)
!
         if(objectset(1,iobject)(19:20).eq.'LE') then
            if(actdist.gt.refdist) then
               dgdxglob(2,iactnode,iobject)=1.0d0
            endif
         else
            if(actdist.lt.refdist) then
               dgdxglob(2,iactnode,iobject)=1.0d0
            endif
         endif      
      enddo
!
      return