COSMOS v7.655  COSMOSv7655
(AirShowerMC)
Zatmos.h
Go to the documentation of this file.
1 ! for length to thickness conversion or v.v
2 !
3  integer maxl2t, maxnodes
4  real*8 LenStep, Hinf
5 ! This is used when making table for 1-dim simulation.
6 ! The slant length for vertical height 0 to 100 km is
7 ! divided by LenStep m steps.
8 ! It can cover the slant length of about 1200 km for cos=0.
9 ! When using the table,
10 ! Non-integral point (i.e., values not in the table) is
11 ! obtained by using interpolation using 3 integral points.
12 ! (see cl2tTblA.f)
13 ! maxnodes: maximum no. of nodes to express the atmosphere.
14  parameter (LenStep=1000.d0, maxl2t=1200.0d3/LenStep,
15  * maxnodes = 40)
16  parameter (Hinf = 1000.d3)
17  integer NumStep
18  real*8 Htop, Hbase, ThickTbl, LenTbl, CosTbl, HeightTbl,
19  * Zsave, ThickSave
20 
21  common /Zatmosc/ ThickTbl(maxl2t), LenTbl(maxl2t), CosTbl(maxl2t),
22  * HeightTbl(maxl2t), Htop, Hbase, Zsave, ThickSave, NumStep
23 
24 
25  type atmosph
26  sequence
27  real*8 z(maxnodes)
28  real*8 zi(maxnodes)
29  real*8 T(maxnodes)
30  real*8 P(maxnodes)
31  real*8 logP(maxnodes)
32  real*8 rho(maxnodes)
33  real*8 logrho(maxnodes)
34  real*8 logrhoi(maxnodes)
35  real*8 H(maxnodes)
36  real*8 a(maxnodes)
37  real*8 b(maxnodes)
38  real*8 cumd(maxnodes)
39  real*8 logcumd(maxnodes)
40  real*8 logcumdi(maxnodes)
41  real*8 d0(maxnodes)
42  real*8 coefh2r(maxnodes,3)
43  real*8 coefh2T(maxnodes,3)
44  real*8 coefh2d(maxnodes,3)
45 ! real*8 coefh2P(maxnodes,3)
46  real*8 coefh2H(maxnodes,3)
47 ! real*8 coefr2h(maxnodes,3)
48  real*8 coefd2r(maxnodes,3)
49  real*8 coefd2h(maxnodes,3)
50 ! real*8 coefP2h(maxnodes,3)
51  integer nodes
52  end type atmosph
53 
54  type(atmosph):: atmos
55 
56  common /Zatmosc2/ atmos
57 
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
real(8), parameter, public m
Definition: cpdgXs.f:13
nodes z
float real
Definition: Zdef.h:2
int integer
Definition: Zdef.h:3
********************block data cblkHeavy ********************integer j data *HeavyG2symbol H
Definition: cblkHeavy.h:7
nodes i
! for length to thickness conversion or v v ! integer maxnodes real Hinf ! This is used when making table for dim simulation ! The slant length for vertical height to km is ! divided by LenStep m steps ! It can cover the slant length of about km for cos
Definition: Zatmos.h:8
block data cblkIncident data *Za1ry *HeightOfInj d3
Definition: cblkIncident.h:7
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
Definition: Zfit.h:15
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
! for length to thickness conversion or v v ! integer maxnodes real * LenStep
Definition: Zatmos.h:3
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos height
Definition: ZavoidUnionMap.h:1
! for length to thickness conversion or v v ! integer maxl2t
Definition: Zatmos.h:3
! type magfield sequence ! Note that position vector where the magnetic field is given ! is not included here ! unit of field strength is in T(1 gauss=10 **-4 T) union map real *8 x
dE dx *! Nuc Int sampling table b
Definition: cblkMuInt.h:130
real(4), save a
Definition: cNRLAtmos.f:20
block data cblkIncident data *Za1ry is
Definition: cblkIncident.h:5