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

Go to the source code of this file.

Functions/Subroutines

program __integratededx.f__
 
subroutine integrate (eventno, x, y, npa)
 

Function/Subroutine Documentation

◆ __integratededx.f__()

program __integratededx.f__ ( )

Definition at line 5 of file integratedEdx.f.

References integrate(), nc, parameter(), true, x, and y.

5  integer np, nc, i
i
nodes i
Definition: cstdatmos0.fNew2.h:93
nc
real cut integer nc
Definition: Zprivate.h:1
Here is the call graph for this function:

◆ integrate()

subroutine integrate ( integer  eventno,
real*8, dimension(np)  x,
real*8, dimension(np)  y,
integer  npa 
)

Definition at line 46 of file integratedEdx.f.

References kcsplcoef(), kcsplinteg(), ktrpzintt2(), and parameter().

46  implicit none
47  integer npa, eventno
48  integer np, nc
49  parameter(np=100, nc=100)
50  real*8 x(np), y(np)
51  real*8 ans, coef(nc, 3)
52  real*8 a, b, s
53  integer i
54 
55 
56  a = x(1)
57  b = x(npa)
58  call ktrpzintt2(y, 1, npa, x, 1, a, b, ans)
59  call kcsplcoef(x, y, npa, coef, nc)
60  call kcsplinteg(x, y, npa, coef, nc, a, b, s)
61  write(*,'(i5, 1p 2g14.3, i4)') eventno, ans, s, npa
kcsplinteg
subroutine kcsplinteg(x, y, n, coef, nc, a, b, s)
Definition: kcsplInteg.f:2
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
kcsplcoef
subroutine kcsplcoef(x, y, n, coef, nc)
Definition: kcsplCoef.f:2
ans
atmos%rho(atmos%nodes) **exp(-(z-atmos%z(atmos%nodes))/Hinf) elseif(z .lt. atmos%z(1)) then ans=atmos%rho(1) **exp((atmos%z(1) -z)/atmos%H(1)) else call kdwhereis(z, atmos%nodes, atmos%z, 1, i) a=atmos%a(i) if(a .ne. 0.d0) then ans=atmos%rho(i) **(1+a *(z-atmos%z(i))/atmos%H(i)) **(-1.0d0-1.d0/a) else ans=*atmos%rho(i) *exp(-(z-atmos%z(i))/atmos%H(i)) endif endif ! zsave=z ! endif cvh2den=ans end ! ---------------------------------- real *8 function cvh2temp(z) implicit none ! vettical height to temperatur(Kelvin) real *8 z ! input. vertical height in m ! output is temperature of the atmospher in Kelvin real *8 ans integer i if(z .gt. atmos%z(atmos%nodes)) then ans=atmos%T(atmos%nodes) elseif(z .lt. atmos%z(1)) then ans=atmos%T(1)+atmos%b(1) *(z - atmos%z(1)) else call kdwhereis(z, atmos%nodes, atmos%z, 1, i) ans=atmos%T(i)+atmos%b(i) *(z-atmos%z(i)) endif cvh2temp=ans end !--------------------------------------------- real *8 function cthick2h(t) implicit none real *8 t ! input. air thickness in kg/m^2 real *8 logt, ans integer i real *8 dod0, fd, a logt=log(t) if(t .ge. atmos%cumd(1)) then ans=atmos%z(1) - *(logt - atmos%logcumd(1)) *atmos%H(1) elseif(t .le. atmos%cumd(atmos%nodes)) then ans=atmos%z(atmos%nodes) - *Hinf *log(t/atmos%cumd(atmos%nodes)) else call kdwhereis(t, atmos%nodes, atmos%cumd, 1, i) ! i is such that X(i) > x >=x(i+1) ans
Definition: cstdatmos0.fNew2.h:82
i
nodes i
Definition: cstdatmos0.fNew2.h:93
ktrpzintt2
subroutine ktrpzintt2(t, intv, n, xt, intvx, a, b, ans)
Definition: ktrpzIntT2.f:5
y
latitude latitude this system is used *****************************************************************! type coord sequence union map real y
Definition: Zcoord.h:25
nrl_atmos::a
real(4), save a
Definition: cNRLAtmos.f:20
nrl_atmos::b
real(4), save b
Definition: cNRLAtmos.f:21
nc
real cut integer nc
Definition: Zprivate.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
Here is the call graph for this function: