ceCent2sph.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cecent2sph (a, bb)

Function/Subroutine Documentation

cecent2sph()

subroutine cecent2sph	(	type(coord)	a,
		type(coord)	bb
	)

Definition at line 2 of file ceCent2sph.f.

References d0.

Referenced by chookenevent(), and ctranscoord2().

!            rectangular coord to spherical coord
      implicit none

#include  "Zglobalc.h"

#include  "Zcoord.h"
      type(coord)::a
      type(coord)::bb
!
      type(coord)::b
      real*8 sint, cost, sinphi, cosphi
!
!      b.radius = sqrt(a.r(1)**2 + a.r(2)**2 + a.r(3)**2)
      b%r(3) = sqrt(a%r(1)**2 + a%r(2)**2 + a%r(3)**2)
!      if(b.radius .gt. 0.) then
      if(b%r(3) .gt. 0.) then
!           sint =sqrt( (a.r(1)/b.radius)**2 + (a.r(2)/b.radius)**2)
           sint =sqrt( (a%r(1)/b%r(3))**2 + (a%r(2)/b%r(3))**2)
!           cost = a.r(3)/b.radius
           cost = a%r(3)/b%r(3)
!           b.theta = atan2(sint, cost)*Todeg
           b%r(1) = atan2(sint, cost)*todeg
           if(sint .ne. 0.d0) then
!              cosphi = a.r(1) / b.radius*sint
!              sinphi = a.r(2) / b.radius *sint
              cosphi = a%r(1) / b%r(3)*sint
              sinphi = a%r(2) / b%r(3) *sint
!              b.phi = atan2(sinphi, cosphi)*Todeg
              b%r(2) = atan2(sinphi, cosphi)*todeg
           else
!              b.phi = 0.d0
              b%r(2) = 0.d0
           endif   
       else
!           b.theta = 0.d0
!           b.phi = 0.d0
           b%r(1) = 0.d0
           b%r(2) = 0.d0
       endif
       b%sys='sph'
       bb = b

Here is the caller graph for this function: