*$ 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