1 !
for length
to thickness conversion or v.v
5 ! This
is used when making table
for 1-dim simulation.
6 ! The slant length
for vertical
height 0
to 100 km
is 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.
13 ! maxnodes: maximum no. of nodes
to express the atmosphere.
18 real*8 Htop, Hbase, ThickTbl, LenTbl, CosTbl, HeightTbl,
22 * HeightTbl(
maxl2t), Htop, Hbase, Zsave, ThickSave, NumStep
33 real*8 logrho(maxnodes)
34 real*8 logrhoi(maxnodes)
39 real*8 logcumd(maxnodes)
40 real*8 logcumdi(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)
56 common /Zatmosc2/ atmos
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
real(8), parameter, public m
********************block data cblkHeavy ********************integer j data *HeavyG2symbol H
! 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
block data cblkIncident data *Za1ry *HeightOfInj d3
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
! for length to thickness conversion or v v ! integer maxnodes real * LenStep
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos height
! for length to thickness conversion or v v ! integer maxl2t
! 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
block data cblkIncident data *Za1ry is