cvec2dcos.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cvec2dcos (vx, vy, dircos)
subroutine	cdispvec2dcos (dx, dy, len, dircos)

Function/Subroutine Documentation

cdispvec2dcos()

subroutine cdispvec2dcos	(	real*8	dx,
		real*8	dy,
		real*8	len,
		type(coord)	dircos
	)

Definition at line 41 of file cvec2dcos.f.

!          len:  output
!        dircos: output
!
      implicit none
!----      include 'Zcoord.h'
#include  "Zcoord.h"
      real*8 dx, dy, len
      type(coord)::dircos
!
      if(dx .eq. 0. .and.  dy .eq. 0.) then
         len = 0.
         dircos%r(1) = 0.
         dircos%r(2) = 0
         dircos%r(3) = 1.
      else
         len = sqrt(dx**2 + dy**2)
         dircos%r(1) = dx/len
         dircos%r(2) = dy/len
         dircos%r(3) = 0.
      endif

cvec2dcos()

subroutine cvec2dcos	(	real*8	vx,
		real*8	vy,
		type(coord)	dircos
	)

Definition at line 5 of file cvec2dcos.f.

References d0.

!        small vector( vx, vy ) is converted into direction cos.
!        valid only for samll vec, which is alw.r(2)s assumed.
      implicit none

#include  "Zcoord.h"
      real*8 vx, vy
      type(coord)::dircos
!
      real*8 sum
!
      if(vx .eq. 0 .and. vy .eq. 0.) then
         dircos%r(1) = 0. 
         dircos%r(2) = 0.
         dircos%r(3) = 1.
      else
         if(abs(vx) .lt. 1.) then
            dircos%r(1) = vx
         else
            dircos%r(1) = 1.
         endif
         if(abs(vy) .lt. 1.) then
            dircos%r(2) = vy
         else
            dircos%r(2) = 1.
         endif
         dircos%r(3) = max(0.d0, 1.d0 - (vx**2 + vy**2)/2)
!               make consistent as direction cos
         sum = sqrt(dircos%r(1)**2 + dircos%r(2)**2 + dircos%r(3)**2)
         dircos%r(1) = dircos%r(1)/sum
         dircos%r(2) = dircos%r(2)/sum
         dircos%r(3) = dircos%r(3)/sum
      endif