det2Exyz.f

Go to the documentation of this file.

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

#include  "BlockData/cblkGene.h"
      implicit none

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




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


!      open(10, file="param")

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

      call cvecprod(detzaxis, detxaxis, detyaxis)
!      write(*,*) ' transformation matrix '
!              (x,y,z,1) = (X,Y,Z, 1) * Matrix
      write(* ,
     * '("# place: long=",f8.2, " lat=",f7.2, " height=",f10.1)')
     *  longitofsite,  latitofsite,  obssites(noofsites)%pos%height
      write(*,'("#  E-xyz values are",3G17.9)') det%x, det%y, det%z
      write(*,'(3g17.9, a)')  detxaxis%x, detxaxis%y, detxaxis%z, "  0"
      write(*,'(3g17.9, a)')  detyaxis%x, detyaxis%y, detyaxis%z, "  0"
      write(*,'(3g17.9, a)')  detzaxis%x, detzaxis%y, detzaxis%z, "  0"
      write(*,'(3g17.9, a)')  det%x, det%y, det%z, "  1"

!               this is for colum vector
!      write(*,'(3g17.9, a)')  DetXaxis.x, DetYaxis.x, DetZaxis.x, "  0"
!      write(*,'(3g17.9, a)')  DetXaxis.y, DetYaxis.y, DetZaxis.y, "  0"
!      write(*,'(3g17.9, a)')  DetXaxis.z, DetYaxis.z, DetZaxis.z, "  0"
!      write(*,'(3g17.9, a)')  det.x, det.y, det.z, "  1"


!      do while(.true.)
!         write(*,*) 'enter x,y,z in det sys'
!         read(*,*)  a.x, a.y, a.z
!         call cdet2xyz(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