COSMOS v7.655  COSMOSv7655
(AirShowerMC)
Functions/Subroutines
psran.f File Reference

Go to the source code of this file.

Functions/Subroutines

real *8 function psran (seed)
 
double precision function gamfun (Y)
 

Function/Subroutine Documentation

◆ gamfun()

double precision function gamfun (   Y)

Definition at line 12 of file psran.f.

References a, d, d0, h, i, n, o, pi, r, t, x, y, and z.

12 ! Gamma function : See Abramowitz, page 257, form. 6.4.40
13 !-----------------------------------------------------------------------
14  IMPLICIT DOUBLE PRECISION(a-h,o-z)
15  double precision
16  + y,r,s,t,afspl,x,
17  + coef(10),pi,zerod,halfd,oned,twod,tend
18  SAVE
19 !
20  DATA coef/8.3333333333333334d-02,-2.7777777777777778d-03,
21  . 7.9365079365079365d-04,-5.9523809523809524d-04,
22  . 8.4175084175084175d-04,-1.9175269175269175d-03,
23  . 6.4102564102564103d-03,-2.9550653594771242d-02,
24  . 0.1796443723688306 ,-0.6962161084529506 /
25  DATA pi/ 3.141592653589793d0/
26  DATA zerod/0.d0/,halfd/0.5d0/,oned/1.d0/,twod/2.d0/,tend/10.d0/
27 !
28  x=y
29  afspl=oned
30  n=int(tend-y)
31  DO 10 i=0,n
32  afspl=afspl*x
33  x=x+oned
34 10 CONTINUE
35  r=(x-halfd)* log(x)-x+halfd* log(twod*pi)
36  s=x
37  t=zerod
38  DO 20 i=1,10
39  t=t+coef(i)/s
40  s=s*x**2
41 20 CONTINUE
42  gamfun = exp(r+t)/afspl
z
nodes z
Definition: cstdatmos0.fNew2.h:93
r
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz r
Definition: ZavoidUnionMap.h:1
i
nodes i
Definition: cstdatmos0.fNew2.h:93
nrl_atmos::h
real(4), dimension(:), allocatable, save h
Definition: cNRLAtmos.f:28
d0
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
o
struct ob o[NpMax]
Definition: Zprivate.h:34
pi
! constants thru Cosmos real * pi
Definition: Zglobalc.h:2
gamfun
double precision function gamfun(Y)
Definition: psran.f:12
y
latitude latitude this system is used *****************************************************************! type coord sequence union map real y
Definition: Zcoord.h:25
d
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
nrl_atmos::a
real(4), save a
Definition: cNRLAtmos.f:20
t
nodes t
Definition: cstdatmos0.fNew2.h:84
n
integer n
Definition: Zcinippxc.h:1
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

◆ psran()

real*8 function psran ( real*8  seed)

Definition at line 3 of file psran.f.

References rndc().

3  implicit none
4  real*8 seed
5  real*8 u
6  call rndc(u)
7  psran = u
rndc
subroutine rndc(u)
Definition: rnd.f:91
psran
real *8 function psran(seed)
Definition: psran.f:3
Here is the call graph for this function: