c(((((((((((((((((((((((((((((((( zgenas (((((((((((((((((((((((((((((((
c          folloing include consts for index of gamma, proton etc
      include 'zsimcn'
             common /zgenas/ height, ehght,
     1  k1ry, e0, elg, a, b, c, d, tm, toff, cos1ry, tcog, zob,
     2 jscint, discle, disclp, deone, deonea(nidsc), aneb, modex,
     3 sage, rmu, rmu2, enff, frcde, efne,  firstz, zhcol, elgpn,
     4 e0pn, ihg1, mass1, nucc,  tcogh, znuc,  aneh,  shvy
c
       dimension tcogh(maxms), znuc(maxms), aneh(maxms), shvy(maxms)
c
       real*8 height, ehght
c
c
             common/ zgenac/
     1 idsc(nidsc), versn
       character*16 idsc
       character*8  versn
c
c height: first collision point height (slant) in m from the
c         observation point.
c  ehght: effective first collision height of the current shower
c         for heavy 1ry. (ne-weighted average of various heights)
c   k1ry: 1 ry particle type (0=gamma, 1=proton)
c     e0: 1 ry energy in tev
c   e0pn: per nucleon energy of 1ry.  for gamma, it is e0
c    elg: log10(e0)
c  elgpn: log10(e0pn)
c     a,b,c,d,tm,toff: coefficients for an air shower generated
c cos1ry: cos of zenith angle of 1ry
c   tcog: slant depth of cog of the current shower (g/cm**2)
c         for heavy, average of tcogh
c  tcogh: tcogh(i) is the tcog of i-th shower in heavy 1ry
c  znuc:  znuc(i) is the collision depth of i-th nucleon in heavy 1ry
c  aneh:  aneh(i) is the ne(b) by the i-th shower at observation depth
c  shvy:  shvy(i) is the age of the i-th shower at observation depth
c   nucc: # of nucleons in heavy 1ry
c  mass1: =nucc
c  ihg1:  heavy group index of the 1ry.  for gamma, undef
c    zob: specified observation depth (vertical g/cm**2)
c jscint: 0---> scintillator structure not yet specifed
c         1---> .1cmfe + 3.5cm scinti
c         2---> .5cmpb + .1cm fe + 3.5cm scinti
c discle: discrimination level in terms of energy loss in scinti
c         (for vertical relativistic particles)
c disclp: discrimination level in terms of equiv. ptcl # in scinti
c         (for vertical relativistic particles)
c  deone: energy loss in the current scinti in mev corresponding to
c            one particle.
c deonea: array to contain energy loss for one particle for each
c         structure of scinti.
c   aneb: ne corresponding to b-approx for the current shower
c         for heavy, sum of aneh
c  modex: 0---> 1 dimensional sim.
c         1---> 4 dimensional sim.
c   sage: age of the shower ;obtained when depob is  called
c         for heavy, ne-weighted average of shvy
c    rmu: moliere unit (m) for the observation depth and current
c         shower
c   rmu2: rmu**2
c   enff: n(de>0)/ne(by b approx)  where n(de>0) is the # of ptcls
c         which give de>0 in scinti
c  frcde: over all fraction of particles with energy deposit >
c         discle (see m in kdig)
c   efne: effective # of particles which give de>discl in scinti.
c firstz: first collision depth (slant one) of the current
c         shower (in g/cm**2). for heavy the first depth where
c         interacting nucleon appeared.
c zhcol:  for heavy, interaction depth of 1ry (slant)
c   idsc: to contain id's for usable structure of scintillator.
c  versn: version # of genas
c
c)))))))))))))))))))))))))))))))) zgenas )))))))))))))))))))))))))))))))
