cmymain.f

Go to the documentation of this file.

!
!     this is not the main program. for C++ interface.
!          this is called from C++ user interface.
    block data cblkelemag


!         parameters for Elemag process.
!   (->   ----------------------------------------------

    real*8  recoilkinemine  !2  Recoil Kinetic Min Energy above which the recoil (=knock-on process)
                                ! is treated. Below this energy, the effect is included as continuous
                                ! energy loss.  Used only if KnockOnRatio $>$ 1.
                                ! See also KnockOnRatio.
        real*8 knockonratio   !2  KnockOnRatio* KEminoObs is used instead of RecoilKineMinE if KnockOnRatio $<$1.
    real*8  x0            !2  Radiation length in kg/m$^2$ for air. Normally the user should not touch this.
    real*8  ecrit         !2  Critical energy in GeV. \newline
                              !   Employed only when calculating air shower size in the hybrid 
                              !  air shower generation.  The value would be dependent on the
                              !  experimental purpose.  The default  value, 81 MeV, is bit too
                              !  small in many applications (The air shower size is overestimated). 
                              !  Comparisons of sizes  by the hybrid method and by the  full Monte 
                              !  Carlo tell that \newline 
                              !  $N_e$ (full 3-D M.C) $ < N_e$ (hybrid AS with $E_c=81$  MeV ) $  < N_e$ (full 1-D M.C)
                              !  $ {\ \lower-1.2pt\vbox{\hbox{\rlap{$<$}\lower5pt\vbox{\hbox{$\sim$}}}}\ }
                              !  N_e$(hybrid AS with $E_c={76}$ MeV)  at around shower maximum.
                              !  Hybrid AS is always essentially 1-D. 

    logical knockon       !2  Obsolete. Don't use this. See RecoilKineMinE
                              !    and KnockonRatio.
    real*8  anihie        !2  If E(positron) $<$ AnihiE, annihilation is  considered. 
    real*8  es            !2  Modified scattering constant. 19.3d-3 GeV
    real*8  maxcomptone   !2  Above this energy, Compton scattering is neglected. 
    real*8  maxphotoe     !2  Above this energy, photoelectric effect is neglected.
    real*8  minphotoprode !1  Below this energy, no photo-prod of hadron.  See also PhotoProd.
    logical photoprod     !1  Switch. if .false., no photo prod. of hadron is considered at all. 
                              !   See also MinPhotoProdE.
        integer moliere       !2  0$\rightarrow$ use Gaussian approx always (with air density change and
                  !    energy loss effect)\newline
                              !   1$\rightarrow$ use Moli\`ere scattering for non-electrons (default)\newline
                              !   2$\rightarrow$ use Moli\`ere scattering for all charged particles.\newline
                              !    If negative, anglular-correlated displacement is made to be 0 since Moli\`ere
                              !    theory cannot give it. (if $>0$, we use Gaussian approximation for correlation). 


!   <-) ----------------------------------------------

    common /zelemagc/  recoilkinemine, knockonratio,
     *  anihie, maxcomptone,
     *  maxphotoe, minphotoprode,   es, x0, ecrit, knockon, 
     *  photoprod, moliere



    data 
     *  recoilkinemine /0.2d-3/,     !
     *  knockonratio  /1.0d0/,       ! 
     *  anihie /30.e-3/,            ! Eposi < 30 MeV, anihilation considered
     *  x0 /365.667/,                ! radiation length of air in kg/m2
     *  ecrit /81.e-3/,             ! critical energy of air in GeV
     *  maxcomptone /1./,       ! compton is considered below 75 MeV
     *  maxphotoe /1.e-3/,            ! above this, PhotoElectric effect neg.
     *  minphotoprode /150.e-3/,     ! below 150 MeV, no gp --> hadrons
     *  es /19.3d-3/,                 ! scattering const. Note, not 21 MeV.
     *  knockon /.true./,            ! knockon is considered. Obsolete
     *  photoprod /.false./,        ! gp--> hadrons not considered.
     *  moliere /2/                 ! Moliere scattering for non-electron 
    end
      block data cblkevhnp

!            hadronic collision parameters 
!   (->  ---------------------------------------------

        real*8  cepic0  !2  Obsolute
        real*8  cekaon  !2  Obsolute
        integer sucint  !2  The number of successive interactions inside A is affected by this parameter.\newline
                        !   If $0\rightarrow$ old formula (before uv3.500) is used, which give rather         
                        !   smaller number ($<Nsuc>$ in Air = 1.7 for 30 mb pp), \newline
                        !   if $1\rightarrow$ new one $<Nsuc>=2.2$ for 30 mb pp). \newline
                        !   Default is 0 (from V5.00 again).
        real*8  ceneuc  !2 \verb| p -> n ; n-> p; p~->n~; n~->p~| prob.
        real*8  mudirp  !2 \verb| DD~| \# enhancement factor. D is only for prompt muon.
        real*8  kpilog  !2  K$_{ch}/\pi_{ch}$=(Kpilog*log(ss+.069)+Kpicns)*exp(-8/s') where ss(GeV**2)= 
                        !   effective s. s'(GeV**2)=s - 4.63.  See also Kpicns.
        real*8  kpicns  !2  See Kpilog. 0.077
        real*8  elund   !2   obsolete (from v6.0)
!                  old def
!        real*8  Elund   !2   E $<$ Elund$\Rightarrow$  Lund M.C is used. in h-n. (IntModel='int1'). 
!                        !    If IntModel='int2' $\Rightarrow$    Gheisha is used. Elund can be as
!                        !    small as 4.99. (Default is 500. For NEXT 4.99).
        real*8  elund2  !2  obsolete (from v6.0)
!              old def
!        real*8  Elund2  !2   Elund $<$ E $<$  Elund2 $\Rightarrow$   ad hoc model is used; Elund2 must
!                        !    be {\ \lower-1.2pt\vbox{\hbox{\rlap{$>$}\lower5pt\vbox{\hbox{$\sim$}}}}c                        !     }100, if Elund2 $<$     Elund3.   Defualt is Elund
         real*8 elund3   !2 obsolete  (from v6.0).
!                old def
!        real*8  Elund3  !2  Elund2$<$ E $<$  Elund3 $\Rightarrow$  New Lund Fritiof is used.
!                        !   Elund3 can be $<$     Elund2. If E $>$  Elund3, ad hoc model is used.
!                        !   Default is Elund.
    real*8  efermi  !2  If Kinetic E $<$ Efermi, Fermi Momentum is considered for Nucleus target.\\

        character*64  intmodel !1 Interaction model description. Usage was changed from v6.0. 
                      !    One may list  code name and upper energy bound  for the code.\newline
                      ! E.g.  IntModel = '"dpmjet3"' ; to specify the dpmjet3 in the entire energy region
                      !       IntModel = '"dpmjet3" 100 "nexus2" to specify dpmjet3 at $<$ 100 GeV and nexus2 
                      !        at E$>$100 GeV \newline
                      !       IntModel = '"nucrin" 5 "fritiof1.6" 500 "adhoc" to specify  Nucrin,
                      !               fritiof1.6, and ad-hoc model in the respective energy region. This 
                      !       corresponds to the old IntModel='int1'. \newline
                      !       IntModel = '"nucrin" 5 "fritiof1.6" 10 "fritiof7.02"  and \newline
                      !       IntModel = '"dpmjet3"' \newline
                      !        are most promissing models that fit the observed 
                      !       data (muons and gamma rays) for which the primary is well known by
                      !       BESS and AMS observations ($<$ 100 GeV).

!              next is old def.
!        character*16  IntModel   !1 interaction model.  If 'int1', Lund Fritiof/Nucrin/adhoc models are used
!                       ! for hadronic interactions.  If 'int2' Gheisha code is used for hadronic interactions,
!                       ! at low energies. (See, Elund?).
!                       !  If 'inclusive', an inclusive treatment is employed as given in Contrib/Inclusive.
!                       !  To make this work actually, INCLUSIVE in cosmos/Zcondc.h must be given a value of 1
!                       !  when making the cosmos library.  At execution time, InclusiveFile must be given
!                       !  the path to the xdist.d file. (normally in Contrib/Inclusive/xdist.d)
!                       !  The treatmen is one dimesional. May be used for neutrino calculations.
        character*100 inclusivefile !2 The path to the inclusive data file, xdist.d. Default is
                       !   "../Contrib/Inclusive/xdist.d"
        real*8 sucpw   !2  In the 2nd, 3rd,.. collision of a nucleon inside a nucleus, the collision is 
                       !  made to be more elastic than normal one. The leading particle spectrum is
                       ! sampled from x**SucPw dx. SucPw should be in 1 to 2.
        real*8  eta2pi0 !2  eta/pi0 ratio. this is used to see the effect due to non-decay of pi0
                      ! at very high energies. Only source of h.e gamma can be eta and LPM may work
                      ! for them. default is 0.2
        integer mullow !2 if 1, QCD predicted multiplicity law is used in the adhoc model else UA5
                       !  parametalization is used. Default is 1. (from v5), 
                       !  0.6135exp(23/18sqrt(2log(roots/0.3))) is QCD jet prediction. 
                       !  7.2roots**0.254 -7 is UA5 data. The branch point is set at roots = 900 GeV. 
                       !  (I have adjusted 0.6135 so that 900 GeV is the b.p)
        integer lundpara !2 To control Lund program. LundPara(1) is set to kfr(7); kfr(7)=1 is for Frititof
                       !  hard scattering. 2 is for Pythia H.S. 2 gives higher multiplicity but shape is 
                       !  strange.  Default is 1. LundPara(2) is set to kfr(12): 1 by for OPAL hard scattering
                       !   parameterization. 2 by DELPHI. Default is 2. (2 gives bit higher PT). LundPara(3)
                       !   $>$ 0 $\Rightarrow$    Pythia message will appear. LundPara(4) $>$ 0 Fritiof
                       !  message; both on ErrorOut. LundPara(5) =0 $\Rightarrow$    All kaons collisions
                       !  are treated as pi- in Fritiof, else they are treated by adhoc model as they are.

!   <-)     -----------------------------
!            next are not input parameters.
        integer nmdls
        parameter(nmdls = 8)
        character*16 regmdls(nmdls)

        common /zevhnp/ cepic0, cekaon, ceneuc, mudirp, kpilog, eta2pi0,
     *         kpicns,  efermi, elund, elund2, elund3, sucpw, mullow, 
     *         lundpara(10), sucint
        common /zevhnc/regmdls, inclusivefile, intmodel 
!                 this should be in block common. 
!        data  Cepic0/0.3d0/, Cekaon/0.d0/, Ceneuc/0.35d0/,
!    *         Mudirp/1.d0/, Kpilog/0.0062/, Kpicns/0.077/,
!    *         Efermi/2.d0/

!            currently usable models
      data regmdls/'dpmjet3', 'nexus2', 'fritiof7.02', 'fritiof1.6',
     *             'gheisha', 'nucrin', 'ad-hoc', 'incdpm3'/
      
       data  
     * cepic0 /0.00d0/ , 
     * cekaon /0.d0/ ,
     * ceneuc /0.35d0/ ,
     * mudirp /1.d0/ , 
     * kpilog /0.0062/ ,
     * kpicns /0.077/ ,

     * elund /500./ , 
     * elund2 /500./,
     * elund3 /500./,

     * sucpw /1.5/,
     * mullow /1/,      ! old one is 0. bef. v.5.0
     * lundpara /1, 2, 0, 0, 1, 5*0/,
     * intmodel /'"dpmjet3" 1.d8'/ ,
     * inclusivefile /'../../Contrib/Inclusive/xdist.d'/ ,
     * efermi /10.e-3/
       data 
     * eta2pi0/0.2/,

     * sucint /1/     ! old one was 1 < v5.0; 0 <=v5.10   1 >= v5.11

      end
!     ********************
      block data cblkheavy
!     ********************

!             ptcl kind code; kindmx is the no. of observable ptcls
!             klast; max ptcl code in the system.
!
        integer  kphoton, kelec, kmuon, kpion,  kkaon, knuc,
     1  kneue, kneumu, kindmx, knnb, kddb, kdmes, krho,
     2  komega, kphi, keta, kgnuc, kalfa, klibe, kcno, khvy, kvhvy,
     3  kiron, khvymax, klast, klambda, ksigma, kgzai, kbomega,
     4  ktriton, klambdac, krare
!            subcode
        integer
     1  regptcl, antip, k0s,  k0l, kneutron, 
     4  kneutronb, kd0, kd0b, kdirectg, kcasg

!
        parameter(
     1  kphoton=1, kelec=2, kmuon=3, kpion=4,  kkaon=5,
     2  knuc=6,
     3  kneue=7, kneumu=8, kgnuc=9, kalfa=10, klibe=11, kcno=12, 
     4  khvy=13, kvhvy=14, kiron=15, kdmes=16, 
!          next line added Nov. 17,'95. 
     5  ktriton=17, klambda=18, ksigma=19, kgzai=20, klambdac=21,
     6  kbomega=22,  kindmx=kbomega, krare = 0,

     7  knnb=kindmx+1, kddb=knnb+1,  krho=kddb+1,
     8  komega=krho+1, kphi=komega+1, keta=kphi+1,
     9  klast=keta, khvymax = kiron)
        parameter(
     1  regptcl=-1, antip=1,
     2  kdirectg=2, kcasg=3,
     3  k0s = 4,  k0l= 5,
     4  kneutron=regptcl,
     5  kneutronb = antip, kd0 =-8, kd0b =-kd0)
!
    integer maxheavymassn, maxheavycharge, maxheavyg
    parameter(maxheavymassn = 56, maxheavycharge = 26, 
     1             maxheavyg = 7)
!       kphoton: gamma ray 
!        kelec: electron, positiron
!        kmuon: muon
!        kpion: pion
!        kkaon: kaon
!        knuc: neucleon
!        kneue: electron neutrino
!       kneumu: muon neutrino
!        kgnuc: general nucleus(A>=2.)
!        kalfa: alpha  (heliunm)
!        klibe: Li, Be, B
!         kcno: C, N, O 
!         khvy: heavy such as, Na/Mg/Si
!        kvhvy: very heavy such as S/Cl/Ar
!        kiron: iron group
!        regptcl: particle index
!        antip: anti-particle index
!        krare:  used to set very rare particle code
!                which might come from imported soft.
!                such as tau. They are neglected in
!                Cosmos.





!
!   (-> ------------------------------------------

    integer charge2heavyg   !2  charge of heavy $\rightarrow$  heavy group index conversion array.
        integer heavyg2massn    !2  heavy group index $\rightarrow$     mass number conversion array.
    integer heavyg2charge   !2  heavy group index $\rightarrow$     charge of heavy conversion array.
        integer heavyg2code     !2  heavy group index $\rightarrow$     particle code conversion array.
        integer code2massn      !2  particle code $\rightarrow$     mass number conversion array.
        integer code2heavyg !2  particle code $\rightarrow$     heavy group index conversion array.
        real*8  fragmenttbl !2  tbl(i,j)=$<$Number$>$  of frag. j when a heavy of heavy group index i
                                !    breaks up at air.
        real*8  ptavnonintenuc  !2  $<$Pt$>$  of non interacting nucleons.
    real*8  ptavfrag        !2  $<$Pt$>$  of heavy fragments.
    character*4 heavyg2symbol !2   heavy group index $\rightarrow$  'Fe' etc conversion array.
     integer howintnuc       !2 If 0, the  number of interacting nucleons among a projectile heavy nucleus is 
                                 !  determined as the number of first collision of each interacting nucleon inside 
                                ! the  nucleus.  If 1, the number is determined as the total number of collisions 
                                !   including successive interactions. Default is 1. (There is uncertaninity in
                                !  interpretation of the formula; value 1 gives larger number of interacting
                                !  nucleons.)


 
!   <-) --------------------------------------
        

    common /zheavyc/
     *   ptavnonintenuc, ptavfrag,
     *   fragmenttbl(maxheavyg, maxheavyg), 
     *   charge2heavyg(maxheavycharge),
     *   heavyg2massn(maxheavyg), heavyg2charge(maxheavyg),
     *   heavyg2code(maxheavyg), code2massn(khvymax),
     *   code2heavyg(khvymax), howintnuc
        common /zheavycc/ heavyg2symbol(maxheavyg)



    
      integer i, j

      data
     *  heavyg2massn /1,   4,    8,    14,    25,    35,     56/ ,
     *  heavyg2charge /1,   2,    4,     7,    12,    17,     26/ ,
     *  heavyg2symbol /'p', 'Alfa',  'L', 'CNO', 'H', 'VH', 'Fe'/ ,
     *  charge2heavyg /1,   2,   3,3,  4,4,4,  5*5,   5*6,   9*7/

      data 
     *  heavyg2code /knuc, kalfa, klibe, kcno, khvy, kvhvy, kiron/    
      data
     * code2heavyg(knuc)/1/ ,
     * code2heavyg(kalfa)/2/ ,
     * code2heavyg(klibe)/3/ ,
     * code2heavyg(kcno)/4/ , 
     * code2heavyg(khvy)/5/ ,
     * code2heavyg(kvhvy)/6/ ,
     * code2heavyg(kiron)/7/
      data
     * code2massn(knuc)/1/ ,
     * code2massn(kalfa)/4/ , 
     * code2massn(klibe)/8/ ,
     * code2massn(kcno) /14/,
     * code2massn(khvy)/25/,
     * code2massn(kvhvy)/35/,
     * code2massn(kiron)/56/

      data  howintnuc/1/
!
      data ( ( fragmenttbl(i,j), j=1,maxheavyg), i=1,maxheavyg)/
     1             1.,    0.,    0.,    0.,    0.,    0.,   0.,
     2             4.,    0.,    0.,    0.,    0.,    0.,   0.,
     3             5.81,  .61,   .11,   0.,    0.,    0.,   0.,
     4             6.62,  .72,   .24,   .17,   0.,    0.,   0.,
     5             10.5,  .77,   .21,   .39,   .16,   0.,   0.,
     6             12.8,  1.17,  .17,   .20,   .42,   .06,  0.,
     7             18.1,  1.71,  .24,   .17,   .22,   .20,  .17/

    data ptavnonintenuc/50.e-3/   ! 50 MeV
    data ptavfrag/50.e-3/     ! 50 MeV
      end
       block data cblkincident
!                parameters for primary angle sampling
!   (-> ---------------------------------------------------

           complex*16  coszenith   !1  Range of cos(zenith angle). Say, (0.5, 1.0). Used when Za2ry is 'is' 
                                   !   If ObsPlane=3 (spherical), real(CosZenith) must be $>$0, and means
                                   !   the zenith angle range at the incident point (not in Exyz system).
                                   !   In that case, azimuth is 0 to 2pi.
           complex*16  azimuth     !1  Range of azimuthal angle in deg. Say, (0, 45). Default is (0,360).
                                   !   Can be such as (300., 390.).  Used when Za1ry is 'is'\newline
                                   !   If ObsPlane=3 (spherical), this is used to show the half opening angle
                                   !   range where the primary injection position is uniformly distributed 
                                   !   on a sphere.  The center of the opening angle is (Latit, Longit, HeightOfInj).
                                   !   In this case, for the upper opening angle,  min( Imag(Azimuth),180.) is used.
           character*4 za1ry       !1  Specify the primary angle sampling method by one of 'is', 'ps' or 'aps'.\newline
                        !   "is" is isotropic. The range is by CosZenith.\newline
                        !   "ps" is for point source (See also SourceDec)\newline
                        !   "aps" is around point source (See also SourceDec and  Ddelta) \newline
                        !   If ObsPlane=3 (spherical), this must be "is".
           real*8  sourcedec    !1  Source declination of point source.(deg)
       real*8  ddelta       !1  SourceDec $\pm$ Ddelta is the region for 'aps' (deg).
           real*8  heightofinj  !1  The vertical height of primary injection point (m).
                                !   If this is $<$ deepest obs. level and zeinth angle of primary is $< 0$, 
                                !   the primary is  assumed to be upgoing even if Reverse =0.
                            !   NOTE: BorderHeightH must be given explicitly in this case.
       real*8  offsetheight !2  The vertical offset  height from  the deepest detector. 
                                ! The  primary is directed to this height above the detector.
                                ! If ObsPlane is  3 (spherical), not used.

!   <-) ----------------------------------------------------

           common /zincident/ azimuth, coszenith, sourcedec, ddelta,
     *      heightofinj,  offsetheight
           common /zincidentc/ za1ry
         data 
     *   azimuth /(0.d0, 360.d0)/,
     *   za1ry /'is'/ ,
     *   coszenith /(0.6d0, 1.d0)/ ,
     *   sourcedec /30./ ,
     *   ddelta /5./,
     *   heightofinj /100.d3/ ,
     *   offsetheight /0./
       end
      block data  cblkmanager


















!   Parameters   needed  for the Launcher.
!
!   (-> ------------------------------------

    integer errorout    !2 Error output logical  dev number.
    character*70  primaryfile  !1  Primary Spectrum data file (full or relative path)
    character*70  cutofffile   !1  Geomagnetic cut-off file
    character*70  contfile     !1  Job continuation information file  (full or relative path).
                                   !   default is "ContInfo".  This will be created when job
                                   !   is finished normally.
    character*70  geomagfile   !2  IGRF or WMM file path which contains geomagnetic field expansion 
                                   !   coefficients.  Their format is the same one as given in their web 
                                   !   page.  If ' ' (default), Cosmos/Data/Geomag/igrf is used.
    character*70  skeletonfile !1   Skeleton information file (full or relative path). created if Job =
                                   !    'skeleton'. Default is 'skeletonParam'.  This is the Namelist data
                                   !     referred by Cosmos automatically  if Job='flesh' is specified. For
                                   !     Job='flesh', you have to modify some part of  this file.
        character*70  dpmfile      !2  control card to specify the dpmjet execution conditions. If ' ',
                               !   Cosmos/Data/DPM/atmos.inp is assumed.  
    character*10  job          !1  What kind of job you are going to do.\newline
                                   !   =' ' (default).  nothing special.\newline
                                   !   ='skeleton'.  Makes skeleton. \newline
                                   !   ='flesh'. Flesh skeleton events.  See manual.\newline
                                   !   ='newskel'   \newline
                                   !   ='newflesh'  see manual. \newline
        character*70  seedfile     !1   File to  contain the initial random numbers for those events to 
                                   !    which you want to flesh. You can create the file by calling
                                   !     cwriteSeed in a user hook routine (say, in chookEnEvent) at 
                                   !     skeleton making time. Default is 'Seed'.  For a normal run with
                                   !      Job=' ', if SeedFile is not ' ',  two integer initial random numbers
                                   !      and the event number are  automatically output on the speicfied disk file.
        integer       seedfiledev  !2   logical device number of SeedFile.
    logical       cont         !1  If T, continuation from a previous job is assumed. Contfile content is used.
    integer       initrn       !1  Initial random number seed. 2 integers. If InitRN(1) $<$ 0, file dev  \# 14
                                   !    is  assumed to have  pairs of IR in each row, and they are read to
                                   !    initialize each event.  This feature is ignored when Job = 'flesh' or 
                                   !    'newflesh'. The \# 14 file should be opened by the user routine
                                   !    (chookBgRun). This is almost debug purpose.\newline
                                   !   If InitRn(2)$<$0, timer, hostname and process number are used for the 
                                   !    initialization.
    integer       eventno      !2  cumulative event number counter.(excluding discarded ones due to cutoff).
    integer       eventsintherun !2  Counter for event number in the run. Internal use.
                                     !          (excluding discarded ones due to cutoff).
    integer       desteventno    !1 2 integers: Final event no. to be generated and events to be generated
                                     !  in the current run.  If negative, their absolute is used and counting 
                                     !  includes discarded ones due to rigidity cutoff.
                                     !  If DestEventNo(2)=0, DestEventNo(1) is used. If it is negative, only
                                     !  DestEventNo(2) is checked to see events in the current run. For the
                                     !  flux calculation, negative ones are better.
    logical       hidden         !1  Make T, if hidden parameters are to be written.
    integer       tempdev      !2  Logical Dev. number for temporary disk use.
    integer       preveventno  !2  The event number already finished.  System use for Cont job.
                                   !        (excluding discarded ones due to cutoff).
    character*8   deadline     !1  The dead line before which the job should terminate.
                                   !   Should be given like '10.11.15' which means the nearest 10th, 11 O'clock,
                                   !   15 min.  Not used if Within has non zero value.  
        integer       within       !1  The job should end within this minutes from now.  Default is 99999.
                                   !   If 0 is given,  DeadLine is used.
        real*8        basetime     !1  Rough cpu time needed for completing one event (say, for protons, or
                                   !   gamma rays) with energy BaseErg.  The cpu time estimation is based on 
                                   !   A * ( E1ry par nucleon )**BasePower / BaseErg * BaseTime, where A is mass number
                                   !  (for nucleus; otherwise 1).
        real*8        baseerg      !2  See BaseTime.  The default is  1000 (GeV).
        real*8        basepower    !2  See BaseTime.   Default is 1.0
        character*100 userhookc    !2  array size is MAX\_USERHOOKC(=5). Usage is left for the user. To get the i-th
                                   !   component, the use may 'call  cqUHookc(i, cv)' in the userHook routine, 
                                   !   where cv is a character variable to receive the data.
        real*8        userhookr    !2  array size is MAX\_USERHOOKR(=10). Usage is left for the user. To get the i-th
                                   !   component, the use may 'call cqUHookr(i, rv)' in the userHook routine,
                                   !   where rv is a real*8 variable to receive the data.
        integer       userhooki    !2  array size is MAX\_USERHOOKI(=10). Usage is left for the user.  To get the i-th
                                   !   component, the use may 'call ccqUHooki(i, iv)' in the userHook rouitne,
                                   !   where iv is an integer varialbe to receive the data.
        character*100 atmosfile    !2  path to the atmospheric data as in 'Cosmos/Data/Atmos/stdatmos2.d'

        character*32  atenv        !2  If this is non blank, an environmental variable with that name is
                                   !   assumed to exist and Cosmos tries to get the value of that env variable.
                                   !   If the value is obtained, the \verb/@/ in \verb/@_/ or \verb/@./
                                   !   expressing a part of a file name is replaced by that value. 
                                   !   (default is blank and in that case the \verb/@/ is replaced by
                                   !    the host name where  the job runs.)

        character*32 sharpenv      !2  If this is non blank, an environmental variable with that name is
                                   !   assumed to exist and Cosmos tries to get the value of that env variable.
                                   !   If the value is obtained, the \verb/#/ in \verb/#_/ or \verb/#./ 
                                   !   expressing a  part of a file name is replaced by that value. 
                                   !   (default is blank and in that case the \verb/#/ is replaced by
                                   !    the process number of the run).

        character*32 percentenv    !2  If this is non blank, an environmental variable with that name is
                                   !   assumed to exist and Cosmos tries to get the value of that env variable.
                                   !   If the value is obtained, the \verb/%/ in \verb/%_/ or \verb/%./ 
                                   !   expressing a  part of a file name is replaced by that value. 
                                   !   (default is blank and in that case the \verb/%/ is replaced by
                                   !    the USER name).


!   <-) -------------------------------------
    common /zmanagerpc/
     *  basetime,  baseerg, basepower, within, userhookr(10),
     *  errorout, cont, initrn(2), userhooki(10),
     *  eventsintherun, desteventno(2), hidden, tempdev, 
     *  preveventno, seedfiledev, eventno


       common /zmanagerpc2/
     * userhookc(5), primaryfile,
     * cutofffile,  job, contfile, atmosfile, geomagfile,
     * skeletonfile, seedfile, dpmfile, deadline, sharpenv,
     * percentenv, atenv
 

       data 
     *  errorout /0/   ! sun4 ; others.  if 0 is ng, give by namelist

       data  
     * desteventno /1,0/ ,
     * hidden /.false./ ,
     * tempdev /11/ ,
     * primaryfile /' '/ ,
     * job /' '/ ,
     * contfile /'ContInfo'/ , 
     * skeletonfile /'SkeletonParam'/ ,
     * seedfile /'Seed'/ ,
     * preveventno /0/ ,
     * eventno /0/, 
     * seedfiledev /22/ ,
     * deadline /' '/ ,
     * within /99999/ , 
     * basetime /10./ ,
     * baseerg /1000./ ,
     * basepower /1./
       data
     * cutofffile /' '/ 
       data
     * geomagfile /' '/
       data
     * userhookc /5*' '/
       data
     * userhooki /10*99999999/
       data
     * userhookr /10*1.d65/
       data
     * atmosfile /' '/
       data
     * dpmfile /' '/
       data
     * percentenv /' '/
       data
     * sharpenv /' '/
       data
     * atenv /' '/
      end
      block data cblkmuint
      integer i
!         muon interaction related variables
!

!          
      real*8 zeff, zeff3  ! Zeff**(1/3)
      real*8 muprvmin, muprdetx, muprde, mupremin,
     *       mupremax, muprdu,  mupremax1
      integer    muprusize, mupresize,  muprtxt
!
      real*8 mubrvmin, mubrdetx, mubrde, mubremin, mubremax,
     *       mubrdu, mubremax1
      integer mubrusize, mubresize, mubrtxt
!
      real*8 munvmin, mundetx, munde, munemin,
     *  munemax, mundu, munemax1
      integer munusize, munesize, muntxt
!
      real*8 munpwtx, munpwdedx0, munpwdedxt

      real*8 muprtx(38), muprdedx0(38), muprdedxt(38)
      real*8 mubrtx(36), mubrdedx0(36), mubrdedxt(36)
      real*8 muntx(34), mundedx0(34), mundedxt(34)

      real*8 muntbl(101, 17)

      real*8  mupa, mura, mupb, muqb, muak, muakm, muakm2,
     *        mupointlike, mushadow, mulogf0


      real*8  munlemin, muprlemin, mubrlemin
      common /muintc/ muntbl,  muprtx, muprdedx0, muprdedxt,
     *    mubrtx, mubrdedx0, mubrdedxt, muntx, mundedx0, mundedxt,
     *    munpwtx, munpwdedx0, munpwdedxt, muprvmin, muprdetx, 
     *    muprde, mupremin, mupremax, muprdu,  mupremax1,
     *    mubrvmin, mubrdetx, mubrde, mubremin, mubremax,
     *    mubrdu, mubremax1, munvmin, mundetx, munde, munemin,
     *    munemax, mundu, munemax1,  munlemin, muprlemin, mubrlemin,
     *    mupa, mura, mupb, muqb, muak, muakm, muakm2,
     *    mupointlike, mushadow, mulogf0, zeff, zeff3,
!
     *    muprusize, mupresize,  muprtxt, mubrusize, mubresize,
     *    mubrtxt, munusize, munesize, muntxt

!          Pair total X-sec. 
          data ( muprtx(i), i=  1,   38)/
     1  0.431132e-02, 0.473175e-02, 0.514887e-02, 0.556010e-02,
     2  0.596286e-02, 0.635461e-02, 0.673285e-02, 0.709522e-02,
     3  0.743955e-02, 0.776396e-02, 0.806695e-02, 0.834742e-02,
     4  0.860479e-02, 0.883891e-02, 0.905012e-02, 0.923915e-02,
     5  0.940706e-02, 0.955517e-02, 0.968498e-02, 0.979806e-02,
     6  0.989601e-02, 0.998043e-02, 0.100528e-01, 0.101146e-01,
     7  0.101672e-01, 0.102116e-01, 0.102491e-01, 0.102805e-01,
     8  0.103068e-01, 0.103286e-01, 0.103467e-01, 0.103616e-01,
     9  0.103739e-01, 0.103839e-01, 0.103920e-01, 0.103987e-01,
     a  0.104040e-01, 0.104083e-01                                      
     * /   

!            Pair  dE/dx(v<vmin)/E /(g/cm2)
          data ( muprdedx0(i), i=  1,   38)/
     1  0.415696e-07, 0.508098e-07, 0.606769e-07, 0.710744e-07,
     2  0.819072e-07, 0.930817e-07, 0.104506e-06, 0.116091e-06,
     3  0.127747e-06, 0.139387e-06, 0.150927e-06, 0.162285e-06,
     4  0.173381e-06, 0.184145e-06, 0.194509e-06, 0.204417e-06,
     5  0.213820e-06, 0.222682e-06, 0.230976e-06, 0.238687e-06,
     6  0.245811e-06, 0.252353e-06, 0.258325e-06, 0.263746e-06,
     7  0.268643e-06, 0.273044e-06, 0.276982e-06, 0.280489e-06,
     8  0.283600e-06, 0.286349e-06, 0.288768e-06, 0.290890e-06,
     9  0.292744e-06, 0.294359e-06, 0.295761e-06, 0.296974e-06,
     a  0.298020e-06, 0.298920e-06
     * /
!
!        Pair; dE/dx/E (v<all)  /(g/cm2)
          data ( muprdedxt(i), i=  1,   38)/
     1  0.664623e-06, 0.725450e-06, 0.786305e-06, 0.846786e-06,
     2  0.906509e-06, 0.965081e-06, 0.102213e-05, 0.107730e-05,
     3  0.113026e-05, 0.118071e-05, 0.122843e-05, 0.127320e-05,
     4  0.131491e-05, 0.135348e-05, 0.138889e-05, 0.142119e-05,
     5  0.145044e-05, 0.147678e-05, 0.150034e-05, 0.152131e-05,
     6  0.153987e-05, 0.155620e-05, 0.157052e-05, 0.158301e-05,
     7  0.159386e-05, 0.160325e-05, 0.161134e-05, 0.161830e-05,
     8  0.162426e-05, 0.162936e-05, 0.163370e-05, 0.163739e-05,
     9  0.164053e-05, 0.164318e-05, 0.164542e-05, 0.164731e-05,
     a  0.164890e-05, 0.165023e-05
     * /
!
!
!       Brem; Total X-sec
          data ( mubrtx(i), i=  1,   36)/
     1  0.348421e-03, 0.353857e-03, 0.358933e-03, 0.363712e-03,
     2  0.368194e-03, 0.372385e-03, 0.376290e-03, 0.379917e-03,
     3  0.383272e-03, 0.386366e-03, 0.389207e-03, 0.391806e-03,
     4  0.394175e-03, 0.396325e-03, 0.398269e-03, 0.400021e-03,
     5  0.401591e-03, 0.402995e-03, 0.404244e-03, 0.405351e-03,
     6  0.406329e-03, 0.407189e-03, 0.407942e-03, 0.408600e-03,
     7  0.409173e-03, 0.409669e-03, 0.410098e-03, 0.410467e-03,
     8  0.410784e-03, 0.411055e-03, 0.411286e-03, 0.411483e-03,
     9  0.411650e-03, 0.411791e-03, 0.411910e-03, 0.412010e-03
     * /

!        Brem:  dE/dx/E (v<vmin)  /(g/cm2)
          data ( mubrdedx0(i), i=  1,   36)/
     1  0.178842e-08, 0.178961e-08, 0.179057e-08, 0.179133e-08,
     2  0.179194e-08, 0.179242e-08, 0.179281e-08, 0.179312e-08,
     3  0.179336e-08, 0.179356e-08, 0.179371e-08, 0.179384e-08,
     4  0.179393e-08, 0.179401e-08, 0.179407e-08, 0.179412e-08,
     5  0.179416e-08, 0.179419e-08, 0.179422e-08, 0.179424e-08,
     6  0.179425e-08, 0.179426e-08, 0.179428e-08, 0.179428e-08,
     7  0.179429e-08, 0.179429e-08, 0.179430e-08, 0.179430e-08,
     8  0.179430e-08, 0.179431e-08, 0.179431e-08, 0.179431e-08,
     9  0.179431e-08, 0.179431e-08, 0.179431e-08, 0.179431e-08
     * /

!       Brem: dE/dx/E(v<all) /(g/cm2)
          data ( mubrdedxt(i), i=  1,   36)/
     1  0.785146e-06, 0.816296e-06, 0.845259e-06, 0.873845e-06,
     2  0.901994e-06, 0.929644e-06, 0.956731e-06, 0.983189e-06,
     3  0.100895e-05, 0.103394e-05, 0.105811e-05, 0.108137e-05,
     4  0.110368e-05, 0.112497e-05, 0.114520e-05, 0.116433e-05,
     5  0.118231e-05, 0.119914e-05, 0.121480e-05, 0.122929e-05,
     6  0.124263e-05, 0.125484e-05, 0.126595e-05, 0.127601e-05,
     7  0.128506e-05, 0.129317e-05, 0.130038e-05, 0.130677e-05,
     8  0.131240e-05, 0.131734e-05, 0.132165e-05, 0.132540e-05,
     9  0.132864e-05, 0.133143e-05, 0.133383e-05, 0.133588e-05
     * /

!       Nuc.Int Total X-sec
          data ( muntx(i), i=  1,   34)/
     1  0.131090e-03, 0.144395e-03, 0.158387e-03, 0.172382e-03,
     2  0.187685e-03, 0.202425e-03, 0.218363e-03, 0.234279e-03,
     3  0.251198e-03, 0.268982e-03, 0.286553e-03, 0.305978e-03,
     4  0.324348e-03, 0.335809e-03, 0.325788e-03, 0.318789e-03,
     5  0.316786e-03, 0.314891e-03, 0.313318e-03, 0.312124e-03,
     6  0.311322e-03, 0.310982e-03, 0.311065e-03, 0.311680e-03,
     7  0.312847e-03, 0.314774e-03, 0.317242e-03, 0.319977e-03,
     8  0.323055e-03, 0.326491e-03, 0.330301e-03, 0.334529e-03,
     9  0.339117e-03, 0.343873e-03
     * /

!        Nuc.Int. dE/dx/E (v<vmin) /(g/cm2)
          data ( mundedx0(i), i=  1,   34)/
     1  0.000000e+00, 0.000000e+00, 0.000000e+00, 0.000000e+00,
     2  0.000000e+00, 0.000000e+00, 0.000000e+00, 0.000000e+00,
     3  0.000000e+00, 0.000000e+00, 0.000000e+00, 0.000000e+00,
     4  0.000000e+00, 0.220165e-09, 0.915062e-09, 0.139798e-08,
     5  0.167378e-08, 0.189541e-08, 0.207196e-08, 0.220734e-08,
     6  0.230651e-08, 0.238182e-08, 0.243526e-08, 0.247029e-08,
     7  0.248771e-08, 0.248822e-08, 0.248032e-08, 0.247220e-08,
     8  0.246348e-08, 0.245385e-08, 0.244290e-08, 0.242993e-08,
     9  0.241769e-08, 0.241067e-08
     * /

!        Nuc.Int; dE/dx/E /(g/cm2)
          data ( mundedxt(i), i=  1,   34)/
     1  0.402333e-06, 0.402855e-06, 0.402707e-06, 0.402319e-06,
     2  0.401983e-06, 0.401883e-06, 0.402184e-06, 0.403040e-06,
     3  0.404575e-06, 0.406782e-06, 0.409625e-06, 0.413063e-06,
     4  0.417039e-06, 0.421500e-06, 0.426368e-06, 0.431588e-06,
     5  0.437116e-06, 0.442897e-06, 0.448886e-06, 0.455043e-06,
     6  0.461337e-06, 0.467737e-06, 0.474219e-06, 0.480764e-06,
     7  0.487353e-06, 0.493973e-06, 0.500611e-06, 0.507257e-06,
     8  0.513904e-06, 0.520544e-06, 0.527173e-06, 0.533785e-06,
     9  0.540378e-06, 0.546948e-06
     * /

!
!      Nuc. Int. sampling table.
!      
          data ( muntbl(i,   1), i=  1,   68)/
     1 -0.172472e+01,-0.145269e+01,-0.125199e+01,-0.111430e+01,
     2 -0.101163e+01,-0.930870e+00,-0.865230e+00,-0.810390e+00,
     3 -0.763540e+00,-0.722800e+00,-0.686860e+00,-0.654790e+00,
     4 -0.625860e+00,-0.599560e+00,-0.575470e+00,-0.553250e+00,
     5 -0.532640e+00,-0.513430e+00,-0.495450e+00,-0.478620e+00,
     6 -0.462850e+00,-0.448010e+00,-0.434010e+00,-0.420760e+00,
     7 -0.408200e+00,-0.396260e+00,-0.384890e+00,-0.374040e+00,
     8 -0.363660e+00,-0.353720e+00,-0.344190e+00,-0.335030e+00,
     9 -0.326220e+00,-0.317730e+00,-0.309550e+00,-0.301640e+00,
     a -0.294000e+00,-0.286610e+00,-0.279440e+00,-0.272500e+00,
     b -0.265760e+00,-0.259210e+00,-0.252850e+00,-0.246650e+00,
     c -0.240630e+00,-0.234750e+00,-0.229020e+00,-0.223430e+00,
     d -0.217980e+00,-0.212650e+00,-0.207430e+00,-0.202330e+00,
     e -0.197340e+00,-0.192450e+00,-0.187660e+00,-0.182970e+00,
     f -0.178360e+00,-0.173830e+00,-0.169390e+00,-0.165030e+00,
     g -0.160740e+00,-0.156510e+00,-0.152360e+00,-0.148270e+00,
     h -0.144240e+00,-0.140260e+00,-0.136350e+00,-0.132480e+00          
     * /   
          data ( muntbl(i,   1), i=  69,   101)/
     1 -0.128660e+00,-0.124890e+00,-0.121160e+00,-0.117480e+00,
     2 -0.113830e+00,-0.110220e+00,-0.106640e+00,-0.103090e+00,
     3 -0.995700e-01,-0.960700e-01,-0.926000e-01,-0.891400e-01,
     4 -0.857000e-01,-0.822700e-01,-0.788500e-01,-0.754400e-01,
     5 -0.720200e-01,-0.686100e-01,-0.651800e-01,-0.617300e-01,
     6 -0.582600e-01,-0.547600e-01,-0.512200e-01,-0.476200e-01,
     7 -0.439500e-01,-0.401800e-01,-0.363100e-01,-0.321800e-01,
     8 -0.279000e-01,-0.225300e-01,-0.167700e-01,-0.132000e-01,
     9 -0.113000e-02                                                    
     * /   
          data ( muntbl(i,   2), i=  1,   68)/
     1 -0.192472e+01,-0.161718e+01,-0.139522e+01,-0.124252e+01,
     2 -0.112845e+01,-0.103907e+01,-0.966380e+00,-0.905490e+00,
     3 -0.853330e+00,-0.807830e+00,-0.767530e+00,-0.731410e+00,
     4 -0.698690e+00,-0.668970e+00,-0.641900e+00,-0.617080e+00,
     5 -0.594210e+00,-0.573010e+00,-0.553280e+00,-0.534850e+00,
     6 -0.517560e+00,-0.501300e+00,-0.485950e+00,-0.471420e+00,
     7 -0.457640e+00,-0.444530e+00,-0.432040e+00,-0.420110e+00,
     8 -0.408700e+00,-0.397770e+00,-0.387270e+00,-0.377180e+00,
     9 -0.367470e+00,-0.358110e+00,-0.349070e+00,-0.340340e+00,
     a -0.331890e+00,-0.323710e+00,-0.315770e+00,-0.308080e+00,
     b -0.300600e+00,-0.293330e+00,-0.286260e+00,-0.279370e+00,
     c -0.272660e+00,-0.266120e+00,-0.259730e+00,-0.253490e+00,
     d -0.247400e+00,-0.241430e+00,-0.235600e+00,-0.229890e+00,
     e -0.224290e+00,-0.218800e+00,-0.213420e+00,-0.208140e+00,
     f -0.202950e+00,-0.197860e+00,-0.192850e+00,-0.187930e+00,
     g -0.183080e+00,-0.178310e+00,-0.173600e+00,-0.168970e+00,
     h -0.164400e+00,-0.159900e+00,-0.155450e+00,-0.151050e+00          
     * /   
          data ( muntbl(i,   2), i=  69,   101)/
     1 -0.146710e+00,-0.142420e+00,-0.138170e+00,-0.133970e+00,
     2 -0.129810e+00,-0.125690e+00,-0.121600e+00,-0.117540e+00,
     3 -0.113510e+00,-0.109510e+00,-0.105530e+00,-0.101570e+00,
     4 -0.976200e-01,-0.936900e-01,-0.897600e-01,-0.858400e-01,
     5 -0.819100e-01,-0.779900e-01,-0.740500e-01,-0.701000e-01,
     6 -0.661200e-01,-0.621000e-01,-0.580400e-01,-0.539200e-01,
     7 -0.497200e-01,-0.454300e-01,-0.409800e-01,-0.363800e-01,
     8 -0.314100e-01,-0.262600e-01,-0.184200e-01,-0.141700e-01,
     9 -0.710000e-03                                                    
     * /   
          data ( muntbl(i,   3), i=  1,   68)/
     1 -0.212472e+01,-0.178196e+01,-0.153838e+01,-0.137024e+01,
     2 -0.124475e+01,-0.114634e+01,-0.106600e+01,-0.998400e+00,
     3 -0.940160e+00,-0.889050e+00,-0.844020e+00,-0.804090e+00,
     4 -0.768300e+00,-0.735960e+00,-0.706500e+00,-0.679500e+00,
     5 -0.654610e+00,-0.631550e+00,-0.610080e+00,-0.590020e+00,
     6 -0.571190e+00,-0.553470e+00,-0.536740e+00,-0.520910e+00,
     7 -0.505870e+00,-0.491570e+00,-0.477930e+00,-0.464890e+00,
     8 -0.452420e+00,-0.440450e+00,-0.428960e+00,-0.417900e+00,
     9 -0.407250e+00,-0.396970e+00,-0.387040e+00,-0.377430e+00,
     a -0.368130e+00,-0.359120e+00,-0.350370e+00,-0.341870e+00,
     b -0.333600e+00,-0.325560e+00,-0.317730e+00,-0.310090e+00,
     c -0.302640e+00,-0.295360e+00,-0.288250e+00,-0.281310e+00,
     d -0.274510e+00,-0.267870e+00,-0.261370e+00,-0.255010e+00,
     e -0.248780e+00,-0.242660e+00,-0.236670e+00,-0.230790e+00,
     f -0.225010e+00,-0.219330e+00,-0.213750e+00,-0.208260e+00,
     g -0.202860e+00,-0.197550e+00,-0.192310e+00,-0.187150e+00,
     h -0.182060e+00,-0.177040e+00,-0.172090e+00,-0.167190e+00          
     * /   
          data ( muntbl(i,   3), i=  69,   101)/
     1 -0.162360e+00,-0.157580e+00,-0.152850e+00,-0.148180e+00,
     2 -0.143540e+00,-0.138960e+00,-0.134410e+00,-0.129890e+00,
     3 -0.125410e+00,-0.120960e+00,-0.116540e+00,-0.112140e+00,
     4 -0.107760e+00,-0.103390e+00,-0.990400e-01,-0.946900e-01,
     5 -0.903400e-01,-0.859900e-01,-0.816200e-01,-0.772400e-01,
     6 -0.728300e-01,-0.683800e-01,-0.638800e-01,-0.593100e-01,
     7 -0.546700e-01,-0.499000e-01,-0.450000e-01,-0.398900e-01,
     8 -0.345400e-01,-0.286700e-01,-0.206700e-01,-0.147500e-01,
     9 -0.450000e-03                                                    
     * /   
          data ( muntbl(i,   4), i=  1,   68)/
     1 -0.232472e+01,-0.194684e+01,-0.168123e+01,-0.149722e+01,
     2 -0.135994e+01,-0.125171e+01,-0.116276e+01,-0.108731e+01,
     3 -0.102278e+01,-0.967050e+00,-0.918210e+00,-0.874890e+00,
     4 -0.836080e+00,-0.800980e+00,-0.769020e+00,-0.739710e+00,
     5 -0.712670e+00,-0.687610e+00,-0.664270e+00,-0.642440e+00,
     6 -0.621950e+00,-0.602640e+00,-0.584400e+00,-0.567110e+00,
     7 -0.550680e+00,-0.535040e+00,-0.520100e+00,-0.505820e+00,
     8 -0.492130e+00,-0.478980e+00,-0.466370e+00,-0.454250e+00,
     9 -0.442580e+00,-0.431330e+00,-0.420480e+00,-0.409990e+00,
     a -0.399850e+00,-0.390020e+00,-0.380490e+00,-0.371250e+00,
     b -0.362270e+00,-0.353540e+00,-0.345040e+00,-0.336760e+00,
     c -0.328690e+00,-0.320820e+00,-0.313140e+00,-0.305630e+00,
     d -0.298300e+00,-0.291120e+00,-0.284090e+00,-0.277210e+00,
     e -0.270470e+00,-0.263850e+00,-0.257360e+00,-0.250990e+00,
     f -0.244730e+00,-0.238580e+00,-0.232530e+00,-0.226580e+00,
     g -0.220720e+00,-0.214950e+00,-0.209260e+00,-0.203660e+00,
     h -0.198130e+00,-0.192670e+00,-0.187290e+00,-0.181960e+00          
     * /   
          data ( muntbl(i,   4), i=  69,   101)/
     1 -0.176700e+00,-0.171500e+00,-0.166360e+00,-0.161260e+00,
     2 -0.156220e+00,-0.151210e+00,-0.146260e+00,-0.141340e+00,
     3 -0.136450e+00,-0.131600e+00,-0.126770e+00,-0.121970e+00,
     4 -0.117190e+00,-0.112420e+00,-0.107660e+00,-0.102910e+00,
     5 -0.981700e-01,-0.934100e-01,-0.886400e-01,-0.838600e-01,
     6 -0.790400e-01,-0.741800e-01,-0.692700e-01,-0.642900e-01,
     7 -0.592200e-01,-0.540300e-01,-0.486900e-01,-0.431300e-01,
     8 -0.373000e-01,-0.309300e-01,-0.235700e-01,-0.151500e-01,
     9 -0.280000e-03                                                    
     * /   
          data ( muntbl(i,   5), i=  1,   68)/
     1 -0.252472e+01,-0.211161e+01,-0.182344e+01,-0.162318e+01,
     2 -0.147320e+01,-0.135396e+01,-0.125521e+01,-0.117290e+01,
     3 -0.110307e+01,-0.104274e+01,-0.989840e+00,-0.942900e+00,
     4 -0.900800e+00,-0.862710e+00,-0.827980e+00,-0.796100e+00,
     5 -0.766670e+00,-0.739340e+00,-0.713860e+00,-0.689990e+00,
     6 -0.667580e+00,-0.646520e+00,-0.626670e+00,-0.607900e+00,
     7 -0.590120e+00,-0.573220e+00,-0.557120e+00,-0.541770e+00,
     8 -0.527090e+00,-0.513020e+00,-0.499530e+00,-0.486560e+00,
     9 -0.474080e+00,-0.462060e+00,-0.450450e+00,-0.439240e+00,
     a -0.428390e+00,-0.417880e+00,-0.407700e+00,-0.397810e+00,
     b -0.388210e+00,-0.378870e+00,-0.369780e+00,-0.360930e+00,
     c -0.352300e+00,-0.343890e+00,-0.335670e+00,-0.327640e+00,
     d -0.319790e+00,-0.312110e+00,-0.304590e+00,-0.297220e+00,
     e -0.290000e+00,-0.282920e+00,-0.275970e+00,-0.269140e+00,
     f -0.262440e+00,-0.255850e+00,-0.249370e+00,-0.242990e+00,
     g -0.236710e+00,-0.230520e+00,-0.224420e+00,-0.218410e+00,
     h -0.212480e+00,-0.206620e+00,-0.200840e+00,-0.195130e+00          
     * /   
          data ( muntbl(i,   5), i=  69,   101)/
     1 -0.189480e+00,-0.183900e+00,-0.178370e+00,-0.172900e+00,
     2 -0.167470e+00,-0.162100e+00,-0.156770e+00,-0.151480e+00,
     3 -0.146230e+00,-0.141010e+00,-0.135830e+00,-0.130660e+00,
     4 -0.125520e+00,-0.120400e+00,-0.115280e+00,-0.110180e+00,
     5 -0.105070e+00,-0.999600e-01,-0.948400e-01,-0.896900e-01,
     6 -0.845200e-01,-0.793000e-01,-0.740300e-01,-0.686800e-01,
     7 -0.632400e-01,-0.576700e-01,-0.519400e-01,-0.459900e-01,
     8 -0.397200e-01,-0.329800e-01,-0.256100e-01,-0.154800e-01,
     9 -0.180000e-03                                                    
     * /   
          data ( muntbl(i,   6), i=  1,   68)/
     1 -0.272472e+01,-0.227609e+01,-0.196461e+01,-0.174754e+01,
     2 -0.158359e+01,-0.145190e+01,-0.134483e+01,-0.125594e+01,
     3 -0.118044e+01,-0.111513e+01,-0.105780e+01,-0.100684e+01,
     4 -0.961060e+00,-0.919560e+00,-0.881650e+00,-0.846920e+00,
     5 -0.814980e+00,-0.785470e+00,-0.758060e+00,-0.732510e+00,
     6 -0.708590e+00,-0.686130e+00,-0.664970e+00,-0.644970e+00,
     7 -0.626030e+00,-0.608040e+00,-0.590910e+00,-0.574580e+00,
     8 -0.558970e+00,-0.544030e+00,-0.529690e+00,-0.515920e+00,
     9 -0.502680e+00,-0.489910e+00,-0.477600e+00,-0.465710e+00,
     a -0.454210e+00,-0.443070e+00,-0.432280e+00,-0.421810e+00,
     b -0.411640e+00,-0.401750e+00,-0.392140e+00,-0.382770e+00,
     c -0.373640e+00,-0.364740e+00,-0.356040e+00,-0.347550e+00,
     d -0.339250e+00,-0.331140e+00,-0.323190e+00,-0.315410e+00,
     e -0.307780e+00,-0.300300e+00,-0.292960e+00,-0.285750e+00,
     f -0.278680e+00,-0.271720e+00,-0.264870e+00,-0.258130e+00,
     g -0.251500e+00,-0.244960e+00,-0.238520e+00,-0.232160e+00,
     h -0.225890e+00,-0.219690e+00,-0.213580e+00,-0.207530e+00          
     * /   
          data ( muntbl(i,   6), i=  69,   101)/
     1 -0.201550e+00,-0.195630e+00,-0.189780e+00,-0.183980e+00,
     2 -0.178230e+00,-0.172530e+00,-0.166870e+00,-0.161260e+00,
     3 -0.155680e+00,-0.150140e+00,-0.144620e+00,-0.139130e+00,
     4 -0.133660e+00,-0.128210e+00,-0.122770e+00,-0.117330e+00,
     5 -0.111890e+00,-0.106440e+00,-0.100980e+00,-0.954900e-01,
     6 -0.899700e-01,-0.844000e-01,-0.787700e-01,-0.730600e-01,
     7 -0.672400e-01,-0.613000e-01,-0.551800e-01,-0.488300e-01,
     8 -0.421600e-01,-0.350400e-01,-0.270700e-01,-0.158900e-01,
     9 -0.110000e-03                                                    
     * /   
          data ( muntbl(i,   7), i=  1,   68)/
     1 -0.292472e+01,-0.244015e+01,-0.210466e+01,-0.186977e+01,
     2 -0.169000e+01,-0.154712e+01,-0.143174e+01,-0.133577e+01,
     3 -0.125410e+01,-0.118330e+01,-0.112098e+01,-0.106546e+01,
     4 -0.101583e+01,-0.971160e+00,-0.930630e+00,-0.893620e+00,
     5 -0.859620e+00,-0.828230e+00,-0.799110e+00,-0.771980e+00,
     6 -0.746620e+00,-0.722810e+00,-0.700410e+00,-0.679260e+00,
     7 -0.659240e+00,-0.640240e+00,-0.622170e+00,-0.604960e+00,
     8 -0.588520e+00,-0.572790e+00,-0.557710e+00,-0.543250e+00,
     9 -0.529340e+00,-0.515950e+00,-0.503050e+00,-0.490590e+00,
     a -0.478550e+00,-0.466900e+00,-0.455620e+00,-0.444660e+00,
     b -0.434030e+00,-0.423690e+00,-0.413630e+00,-0.403840e+00,
     c -0.394290e+00,-0.384970e+00,-0.375880e+00,-0.366990e+00,
     d -0.358310e+00,-0.349800e+00,-0.341480e+00,-0.333320e+00,
     e -0.325330e+00,-0.317480e+00,-0.309780e+00,-0.302210e+00,
     f -0.294770e+00,-0.287460e+00,-0.280260e+00,-0.273180e+00,
     g -0.266200e+00,-0.259320e+00,-0.252530e+00,-0.245840e+00,
     h -0.239230e+00,-0.232700e+00,-0.226250e+00,-0.219870e+00          
     * /   
          data ( muntbl(i,   7), i=  69,   101)/
     1 -0.213560e+00,-0.207320e+00,-0.201130e+00,-0.195010e+00,
     2 -0.188930e+00,-0.182910e+00,-0.176930e+00,-0.170990e+00,
     3 -0.165080e+00,-0.159210e+00,-0.153370e+00,-0.147550e+00,
     4 -0.141750e+00,-0.135970e+00,-0.130200e+00,-0.124430e+00,
     5 -0.118650e+00,-0.112870e+00,-0.107060e+00,-0.101230e+00,
     6 -0.953700e-01,-0.894500e-01,-0.834600e-01,-0.774000e-01,
     7 -0.712200e-01,-0.649000e-01,-0.584000e-01,-0.516500e-01,
     8 -0.445800e-01,-0.370000e-01,-0.285400e-01,-0.164300e-01,
     9 -0.700000e-04                                                    
     * /   
          data ( muntbl(i,   8), i=  1,   68)/
     1 -0.300000e+01,-0.246289e+01,-0.211859e+01,-0.187312e+01,
     2 -0.168713e+01,-0.154212e+01,-0.142447e+01,-0.132610e+01,
     3 -0.124206e+01,-0.116974e+01,-0.110669e+01,-0.105107e+01,
     4 -0.100149e+01,-0.956920e+00,-0.916540e+00,-0.879710e+00,
     5 -0.845940e+00,-0.814790e+00,-0.785940e+00,-0.759110e+00,
     6 -0.734050e+00,-0.710570e+00,-0.688510e+00,-0.667700e+00,
     7 -0.648020e+00,-0.629360e+00,-0.611610e+00,-0.594710e+00,
     8 -0.578560e+00,-0.563120e+00,-0.548320e+00,-0.534110e+00,
     9 -0.520440e+00,-0.507280e+00,-0.494600e+00,-0.482340e+00,
     a -0.470500e+00,-0.459030e+00,-0.447920e+00,-0.437140e+00,
     b -0.426670e+00,-0.416490e+00,-0.406590e+00,-0.396950e+00,
     c -0.387550e+00,-0.378380e+00,-0.369430e+00,-0.360680e+00,
     d -0.352130e+00,-0.343760e+00,-0.335560e+00,-0.327540e+00,
     e -0.319660e+00,-0.311940e+00,-0.304360e+00,-0.296910e+00,
     f -0.289590e+00,-0.282400e+00,-0.275310e+00,-0.268340e+00,
     g -0.261470e+00,-0.254700e+00,-0.248030e+00,-0.241440e+00,
     h -0.234940e+00,-0.228520e+00,-0.222180e+00,-0.215910e+00          
     * /   
          data ( muntbl(i,   8), i=  69,   101)/
     1 -0.209710e+00,-0.203570e+00,-0.197490e+00,-0.191460e+00,
     2 -0.185490e+00,-0.179570e+00,-0.173690e+00,-0.167860e+00,
     3 -0.162060e+00,-0.156290e+00,-0.150550e+00,-0.144840e+00,
     4 -0.139140e+00,-0.133460e+00,-0.127790e+00,-0.122130e+00,
     5 -0.116460e+00,-0.110780e+00,-0.105090e+00,-0.993600e-01,
     6 -0.936100e-01,-0.878000e-01,-0.819200e-01,-0.759700e-01,
     7 -0.699100e-01,-0.637100e-01,-0.573300e-01,-0.507100e-01,
     8 -0.437600e-01,-0.363100e-01,-0.277500e-01,-0.187300e-01,
     9 -0.400000e-04                                                    
     * /   
          data ( muntbl(i,   9), i=  1,   68)/
     1 -0.300000e+01,-0.243876e+01,-0.208535e+01,-0.183439e+01,
     2 -0.164808e+01,-0.150166e+01,-0.138275e+01,-0.128472e+01,
     3 -0.120212e+01,-0.113132e+01,-0.106974e+01,-0.101555e+01,
     4 -0.967370e+00,-0.924170e+00,-0.885140e+00,-0.849610e+00,
     5 -0.817050e+00,-0.787040e+00,-0.759240e+00,-0.733370e+00,
     6 -0.709210e+00,-0.686550e+00,-0.665230e+00,-0.645120e+00,
     7 -0.626090e+00,-0.608030e+00,-0.590860e+00,-0.574500e+00,
     8 -0.558870e+00,-0.543920e+00,-0.529590e+00,-0.515830e+00,
     9 -0.502600e+00,-0.489850e+00,-0.477560e+00,-0.465690e+00,
     a -0.454220e+00,-0.443110e+00,-0.432340e+00,-0.421900e+00,
     b -0.411760e+00,-0.401900e+00,-0.392310e+00,-0.382960e+00,
     c -0.373860e+00,-0.364980e+00,-0.356310e+00,-0.347830e+00,
     d -0.339550e+00,-0.331450e+00,-0.323520e+00,-0.315740e+00,
     e -0.308120e+00,-0.300650e+00,-0.293310e+00,-0.286110e+00,
     f -0.279030e+00,-0.272060e+00,-0.265220e+00,-0.258470e+00,
     g -0.251830e+00,-0.245290e+00,-0.238840e+00,-0.232480e+00,
     h -0.226200e+00,-0.220000e+00,-0.213880e+00,-0.207820e+00          
     * /   
          data ( muntbl(i,   9), i=  69,   101)/
     1 -0.201830e+00,-0.195910e+00,-0.190050e+00,-0.184240e+00,
     2 -0.178480e+00,-0.172770e+00,-0.167100e+00,-0.161480e+00,
     3 -0.155890e+00,-0.150340e+00,-0.144810e+00,-0.139310e+00,
     4 -0.133830e+00,-0.128360e+00,-0.122910e+00,-0.117460e+00,
     5 -0.112010e+00,-0.106550e+00,-0.101070e+00,-0.955700e-01,
     6 -0.900400e-01,-0.844600e-01,-0.788100e-01,-0.730900e-01,
     7 -0.672600e-01,-0.613100e-01,-0.551800e-01,-0.488200e-01,
     8 -0.421400e-01,-0.349700e-01,-0.265800e-01,-0.181800e-01,
     9 -0.300000e-04                                                    
     * /   
          data ( muntbl(i,   10), i=  1,   68)/
     1 -0.300000e+01,-0.242609e+01,-0.205915e+01,-0.180517e+01,
     2 -0.161495e+01,-0.146767e+01,-0.135004e+01,-0.125346e+01,
     3 -0.117244e+01,-0.110329e+01,-0.104339e+01,-0.990760e+00,
     4 -0.943970e+00,-0.901990e+00,-0.864000e+00,-0.829390e+00,
     5 -0.797650e+00,-0.768380e+00,-0.741270e+00,-0.716030e+00,
     6 -0.692440e+00,-0.670320e+00,-0.649500e+00,-0.629860e+00,
     7 -0.611260e+00,-0.593620e+00,-0.576840e+00,-0.560840e+00,
     8 -0.545570e+00,-0.530950e+00,-0.516930e+00,-0.503480e+00,
     9 -0.490540e+00,-0.478070e+00,-0.466050e+00,-0.454450e+00,
     a -0.443220e+00,-0.432360e+00,-0.421820e+00,-0.411610e+00,
     b -0.401690e+00,-0.392040e+00,-0.382660e+00,-0.373520e+00,
     c -0.364620e+00,-0.355930e+00,-0.347450e+00,-0.339170e+00,
     d -0.331070e+00,-0.323150e+00,-0.315390e+00,-0.307790e+00,
     e -0.300340e+00,-0.293030e+00,-0.285860e+00,-0.278820e+00,
     f -0.271900e+00,-0.265100e+00,-0.258410e+00,-0.251830e+00,
     g -0.245340e+00,-0.238950e+00,-0.232650e+00,-0.226440e+00,
     h -0.220310e+00,-0.214260e+00,-0.208290e+00,-0.202380e+00          
     * /   
          data ( muntbl(i,   10), i=  69,   101)/
     1 -0.196540e+00,-0.190760e+00,-0.185040e+00,-0.179370e+00,
     2 -0.173760e+00,-0.168190e+00,-0.162670e+00,-0.157190e+00,
     3 -0.151750e+00,-0.146340e+00,-0.140950e+00,-0.135600e+00,
     4 -0.130260e+00,-0.124940e+00,-0.119630e+00,-0.114320e+00,
     5 -0.109010e+00,-0.103700e+00,-0.983700e-01,-0.930200e-01,
     6 -0.876300e-01,-0.822000e-01,-0.767200e-01,-0.711500e-01,
     7 -0.654800e-01,-0.596900e-01,-0.537300e-01,-0.475400e-01,
     8 -0.410500e-01,-0.340600e-01,-0.258200e-01,-0.178000e-01,
     9 -0.200000e-04                                                    
     * /   
          data ( muntbl(i,   11), i=  1,   68)/
     1 -0.300000e+01,-0.240477e+01,-0.203076e+01,-0.177108e+01,
     2 -0.158100e+01,-0.143556e+01,-0.132020e+01,-0.122609e+01,
     3 -0.114729e+01,-0.107998e+01,-0.102156e+01,-0.970180e+00,
     4 -0.924480e+00,-0.883450e+00,-0.846310e+00,-0.812460e+00,
     5 -0.781400e+00,-0.752760e+00,-0.726210e+00,-0.701490e+00,
     6 -0.678380e+00,-0.656710e+00,-0.636310e+00,-0.617050e+00,
     7 -0.598820e+00,-0.581530e+00,-0.565070e+00,-0.549380e+00,
     8 -0.534400e+00,-0.520060e+00,-0.506310e+00,-0.493110e+00,
     9 -0.480420e+00,-0.468190e+00,-0.456390e+00,-0.445000e+00,
     a -0.433990e+00,-0.423330e+00,-0.412990e+00,-0.402970e+00,
     b -0.393230e+00,-0.383770e+00,-0.374560e+00,-0.365600e+00,
     c -0.356860e+00,-0.348340e+00,-0.340020e+00,-0.331890e+00,
     d -0.323950e+00,-0.316170e+00,-0.308570e+00,-0.301110e+00,
     e -0.293810e+00,-0.286640e+00,-0.279610e+00,-0.272710e+00,
     f -0.265930e+00,-0.259260e+00,-0.252700e+00,-0.246250e+00,
     g -0.239890e+00,-0.233630e+00,-0.227460e+00,-0.221380e+00,
     h -0.215380e+00,-0.209450e+00,-0.203600e+00,-0.197810e+00          
     * /   
          data ( muntbl(i,   11), i=  69,   101)/
     1 -0.192090e+00,-0.186440e+00,-0.180840e+00,-0.175300e+00,
     2 -0.169800e+00,-0.164360e+00,-0.158960e+00,-0.153600e+00,
     3 -0.148270e+00,-0.142980e+00,-0.137720e+00,-0.132480e+00,
     4 -0.127260e+00,-0.122060e+00,-0.116870e+00,-0.111690e+00,
     5 -0.106500e+00,-0.101310e+00,-0.961100e-01,-0.908800e-01,
     6 -0.856200e-01,-0.803200e-01,-0.749600e-01,-0.695200e-01,
     7 -0.639900e-01,-0.583300e-01,-0.525100e-01,-0.464600e-01,
     8 -0.401200e-01,-0.333000e-01,-0.252000e-01,-0.174700e-01,
     9 -0.100000e-04                                                    
     * /   
          data ( muntbl(i,   12), i=  1,   68)/
     1 -0.300000e+01,-0.238211e+01,-0.199913e+01,-0.173882e+01,
     2 -0.155119e+01,-0.140901e+01,-0.129651e+01,-0.120452e+01,
     3 -0.112740e+01,-0.106148e+01,-0.100422e+01,-0.953830e+00,
     4 -0.909000e+00,-0.868730e+00,-0.832260e+00,-0.799010e+00,
     5 -0.768490e+00,-0.740340e+00,-0.714240e+00,-0.689940e+00,
     6 -0.667210e+00,-0.645890e+00,-0.625820e+00,-0.606880e+00,
     7 -0.588940e+00,-0.571910e+00,-0.555720e+00,-0.540280e+00,
     8 -0.525520e+00,-0.511410e+00,-0.497870e+00,-0.484870e+00,
     9 -0.472370e+00,-0.460330e+00,-0.448710e+00,-0.437500e+00,
     a -0.426650e+00,-0.416150e+00,-0.405970e+00,-0.396100e+00,
     b -0.386510e+00,-0.377190e+00,-0.368130e+00,-0.359300e+00,
     c -0.350690e+00,-0.342300e+00,-0.334110e+00,-0.326110e+00,
     d -0.318290e+00,-0.310630e+00,-0.303140e+00,-0.295810e+00,
     e -0.288620e+00,-0.281570e+00,-0.274650e+00,-0.267850e+00,
     f -0.261180e+00,-0.254620e+00,-0.248170e+00,-0.241820e+00,
     g -0.235570e+00,-0.229410e+00,-0.223340e+00,-0.217360e+00,
     h -0.211450e+00,-0.205630e+00,-0.199870e+00,-0.194190e+00          
     * /   
          data ( muntbl(i,   12), i=  69,   101)/
     1 -0.188570e+00,-0.183010e+00,-0.177510e+00,-0.172060e+00,
     2 -0.166660e+00,-0.161310e+00,-0.156010e+00,-0.150740e+00,
     3 -0.145510e+00,-0.140320e+00,-0.135150e+00,-0.130010e+00,
     4 -0.124880e+00,-0.119780e+00,-0.114680e+00,-0.109600e+00,
     5 -0.104510e+00,-0.994100e-01,-0.943100e-01,-0.891800e-01,
     6 -0.840200e-01,-0.788100e-01,-0.735600e-01,-0.682200e-01,
     7 -0.628000e-01,-0.572400e-01,-0.515300e-01,-0.456000e-01,
     8 -0.393800e-01,-0.327000e-01,-0.247300e-01,-0.172000e-01,
     9 -0.100000e-04                                                    
     * /   
          data ( muntbl(i,   13), i=  1,   68)/
     1 -0.300000e+01,-0.235981e+01,-0.197039e+01,-0.171262e+01,
     2 -0.152890e+01,-0.138951e+01,-0.127903e+01,-0.118860e+01,
     3 -0.111274e+01,-0.104785e+01,-0.991460e+00,-0.941810e+00,
     4 -0.897620e+00,-0.857910e+00,-0.821940e+00,-0.789130e+00,
     5 -0.759020e+00,-0.731230e+00,-0.705470e+00,-0.681470e+00,
     6 -0.659030e+00,-0.637970e+00,-0.618140e+00,-0.599420e+00,
     7 -0.581700e+00,-0.564870e+00,-0.548870e+00,-0.533610e+00,
     8 -0.519030e+00,-0.505070e+00,-0.491690e+00,-0.478840e+00,
     9 -0.466480e+00,-0.454580e+00,-0.443090e+00,-0.432000e+00,
     a -0.421280e+00,-0.410900e+00,-0.400840e+00,-0.391080e+00,
     b -0.381600e+00,-0.372380e+00,-0.363420e+00,-0.354690e+00,
     c -0.346180e+00,-0.337890e+00,-0.329790e+00,-0.321880e+00,
     d -0.314150e+00,-0.306580e+00,-0.299180e+00,-0.291930e+00,
     e -0.284820e+00,-0.277850e+00,-0.271020e+00,-0.264300e+00,
     f -0.257710e+00,-0.251230e+00,-0.244850e+00,-0.238580e+00,
     g -0.232400e+00,-0.226320e+00,-0.220330e+00,-0.214420e+00,
     h -0.208590e+00,-0.202830e+00,-0.197150e+00,-0.191540e+00          
     * /   
          data ( muntbl(i,   13), i=  69,   101)/
     1 -0.185990e+00,-0.180500e+00,-0.175070e+00,-0.169690e+00,
     2 -0.164360e+00,-0.159080e+00,-0.153850e+00,-0.148650e+00,
     3 -0.143490e+00,-0.138370e+00,-0.133270e+00,-0.128190e+00,
     4 -0.123140e+00,-0.118110e+00,-0.113080e+00,-0.108060e+00,
     5 -0.103050e+00,-0.980200e-01,-0.929900e-01,-0.879300e-01,
     6 -0.828400e-01,-0.777100e-01,-0.725300e-01,-0.672700e-01,
     7 -0.619200e-01,-0.564400e-01,-0.508100e-01,-0.449700e-01,
     8 -0.388200e-01,-0.322500e-01,-0.243800e-01,-0.169900e-01,
     9  0.000000e+00                                                    
     * /   
          data ( muntbl(i,   14), i=  1,   68)/
     1 -0.300000e+01,-0.234033e+01,-0.195058e+01,-0.169674e+01,
     2 -0.151557e+01,-0.137792e+01,-0.126871e+01,-0.117927e+01,
     3 -0.110419e+01,-0.103994e+01,-0.984080e+00,-0.934890e+00,
     4 -0.891080e+00,-0.851710e+00,-0.816040e+00,-0.783500e+00,
     5 -0.753630e+00,-0.726060e+00,-0.700490e+00,-0.676670e+00,
     6 -0.654390e+00,-0.633490e+00,-0.613810e+00,-0.595220e+00,
     7 -0.577620e+00,-0.560910e+00,-0.545020e+00,-0.529860e+00,
     8 -0.515380e+00,-0.501520e+00,-0.488230e+00,-0.475460e+00,
     9 -0.463180e+00,-0.451360e+00,-0.439950e+00,-0.428930e+00,
     a -0.418280e+00,-0.407960e+00,-0.397970e+00,-0.388270e+00,
     b -0.378850e+00,-0.369700e+00,-0.360790e+00,-0.352120e+00,
     c -0.343670e+00,-0.335430e+00,-0.327380e+00,-0.319520e+00,
     d -0.311840e+00,-0.304320e+00,-0.296970e+00,-0.289770e+00,
     e -0.282710e+00,-0.275780e+00,-0.268990e+00,-0.262320e+00,
     f -0.255770e+00,-0.249330e+00,-0.243000e+00,-0.236770e+00,
     g -0.230640e+00,-0.224600e+00,-0.218640e+00,-0.212770e+00,
     h -0.206980e+00,-0.201270e+00,-0.195630e+00,-0.190050e+00          
     * /   
          data ( muntbl(i,   14), i=  69,   101)/
     1 -0.184540e+00,-0.179090e+00,-0.173700e+00,-0.168360e+00,
     2 -0.163080e+00,-0.157830e+00,-0.152640e+00,-0.147480e+00,
     3 -0.142360e+00,-0.137270e+00,-0.132210e+00,-0.127170e+00,
     4 -0.122160e+00,-0.117160e+00,-0.112180e+00,-0.107200e+00,
     5 -0.102220e+00,-0.972300e-01,-0.922400e-01,-0.872200e-01,
     6 -0.821700e-01,-0.770800e-01,-0.719400e-01,-0.667200e-01,
     7 -0.614100e-01,-0.559800e-01,-0.504000e-01,-0.446000e-01,
     8 -0.385000e-01,-0.319800e-01,-0.241800e-01,-0.168600e-01,
     9  0.000000e+00                                                    
     * /   
          data ( muntbl(i,   15), i=  1,   68)/
     1 -0.300000e+01,-0.232414e+01,-0.193753e+01,-0.168662e+01,
     2 -0.150721e+01,-0.137075e+01,-0.126241e+01,-0.117362e+01,
     3 -0.109906e+01,-0.103522e+01,-0.979710e+00,-0.930800e+00,
     4 -0.887250e+00,-0.848090e+00,-0.812610e+00,-0.780240e+00,
     5 -0.750520e+00,-0.723080e+00,-0.697630e+00,-0.673920e+00,
     6 -0.651740e+00,-0.630930e+00,-0.611340e+00,-0.592830e+00,
     7 -0.575300e+00,-0.558670e+00,-0.542840e+00,-0.527740e+00,
     8 -0.513320e+00,-0.499510e+00,-0.486270e+00,-0.473560e+00,
     9 -0.461330e+00,-0.449550e+00,-0.438180e+00,-0.427200e+00,
     a -0.416590e+00,-0.406310e+00,-0.396360e+00,-0.386690e+00,
     b -0.377310e+00,-0.368190e+00,-0.359320e+00,-0.350680e+00,
     c -0.342260e+00,-0.334050e+00,-0.326030e+00,-0.318200e+00,
     d -0.310550e+00,-0.303060e+00,-0.295730e+00,-0.288550e+00,
     e -0.281520e+00,-0.274620e+00,-0.267860e+00,-0.261210e+00,
     f -0.254680e+00,-0.248270e+00,-0.241960e+00,-0.235750e+00,
     g -0.229640e+00,-0.223630e+00,-0.217700e+00,-0.211850e+00,
     h -0.206080e+00,-0.200390e+00,-0.194770e+00,-0.189220e+00          
     * /   
          data ( muntbl(i,   15), i=  69,   101)/
     1 -0.183730e+00,-0.178300e+00,-0.172930e+00,-0.167610e+00,
     2 -0.162350e+00,-0.157130e+00,-0.151950e+00,-0.146810e+00,
     3 -0.141710e+00,-0.136640e+00,-0.131610e+00,-0.126590e+00,
     4 -0.121600e+00,-0.116620e+00,-0.111660e+00,-0.106700e+00,
     5 -0.101740e+00,-0.967800e-01,-0.918000e-01,-0.868100e-01,
     6 -0.817800e-01,-0.767100e-01,-0.715900e-01,-0.664000e-01,
     7 -0.611100e-01,-0.557100e-01,-0.501500e-01,-0.443700e-01,
     8 -0.383000e-01,-0.318200e-01,-0.240600e-01,-0.167800e-01,
     9  0.000000e+00                                                    
     * /   
          data ( muntbl(i,   16), i=  1,   68)/
     1 -0.300000e+01,-0.231337e+01,-0.193063e+01,-0.168158e+01,
     2 -0.150326e+01,-0.136752e+01,-0.125968e+01,-0.117126e+01,
     3 -0.109697e+01,-0.103336e+01,-0.978030e+00,-0.929270e+00,
     4 -0.885830e+00,-0.846780e+00,-0.811390e+00,-0.779100e+00,
     5 -0.749440e+00,-0.722060e+00,-0.696670e+00,-0.673010e+00,
     6 -0.650870e+00,-0.630100e+00,-0.610540e+00,-0.592070e+00,
     7 -0.574570e+00,-0.557960e+00,-0.542150e+00,-0.527080e+00,
     8 -0.512680e+00,-0.498900e+00,-0.485680e+00,-0.472980e+00,
     9 -0.460770e+00,-0.449000e+00,-0.437650e+00,-0.426690e+00,
     a -0.416090e+00,-0.405820e+00,-0.395880e+00,-0.386230e+00,
     b -0.376860e+00,-0.367750e+00,-0.358880e+00,-0.350250e+00,
     c -0.341840e+00,-0.333640e+00,-0.325630e+00,-0.317810e+00,
     d -0.310170e+00,-0.302690e+00,-0.295370e+00,-0.288200e+00,
     e -0.281170e+00,-0.274280e+00,-0.267520e+00,-0.260880e+00,
     f -0.254360e+00,-0.247950e+00,-0.241650e+00,-0.235450e+00,
     g -0.229350e+00,-0.223340e+00,-0.217410e+00,-0.211570e+00,
     h -0.205810e+00,-0.200120e+00,-0.194510e+00,-0.188960e+00          
     * /   
          data ( muntbl(i,   16), i=  69,   101)/
     1 -0.183480e+00,-0.178060e+00,-0.172690e+00,-0.167380e+00,
     2 -0.162120e+00,-0.156900e+00,-0.151730e+00,-0.146600e+00,
     3 -0.141510e+00,-0.136440e+00,-0.131410e+00,-0.126400e+00,
     4 -0.121410e+00,-0.116440e+00,-0.111480e+00,-0.106530e+00,
     5 -0.101580e+00,-0.966200e-01,-0.916500e-01,-0.866600e-01,
     6 -0.816400e-01,-0.765800e-01,-0.714600e-01,-0.662800e-01,
     7 -0.610000e-01,-0.556000e-01,-0.500500e-01,-0.442900e-01,
     8 -0.382200e-01,-0.317500e-01,-0.240000e-01,-0.167400e-01,
     9  0.000000e+00                                                    
     * /   
          data ( muntbl(i,   17), i=  1,   68)/
     1 -0.300000e+01,-0.231061e+01,-0.192993e+01,-0.168175e+01,
     2 -0.150387e+01,-0.136836e+01,-0.126066e+01,-0.117231e+01,
     3 -0.109807e+01,-0.103448e+01,-0.979160e+00,-0.930400e+00,
     4 -0.886960e+00,-0.847900e+00,-0.812490e+00,-0.780180e+00,
     5 -0.750510e+00,-0.723120e+00,-0.697700e+00,-0.674030e+00,
     6 -0.651880e+00,-0.631090e+00,-0.611510e+00,-0.593020e+00,
     7 -0.575500e+00,-0.558880e+00,-0.543050e+00,-0.527970e+00,
     8 -0.513550e+00,-0.499750e+00,-0.486510e+00,-0.473800e+00,
     9 -0.461570e+00,-0.449790e+00,-0.438430e+00,-0.427450e+00,
     a -0.416830e+00,-0.406560e+00,-0.396600e+00,-0.386930e+00,
     b -0.377550e+00,-0.368420e+00,-0.359550e+00,-0.350900e+00,
     c -0.342480e+00,-0.334260e+00,-0.326240e+00,-0.318410e+00,
     d -0.310750e+00,-0.303260e+00,-0.295930e+00,-0.288750e+00,
     e -0.281710e+00,-0.274810e+00,-0.268030e+00,-0.261380e+00,
     f -0.254850e+00,-0.248430e+00,-0.242120e+00,-0.235910e+00,
     g -0.229790e+00,-0.223770e+00,-0.217830e+00,-0.211980e+00,
     h -0.206200e+00,-0.200510e+00,-0.194880e+00,-0.189320e+00          
     * /   
          data ( muntbl(i,   17), i=  69,   101)/
     1 -0.183830e+00,-0.178390e+00,-0.173020e+00,-0.167690e+00,
     2 -0.162420e+00,-0.157200e+00,-0.152010e+00,-0.146870e+00,
     3 -0.141770e+00,-0.136690e+00,-0.131650e+00,-0.126630e+00,
     4 -0.121630e+00,-0.116650e+00,-0.111680e+00,-0.106710e+00,
     5 -0.101750e+00,-0.967800e-01,-0.918000e-01,-0.868000e-01,
     6 -0.817700e-01,-0.767000e-01,-0.715700e-01,-0.663700e-01,
     7 -0.610800e-01,-0.556800e-01,-0.501100e-01,-0.443300e-01,
     8 -0.382600e-01,-0.317800e-01,-0.240200e-01,-0.167400e-01,
     9  0.000000e+00                                                    
     * /   

      end
!     *******************
      block data cblkobs
!     *******************
      implicit none

















!       Zobs.h     header file for observation sites definition
!
           integer maxnoofsites, maxnoofassites, horizontal,
     *     perpendicular, notused, spherical
       parameter(

     *    maxnoofsites = 20,


     *    maxnoofassites=20,

     *    notused = 0,           ! detector plane is not used
     *    horizontal = 1,        ! detector is horizontal
     *    perpendicular = 2,      ! detector is pependicular to 1ry.
     *    spherical = 3          ! detector is cocentric sphere as the earth
     *            )

!       Zobsvp.h---parameters to be given by input.
!       This must be preceded by Zobs.h

!   (-> ---------------------------------------------------

    real*8  heightlist  !1  Height of observation levels in m. This is  made from DepthList internally. 
                            ! I.e., this one is usually not an input. However, if the DepthList values are 
                            ! negative, this is used as input and corresponding DepthList is computed internally.
        real*8  depthlist   !1  Depth List of Observation level in kg/m$^2$. If $< 0$, HeightList has priority. 
                            !  (See HeightList)
        real*8  asheightlist    !1  This is HeightList for Air Shower observ.  Used only if Generate contains
                            !  "as". See  HeightList.
        real*8  asdepthlist     !1  This is DepthList for AS observation.  Used only if Generate contains 
                            ! "as". See DepthList.
        real*8  latitofsite     !1  Latitude of the deepest observation level in degree.  East is positive.
        real*8  longitofsite    !1  Longitude of the deepest observation level in degree.  North is positive.

        real*8  dtgmt           !1  Difference of the local time of the observation place from GMT (hour).
    real*8  yearofgeomag    !1  Like 1999.5. Year when Geomagnetic field is to be calculated.
    integer obsplane        !1    How to observe particles. \newline
                                !    0$ \Rightarrow $ no detector plane is used for observation. BorderHeightL
                                !    and BorderHeightH are used to detect particles. This is for, say, neutrino
                                !    observation. See BorderHeight{L,H}. However, the primary is directed to
                                !    the deepest depth.  \newline
                                !    1,-1$ \Rightarrow $ detector at the observation place is horizontal. Note 
                                !    that the horizontal means not tangential plane, but rather a spherical surface \newline
                                !    2,-2$ \Rightarrow $ detector is perpendicular to the primary. \newline
                                !    3$ \Rightarrow $ spherical observation. See text. \newline
                                !    For ObsPlane={1,2}, the user observation routine will receive coordinate values in
                                !    the corresponding detector system. However, if it is 0, 3 or negative, Exyz values
                                !    are obtained.
        integer noofsites2    !2   No of Sites for particle observation; not to be touched; for skeleton/flesh use.
    real*8 xaxisfromsouth   !2 Angle between the horizontal detector X-axis and the south(deg). + is counter
                                ! clockwise.  If $|$XaxisFromSouth$| > 360$, it is computed so that the direction is
                                ! to the magnetic east at the deepest observation point. Default is 361.
!   <-) --------------------------------------------

   

        common /zobsc/
     *   heightlist(maxnoofsites),
     *   depthlist(maxnoofsites),
     *   asheightlist(maxnoofassites),
     *   asdepthlist(maxnoofassites),
     *   latitofsite, longitofsite, dtgmt,
     *   xaxisfromsouth, yearofgeomag,
     *   obsplane, noofsites2

          
       data 
     * depthlist /maxnoofsites*0./ ,
     * heightlist /maxnoofsites*0./ ,
     * asdepthlist /maxnoofassites*0./ ,
     * asheightlist /maxnoofassites*0./ ,
     * dtgmt /8.0/ ,
     * latitofsite /30.11/ ,
     * longitofsite /90.53/ ,
     * obsplane /1/ ,   
     * xaxisfromsouth /361.0/,
     * yearofgeomag /2000.5/
      end
!     *************************
      block data  cblktracking
!     *************************
      implicit none

!           Parameters used for Tracking.
!   (-> --------------------------------------------

    logical  exactthick   !2  If T, a given length is converted into thickness with best accuracy even for very 
                              !  inclined trajectory by using numerical integration.
        logical incmuonpolari  !1  if T, consider muon polarization
        logical freec          !1  if F, the first interaction point is forced to be the injection point else
                               !   the interaction poin is randomly sampled.
    integer onedim         !1  If 0, 3 dimensional simulation. if $\ge$1, one
                               !  dimensional simulation is performed.  \newline
                               !  1: onedim without use of table. \newline
                               !  2: table is used for thickness $ \leftrightarrow$ length conversion. if cos $<$ .5 \newline
                               !  3: table is always used for any angle.
                               !  ( for height $>$ 30 km, table is not used in any case). 
        real*8  lamordiv       !2  In the geomagnetic field, a charged particle can travel almost streight 
                               !  in (Lamor Radius)/LamorDiv.  Default is 5. For AMS like tracking 20 may be needed.
    real*8  truncc         !2 coeff. for truncating path.
    real*8  truncn         !2 coeff. for truncating path.
    real*8  truncx         !2 coeff. for truncating path.
    real*8  keminobs(8)      !1  The min kinetic energy of particles for observation.
    real*8  keminobs2(8)     !2  Don't touch this. skeleton/flesh use.
    real*8  ratiotoe0      !2 In the A.S generation,  hadronic interactions are followed down to at
                               !  least  RatioToE0 * E0/nucleon energy.
        real*8  waitratio      !1  Wait A.S generation until the electron energy, Ee, becomes $<$ WaitRatio* E0. 
                               !   This many be 1.0 for hadron origin case. But for gamma/electron primary, 
                               !   this should be as low as 0.01 to enjoy full fluctuation.
        integer endlevel !2 Used for skeleton/flesh-out job. In a normal job, system default value 0 is reset by
            ! the system to be the max number of observation levels. (=NoOfSites).  Its real use is in such a
            ! skeleton/flesh-out job that one first follows the particles up to some high depth and later chooses
            ! events and flesh them out to deeper depths. In such a skeleton-making job, the user must give the
            ! depth list which is used flesh-out job, too.  In the skeleton job, particle tracking is terminated
            ! at the level specified by EndLevel.  In such a flesh-out job, the user must give a larger value 
            ! or 0 to EndLevel 
        integer endlevel2 !2  Don't worry. This  is system use.
    integer trace   !1  Flag for trace information output.\newline
            !   0 $\rightarrow$ no trace information is output.\newline
            !  $<$10$\rightarrow$ x, y, z in the primary system(say, 1)\newline
            !  $<$20 $\rightarrow$ x, y, in the primary sys. z in kg/m$^2$.(say,11)\newline
            !  $<$30 $\rightarrow$ x, y, z in the detector system\newline
            !  $<$40 $\rightarrow$ x, y, in the detector system. z in kg/m2\newline 
            !  $<$50 $\rightarrow$ x, y, z in 'xyz' system.\newline
            !  $<$60 $\rightarrow$ x,y, in 'xyz' and z in kg/m2\newline
            !  61-100 $\rightarrow$  for Cherenkov observation. For Coord system,  subtract 60.\newline
            !   if the value is even,  binary output is made on TraceDev.\newline
            !   if the last digit is 1 or 2, trace is always taken. if the last digit is 3 or 4, trace is taken
            !   only if the particle is located below the heighest observation depeth. 
            !  $>$ 101  $\rightarrow$ subtract 100 and apply the above, but chookTrace or chookCeren are used.\newline
            !  Primary system:  Origin is the deepest detector.  Z-axis is the primary direction. 
            !                   X-axis is Z x Vertical axis.  X-Y plane is orthogonal to the primary.\newline
            !  Detector system: origin is the deepest detector. Z-axis is the vertical one.  X-axis is 
            !                   directed to the magnetic east.  X-Y plane is horizontal.\newline
            !  z in kg/m$^2$ :     Vertical depth in kg/m$^2$  above the  deepest detector to the particle.
    integer tracedev      !2  Logical dev \# for TraceDir/trace1,2,.... 
        character*70  tracedir !1 Directory.  Default Trace information is put TraceDir/trace1, 2,..
                               ! for event 1, 2, ... The directory should exist.  Default is ' ' and in this case
                               ! /tmp/YourLoginName/ is employed. The last "/" should not be given.
                               ! *** NOTE that default Cherenkov output is made only using TraceDev,
                               !    TraceDir is not used.  You have to open the disk file at chookbgRun
                               !    It can by binary or ascii file depending on Trace value.
    logical thinsampling  !1  if F, thinsampling is not tried. if T, alla Hillas thinning. Don't use with
                              !    the  skeleton/flesh method 
        real*8  ethinratio    !2  if ThinsSamplig != F, thin sampling is performed if the energy of a particle is
                              !   $<$ EthinRatio * PrimaryEnergy(/nucleon) (=Ethin) ( EtinRatio$>$ 0).
                              !  If EthinRatio $<$ 0, Ethin will be |EthinRatio| (GeV). 
        logical timestructure !1  If T,  time information is computed
    integer howgeomag     !2  if 1, no magnetic field until first coll. \newline
                              !      2, mag.f always exists. If Reverse not=0, use this. \newline
                              !     11, same as 1 but mag.f is const. \newline
                              !     12, same as 2 but mag.f is const. \newline
                              !     21, same as 1 but mag.f is const. \newline
                              !     22, same as 2 but mag.f is const.  \newline
                              !     31, same as 1 but mag.f is dependent on the position.  \newline
                              !  const value is the one at deepest observation plane. for 11,12  or should be given by
                              !  MagN, MagE, MagD for 21, 22. For normal applications,  11 is good.
                              !  If no magnetic field is applied, energy loss by dE/dx is considered.(bef.4.92,
                              !  and aft. 5.14)
        real*8  magn          !2 See HowGeomag (in Tesla)
    real*8  mage          !2 See HowGeomag (in Tesla)
        real*8  magd          !2 See HowGeomag (in Tesla)

        real*8  magchgdist    !2 Distance where mag. can be seen  as const.(m) at sea level
        integer userungekutta !2 How to calculate deflection by the geomagnetic field. Let L be the distance
                              !   the  particle travels. \newline
                              ! 0$\rightarrow$Don't use RungeKutta method. Use the solution assuming the constant B, which
                              !     is exact if B is const.  Since the particle path is made  short, this is
                              !     enough for normal cases where particles are inside the atmosphere.(default) \newline
                              !     In every case below, if the particle height is $<$ 30km  
                              !    (= cheight in ccomPathEnd.f), the same method as 0 is used. \newline
                              ! 1$\rightarrow$ Use the Euler method.  Time needed is 20\% more than the 0 case.
                              !      As B, use the value at L/2 point obtained by using the current direction. \newline
                              ! 2$\rightarrow$ mixture of 1 and Runge-Kutta-Gill method. If gradient of B is large, RKG is
                              !      employed. This needs $\sim$4 times more cpu time than case of 1 when making a
                              !      cutoff table.  The step size of RKG is $\sim$1/10 of the Lamore radius. \newline
                              ! 3$\rightarrow$ The same as 2 but use the Runge-Kutta-Fehlberg method instead of RKG.
                              !      Step size is automattically adjusted ($\sim$1/20 $\sim$1/30 of Lamor radius) \newline
                              ! 4$\rightarrow$ As a middle point, use the point obtained by assuming the constant B at
                              !      initial point. If grad B is still large, use RKG. \newline
                              ! 5$\rightarrow$ The same as 4 but us RKF instead of RKG. \newline
                              ! 6$\rightarrow$ Use always RKG \newline
                              ! 7$\rightarrow$ Use always RKF.  This takes very long time.(50 times of 0). \newline
    real*8  borderheighth !2 If a particle goes higher than this, discard it.  This should be larger than 
                              !  HeightOfInj or 0.  
                              !  If 0, it is adjusted to be the same as HeightOfInj. NOTE: For upgoin primary cases, you have
                              !  to set this one explicitly. 
    real*8  borderheightl !2 If a particle reaches this hight, call observation routine. No further tracking is done. 
                              !  This is for neutrino observation.  See ObsPlane.
        real*8  backanglimit  !2 If the cosine of the angle between a particle and the primary becomes smaller than
                              !  this value, the particle is discarded. See also BorderHeighH. If you give a value 
                              !  less than -1.0, such rejection will never happen.   Default is -1.0
    character*16 generate !1 specify what should be generated \newline
                              !   1)  Electro-magnetic cascade(em), \newline
                              !   2)  one dimensional  hybrid AS(as/qas) and/or \newline
                              !   3)  AS Lateral distribution(lat). \newline
                              !    If Generate= ' ', hadronic cascade shower is generated. \newline
                              !  For example, you may give as follows: \newline
                              !    Generate='em,as' or 'em/as' (order/case/separator insensitive) is to  generate EM-cascade and AS. \newline
                              !    Generate='as' will  generate AS with some  adequate EM cascade (EM cascade is automatically generated
                              !    so that hybrid A.S can be observed, but the minimum energy in EM cascade is independent of KEminObs). \newline
                              !   If 'qas' is given, quick generation of AS for heavy primaries is tried. See  chookASbyH.f 

        character*16 generate2 !2  don't touch this.  for skeleton/flesh use.

    integer magbrem       !2 If 0, no magnetic bremsstrahlung is considered. \newline
                              ! if 1 and Ee $>$ MagBremEmin, energy loss due to magnetic brems is considered \newline
                              ! if 2 and Ee $>$ MagBremEmin, real sampling of gamma is performed. \newline
                              ! (note, actually upsilon is referred further).
                              ! if generate='as' with really high energy primaries, WaitRatio
                              ! must be made small so that WaitRatio*E0 $\sim$ MagBremEmin 
        integer magpair       !2 If 0, no magnetic pair creation is considered. \newline
                              !  if 1 and Eg > MagPairEmin, real sampling is tried.
                              ! (note, actually upsilon is referred further).  To see these magnetic effects,
                              !  HowGeoMag=2 and HightOfInj $\sim$ 5000 km are desirable.

        logical lpmeffect     !1 If t, the LPM effect is  considered when Ee $>$ LpmBremEmin for electrons and
                              !        Eg $>$ LpmPairEmin for gamma rays.

        real*8  magbrememin   !2  E $>$ this, magnetic bremsstrahlung by electrons may be considered. However, if
                              !   MagBrem = 0, not considered at all \newline
                              !   MagBrem = 1, total energy loss due to brems is considered. \newline
                              !   MagBrem = 2, gamma energy is sampled actually. \newline
                              !   If upsilon (Ee/m * B/Bcr) is small, the effective treatment will be
                              !   the same as MagBrem = 0 case.
        real*8  magpairemin   !2  E $>$ this, magnetic pair creation by gamma may be considered. However, if
                              !   MagPair = 0, not considered at all. \newline
                              !   MagPair = 1, pair creation is sampled.  \newline
                              !   However, again, actual occurrence will be dependent on the angle between
                              !   B and photon direction.
        real*8  upsilonmin    !2  Magnetic bremsstralhung is considered only if upsilon $>$ UpsilonMin.
        real*8  lpmbrememin   !2  The LPM effect is taken into account for bremsstrahlung when LpmEffect is .true. 
                              !   and the electron energy is higher than this.
        real*8  lpmpairemin   !2  The LPM effect is taken into account for pair creation when LpmEffect is .true.
                              !    and the gamma energy is higher than this.
        integer reverse       !2  0$\rightarrow$ Normal tracking. \newline
                              !   1$\rightarrow$ incident is tracked to a direction opposite to the given one.
                              !       the incident is charge-conjugated.
                              !       All interactions are ignored. (Use when to make cut-off table or to see
                              !       a given particle (say, observed anti proton) can go out of Earth. \newline
                              !   2$\rightarrow$ same as 1 but energy gain (not loss) is taken into account
                              !   TimeStructure should be T if Reverse != 0.  See BackAnglLimit.

        real*8 pathlimit      !2  If the sum of (path/beta) of a particle  exceeds this, it is judged as dead.
                              !   (to avoid infinite cyclotron loop).  However, for normal applications,
                              !   this will not be effective because of BackAnglLimit. See Reverse.
                              !   TimeStructure should be T if Reverse != 0 and PathLimit is to be effective.

       integer muni           !2 0$\rightarrow$ nuclear interaction of muon is completely neglected \newline
                              !  1$\rightarrow$ energy loss by n.i is subsumed in dE/dx of muons as a continuous energy loss.  Let v=
                              !     Etransfer/Emu,  the loss here is Int(vc:vmax) of (Emu vdsigma/dv).  (vc $\sim$0, vmax$\sim$1). \newline
                              !  2$\rightarrow$ (Default value). similar to 1 but as the continuous loss only v $<$ vmin=10$^{-3}$ of
                              !     fractional muon energy is subsumed (Int(vc: vmin) of (Emu vdsigma/dv)).  The portion
                              !     of loss by v$>$vmin is treated as a stocastic  process.  However, the product from the
                              !     n.i itself is neglected \newline
                              !  3$\rightarrow$ the same as 2, but the n.i is explicitly included to produce a number of particles.  
                              !     The n.i is treated as a photo-nucleus interaction.
      integer mubr            !2  parameter similar to MuNI but for bremsstrahlung by muons.
       integer mupr           !2  parameter similar to MuNI but for pair creation by muons.

!    <-)    ----------------------------------------------

    common /cztracp/ truncc, truncn, truncx, 
     *  keminobs, keminobs2, ratiotoe0, pathlimit,
     *  waitratio,  ethinratio, backanglimit, lamordiv,
     *  borderheighth,  magn, mage, magd, magchgdist,
     *  borderheightl,  muni, mubr, mupr,
     *  magbrememin, magpairemin, upsilonmin, lpmbrememin, 
     *  lpmpairemin, userungekutta,
     *  thinsampling, timestructure, howgeomag,
     *  trace, tracedev,  exactthick, onedim, reverse,
     *  freec,  incmuonpolari, magbrem, magpair, lpmeffect,
     *  endlevel, endlevel2
 
        common /cztrackpc/ generate, generate2, tracedir

      data 
     *  exactthick /.false./ ,
     *  freec /.true./ ,
     *  tracedev /21/ , 
     *  tracedir /' '/ ,
     *  borderheighth /0./ ,
     *  borderheightl /0./ ,
     *  howgeomag /11/ ,
     *  magn /0./ ,
     *  mage /0./ ,
     *  magd /0./ ,
     *  ratiotoe0 /1.e-5/ ,
     *  userungekutta /0/ ,
     *  magchgdist /20.e3/ ,
     *  timestructure /.true./ ,
     *  trace /0/ ,
     *  waitratio /1./ ,    
     *  truncc/5.0/ ,
     *  truncn /1.e-2/ ,
     *  truncx /2./

      data
     * incmuonpolari /.true./ ,
     * keminobs /2*500d-6,6*50d-3/ ,
     * thinsampling /.false./,
     * ethinratio /3.e-5/ ,
     * generate /'em'/ , 
     * generate2 /' '/ ,
     * backanglimit  /-1.0/ ,
     * onedim /0/,
     * magbrem /2/,
     * magpair /1/,
     * lpmeffect /.true./,
     * magpairemin /2.e10/,         ! Eg > this ==> magnetic pair considered.> 3x10^19 eV
     * magbrememin /3.e8/,          ! Ee > this ==> magnetic brems considered > 10^18 eV
     * upsilonmin  /3.e-3/,         ! if upsilon is < this, no magnetic brem
     * lpmpairemin /1.e9/,         ! Eg > LpmPairEmin && LpmEffect ==> LPM effect for g
     * lpmbrememin /1.e8/           ! Ee > LpmBremEmin && LpmEffect ==> LPM effect for e

       data
     * reverse /0/,
     * pathlimit /13000d4/,         !  ~ 20 x Eradius
     * endlevel/0/,
     * muni /2/, 
     * mubr /2/,
     * mupr /2/,
     * lamordiv /5./
      end
!           block data for Zxsectionp
      block data cblkxsec
!----      include 'Zxsectionp.h'
!           common parameters for tracking
!
!   (->  -------------------------------------

        real*8  deltpp  !2  p-p xsection increases as $E^{Deltpp}$(E$>$ 100GeV)
        real*8  deltpip !2  pi-p xsection increases as $E^{Deltpip}$ (E$>$ 100GeV)
        real*8  deltkp  !2  k-p xsection increases as $E^{Deltkp}$ (E$>$ 100GeV)
        real*8  increasexsec !2   how the xsection increases. 1.0$\rightarrow$  power of E

!   <-)  -----------------------------------

        common /zxsectionp/
     *  deltpp, deltpip, deltkp, increasexsec
!

!
      data 
     *  deltpp /0.08/ ,
     *  deltpip /0.08/ ,
     *  deltkp /0.08/ , 
     *  increasexsec /1.0/
      end
      block data cblkdedx

      common /zdedxair/ stha, sthb, sthc, sthx0, sthx1, sthsa,
     *   wlg0, w0, betasq, tmax, norm, w0inmev,
     *   knckon, jdef

    real*8  stha, sthb, sthc, sthx0, sthx1, sthsa
    real*8  wlg0, w0, betasq, tmax, norm, w0inmev
    integer jdef
        logical knckon

!

      data jdef/0/
!
!        table of sternheimer's consts. by p.r.b vol.3 (1971)3681
!     for air
!                        
      data stha/7.68e-02/
      data sthb/ 17.86/
      data sthc/ -10.79/
      data sthx0/1.809/
      data sthx1/4.0/
      data sthsa/ 0.234/
      end

      block data  cblkspecial

      integer modifyx
      real*8  modifyxpw1, modifyxpw2
      common /zspecial/  modifyxpw1, modifyxpw2, modifyx


       data  
     * modifyx /0/ ,
     * modifyxpw1 /0.5/ ,
     * modifyxpw2 /0.5/
      end

      subroutine cmymain
!         this list block data names as external names so that the
!       block data is surely enabled.





!   dpmjet cannot be used on NEXTSTEP, so
!   you have to make the next 0. 



!            make DEBUG > 0 depending on the debug purpose. 


!
!   choose:    Old atmosphere or new segmented atmosphere
!            define 
!               old atmosphere --> 0
!           or  new with c-spline
!               new atmosphere --> 1
!           or  new with linear interp.
!               new atmosphere --> 2


!     if you want to put a lable on each particle to identify that
!     the one and the same particle crosses a given observation
!     plane more than once, make this 1 or 2.  Then the same particle
!     will have the  same label number in track record.
!     ( aTrack.label ).  If this is 0, aTrack.lable record dose not
!     exists. 
!     If 1; after any interaction (except for continuous energy
!     loss by dE/dx and deflection by B or scattering), label is
!     changed.
!     If 2: For knockon and Bremstrahlung, the survival particle
!     will have the same label. In the case of Moller scattring
!     higher enregy electrons are regarded as the survival one.
!

!     if you want to have a detailed info. for particle tracking
!     make the below >=1.  The user observation routine is called
!     with the following id  on the following  conditions:
!              chookobs(a, id)
!     1)  if it is >=1,  a particle is going to interact at a point given in
!         the track information, id=4
!     2)  if it is >=1,  a particle is going to die, id=5
!     3)  if it is >=2,  a particle is being discarded due to the large
!          angle (cos(angle relative to the parent) > BackAngLimit). id=6
!     4)  if it is >=3,  a particle makes a step. id=7
!        







      external cblkelemag
      external cblkdedx
      external cblkheavy 
      external cblkmanager 

      external cblkxsec
      external cblkevhnp
      external cblkincident
      external cblkobs
      external cblktracking
      external cblkmuint
!              For Lund related block common
      external blkhd1, blkhd2, blkhd3, blkhd4
      external blkdc1, blkdc2
      external ludatac
      external luedatc
      external luhdatc
      external cblkspeical

         call cmanager   ! call Manager for Cosmos Simulation
      end