xyz2det.f

Go to the documentation of this file.

!      convert a given (x,y,z) in the Exyz system into
!     the detector  system.
!
!   input.  latitude, longitude, height(depth)  from param.
!           (name 'param' is fixed)
!   (x,y,z) is from stdin.
!   output.  (X,Y,Z) in det system  on stdout.
!

#include  "BlockData/cblkGene.h"
      implicit none

#include  "Zglobalc.h"
#include  "Ztrack.h"
#include  "Ztrackp.h"
#include  "Ztrackv.h"
#include  "Zobs.h"
#include  "Zobsv.h"




      type(coord):: a, b, det, detyaxis


      open(10, file="param")

      call creadparam(10)
      call cbeginrun
      det = obssites(noofsites)%pos%xyz

      call cvecprod(detzaxis, detxaxis, detyaxis)
      write(*,*) ' transformation matrix '
      write(*,'(4g17.9)')  detxaxis%x, detyaxis%x, detzaxis%x, 0
      write(*,'(4g17.9)')  detxaxis%y, detyaxis%y, detzaxis%y, 0
      write(*,'(4g17.9)')  detxaxis%z, detyaxis%z, detzaxis%z, 0
      write(*,'(4g17.9)')  det%x, det%y, det%z, 1 

      do while(.true.)
         write(*,*) 'enter x,y,z in Exyz sys'
         read(*,*)  a%x, a%y, a%z
         call cxyz2det(det, a, b)
         write(*,*) ' input=', a%x, a%y, a%z
         write(*,*) ' output=', b%x, b%y, b%z
      enddo
      end
      subroutine chooktrace
      end
      subroutine chookceren
      end
      subroutine chookcerens
      end
      subroutine chookcerene
      end
      subroutine chookbgrun
      end
!
       temp%r(1) = dir1%r(1) * xax%r(1) + dir1%r(2) *xax%r(2) +
     *  dir1%r(3)  *xax%r(3)
       temp%r(2) = dir1%r(1) * yvec%r(1)+ dir1%r(2) *yvec%r(2) + 
     *  dir1%r(3) *yvec%r(3)
       temp%r(3) = dir1%r(1) * zax%r(1) + dir1%r(2) *zax%r(2) +
     *    dir1%r(3)  *zax%r(3)
       dir2 = temp