A production height table is formatted as

 I1 I2 I3 I4     p1      p2      p3      p4      p5      p6 ...  p21
  En1           h11     h12     h13     h14     h15     h16 ... h121
  En2           h21     h22     h23     h24     h25     h26 ... h221
  En3           h31     h32     h33     h34     h35     h36 ... h421
  .................

A sample for actual table is as follows (nu_mu for Kamioka)
  1  1  1  1   0.01    0.05    0.10    0.15    0.20    0.25 ...    0.99
  1.059E-01  2187.4  5379.9  7281.8  8634.7  9736.8 10715.2 ... 44116.4
  1.189E-01  2127.7  5352.3  7262.9  8607.7  9706.5 10685.0 ... 44098.6
  ------------------

In the first line,
I1 stands for the experimental site, and I2 for the flavor of the
neutrino. However, those are already shown as the file name.

I3 is for the zenith angle bin, and runs from 1 to 20 corresponding to
1 > cos(Theta_zenith) > 0.9, 0.9 > cos(Theta_zenith) > 0.8 ...
-0.9 > cos(Theta_zenith) > -1 .

I4 is for the azimuth angle bin, and runs from 1 to 12 corresponding to
0 < phai_azimuth < 30, 30 < phai_azimuth < 60, .... 330 < phai_azimuth <
360 in degree.

  Note, the zenith angle and azimuth angle are that for arrival
  direction, and are the same as the flux table including the binnings.
  Then, phai_azimuth = 0 stands for the South direction and phai_azimuth =
  90 degree stands for the East direction.

Also in the first line, p1, p2, .... p21 stand for the probabilities,
used in the following lines.

In the following lines the first column show the neutrino energy in GeV,
and the value in the 2nd column (h11) is the height that the neutrino is
produced below the in the probablity of p1. The value in 3rd column
(h12) is the height that the neutrino is produced below p2, and so on.
After the neutrino energy reaches 1.059E+4 GeV. New table for next
azimuth bin or next zenith angle bin starts.