csPrimAng.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	csprimang (dir)
subroutine	csprimisoang (dir)
subroutine	csprimpsang (dir)
subroutine	csprimapsang (dir)
subroutine	cinisprimang

Function/Subroutine Documentation

cinisprimang()

subroutine cinisprimang

(

)

Definition at line 150 of file csPrimAng.f.

References cerrormsg(), cos, coszenith, d0, kdhtoh(), kdzth2(), and pi.

Referenced by cbeginrun().

      implicit none
#include  "Zglobalc.h"
#include  "Ztrack.h"
#include  "Zincidentp.h"
#include  "Zincidentv.h"
      
      real*8  smx, smn, w3min, hmaxd, hmind, hmax, hmin, 
     *        wh1, wh2, wh3, w3max
      integer icon, jcon
      external cblkincident
      character*200  msg
!
!               used if arround point source
           smx=pi/2 -(sourcedec+ddelta)*torad
           smn=pi/2 -(sourcedec-ddelta)*torad
           cspsmx=cos(smx)
           cspsmn=cos(smn)
           if(za1ry .eq. 'ps' .or. za1ry .eq. 'aps') then
               write(msg,*) '1ry from point source is specified'
               call cerrormsg(msg, 1)
               w3min=real(coszenith)
!                  from declination and zenith angle, get time
!                  from meridian
               call kdzth2(sourcedec, w3min, hmaxd, icon)
               hmax=hmaxd
               if(icon .le. 1) then
                   w3max = imag_p(coszenith)
                   call kdzth2(sourcedec, w3max, hmind, jcon)
                   hmin=hmind
                   if(jcon .eq. 2) then
                       hmin=0.
                       call kdhtoh(sourcedec, 0.d0, wh1, wh2, wh3)
                       w3max=wh3
                   endif
                   obsvhour=cmplx(hmin, hmax, 8)
                   if(hmin .gt. hmax) then
                       write(msg,*) ' invalid CosZenith=',coszenith
                       call cerrormsg(msg, 0)
                   endif
                   if(icon .eq. 1) then
                       write(msg,*) ' max observable zenith is',
     *                 ' smaller than the input'
                       call cerrormsg(msg, 1)
                   endif
                   write(msg,*)' range of  hour =',hmin,' to ',
     *             hmax,'  which corresponds to cos(zenith)=',
     *             w3min, w3max
                   call cerrormsg(msg, 1)
               elseif(icon .eq. 2) then
!                       no observation range
                   write(msg,*) ' range of cos(zenith)=', coszenith,
     *             ' but no such angle is realizable'
                   call cerrormsg(msg, 0)
               endif
           endif

Here is the call graph for this function:

Here is the caller graph for this function:

csprimang()

subroutine csprimang

(

type(coord)

dir

)

Definition at line 14 of file csPrimAng.f.

References cerrormsg(), cos, csprimapsang(), csprimisoang(), csprimpsang(), d0, kcossn(), ksamplin(), rndc(), and spherical.

Referenced by __testsampprim.f__(), and cmkincident().

      implicit none
!          
#include  "Zglobalc.h"
#include  "Zcoord.h"
#include  "Zincidentp.h"
#include  "Zobs.h"
#include  "Zobsp.h"

      type(coord)::dir
      character*50 msg
      real*8 u,  cs, sn, sint, fai, za

      if(za1ry .eq. 'is') then
         call csprimisoang(dir)
      elseif(za1ry .eq. 'ps') then
!                  from point source
         call csprimpsang(dir)
      elseif(za1ry .eq. 'aps') then
!                around point source
         call csprimapsang(dir)
      elseif(za1ry .eq. 'cos' ) then
!          cos dcos
         call ksamplin
     *     (1.0d0, 0.d0,  real(CosZenith), imag(coszenith), za)
         dir%r(3) =  - za   ! going down is negative
         if(obsplane .eq.  spherical) then
!              actually setting here is dummy.
            call kcossn(cs, sn)
         else
            call rndc(u)
            fai = (imag_p(azimuth)- real(azimuth)) *u +
     *           real(azimuth)
            cs = cos(torad*fai)
            sn = sin(torad*fai)
         endif
         sint = sqrt(1.d0-dir%r(3)**2)
         dir%r(1) = - sint*cs   ! - is needed for going down ptcl.
         dir%r(2) = - sint*sn
         dir%sys = 'det'     !   fof ObsPlane=spherical, this is  reset later
      else
         write(msg,*) 'strange Za1ry=',za1ry
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

csprimapsang()

subroutine csprimapsang

(

type(coord)

dir

)

Definition at line 119 of file csPrimAng.f.

References kdhtoh(), khtad(), and rndc().

Referenced by csprimang().

      implicit none
#include  "Zglobalc.h"
#include  "Ztrack.h"
#include  "Zincidentp.h"
#include  "Zincidentv.h"
      type(coord)::dir

      real*8  u, w3p, decx, h, w1, w2, w3
!            sample angle in horizontal system
          call rndc(u)
          w3p=(cspsmx -cspsmn)*u + cspsmn
          decx=acos(w3p)*todeg
!c          if(decx .gt. 90.) then ! bug found Oct.15. by Kawata.
             decx=90.-decx
!c          endif
          call rndc(u)
          h=( imag_p(obsvhour)-real(obsvhour))*u + real(obsvhour)
          call rndc(u)
          if(u .lt. .5) then
                h=-h
          endif
          call kdhtoh(decx, h, w1, w2, w3)
          w3 = - w3   ! note vector direction is opposit.
          w2 = - w2
          w1 = - w1
          call khtad(w1, w2, w3, dir%r(1), dir%r(2), dir%r(3))
          dir%sys = 'det'

Here is the call graph for this function:

Here is the caller graph for this function:

csprimisoang()

subroutine csprimisoang

(

type(coord)

dir

)

Definition at line 68 of file csPrimAng.f.

References cos, coszenith, d0, and rndc().

Referenced by csprimang().

      implicit none
!
#include  "Zglobalc.h"
#include  "Zcoord.h"
#include  "Zincidentp.h"
      type(coord)::dir
!
      real*8 u, cs, sn, sint, fai
          call rndc(u)
          dir%r(3) =  -(  (imag_p(coszenith)- real(CosZenith) )*u +
     *             real(CosZenith) )    ! going down is negative 
          call rndc(u)
          fai = (imag_p(azimuth)- real(azimuth)) *u + 
     *           real(azimuth)
          cs = cos(torad*fai)
          sn = sin(torad*fai)
!ccc          call kcossn(cs, sn)
          sint = sqrt(1.d0-dir%r(3)**2)
          dir%r(1) = - sint*cs  ! - is needed for going down ptcl.
          dir%r(2) = - sint*sn
          dir%sys = 'det'

Here is the call graph for this function:

Here is the caller graph for this function:

csprimpsang()

subroutine csprimpsang

(

type(coord)

dir

)

Definition at line 93 of file csPrimAng.f.

References kdhtoh(), khtad(), and rndc().

Referenced by csprimang().

      implicit none

#include  "Ztrack.h"
#include  "Zincidentp.h"
#include  "Zincidentv.h"
      type(coord)::dir
!
      real*8 u, h, w1, w2, w3
          call rndc(u)
          h=(imag_p(obsvhour)-real(obsvhour))*u + real(obsvhour)
          call rndc(u)
          if(u .lt. .5) then
                h=-h
          endif
!             for source SourceDec at hour h, get horizontal vector
          call kdhtoh(sourcedec, h, w1, w2, w3)
          w3 = - w3     ! note  vector def. is opposit
          w2 = - w2
          w1 = - w1
!             convert to detector system
          call khtad(w1, w2, w3, dir%r(1), dir%r(2), dir%r(3))
          dir%sys = 'det'

Here is the call graph for this function:

Here is the caller graph for this function: