COSMOS v7.655  COSMOSv7655
(AirShowerMC)
test.f
Go to the documentation of this file.
1 ! generate isotropic point on a sphere and extract
2 ! vertical ones. What dose it look like ?
3 ! should be uniform. in the circle
4  implicit none
5  real*8 cs, sn, sint, cost, x, y, z, cosx, u
6 
7  integer i
8 
9  do i = 1, 1000000
10 
11  call kcossn(cs, sn)
12  call rndc(cost)
13  sint = sqrt(1.d0- cost**2)
14 
15  x = sint* cs
16  y = sint* sn
17  z = cost
18 
19  call rndc(cosx)
20  if(cosx .gt. 0.99d0) then
21  call rndc(u)
22  if(cost .gt. u) then
23  write(*,*) sngl(x), sngl(y)
24  endif
25  endif
26  enddo
27  end
z
nodes z
Definition: cstdatmos0.fNew2.h:93
i
nodes i
Definition: cstdatmos0.fNew2.h:93
rndc
subroutine rndc(u)
Definition: rnd.f:91
d0
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
y
latitude latitude this system is used *****************************************************************! type coord sequence union map real y
Definition: Zcoord.h:25
kcossn
subroutine kcossn(cs, sn)
Definition: kcossn.f:13
x
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real * x
Definition: Zptcl.h:21