intersectionpoint.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	intersectionpoint (pa, pb, xcp, t, xinters)

Function/Subroutine Documentation

intersectionpoint()

subroutine intersectionpoint	(	real*8, dimension(*)	pa,
		real*8, dimension(*)	pb,
		real*8, dimension(*)	xcp,
		real*8	t,
		real*8, dimension(*)	xinters
	)

!     
      implicit none
!
      integer k
!     
      real*8 pa(*),pb(*),xcp(*),t,xinters(*),diff,pab(3),
     &     eplane,tnull
!     
      do k=1,3
         pab(k)=pb(k)-pa(k)
      enddo
      diff=0.0
      if(abs(eplane(pab,xcp,0)).lt.1.d-13)then
         write(*,*) 'SH: IP no intersection point can be found'
         write(*,*) 'SH: IP pab paralell to plane! '
         call exit(201)
      else 
         tnull=0.0
         diff=-eplane(pa, xcp,t)/eplane(pab,xcp,tnull)
      endif  
      do k=1,3
         xinters(k)=pa(k)+diff*pab(k)
      enddo
      return