*$ CREATE SOURCE.FOR
*COPY SOURCE
*
*=== source ===========================================================*
*
      SUBROUTINE SOURCE ( NOMORE )

      INCLUDE '(DBLPRC)'
      INCLUDE '(DIMPAR)'
      INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
* Version modified for FFs generation                                  *
      DIMENSION ASYM(0:235)
      DIMENSION ANSYM(0:235)
      DIMENSION ADIST(0:235)
      DIMENSION CUMUL(0:235)
      DIMENSION ENEDGE (0:1099)
      DIMENSION FLUX (0:1099)
      DIMENSION TEGR(0:1099)
*                                                                      *
*----------------------------------------------------------------------*
*
      INCLUDE '(BEAMCM)'
      INCLUDE '(FHEAVY)'
      INCLUDE '(FLKSTK)'
      INCLUDE '(IOIOCM)'
      INCLUDE '(LTCLCM)'
      INCLUDE '(PAPROP)'
      INCLUDE '(SOURCM)'
      INCLUDE '(SUMCOU)'
*
      LOGICAL LFIRST
*
      SAVE LFIRST
      DATA LFIRST / .TRUE. /
* Neutron energy group boundaries in MeV
      DATA ENEDGE /
     &   9.885531e-09,
     &   1.011579e-08,
     &   1.035142e-08,
     &   1.059254e-08,
     &   1.083927e-08,
     &   1.109175e-08,
     &   1.135011e-08,
     &   1.161449e-08,
     &   1.188502e-08,
     &   1.216186e-08,
     &   1.244515e-08,
     &   1.273503e-08,
     &   1.303167e-08,
     &   1.333521e-08,
     &   1.364583e-08,
     &   1.396368e-08,
     &   1.428894e-08,
     &   1.462177e-08,
     &   1.496236e-08,
     &   1.531087e-08,
     &   1.566751e-08,
     &   1.603245e-08,
     &   1.640590e-08,
     &   1.678804e-08,
     &   1.717908e-08,
     &   1.757924e-08,
     &   1.798871e-08,
     &   1.840772e-08,
     &   1.883649e-08,
     &   1.927525e-08,
     &   1.972423e-08,
     &   2.018366e-08,
     &   2.065380e-08,
     &   2.113489e-08,
     &   2.162719e-08,
     &   2.213095e-08,
     &   2.264644e-08,
     &   2.317395e-08,
     &   2.371374e-08,
     &   2.426610e-08,
     &   2.483133e-08,
     &   2.540973e-08,
     &   2.600160e-08,
     &   2.660725e-08,
     &   2.722701e-08,
     &   2.786121e-08,
     &   2.851018e-08,
     &   2.917427e-08,
     &   2.985383e-08,
     &   3.054921e-08,
     &   3.126079e-08,
     &   3.198895e-08,
     &   3.273407e-08,
     &   3.349654e-08,
     &   3.427678e-08,
     &   3.507519e-08,
     &   3.589219e-08,
     &   3.672823e-08,
     &   3.758374e-08,
     &   3.845918e-08,
     &   3.935501e-08,
     &   4.027170e-08,
     &   4.120975e-08,
     &   4.216965e-08,
     &   4.315191e-08,
     &   4.415704e-08,
     &   4.518559e-08,
     &   4.623810e-08,
     &   4.731513e-08,
     &   4.841724e-08,
     &   4.954502e-08,
     &   5.069907e-08,
     &   5.188000e-08,
     &   5.308844e-08,
     &   5.432503e-08,
     &   5.559043e-08,
     &   5.688529e-08,
     &   5.821032e-08,
     &   5.956621e-08,
     &   6.095369e-08,
     &   6.237348e-08,
     &   6.382635e-08,
     &   6.531306e-08,
     &   6.683439e-08,
     &   6.839116e-08,
     &   6.998420e-08,
     &   7.161434e-08,
     &   7.328245e-08,
     &   7.498942e-08,
     &   7.673615e-08,
     &   7.852356e-08,
     &   8.035261e-08,
     &   8.222426e-08,
     &   8.413951e-08,
     &   8.609938e-08,
     &   8.810489e-08,
     &   9.015711e-08,
     &   9.225714e-08,
     &   9.440609e-08,
     &   9.660509e-08,
     &   9.885531e-08,
     &   1.011579e-07,
     &   1.035142e-07,
     &   1.059254e-07,
     &   1.083927e-07,
     &   1.109175e-07,
     &   1.135011e-07,
     &   1.161449e-07,
     &   1.188502e-07,
     &   1.216186e-07,
     &   1.244515e-07,
     &   1.273503e-07,
     &   1.303167e-07,
     &   1.333521e-07,
     &   1.364583e-07,
     &   1.396368e-07,
     &   1.428894e-07,
     &   1.462177e-07,
     &   1.496236e-07,
     &   1.531087e-07,
     &   1.566751e-07,
     &   1.603245e-07,
     &   1.640590e-07,
     &   1.678804e-07,
     &   1.717908e-07,
     &   1.757924e-07,
     &   1.798871e-07,
     &   1.840772e-07,
     &   1.883649e-07,
     &   1.927525e-07,
     &   1.972423e-07,
     &   2.018366e-07,
     &   2.065380e-07,
     &   2.113489e-07,
     &   2.162719e-07,
     &   2.213095e-07,
     &   2.264644e-07,
     &   2.317395e-07,
     &   2.371374e-07,
     &   2.426610e-07,
     &   2.483133e-07,
     &   2.540973e-07,
     &   2.600160e-07,
     &   2.660725e-07,
     &   2.722701e-07,
     &   2.786121e-07,
     &   2.851018e-07,
     &   2.917427e-07,
     &   2.985383e-07,
     &   3.054921e-07,
     &   3.126079e-07,
     &   3.198895e-07,
     &   3.273407e-07,
     &   3.349654e-07,
     &   3.427678e-07,
     &   3.507519e-07,
     &   3.589219e-07,
     &   3.672823e-07,
     &   3.758374e-07,
     &   3.845918e-07,
     &   3.935501e-07,
     &   4.027170e-07,
     &   4.120975e-07,
     &   4.216965e-07,
     &   4.315191e-07,
     &   4.415704e-07,
     &   4.518559e-07,
     &   4.623810e-07,
     &   4.731513e-07,
     &   4.841724e-07,
     &   4.954502e-07,
     &   5.069907e-07,
     &   5.188000e-07,
     &   5.308844e-07,
     &   5.432503e-07,
     &   5.559043e-07,
     &   5.688529e-07,
     &   5.821032e-07,
     &   5.956621e-07,
     &   6.095369e-07,
     &   6.237348e-07,
     &   6.382635e-07,
     &   6.531306e-07,
     &   6.683439e-07,
     &   6.839116e-07,
     &   6.998420e-07,
     &   7.161434e-07,
     &   7.328245e-07,
     &   7.498942e-07,
     &   7.673615e-07,
     &   7.852356e-07,
     &   8.035261e-07,
     &   8.222426e-07,
     &   8.413951e-07,
     &   8.609938e-07,
     &   8.810489e-07,
     &   9.015711e-07,
     &   9.225714e-07,
     &   9.440609e-07,
     &   9.660509e-07,
     &   9.885531e-07,
     &   1.011579e-06,
     &   1.035142e-06,
     &   1.059254e-06,
     &   1.083927e-06,
     &   1.109175e-06,
     &   1.135011e-06,
     &   1.161449e-06,
     &   1.188502e-06,
     &   1.216186e-06,
     &   1.244515e-06,
     &   1.273503e-06,
     &   1.303167e-06,
     &   1.333521e-06,
     &   1.364583e-06,
     &   1.396368e-06,
     &   1.428894e-06,
     &   1.462177e-06,
     &   1.496236e-06,
     &   1.531087e-06,
     &   1.566751e-06,
     &   1.603245e-06,
     &   1.640590e-06,
     &   1.678804e-06,
     &   1.717908e-06,
     &   1.757924e-06,
     &   1.798871e-06,
     &   1.840772e-06,
     &   1.883649e-06,
     &   1.927525e-06,
     &   1.972423e-06,
     &   2.018366e-06,
     &   2.065380e-06,
     &   2.113489e-06,
     &   2.162719e-06,
     &   2.213095e-06,
     &   2.264644e-06,
     &   2.317395e-06,
     &   2.371374e-06,
     &   2.426610e-06,
     &   2.483133e-06,
     &   2.540973e-06,
     &   2.600160e-06,
     &   2.660725e-06,
     &   2.722701e-06,
     &   2.786121e-06,
     &   2.851018e-06,
     &   2.917427e-06,
     &   2.985383e-06,
     &   3.054921e-06,
     &   3.126079e-06,
     &   3.198895e-06,
     &   3.273407e-06,
     &   3.349654e-06,
     &   3.427678e-06,
     &   3.507519e-06,
     &   3.589219e-06,
     &   3.672823e-06,
     &   3.758374e-06,
     &   3.845918e-06,
     &   3.935501e-06,
     &   4.027170e-06,
     &   4.120975e-06,
     &   4.216965e-06,
     &   4.315191e-06,
     &   4.415704e-06,
     &   4.518559e-06,
     &   4.623810e-06,
     &   4.731513e-06,
     &   4.841724e-06,
     &   4.954502e-06,
     &   5.069907e-06,
     &   5.188000e-06,
     &   5.308844e-06,
     &   5.432503e-06,
     &   5.559043e-06,
     &   5.688529e-06,
     &   5.821032e-06,
     &   5.956621e-06,
     &   6.095369e-06,
     &   6.237348e-06,
     &   6.382635e-06,
     &   6.531306e-06,
     &   6.683439e-06,
     &   6.839116e-06,
     &   6.998420e-06,
     &   7.161434e-06,
     &   7.328245e-06,
     &   7.498942e-06,
     &   7.673615e-06,
     &   7.852356e-06,
     &   8.035261e-06,
     &   8.222426e-06,
     &   8.413951e-06,
     &   8.609938e-06,
     &   8.810489e-06,
     &   9.015711e-06,
     &   9.225714e-06,
     &   9.440609e-06,
     &   9.660509e-06,
     &   9.885531e-06,
     &   1.011579e-05,
     &   1.035142e-05,
     &   1.059254e-05,
     &   1.083927e-05,
     &   1.109175e-05,
     &   1.135011e-05,
     &   1.161449e-05,
     &   1.188502e-05,
     &   1.216186e-05,
     &   1.244515e-05,
     &   1.273503e-05,
     &   1.303167e-05,
     &   1.333521e-05,
     &   1.364583e-05,
     &   1.396368e-05,
     &   1.428894e-05,
     &   1.462177e-05,
     &   1.496236e-05,
     &   1.531087e-05,
     &   1.566751e-05,
     &   1.603245e-05,
     &   1.640590e-05,
     &   1.678804e-05,
     &   1.717908e-05,
     &   1.757924e-05,
     &   1.798871e-05,
     &   1.840772e-05,
     &   1.883649e-05,
     &   1.927525e-05,
     &   1.972423e-05,
     &   2.018366e-05,
     &   2.065380e-05,
     &   2.113489e-05,
     &   2.162719e-05,
     &   2.213095e-05,
     &   2.264644e-05,
     &   2.317395e-05,
     &   2.371374e-05,
     &   2.426610e-05,
     &   2.483133e-05,
     &   2.540973e-05,
     &   2.600160e-05,
     &   2.660725e-05,
     &   2.722701e-05,
     &   2.786121e-05,
     &   2.851018e-05,
     &   2.917427e-05,
     &   2.985383e-05,
     &   3.054921e-05,
     &   3.126079e-05,
     &   3.198895e-05,
     &   3.273407e-05,
     &   3.349654e-05,
     &   3.427678e-05,
     &   3.507519e-05,
     &   3.589219e-05,
     &   3.672823e-05,
     &   3.758374e-05,
     &   3.845918e-05,
     &   3.935501e-05,
     &   4.027170e-05,
     &   4.120975e-05,
     &   4.216965e-05,
     &   4.315191e-05,
     &   4.415704e-05,
     &   4.518559e-05,
     &   4.623810e-05,
     &   4.731513e-05,
     &   4.841724e-05,
     &   4.954502e-05,
     &   5.069907e-05,
     &   5.188000e-05,
     &   5.308844e-05,
     &   5.432503e-05,
     &   5.559043e-05,
     &   5.688529e-05,
     &   5.821032e-05,
     &   5.956621e-05,
     &   6.095369e-05,
     &   6.237348e-05,
     &   6.382635e-05,
     &   6.531306e-05,
     &   6.683439e-05,
     &   6.839116e-05,
     &   6.998420e-05,
     &   7.161434e-05,
     &   7.328245e-05,
     &   7.498942e-05,
     &   7.673615e-05,
     &   7.852356e-05,
     &   8.035261e-05,
     &   8.222426e-05,
     &   8.413951e-05,
     &   8.609938e-05,
     &   8.810489e-05,
     &   9.015711e-05,
     &   9.225714e-05,
     &   9.440609e-05,
     &   9.660509e-05,
     &   9.885531e-05,
     &   1.011579e-04,
     &   1.035142e-04,
     &   1.059254e-04,
     &   1.083927e-04,
     &   1.109175e-04,
     &   1.135011e-04,
     &   1.161449e-04,
     &   1.188502e-04,
     &   1.216186e-04,
     &   1.244515e-04,
     &   1.273503e-04,
     &   1.303167e-04,
     &   1.333521e-04,
     &   1.364583e-04,
     &   1.396368e-04,
     &   1.428894e-04,
     &   1.462177e-04,
     &   1.496236e-04,
     &   1.531087e-04,
     &   1.566751e-04,
     &   1.603245e-04,
     &   1.640590e-04,
     &   1.678804e-04,
     &   1.717908e-04,
     &   1.757924e-04,
     &   1.798871e-04,
     &   1.840772e-04,
     &   1.883649e-04,
     &   1.927525e-04,
     &   1.972423e-04,
     &   2.018366e-04,
     &   2.065380e-04,
     &   2.113489e-04,
     &   2.162719e-04,
     &   2.213095e-04,
     &   2.264644e-04,
     &   2.317395e-04,
     &   2.371374e-04,
     &   2.426610e-04,
     &   2.483133e-04,
     &   2.540973e-04,
     &   2.600160e-04,
     &   2.660725e-04,
     &   2.722701e-04,
     &   2.786121e-04,
     &   2.851018e-04,
     &   2.917427e-04,
     &   2.985383e-04,
     &   3.054921e-04,
     &   3.126079e-04,
     &   3.198895e-04,
     &   3.273407e-04,
     &   3.349654e-04,
     &   3.427678e-04,
     &   3.507519e-04,
     &   3.589219e-04,
     &   3.672823e-04,
     &   3.758374e-04,
     &   3.845918e-04,
     &   3.935501e-04,
     &   4.027170e-04,
     &   4.120975e-04,
     &   4.216965e-04,
     &   4.315191e-04,
     &   4.415704e-04,
     &   4.518559e-04,
     &   4.623810e-04,
     &   4.731513e-04,
     &   4.841724e-04,
     &   4.954502e-04,
     &   5.069907e-04,
     &   5.188000e-04,
     &   5.308844e-04,
     &   5.432503e-04,
     &   5.559043e-04,
     &   5.688529e-04,
     &   5.821032e-04,
     &   5.956621e-04,
     &   6.095369e-04,
     &   6.237348e-04,
     &   6.382635e-04,
     &   6.531306e-04,
     &   6.683439e-04,
     &   6.839116e-04,
     &   6.998420e-04,
     &   7.161434e-04,
     &   7.328245e-04,
     &   7.498942e-04,
     &   7.673615e-04,
     &   7.852356e-04,
     &   8.035261e-04,
     &   8.222426e-04,
     &   8.413951e-04,
     &   8.609938e-04,
     &   8.810489e-04,
     &   9.015711e-04,
     &   9.225714e-04,
     &   9.440609e-04,
     &   9.660509e-04,
     &   9.885531e-04,
     &   1.011579e-03,
     &   1.035142e-03,
     &   1.059254e-03,
     &   1.083927e-03,
     &   1.109175e-03,
     &   1.135011e-03,
     &   1.161449e-03,
     &   1.188502e-03,
     &   1.216186e-03,
     &   1.244515e-03,
     &   1.273503e-03,
     &   1.303167e-03,
     &   1.333521e-03,
     &   1.364583e-03,
     &   1.396368e-03,
     &   1.428894e-03,
     &   1.462177e-03,
     &   1.496236e-03,
     &   1.531087e-03,
     &   1.566751e-03,
     &   1.603245e-03,
     &   1.640590e-03,
     &   1.678804e-03,
     &   1.717908e-03,
     &   1.757924e-03,
     &   1.798871e-03,
     &   1.840772e-03,
     &   1.883649e-03,
     &   1.927525e-03,
     &   1.972423e-03,
     &   2.018366e-03,
     &   2.065380e-03,
     &   2.113489e-03,
     &   2.162719e-03,
     &   2.213095e-03,
     &   2.264644e-03,
     &   2.317395e-03,
     &   2.371374e-03,
     &   2.426610e-03,
     &   2.483133e-03,
     &   2.540973e-03,
     &   2.600160e-03,
     &   2.660725e-03,
     &   2.722701e-03,
     &   2.786121e-03,
     &   2.851018e-03,
     &   2.917427e-03,
     &   2.985383e-03,
     &   3.054921e-03,
     &   3.126079e-03,
     &   3.198895e-03,
     &   3.273407e-03,
     &   3.349654e-03,
     &   3.427678e-03,
     &   3.507519e-03,
     &   3.589219e-03,
     &   3.672823e-03,
     &   3.758374e-03,
     &   3.845918e-03,
     &   3.935501e-03,
     &   4.027170e-03,
     &   4.120975e-03,
     &   4.216965e-03,
     &   4.315191e-03,
     &   4.415704e-03,
     &   4.518559e-03,
     &   4.623810e-03,
     &   4.731513e-03,
     &   4.841724e-03,
     &   4.954502e-03,
     &   5.069907e-03,
     &   5.188000e-03,
     &   5.308844e-03,
     &   5.432503e-03,
     &   5.559043e-03,
     &   5.688529e-03,
     &   5.821032e-03,
     &   5.956621e-03,
     &   6.095369e-03,
     &   6.237348e-03,
     &   6.382635e-03,
     &   6.531306e-03,
     &   6.683439e-03,
     &   6.839116e-03,
     &   6.998420e-03,
     &   7.161434e-03,
     &   7.328245e-03,
     &   7.498942e-03,
     &   7.673615e-03,
     &   7.852356e-03,
     &   8.035261e-03,
     &   8.222426e-03,
     &   8.413951e-03,
     &   8.609938e-03,
     &   8.810489e-03,
     &   9.015711e-03,
     &   9.225714e-03,
     &   9.440609e-03,
     &   9.660509e-03,
     &   9.885531e-03,
     &   1.011579e-02,
     &   1.035142e-02,
     &   1.059254e-02,
     &   1.083927e-02,
     &   1.109175e-02,
     &   1.135011e-02,
     &   1.161449e-02,
     &   1.188502e-02,
     &   1.216186e-02,
     &   1.244515e-02,
     &   1.273503e-02,
     &   1.303167e-02,
     &   1.333521e-02,
     &   1.364583e-02,
     &   1.396368e-02,
     &   1.428894e-02,
     &   1.462177e-02,
     &   1.496236e-02,
     &   1.531087e-02,
     &   1.566751e-02,
     &   1.603245e-02,
     &   1.640590e-02,
     &   1.678804e-02,
     &   1.717908e-02,
     &   1.757924e-02,
     &   1.798871e-02,
     &   1.840772e-02,
     &   1.883649e-02,
     &   1.927525e-02,
     &   1.972423e-02,
     &   2.018366e-02,
     &   2.065380e-02,
     &   2.113489e-02,
     &   2.162719e-02,
     &   2.213095e-02,
     &   2.264644e-02,
     &   2.317395e-02,
     &   2.371374e-02,
     &   2.426610e-02,
     &   2.483133e-02,
     &   2.540973e-02,
     &   2.600160e-02,
     &   2.660725e-02,
     &   2.722701e-02,
     &   2.786121e-02,
     &   2.851018e-02,
     &   2.917427e-02,
     &   2.985383e-02,
     &   3.054921e-02,
     &   3.126079e-02,
     &   3.198895e-02,
     &   3.273407e-02,
     &   3.349654e-02,
     &   3.427678e-02,
     &   3.507519e-02,
     &   3.589219e-02,
     &   3.672823e-02,
     &   3.758374e-02,
     &   3.845918e-02,
     &   3.935501e-02,
     &   4.027170e-02,
     &   4.120975e-02,
     &   4.216965e-02,
     &   4.315191e-02,
     &   4.415704e-02,
     &   4.518559e-02,
     &   4.623810e-02,
     &   4.731513e-02,
     &   4.841724e-02,
     &   4.954502e-02,
     &   5.069907e-02,
     &   5.188000e-02,
     &   5.308844e-02,
     &   5.432503e-02,
     &   5.559043e-02,
     &   5.688529e-02,
     &   5.821032e-02,
     &   5.956621e-02,
     &   6.095369e-02,
     &   6.237348e-02,
     &   6.382635e-02,
     &   6.531306e-02,
     &   6.683439e-02,
     &   6.839116e-02,
     &   6.998420e-02,
     &   7.161434e-02,
     &   7.328245e-02,
     &   7.498942e-02,
     &   7.673615e-02,
     &   7.852356e-02,
     &   8.035261e-02,
     &   8.222426e-02,
     &   8.413951e-02,
     &   8.609938e-02,
     &   8.810489e-02,
     &   9.015711e-02,
     &   9.225714e-02,
     &   9.440609e-02,
     &   9.660509e-02,
     &   9.885531e-02,
     &   1.011579e-01,
     &   1.035142e-01,
     &   1.059254e-01,
     &   1.083927e-01,
     &   1.109175e-01,
     &   1.135011e-01,
     &   1.161449e-01,
     &   1.188502e-01,
     &   1.216186e-01,
     &   1.244515e-01,
     &   1.273503e-01,
     &   1.303167e-01,
     &   1.333521e-01,
     &   1.364583e-01,
     &   1.396368e-01,
     &   1.428894e-01,
     &   1.462177e-01,
     &   1.496236e-01,
     &   1.531087e-01,
     &   1.566751e-01,
     &   1.603245e-01,
     &   1.640590e-01,
     &   1.678804e-01,
     &   1.717908e-01,
     &   1.757924e-01,
     &   1.798871e-01,
     &   1.840772e-01,
     &   1.883649e-01,
     &   1.927525e-01,
     &   1.972423e-01,
     &   2.018366e-01,
     &   2.065380e-01,
     &   2.113489e-01,
     &   2.162719e-01,
     &   2.213095e-01,
     &   2.264644e-01,
     &   2.317395e-01,
     &   2.371374e-01,
     &   2.426610e-01,
     &   2.483133e-01,
     &   2.540973e-01,
     &   2.600160e-01,
     &   2.660725e-01,
     &   2.722701e-01,
     &   2.786121e-01,
     &   2.851018e-01,
     &   2.917427e-01,
     &   2.985383e-01,
     &   3.054921e-01,
     &   3.126079e-01,
     &   3.198895e-01,
     &   3.273407e-01,
     &   3.349654e-01,
     &   3.427678e-01,
     &   3.507519e-01,
     &   3.589219e-01,
     &   3.672823e-01,
     &   3.758374e-01,
     &   3.845918e-01,
     &   3.935501e-01,
     &   4.027170e-01,
     &   4.120975e-01,
     &   4.216965e-01,
     &   4.315191e-01,
     &   4.415704e-01,
     &   4.518559e-01,
     &   4.623810e-01,
     &   4.731513e-01,
     &   4.841724e-01,
     &   4.954502e-01,
     &   5.069907e-01,
     &   5.188000e-01,
     &   5.308844e-01,
     &   5.432503e-01,
     &   5.559043e-01,
     &   5.688529e-01,
     &   5.821032e-01,
     &   5.956621e-01,
     &   6.095369e-01,
     &   6.237348e-01,
     &   6.382635e-01,
     &   6.531306e-01,
     &   6.683439e-01,
     &   6.839116e-01,
     &   6.998420e-01,
     &   7.161434e-01,
     &   7.328245e-01,
     &   7.498942e-01,
     &   7.673615e-01,
     &   7.852356e-01,
     &   8.035261e-01,
     &   8.222426e-01,
     &   8.413951e-01,
     &   8.609938e-01,
     &   8.810489e-01,
     &   9.015711e-01,
     &   9.225714e-01,
     &   9.440609e-01,
     &   9.660509e-01,
     &   9.885531e-01,
     &   1.011579e+00,
     &   1.035142e+00,
     &   1.059254e+00,
     &   1.083927e+00,
     &   1.109175e+00,
     &   1.135011e+00,
     &   1.161449e+00,
     &   1.188502e+00,
     &   1.216186e+00,
     &   1.244515e+00,
     &   1.273503e+00,
     &   1.303167e+00,
     &   1.333521e+00,
     &   1.364583e+00,
     &   1.396368e+00,
     &   1.428894e+00,
     &   1.462177e+00,
     &   1.496236e+00,
     &   1.531087e+00,
     &   1.566751e+00,
     &   1.603245e+00,
     &   1.640590e+00,
     &   1.678804e+00,
     &   1.717908e+00,
     &   1.757924e+00,
     &   1.798871e+00,
     &   1.840772e+00,
     &   1.883649e+00,
     &   1.927525e+00,
     &   1.972423e+00,
     &   2.018366e+00,
     &   2.065380e+00,
     &   2.113489e+00,
     &   2.162719e+00,
     &   2.213095e+00,
     &   2.264644e+00,
     &   2.317395e+00,
     &   2.371374e+00,
     &   2.426610e+00,
     &   2.483133e+00,
     &   2.540973e+00,
     &   2.600160e+00,
     &   2.660725e+00,
     &   2.722701e+00,
     &   2.786121e+00,
     &   2.851018e+00,
     &   2.917427e+00,
     &   2.985383e+00,
     &   3.054921e+00,
     &   3.126079e+00,
     &   3.198895e+00,
     &   3.273407e+00,
     &   3.349654e+00,
     &   3.427678e+00,
     &   3.507519e+00,
     &   3.589219e+00,
     &   3.672823e+00,
     &   3.758374e+00,
     &   3.845918e+00,
     &   3.935501e+00,
     &   4.027170e+00,
     &   4.120975e+00,
     &   4.216965e+00,
     &   4.315191e+00,
     &   4.415704e+00,
     &   4.518559e+00,
     &   4.623810e+00,
     &   4.731513e+00,
     &   4.841724e+00,
     &   4.954502e+00,
     &   5.069907e+00,
     &   5.188000e+00,
     &   5.308844e+00,
     &   5.432503e+00,
     &   5.559043e+00,
     &   5.688529e+00,
     &   5.821032e+00,
     &   5.956621e+00,
     &   6.095369e+00,
     &   6.237348e+00,
     &   6.382635e+00,
     &   6.531306e+00,
     &   6.683439e+00,
     &   6.839116e+00,
     &   6.998420e+00,
     &   7.161434e+00,
     &   7.328245e+00,
     &   7.498942e+00,
     &   7.673615e+00,
     &   7.852356e+00,
     &   8.035261e+00,
     &   8.222426e+00,
     &   8.413951e+00,
     &   8.609938e+00,
     &   8.810489e+00,
     &   9.015711e+00,
     &   9.225714e+00,
     &   9.440609e+00,
     &   9.660509e+00,
     &   9.885531e+00,
     &   1.011579e+01,
     &   1.035142e+01,
     &   1.059254e+01,
     &   1.083927e+01,
     &   1.109175e+01,
     &   1.135011e+01,
     &   1.161449e+01,
     &   1.188502e+01,
     &   1.216186e+01,
     &   1.244515e+01,
     &   1.273503e+01,
     &   1.303167e+01,
     &   1.333521e+01,
     &   1.364583e+01,
     &   1.396368e+01,
     &   1.428894e+01,
     &   1.462177e+01,
     &   1.496236e+01,
     &   1.531087e+01,
     &   1.566751e+01,
     &   1.603245e+01,
     &   1.640590e+01,
     &   1.678804e+01,
     &   1.717908e+01,
     &   1.757924e+01,
     &   1.798871e+01,
     &   1.840772e+01,
     &   1.883649e+01,
     &   1.927525e+01,
     &   1.972423e+01,
     &   2.018366e+01,
     &   2.065380e+01,
     &   2.113489e+01,
     &   2.162719e+01,
     &   2.213095e+01,
     &   2.264644e+01,
     &   2.317395e+01,
     &   2.371374e+01,
     &   2.426610e+01,
     &   2.483133e+01,
     &   2.540973e+01,
     &   2.600160e+01,
     &   2.660725e+01,
     &   2.722701e+01,
     &   2.786121e+01,
     &   2.851018e+01,
     &   2.917427e+01,
     &   2.985383e+01,
     &   3.054921e+01,
     &   3.126079e+01,
     &   3.198895e+01,
     &   3.273407e+01,
     &   3.349654e+01,
     &   3.427678e+01,
     &   3.507519e+01,
     &   3.589219e+01,
     &   3.672823e+01,
     &   3.758374e+01,
     &   3.845918e+01,
     &   3.935501e+01,
     &   4.027170e+01,
     &   4.120975e+01,
     &   4.216965e+01,
     &   4.315191e+01,
     &   4.415704e+01,
     &   4.518559e+01,
     &   4.623810e+01,
     &   4.731513e+01,
     &   4.841724e+01,
     &   4.954502e+01,
     &   5.069907e+01,
     &   5.188000e+01,
     &   5.308844e+01,
     &   5.432503e+01,
     &   5.559043e+01,
     &   5.688529e+01,
     &   5.821032e+01,
     &   5.956621e+01,
     &   6.095369e+01,
     &   6.237348e+01,
     &   6.382635e+01,
     &   6.531306e+01,
     &   6.683439e+01,
     &   6.839116e+01,
     &   6.998420e+01,
     &   7.161434e+01,
     &   7.328245e+01,
     &   7.498942e+01,
     &   7.673615e+01,
     &   7.852356e+01,
     &   8.035261e+01,
     &   8.222426e+01,
     &   8.413951e+01,
     &   8.609938e+01,
     &   8.810489e+01,
     &   9.015711e+01,
     &   9.225714e+01,
     &   9.440609e+01,
     &   9.660509e+01,
     &   9.885531e+01,
     &   1.011579e+02,
     &   1.035142e+02,
     &   1.059254e+02,
     &   1.083927e+02,
     &   1.109175e+02,
     &   1.135011e+02,
     &   1.161449e+02,
     &   1.188502e+02,
     &   1.216186e+02,
     &   1.244515e+02,
     &   1.273503e+02,
     &   1.303167e+02,
     &   1.333521e+02,
     &   1.364583e+02,
     &   1.396368e+02,
     &   1.428894e+02,
     &   1.462177e+02,
     &   1.496236e+02,
     &   1.531087e+02,
     &   1.566751e+02,
     &   1.603245e+02,
     &   1.640590e+02,
     &   1.678804e+02,
     &   1.717908e+02,
     &   1.757924e+02,
     &   1.798871e+02,
     &   1.840772e+02,
     &   1.883649e+02,
     &   1.927525e+02,
     &   1.972423e+02,
     &   2.018366e+02,
     &   2.065380e+02,
     &   2.113489e+02,
     &   2.162719e+02,
     &   2.213095e+02,
     &   2.264644e+02,
     &   2.317395e+02,
     &   2.371374e+02,
     &   2.426610e+02,
     &   2.483133e+02,
     &   2.540973e+02,
     &   2.600160e+02,
     &   2.660725e+02,
     &   2.722701e+02,
     &   2.786121e+02,
     &   2.851018e+02,
     &   2.917427e+02,
     &   2.985383e+02,
     &   3.054921e+02,
     &   3.126079e+02,
     &   3.198895e+02,
     &   3.273407e+02,
     &   3.349654e+02,
     &   3.427678e+02,
     &   3.507519e+02,
     &   3.589219e+02,
     &   3.672823e+02,
     &   3.758374e+02,
     &   3.845918e+02,
     &   3.935501e+02,
     &   4.027170e+02,
     &   4.120975e+02,
     &   4.216965e+02,
     &   4.315191e+02,
     &   4.415704e+02,
     &   4.518559e+02,
     &   4.623810e+02,
     &   4.731513e+02,
     &   4.841724e+02,
     &   4.954502e+02,
     &   5.069907e+02,
     &   5.188000e+02,
     &   5.308844e+02,
     &   5.432503e+02,
     &   5.559043e+02,
     &   5.688529e+02,
     &   5.821032e+02,
     &   5.956621e+02,
     &   6.095369e+02,
     &   6.237348e+02,
     &   6.382635e+02,
     &   6.531306e+02,
     &   6.683439e+02,
     &   6.839116e+02,
     &   6.998420e+02,
     &   7.161434e+02,
     &   7.328245e+02,
     &   7.498942e+02,
     &   7.673615e+02,
     &   7.852356e+02,
     &   8.035261e+02,
     &   8.222426e+02,
     &   8.413951e+02,
     &   8.609938e+02,
     &   8.810489e+02,
     &   9.015711e+02,
     &   9.225714e+02,
     &   9.440609e+02,
     &   9.660509e+02 /
* Flux 
*----------------------------------------------------------------------*
      DATA FLUX/
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   0.000000e+00,
     &   4.185248e+05,
     &   8.346601e+05,
     &   8.591374e+05,
     &   8.821729e+05,
     &   9.098548e+05,
     &   9.224432e+05,
     &   9.414869e+05,
     &   9.601676e+05,
     &   9.847459e+05,
     &   1.014947e+06,
     &   1.042084e+06,
     &   1.061698e+06,
     &   1.082982e+06,
     &   1.100382e+06,
     &   1.133115e+06,
     &   1.152028e+06,
     &   1.170691e+06,
     &   1.190097e+06,
     &   1.212990e+06,
     &   1.226195e+06,
     &   1.251380e+06,
     &   1.263987e+06,
     &   1.285456e+06,
     &   1.297969e+06,
     &   1.314806e+06,
     &   1.331022e+06,
     &   1.343654e+06,
     &   1.349398e+06,
     &   1.354648e+06,
     &   1.365245e+06,
     &   1.374926e+06,
     &   1.377800e+06,
     &   1.380941e+06,
     &   1.380295e+06,
     &   1.391406e+06,
     &   1.384421e+06,
     &   1.388290e+06,
     &   1.382669e+06,
     &   1.364683e+06,
     &   1.362339e+06,
     &   1.349702e+06,
     &   1.337455e+06,
     &   1.328654e+06,
     &   1.318772e+06,
     &   1.313457e+06,
     &   1.287248e+06,
     &   1.272208e+06,
     &   1.248578e+06,
     &   1.224827e+06,
     &   1.199892e+06,
     &   1.175833e+06,
     &   1.159690e+06,
     &   1.123738e+06,
     &   1.095486e+06,
     &   1.067712e+06,
     &   1.041297e+06,
     &   1.012985e+06,
     &   9.804644e+05,
     &   9.578788e+05,
     &   9.244053e+05,
     &   8.950790e+05,
     &   8.613314e+05,
     &   8.354486e+05,
     &   8.038064e+05,
     &   7.746850e+05,
     &   7.393953e+05,
     &   7.141926e+05,
     &   6.821736e+05,
     &   6.600359e+05,
     &   6.251809e+05,
     &   6.045391e+05,
     &   5.725173e+05,
     &   5.490689e+05,
     &   5.213306e+05,
     &   4.982834e+05,
     &   4.774998e+05,
     &   4.591778e+05,
     &   4.333620e+05,
     &   4.234711e+05,
     &   3.979308e+05,
     &   3.869187e+05,
     &   3.667044e+05,
     &   3.544773e+05,
     &   3.379223e+05,
     &   3.260191e+05,
     &   3.124472e+05,
     &   3.043739e+05,
     &   2.962404e+05,
     &   2.857787e+05,
     &   2.695668e+05,
     &   2.689732e+05,
     &   2.613273e+05,
     &   2.560571e+05,
     &   2.521958e+05,
     &   2.463532e+05,
     &   2.410197e+05,
     &   2.399670e+05,
     &   2.357220e+05,
     &   2.354016e+05,
     &   2.326403e+05,
     &   2.298855e+05,
     &   2.295310e+05,
     &   2.282441e+05,
     &   2.277415e+05,
     &   2.292320e+05,
     &   2.266498e+05,
     &   2.259446e+05,
     &   2.202129e+05,
     &   2.255865e+05,
     &   2.226088e+05,
     &   2.244657e+05,
     &   2.210812e+05,
     &   2.239685e+05,
     &   2.210413e+05,
     &   2.242192e+05,
     &   2.225634e+05,
     &   2.178596e+05,
     &   2.227518e+05,
     &   2.205715e+05,
     &   2.249983e+05,
     &   2.217995e+05,
     &   2.206701e+05,
     &   2.213428e+05,
     &   2.220526e+05,
     &   2.177940e+05,
     &   2.241261e+05,
     &   2.186185e+05,
     &   2.137404e+05,
     &   2.219040e+05,
     &   2.217943e+05,
     &   2.154935e+05,
     &   2.181336e+05,
     &   2.128390e+05,
     &   2.152177e+05,
     &   2.137602e+05,
     &   2.092125e+05,
     &   2.161979e+05,
     &   2.119642e+05,
     &   2.137131e+05,
     &   2.099175e+05,
     &   2.124202e+05,
     &   2.181669e+05,
     &   2.136457e+05,
     &   2.142440e+05,
     &   2.116618e+05,
     &   2.100304e+05,
     &   2.088772e+05,
     &   2.115143e+05,
     &   2.061360e+05,
     &   2.060603e+05,
     &   2.122702e+05,
     &   2.122231e+05,
     &   2.068564e+05,
     &   2.106929e+05,
     &   2.092237e+05,
     &   2.103810e+05,
     &   2.091946e+05,
     &   2.103226e+05,
     &   2.090735e+05,
     &   2.073633e+05,
     &   2.109546e+05,
     &   2.098967e+05,
     &   2.125774e+05,
     &   2.053188e+05,
     &   2.035125e+05,
     &   1.949585e+05,
     &   1.989696e+05,
     &   2.149785e+05,
     &   2.207410e+05,
     &   2.128099e+05,
     &   2.196114e+05,
     &   2.350130e+05,
     &   2.218104e+05,
     &   2.129364e+05,
     &   2.135439e+05,
     &   2.142611e+05,
     &   2.071583e+05,
     &   2.121385e+05,
     &   2.105370e+05,
     &   2.168169e+05,
     &   2.160417e+05,
     &   2.172339e+05,
     &   2.123925e+05,
     &   2.166062e+05,
     &   2.036298e+05,
     &   2.088495e+05,
     &   2.219941e+05,
     &   2.107570e+05,
     &   2.013866e+05,
     &   2.083434e+05,
     &   2.007432e+05,
     &   1.923108e+05,
     &   1.975618e+05,
     &   2.077128e+05,
     &   2.075472e+05,
     &   2.029598e+05,
     &   2.062783e+05,
     &   2.001948e+05,
     &   2.028408e+05,
     &   2.073575e+05,
     &   1.842694e+05,
     &   2.086808e+05,
     &   2.176734e+05,
     &   2.028693e+05,
     &   2.019262e+05,
     &   2.053126e+05,
     &   2.080817e+05,
     &   2.011997e+05,
     &   2.174922e+05,
     &   1.908998e+05,
     &   2.001953e+05,
     &   2.150671e+05,
     &   2.156811e+05,
     &   2.128461e+05,
     &   2.434477e+05,
     &   2.922631e+05,
     &   3.184095e+05,
     &   2.302283e+05,
     &   2.146085e+05,
     &   1.980859e+05,
     &   1.981225e+05,
     &   2.152058e+05,
     &   1.977064e+05,
     &   1.958248e+05,
     &   2.023029e+05,
     &   2.168032e+05,
     &   1.890375e+05,
     &   1.891562e+05,
     &   1.898486e+05,
     &   2.046734e+05,
     &   2.053603e+05,
     &   1.968216e+05,
     &   2.071469e+05,
     &   1.991346e+05,
     &   2.172040e+05,
     &   2.049796e+05,
     &   2.081311e+05,
     &   2.155971e+05,
     &   2.657664e+05,
     &   2.670185e+05,
     &   2.451573e+05,
     &   2.273562e+05,
     &   2.108210e+05,
     &   2.080627e+05,
     &   2.171487e+05,
     &   2.241107e+05,
     &   1.770060e+05,
     &   2.061205e+05,
     &   1.978409e+05,
     &   1.992734e+05,
     &   1.990852e+05,
     &   2.138651e+05,
     &   2.102188e+05,
     &   2.189098e+05,
     &   2.220548e+05,
     &   2.153634e+05,
     &   2.145701e+05,
     &   2.255854e+05,
     &   2.173148e+05,
     &   2.075960e+05,
     &   2.152141e+05,
     &   2.128928e+05,
     &   2.197683e+05,
     &   2.236416e+05,
     &   2.011529e+05,
     &   2.151232e+05,
     &   2.248489e+05,
     &   2.249814e+05,
     &   2.318220e+05,
     &   2.227958e+05,
     &   2.080191e+05,
     &   2.124356e+05,
     &   2.089477e+05,
     &   2.290812e+05,
     &   2.422479e+05,
     &   2.320837e+05,
     &   2.096925e+05,
     &   2.350609e+05,
     &   2.218233e+05,
     &   2.175916e+05,
     &   2.026321e+05,
     &   2.327389e+05,
     &   2.190637e+05,
     &   2.147678e+05,
     &   2.114962e+05,
     &   2.234103e+05,
     &   2.207616e+05,
     &   2.470473e+05,
     &   2.573350e+05,
     &   2.289897e+05,
     &   2.535524e+05,
     &   2.286025e+05,
     &   2.266169e+05,
     &   2.309979e+05,
     &   2.169859e+05,
     &   2.148126e+05,
     &   2.162785e+05,
     &   2.125330e+05,
     &   2.085513e+05,
     &   2.082456e+05,
     &   2.271908e+05,
     &   2.339613e+05,
     &   2.301668e+05,
     &   2.602111e+05,
     &   2.353801e+05,
     &   2.213622e+05,
     &   2.170628e+05,
     &   2.150211e+05,
     &   2.522061e+05,
     &   2.159688e+05,
     &   2.255522e+05,
     &   2.110847e+05,
     &   2.241863e+05,
     &   2.692415e+05,
     &   2.560880e+05,
     &   2.197127e+05,
     &   2.177472e+05,
     &   2.197635e+05,
     &   2.355513e+05,
     &   2.308828e+05,
     &   2.474928e+05,
     &   2.123993e+05,
     &   1.982859e+05,
     &   2.435011e+05,
     &   2.174241e+05,
     &   2.173824e+05,
     &   2.337451e+05,
     &   2.205935e+05,
     &   2.148278e+05,
     &   2.223786e+05,
     &   2.534614e+05,
     &   2.216740e+05,
     &   2.147439e+05,
     &   2.211671e+05,
     &   2.364271e+05,
     &   2.327104e+05,
     &   2.441189e+05,
     &   2.222182e+05,
     &   2.590299e+05,
     &   2.287148e+05,
     &   2.287913e+05,
     &   2.238847e+05,
     &   2.310017e+05,
     &   2.189929e+05,
     &   2.228395e+05,
     &   2.438794e+05,
     &   2.362620e+05,
     &   2.301266e+05,
     &   2.307189e+05,
     &   2.194234e+05,
     &   2.255697e+05,
     &   2.270590e+05,
     &   2.331393e+05,
     &   2.188907e+05,
     &   2.370612e+05,
     &   2.260237e+05,
     &   2.406175e+05,
     &   2.328305e+05,
     &   2.303307e+05,
     &   2.276853e+05,
     &   2.338142e+05,
     &   2.163234e+05,
     &   2.418659e+05,
     &   2.270508e+05,
     &   2.275633e+05,
     &   2.273890e+05,
     &   2.338829e+05,
     &   2.068271e+05,
     &   2.061681e+05,
     &   2.217071e+05,
     &   2.180983e+05,
     &   2.348543e+05,
     &   2.367062e+05,
     &   2.345670e+05,
     &   2.216382e+05,
     &   2.486776e+05,
     &   2.443730e+05,
     &   2.481804e+05,
     &   2.382928e+05,
     &   2.276090e+05,
     &   2.491266e+05,
     &   2.261721e+05,
     &   2.309070e+05,
     &   2.291059e+05,
     &   2.296551e+05,
     &   2.310916e+05,
     &   2.292008e+05,
     &   2.299934e+05,
     &   2.245731e+05,
     &   2.373962e+05,
     &   2.212671e+05,
     &   2.288620e+05,
     &   2.409798e+05,
     &   2.268266e+05,
     &   2.289419e+05,
     &   2.240131e+05,
     &   2.255466e+05,
     &   2.180291e+05,
     &   2.166711e+05,
     &   2.291622e+05,
     &   2.173868e+05,
     &   2.293910e+05,
     &   2.260684e+05,
     &   2.271568e+05,
     &   2.374246e+05,
     &   2.128273e+05,
     &   2.239152e+05,
     &   2.050905e+05,
     &   1.753777e+05,
     &   1.405254e+05,
     &   1.471528e+05,
     &   1.816299e+05,
     &   1.965704e+05,
     &   2.323932e+05,
     &   2.147508e+05,
     &   2.295783e+05,
     &   2.315369e+05,
     &   2.245020e+05,
     &   2.430652e+05,
     &   2.252297e+05,
     &   2.272121e+05,
     &   2.347236e+05,
     &   2.429399e+05,
     &   2.270091e+05,
     &   2.629396e+05,
     &   2.420301e+05,
     &   2.320503e+05,
     &   2.393105e+05,
     &   2.311245e+05,
     &   2.397980e+05,
     &   2.379296e+05,
     &   2.507702e+05,
     &   2.403597e+05,
     &   2.344656e+05,
     &   2.364810e+05,
     &   2.286034e+05,
     &   2.388813e+05,
     &   2.336096e+05,
     &   2.467162e+05,
     &   2.611717e+05,
     &   2.242492e+05,
     &   2.386721e+05,
     &   2.414772e+05,
     &   2.492089e+05,
     &   2.411433e+05,
     &   2.265188e+05,
     &   2.466112e+05,
     &   2.272369e+05,
     &   2.157473e+05,
     &   2.376559e+05,
     &   2.501900e+05,
     &   2.692589e+05,
     &   2.495516e+05,
     &   2.500352e+05,
     &   2.522467e+05,
     &   2.385323e+05,
     &   2.693454e+05,
     &   2.413054e+05,
     &   2.506995e+05,
     &   2.364481e+05,
     &   2.318431e+05,
     &   2.154069e+05,
     &   2.411099e+05,
     &   2.452092e+05,
     &   2.662934e+05,
     &   2.351264e+05,
     &   2.207560e+05,
     &   2.546877e+05,
     &   2.410834e+05,
     &   2.472111e+05,
     &   2.317727e+05,
     &   2.658618e+05,
     &   2.469725e+05,
     &   2.327264e+05,
     &   2.354185e+05,
     &   2.624753e+05,
     &   2.392705e+05,
     &   2.530969e+05,
     &   2.364646e+05,
     &   2.485427e+05,
     &   2.537232e+05,
     &   2.441409e+05,
     &   2.359693e+05,
     &   2.416292e+05,
     &   2.436685e+05,
     &   2.337767e+05,
     &   2.446672e+05,
     &   2.542168e+05,
     &   2.220403e+05,
     &   2.348450e+05,
     &   2.542109e+05,
     &   2.515979e+05,
     &   2.586752e+05,
     &   2.335432e+05,
     &   1.938093e+05,
     &   2.483029e+05,
     &   2.250437e+05,
     &   2.635393e+05,
     &   2.372322e+05,
     &   2.092300e+05,
     &   2.113146e+05,
     &   2.664839e+05,
     &   2.388176e+05,
     &   2.363810e+05,
     &   2.336326e+05,
     &   2.577475e+05,
     &   2.382547e+05,
     &   2.367262e+05,
     &   2.940045e+05,
     &   1.903993e+05,
     &   2.361094e+05,
     &   2.534476e+05,
     &   2.612722e+05,
     &   2.343528e+05,
     &   2.140559e+05,
     &   2.422485e+05,
     &   2.564255e+05,
     &   2.206072e+05,
     &   2.534623e+05,
     &   2.689700e+05,
     &   2.609573e+05,
     &   2.553056e+05,
     &   2.646958e+05,
     &   2.546213e+05,
     &   2.453193e+05,
     &   2.338023e+05,
     &   2.429253e+05,
     &   2.727945e+05,
     &   2.482522e+05,
     &   2.510831e+05,
     &   2.204783e+05,
     &   2.676603e+05,
     &   2.375120e+05,
     &   2.476273e+05,
     &   1.899422e+05,
     &   2.626766e+05,
     &   2.478035e+05,
     &   2.618030e+05,
     &   2.453952e+05,
     &   2.285365e+05,
     &   2.415749e+05,
     &   2.902999e+05,
     &   2.442543e+05,
     &   2.647917e+05,
     &   2.484076e+05,
     &   2.238272e+05,
     &   2.795363e+05,
     &   2.768267e+05,
     &   2.323742e+05,
     &   2.659097e+05,
     &   2.713711e+05,
     &   2.684776e+05,
     &   2.746995e+05,
     &   2.713295e+05,
     &   2.947338e+05,
     &   2.657728e+05,
     &   3.053282e+05,
     &   2.588394e+05,
     &   2.617177e+05,
     &   2.635238e+05,
     &   2.389034e+05,
     &   2.621284e+05,
     &   2.859589e+05,
     &   2.589300e+05,
     &   2.621372e+05,
     &   2.640006e+05,
     &   2.806105e+05,
     &   2.425422e+05,
     &   2.680405e+05,
     &   2.877137e+05,
     &   3.056120e+05,
     &   2.607395e+05,
     &   2.602379e+05,
     &   2.780353e+05,
     &   2.952448e+05,
     &   2.563302e+05,
     &   2.368932e+05,
     &   2.506262e+05,
     &   2.536372e+05,
     &   2.888792e+05,
     &   2.774686e+05,
     &   2.614717e+05,
     &   3.162446e+05,
     &   3.135645e+05,
     &   2.685958e+05,
     &   2.929958e+05,
     &   2.939644e+05,
     &   3.302479e+05,
     &   2.552742e+05,
     &   2.659992e+05,
     &   2.941809e+05,
     &   3.619655e+05,
     &   2.757534e+05,
     &   3.032413e+05,
     &   3.170813e+05,
     &   3.074532e+05,
     &   3.269161e+05,
     &   3.737800e+05,
     &   3.490778e+05,
     &   3.055630e+05,
     &   2.884601e+05,
     &   2.290889e+05,
     &   2.865586e+05,
     &   2.773170e+05,
     &   2.960482e+05,
     &   3.384271e+05,
     &   3.650086e+05,
     &   3.120484e+05,
     &   2.903981e+05,
     &   2.774383e+05,
     &   2.139950e+05,
     &   2.142932e+05,
     &   2.243593e+05,
     &   2.254240e+05,
     &   2.555569e+05,
     &   3.047714e+05,
     &   3.136565e+05,
     &   3.541802e+05,
     &   3.223342e+05,
     &   3.236614e+05,
     &   3.267683e+05,
     &   3.723630e+05,
     &   3.065080e+05,
     &   3.545644e+05,
     &   3.632363e+05,
     &   3.274078e+05,
     &   3.772369e+05,
     &   3.413182e+05,
     &   3.953847e+05,
     &   3.944385e+05,
     &   4.136842e+05,
     &   4.324584e+05,
     &   3.927315e+05,
     &   3.835275e+05,
     &   3.912660e+05,
     &   3.884568e+05,
     &   4.196135e+05,
     &   4.449969e+05,
     &   4.151548e+05,
     &   4.301767e+05,
     &   4.295379e+05,
     &   3.925692e+05,
     &   4.566231e+05,
     &   4.286886e+05,
     &   4.461777e+05,
     &   4.459970e+05,
     &   4.233813e+05,
     &   3.653368e+05,
     &   3.794878e+05,
     &   3.340632e+05,
     &   3.265419e+05,
     &   3.609951e+05,
     &   3.480886e+05,
     &   2.939712e+05,
     &   4.143422e+05,
     &   4.140288e+05,
     &   4.213246e+05,
     &   4.723214e+05,
     &   4.572045e+05,
     &   5.204847e+05,
     &   4.978997e+05,
     &   5.039750e+05,
     &   5.157436e+05,
     &   5.777926e+05,
     &   5.296664e+05,
     &   4.625469e+05,
     &   6.395889e+05,
     &   6.089052e+05,
     &   5.568618e+05,
     &   6.638622e+05,
     &   5.788851e+05,
     &   5.666222e+05,
     &   5.087365e+05,
     &   4.803344e+05,
     &   5.228159e+05,
     &   5.397148e+05,
     &   5.382583e+05,
     &   5.215486e+05,
     &   5.540928e+05,
     &   6.573626e+05,
     &   6.605966e+05,
     &   6.748489e+05,
     &   6.743839e+05,
     &   7.069578e+05,
     &   6.754547e+05,
     &   6.789572e+05,
     &   7.374771e+05,
     &   7.373871e+05,
     &   6.290074e+05,
     &   7.333486e+05,
     &   6.856844e+05,
     &   7.223512e+05,
     &   7.251059e+05,
     &   7.774475e+05,
     &   7.778654e+05,
     &   7.877384e+05,
     &   9.043667e+05,
     &   8.987544e+05,
     &   7.959597e+05,
     &   9.743281e+05,
     &   9.715410e+05,
     &   9.152392e+05,
     &   7.975610e+05,
     &   8.976858e+05,
     &   9.533612e+05,
     &   9.711022e+05,
     &   1.010972e+06,
     &   1.036101e+06,
     &   9.987400e+05,
     &   1.113820e+06,
     &   1.087490e+06,
     &   1.110872e+06,
     &   1.127596e+06,
     &   8.981161e+05,
     &   1.027670e+06,
     &   1.130403e+06,
     &   1.065488e+06,
     &   1.071767e+06,
     &   9.765247e+05,
     &   8.697885e+05,
     &   8.555231e+05,
     &   8.499964e+05,
     &   1.076183e+06,
     &   1.234454e+06,
     &   1.315493e+06,
     &   1.535147e+06,
     &   1.651925e+06,
     &   1.816908e+06,
     &   1.522120e+06,
     &   1.413080e+06,
     &   1.719761e+06,
     &   1.785883e+06,
     &   1.641584e+06,
     &   1.664552e+06,
     &   1.817229e+06,
     &   1.811736e+06,
     &   1.796624e+06,
     &   1.761726e+06,
     &   1.896275e+06,
     &   1.993623e+06,
     &   1.879910e+06,
     &   1.959659e+06,
     &   1.841827e+06,
     &   1.881759e+06,
     &   2.168161e+06,
     &   1.944512e+06,
     &   1.880834e+06,
     &   1.818204e+06,
     &   1.654891e+06,
     &   1.818524e+06,
     &   2.010541e+06,
     &   1.969779e+06,
     &   1.975429e+06,
     &   1.882081e+06,
     &   1.760455e+06,
     &   1.663075e+06,
     &   1.432198e+06,
     &   1.461425e+06,
     &   1.670035e+06,
     &   1.858379e+06,
     &   2.078592e+06,
     &   2.168417e+06,
     &   2.153011e+06,
     &   1.934614e+06,
     &   2.271377e+06,
     &   2.289134e+06,
     &   2.424611e+06,
     &   2.062967e+06,
     &   1.915146e+06,
     &   1.810309e+06,
     &   2.290952e+06,
     &   2.199903e+06,
     &   2.251468e+06,
     &   2.107310e+06,
     &   2.135050e+06,
     &   2.376991e+06,
     &   2.337894e+06,
     &   2.580318e+06,
     &   2.023119e+06,
     &   1.991477e+06,
     &   2.116253e+06,
     &   1.868880e+06,
     &   2.600121e+06,
     &   2.075041e+06,
     &   2.315912e+06,
     &   2.026302e+06,
     &   2.010958e+06,
     &   2.053581e+06,
     &   2.046997e+06,
     &   1.948828e+06,
     &   2.104972e+06,
     &   2.005781e+06,
     &   1.783363e+06,
     &   2.122030e+06,
     &   2.205766e+06,
     &   1.934133e+06,
     &   1.936410e+06,
     &   1.645268e+06,
     &   1.558050e+06,
     &   1.396694e+06,
     &   1.559050e+06,
     &   1.596779e+06,
     &   1.133139e+06,
     &   1.456160e+06,
     &   1.486070e+06,
     &   1.327474e+06,
     &   1.269854e+06,
     &   1.032856e+06,
     &   1.227260e+06,
     &   2.351807e+06,
     &   1.066892e+06,
     &   9.383027e+05,
     &   8.629653e+05,
     &   8.720277e+05,
     &   9.101274e+05,
     &   8.858571e+05,
     &   8.559527e+05,
     &   8.543814e+05,
     &   7.972411e+05,
     &   8.233788e+05,
     &   7.422031e+05,
     &   7.755671e+05,
     &   7.027527e+05,
     &   7.036158e+05,
     &   6.794818e+05,
     &   6.965229e+05,
     &   6.624973e+05,
     &   6.562693e+05,
     &   6.764554e+05,
     &   6.910052e+05,
     &   6.184386e+05,
     &   6.242189e+05,
     &   7.138040e+05,
     &   7.294442e+05,
     &   6.739023e+05,
     &   5.108970e+05,
     &   5.897131e+05,
     &   5.898056e+05,
     &   5.790124e+05,
     &   6.109988e+05,
     &   5.879024e+05,
     &   5.471592e+05,
     &   5.027038e+05,
     &   5.258947e+05,
     &   6.086648e+05,
     &   5.608282e+05,
     &   5.660968e+05,
     &   4.834525e+05,
     &   4.624152e+05,
     &   4.635418e+05,
     &   4.876818e+05,
     &   5.052717e+05,
     &   4.759797e+05,
     &   4.652638e+05,
     &   4.290187e+05,
     &   4.640515e+05,
     &   5.038368e+05,
     &   4.735362e+05,
     &   4.590212e+05,
     &   4.246633e+05,
     &   4.153149e+05,
     &   3.946339e+05,
     &   3.923516e+05,
     &   4.384874e+05,
     &   4.203683e+05,
     &   3.956254e+05,
     &   3.572838e+05,
     &   3.751190e+05,
     &   3.703819e+05,
     &   3.677620e+05,
     &   3.568768e+05,
     &   3.418222e+05,
     &   3.674442e+05,
     &   3.965760e+05,
     &   3.903662e+05,
     &   3.422807e+05,
     &   3.390420e+05,
     &   3.367579e+05,
     &   3.409624e+05,
     &   4.105005e+05,
     &   3.817958e+05,
     &   3.584572e+05,
     &   3.469658e+05,
     &   3.793839e+05,
     &   3.670681e+05,
     &   3.467344e+05,
     &   3.831094e+05,
     &   3.887991e+05,
     &   3.573912e+05,
     &   3.574066e+05,
     &   3.831616e+05,
     &   3.584989e+05,
     &   3.585900e+05,
     &   3.730100e+05,
     &   3.345992e+05,
     &   3.399287e+05,
     &   3.347308e+05,
     &   3.642304e+05,
     &   3.769076e+05,
     &   3.467138e+05,
     &   3.991438e+05,
     &   3.602398e+05,
     &   4.133876e+05,
     &   4.349365e+05,
     &   4.182949e+05,
     &   3.694419e+05,
     &   4.067498e+05,
     &   3.658472e+05,
     &   3.435403e+05,
     &   3.474914e+05,
     &   3.737446e+05,
     &   4.263032e+05,
     &   3.682348e+05,
     &   3.925871e+05,
     &   4.633817e+05,
     &   4.629837e+05,
     &   4.309462e+05,
     &   4.443698e+05,
     &   4.257695e+05,
     &   4.525926e+05,
     &   4.474688e+05,
     &   3.903425e+05,
     &   3.720991e+05,
     &   4.218082e+05,
     &   4.392820e+05,
     &   3.803293e+05,
     &   3.917582e+05,
     &   3.805586e+05,
     &   3.958643e+05,
     &   3.954733e+05,
     &   3.697002e+05,
     &   4.344884e+05,
     &   3.963692e+05,
     &   4.172907e+05,
     &   4.259444e+05,
     &   4.575718e+05,
     &   4.323257e+05,
     &   4.771712e+05,
     &   4.683998e+05,
     &   4.249831e+05,
     &   4.431483e+05,
     &   5.155095e+05,
     &   5.007351e+05,
     &   4.371650e+05,
     &   4.839679e+05,
     &   4.365660e+05,
     &   4.526877e+05,
     &   4.208771e+05,
     &   4.737950e+05,
     &   4.851971e+05,
     &   4.500327e+05,
     &   3.944374e+05,
     &   4.332980e+05,
     &   4.864081e+05,
     &   4.669513e+05,
     &   4.613579e+05,
     &   4.395830e+05,
     &   4.670578e+05,
     &   4.768357e+05,
     &   4.820442e+05,
     &   4.351215e+05,
     &   4.636038e+05,
     &   4.302619e+05,
     &   4.470690e+05,
     &   4.555834e+05,
     &   4.686113e+05,
     &   4.792960e+05,
     &   4.740601e+05,
     &   4.788443e+05,
     &   4.324882e+05,
     &   4.077432e+05,
     &   4.108872e+05,
     &   4.367994e+05,
     &   4.613690e+05,
     &   4.325631e+05,
     &   4.191638e+05,
     &   4.163018e+05,
     &   3.733343e+05,
     &   3.974018e+05,
     &   4.037606e+05,
     &   3.455294e+05,
     &   3.925618e+05,
     &   3.817015e+05,
     &   3.795653e+05,
     &   3.811198e+05,
     &   3.128170e+05,
     &   3.526070e+05,
     &   3.898050e+05,
     &   3.520024e+05,
     &   3.731558e+05,
     &   3.495723e+05,
     &   3.326037e+05,
     &   3.298362e+05,
     &   3.442338e+05,
     &   3.218719e+05,
     &   3.294875e+05,
     &   2.869419e+05,
     &   2.925400e+05,
     &   3.200465e+05,
     &   2.766477e+05,
     &   2.876634e+05,
     &   2.804648e+05,
     &   3.198658e+05,
     &   2.585200e+05,
     &   3.036476e+05,
     &   3.038575e+05,
     &   2.499644e+05,
     &   2.669444e+05,
     &   2.508986e+05,
     &   2.397752e+05,
     &   2.237389e+05,
     &   2.279498e+05,
     &   2.182148e+05,
     &   2.068307e+05,
     &   2.035699e+05,
     &   2.069498e+05,
     &   1.962502e+05,
     &   2.124918e+05,
     &   1.856580e+05,
     &   1.822916e+05,
     &   1.848643e+05,
     &   1.547899e+05,
     &   1.492604e+05,
     &   1.427259e+05,
     &   1.402433e+05,
     &   1.422935e+05,
     &   1.145369e+05,
     &   1.108569e+05,
     &   1.041398e+05,
     &   8.401086e+04,
     &   9.144321e+04,
     &   7.859905e+04,
     &   6.756712e+04,
     &   6.229154e+04,
     &   5.509066e+04,
     &   5.502643e+04,
     &   4.999515e+04,
     &   4.402821e+04,
     &   3.941741e+04,
     &   3.906271e+04,
     &   3.631787e+04,
     &   3.005092e+04,
     &   2.911373e+04,
     &   2.375045e+04,
     &   2.224273e+04,
     &   2.551665e+04,
     &   1.685478e+04,
     &   1.647112e+04,
     &   1.631669e+04,
     &   1.228216e+04,
     &   1.368716e+04,
     &   1.378717e+04,
     &   1.005472e+04,
     &   7.664153e+03,
     &   8.449919e+03,
     &   5.442496e+03  /
*======================================================================*
*                                                                      *
*                 BASIC VERSION                                        *
*                                                                      *
*======================================================================*
      NOMORE = 0
*  +-------------------------------------------------------------------*
*  |  First call initializations:
      IF ( LFIRST ) THEN
*  |  *** The following 3 cards are mandatory ***
         TKESUM = ZERZER
         LFIRST = .FALSE.
         LUSSRC = .TRUE.
* Calculate the cumulative function 
         TEGR(0)=FLUX(0)
	 DO I = 1, 1099
	      TEGR(I) = TEGR(I-1) + FLUX(I)
c If you want to use a simple trapezium integratio, decomment the
c following two lines	      
c	      TEGR(I) = TEGR(I-1) + HLFHLF*(FLUX(I)+FLUX(I-1)) *
c     &   (I-(I-1))
	 END DO	 
* Normalize the cumulative function to 1
          SOMMA = TEGR(1099)
	  DO I =0, 1099
	      TEGR(I)=TEGR(I)/SOMMA
	  END DO
*  |  *** User initialization ***
      END IF
*  |
*  +-------------------------------------------------------------------*
*  Push one source particle to the stack. Note that you could as well
*  push many but this way we reserve a maximum amount of space in the
*  stack for the secondaries to be generated
* Npflka is the stack counter: of course any time source is called it
* must be =0
      NPFLKA = NPFLKA + 1
* Wt is the weight of the particle
      WTFLK  (NPFLKA) = ONEONE
      WEIPRI = WEIPRI + WTFLK (NPFLKA)
* Particle type (1=proton.....). Ijbeam is the type set by the BEAM
* card
*  +-------------------------------------------------------------------*
*  Insert here the excitation energy of the compound nucleus
*  The excitation energy U=E(235U)+E(n)+T(n)-E(236)
*  with E(235U) = mass_energy of 235U
*       E(n) = mass_energy of neutron
*       T(n) = kinietic energy of the incident spallation neutron 
*       E(236U) = mass_energy of 236U
*  Keep in mind that E(235U)+E(n)-E(236U) = 6.5MeV
*
*
*
*  The kinetic energy of the spallation neutron si sampled from the 
*  neutron flux spectrum at the n_TOF facility
*
*
*
      XYZ= FLRNDM(XYZW)
      DO 10 I=0, 1098
           IF (TEGR(I) .GE. XYZ) THEN
	       TNEUTENE = ENEDGE(I) + (ENEDGE(I+1) -ENEDGE(I))*
     &  (XYZ-TEGR(I))/(TEGR(I+1) - TEGR(I))
               GO TO 20
	   END IF 
 10   CONTINUE
 20   CONTINUE
C DA MODIFICARE PER RENDERE PIU' CONTINUA L' ENERGIA!!!
      EMIN = 6.5
      U = EMIN + TNEUTENE
      IF (U .LE. 16.25) THEN 
         WA = EXP(0.538*U-9.9564)
      END IF
      IF (U .GT. 16.25) THEN
         WA = EXP(0.09197*U-2.7003)
      END IF
*
*
*
*
*
*
*
*
*
*
*  +-------------------------------------------------------------------*
*  |  (Radioactive) isotope:
      IF ( IJBEAM .EQ. -2 .AND. LRDBEA ) THEN
         IARES  = IPROA
         IZRES  = IPROZ
         IISRES = IPROM
         CALL STISBM ( IARES, IZRES, IISRES )
         IJHION = IPROZ  * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID  = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
*  |
*  +-------------------------------------------------------------------*
*  |  Heavy ion:
      ELSE IF ( IJBEAM .EQ. -2 ) THEN
C         ASYM(0) = ZERZER 
C	 ANSYM(0) = ZERZER
C	 ADIST(0) = ZERZER
	 CUMUL(0) = ZERZER
	 SIGMA = EXP(0.00553*U+2.1386)
	 SIGMA1 = 0.5 *5.6 
	 SIGMA2 = 5.6 
	 A1 = 141.0
	 A2 = 92.0
	 AINI= 235.0
	 AMEZ=235./2
	 AMEZZI =INT(AMEZ)
	 AHS1=(A1+A2)/2
	 AHS=INT(AH1S)
*	 AMEZZI=235./2
*	 AHS=(A1+A2)/2
* Take A/2=120 FOR 241Am
         DO I=0, 235
	    ASYM(I) = EXP ( -((I-AMEZZI)**2)/(2*(SIGMA**2)))
	    ANSYM(I) = EXP(-((I-A2)**2)/(2*SIGMA2**2)) +
     &   EXP(-((I-(AINI-A2))**2)/(2*SIGMA2**2)) + HLFHLF *
     &  (EXP(-((I-A1)**2)/(2*SIGMA1**2)) +EXP(-((I-(AINI-A1))
     &   **2)/(2*SIGMA1**2)) )	    
*	    WRITE (lnuout,*) ASYM(I), ANSYM(I)   
	 END DO        
	 W = ((WA - ANSYM(AMEZZI))/(1-WA*(ASYM(AHS))))
         DO I=0,235
	     ADIST(I)= ANSYM(I) + W*ANSYM(I)
	 END DO    
* Calculate the cumulative function (simple trapezium integration)	 
	 DO I = 1, 235
	     CUMUL(I) = CUMUL(I-1) + HLFHLF*(ADIST(I)+ADIST(I-1))
	 END DO	 
	 SUM = CUMUL(235)
	 DO I = 1, 235
	    CUMUL(I)=CUMUL(I)/SUM
	 END DO
* Sample from the mass distribution
         XI = FLRNDM(DUMMY)
*	 WRITE(lnuout,*) XI
	 DO 1 I=1, 234
	    IF (CUMUL(I) .GE. XI) THEN
	        SAMPLE = I+((XI-CUMUL(I))/(CUMUL(I+1)-CUMUL(I)))
		GO TO 2
	    END IF
  1      CONTINUE
  2      CONTINUE
         MASS = INT(SAMPLE)
         IPROA = MASS	 
	 IF (MASS .GE. 134) THEN
	     DELTAZ= -0.45
	 END IF
	 IF ((MASS .GT. 101).AND.(MASS .LT. 134)) THEN
	    DELTAZ = -0.45*(MASS -117)/(134-117)
	 END IF
	 IF (MASS .LE. 101) THEN
	    DELTAZ = 0.45
	 END IF
	 ZPEAK = MASS * 92/235 + DELTAZ
	 CALL FLNRRN(RGAUSS)
*         IF (MASS.LE.117) THEN
*	 ZPEAK=38.0
*	 Z = ZPEAK -2*RGAUSS
*	 END IF
*	 IF (MASS.GT.117) THEN
*	 ZPEAK= 55.0
*	 Z = ZPEAK -1.66*RGAUSS
*	 END IF 
	 CALL FLNRRN(RGAUSS)
	 Z = ZPEAK - 0.6*RGAUSS
	 Z = INT(Z)	 
	 IPROZ = Z
         IJHION = IPROZ  * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID  = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
         ILOFLK (NPFLKA) = IJHION
*  |  Flag this is prompt radiation
         LRADDC (NPFLKA) = .FALSE.
*  |
*  +-------------------------------------------------------------------*
*  |  Normal hadron:
      ELSE
         IONID = IJBEAM
         ILOFLK (NPFLKA) = IJBEAM
*  |  Flag this is prompt radiation
         LRADDC (NPFLKA) = .FALSE.
      END IF
*  |
*  +-------------------------------------------------------------------*
* From this point .....
* Particle generation (1 for primaries)
      LOFLK  (NPFLKA) = 1
* User dependent flag:
      LOUSE  (NPFLKA) = 0
* User dependent spare variables:
      DO 100 ISPR = 1, MKBMX1
         SPAREK (ISPR,NPFLKA) = ZERZER
 100  CONTINUE
* User dependent spare flags:
      DO 200 ISPR = 1, MKBMX2
         ISPARK (ISPR,NPFLKA) = 0
 200  CONTINUE
* Save the track number of the stack particle:
      ISPARK (MKBMX2,NPFLKA) = NPFLKA
      NPARMA = NPARMA + 1
      NUMPAR (NPFLKA) = NPARMA
      NEVENT (NPFLKA) = 0
      DFNEAR (NPFLKA) = +ZERZER
* ... to this point: don't change anything
* Particle age (s)
      AGESTK (NPFLKA) = +ZERZER
      AKNSHR (NPFLKA) = -TWOTWO
* Group number for "low" energy neutrons, set to 0 anyway
      IGROUP (NPFLKA) = 0
*------------------------------------------------------------*
*                                                            *
*       Kinetic energy of the FFs (GeV/nucleons)             *
*                                                            *
*------------------------------------------------------------*
*
*
* The total kinetic energy of the fragments is about 80% * (Q +T)
* You have to take into account that fission in this case is
* induced by neutrons with energy up to 1 GeV!!!
* scrivi anche il Q della reazione...
*      TMAX = TNEUTENE + 200.0
      TMAX = TNEUTENE *0.8+ 167.37
*      TMAX = 167.37
      TTOT=TMAX*(AINI-MASS)/AINI 
      CALL FLNRRN(RGAUSS) 
      IF (MASS.LE.117) THEN 
          ENE = TTOT - RGAUSS* 5.33        
      END IF
      IF (MASS.GT.117) THEN 
          ENE = TTOT - RGAUSS* 6.66        
      END IF      
*      DELTAENE = 20/10000.0
*      DELTAGAUSS =6.0/10000.0
*      CALL FLNRRN(RGAUSS)
*      WRITE(lnuout,*) RGAUSS
*      REAL TTR=RGAUSS      
*      DO 8 I = 0, 9999
*           IF (TTR .LE. (-3.0 + I*DELTAGAUSS)) THEN
*	       ENE = (TTOT-30) + I * DELTAENE
*	       GO TO 9
*	  END IF
*8     CONTINUE
*9     CONTINUE
      ENE = ENE *0.001
      TKEFLK (NPFLKA) = ENE
      WRITE (LNUOUT,*) MASS,'    ',Z,'    ',TNEUTENE,'           ',ENE
* Particle momentum
*      PMOFLK (NPFLKA) = PBEAM
      PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
     &                       + TWOTWO * AM (IONID) ) )
*
*
*
* Cosines (tx,ty,tz)
* You want a 4PI isotropi angular distribution
       COSTH = ONEONE -2.D0 * FLRNDM (RRRR)
       IF (COSTH.GT.ONEONE) COSTH = ONEONE
       IF (COSTH.LT.-ONEONE) COSTH = -ONEONE
       ANG = TWOPIP * FLRNDM (SSSS)
       UBEAM=SQRT( ONEONE - COSTH * COSTH) * COS(ANG)
       VBEAM=SQRT( ONEONE - COSTH * COSTH) * SIN(ANG)
       WBEAM=COSTH              
      TXFLK  (NPFLKA) = UBEAM
      TYFLK  (NPFLKA) = VBEAM
      TZFLK  (NPFLKA) = WBEAM
*     TZFLK  (NPFLKA) = SQRT ( ONEONE - TXFLK (NPFLKA)**2
*    &                       - TYFLK (NPFLKA)**2 )
* Polarization cosines:
      TXPOL  (NPFLKA) = -TWOTWO
      TYPOL  (NPFLKA) = +ZERZER
      TZPOL  (NPFLKA) = +ZERZER
* Particle coordinates
 99   X1 = 0.D0
      X2 = 4.D0   
      SEGNOX = FLRNDM (SSSS)
      IF (SEGNOX.LT.0.5)   XFLK (NPFLKA) = -( X2 - X1 ) * FLRNDM(XXXX)
      IF (SEGNOX.GE.0.5)   XFLK (NPFLKA) = +( X2 - X1 ) * FLRNDM(XXXX)
      Y1 = 0.D0
      Y2 = 4.D0   
      SEGNOY = FLRNDM (SSZZ)
      IF (SEGNOY.LT.0.5)   YFLK (NPFLKA) = -( Y2 - Y1 ) * FLRNDM(XXXX)
      IF (SEGNOY.GE.0.5)   YFLK (NPFLKA) = +( Y2 - Y1 ) * FLRNDM(XXXX)
      CHECK = YFLK (NPFLKA)**2 + XFLK (NPFLKA)**2
      IF (CHECK.GT.16.0) GO TO 99
      Z1 = 0.D0
      Z2 = 4.2633E-5
      ZFLK (NPFLKA) = ( Z2 -Z1 ) * FLRNDM(ZZZZ)
*  Calculate the total kinetic energy of the primaries: don't change
      IF ( ILOFLK (NPFLKA) .EQ. -2 .OR. ILOFLK (NPFLKA) .GT. 100000 )
     &   THEN
         TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
      ELSE IF ( ILOFLK (NPFLKA) .NE. 0 ) THEN
         TKESUM = TKESUM + ( TKEFLK (NPFLKA) + AMDISC (ILOFLK(NPFLKA)) )
     &          * WTFLK (NPFLKA)
      ELSE
         TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
      END IF
      RADDLY (NPFLKA) = ZERZER
*  Here we ask for the region number of the hitting point.
*     NREG (NPFLKA) = ...
*  The following line makes the starting region search much more
*  robust if particles are starting very close to a boundary:
      CALL GEOCRS ( TXFLK (NPFLKA), TYFLK (NPFLKA), TZFLK (NPFLKA) )
      CALL GEOREG ( XFLK  (NPFLKA), YFLK  (NPFLKA), ZFLK  (NPFLKA),
     &              NRGFLK(NPFLKA), IDISC )
*  Do not change these cards:
      CALL GEOHSM ( NHSPNT (NPFLKA), 1, -11, MLATTC )
      NLATTC (NPFLKA) = MLATTC
      CMPATH (NPFLKA) = ZERZER
      CALL SOEVSV
      RETURN
*=== End of subroutine Source =========================================*
      END

