*$ CREATE SOURCE.FOR *COPY SOURCE * *=== source ===========================================================* SUBROUTINE SOURCE ( NOMORE ) INCLUDE '(DBLPRC)' INCLUDE '(DIMPAR)' INCLUDE '(IOUNIT)' * * * Used DIMENSION ENEDGE (0:799) DIMENSION FLUX (0:799) DIMENSION TEGR (0:799) DIMENSION RATIO (0:799) * Not used DIMENSION CSEDGE (0:4000) DIMENSION PANGLE (0:4000) DIMENSION TRIAL (0:4000) * * *----------------------------------------------------------------------* * INCLUDE '(BEAMCM)' INCLUDE '(FHEAVY)' INCLUDE '(FLKSTK)' INCLUDE '(IOIOCM)' INCLUDE '(LTCLCM)' INCLUDE '(PAPROP)' INCLUDE '(SOURCM)' INCLUDE '(SUMCOU)' INCLUDE '(CASLIM)' * LOGICAL LFIRST * SAVE LFIRST DATA LFIRST / .TRUE. / * Neutron energy group boundaries in MeV DATA ENEDGE / & 1.011646e-08, & 1.035211e-08, & 1.059324e-08, & 1.083999e-08, & 1.109248e-08, & 1.135086e-08, & 1.161526e-08, & 1.188581e-08, & 1.216267e-08, & 1.244597e-08, & 1.273587e-08, & 1.303253e-08, & 1.333610e-08, & 1.364674e-08, & 1.396461e-08, & 1.428989e-08, & 1.462274e-08, & 1.496335e-08, & 1.531189e-08, & 1.566855e-08, & 1.603352e-08, & 1.640698e-08, & 1.678915e-08, & 1.718022e-08, & 1.758040e-08, & 1.798990e-08, & 1.840894e-08, & 1.883774e-08, & 1.927653e-08, & 1.972553e-08, & 2.018500e-08, & 2.065517e-08, & 2.113629e-08, & 2.162862e-08, & 2.213241e-08, & 2.264794e-08, & 2.317548e-08, & 2.371531e-08, & 2.426771e-08, & 2.483298e-08, & 2.541141e-08, & 2.600332e-08, & 2.660901e-08, & 2.722882e-08, & 2.786306e-08, & 2.851207e-08, & 2.917620e-08, & 2.985580e-08, & 3.055124e-08, & 3.126286e-08, & 3.199107e-08, & 3.273624e-08, & 3.349876e-08, & 3.427905e-08, & 3.507751e-08, & 3.589457e-08, & 3.673066e-08, & 3.758623e-08, & 3.846173e-08, & 3.935762e-08, & 4.027437e-08, & 4.121248e-08, & 4.217244e-08, & 4.315477e-08, & 4.415997e-08, & 4.518859e-08, & 4.624117e-08, & 4.731826e-08, & 4.842045e-08, & 4.954830e-08, & 5.070243e-08, & 5.188344e-08, & 5.309196e-08, & 5.432864e-08, & 5.559411e-08, & 5.688906e-08, & 5.821418e-08, & 5.957016e-08, & 6.095773e-08, & 6.237762e-08, & 6.383058e-08, & 6.531738e-08, & 6.683882e-08, & 6.839570e-08, & 6.998884e-08, & 7.161909e-08, & 7.328731e-08, & 7.499439e-08, & 7.674123e-08, & 7.852877e-08, & 8.035794e-08, & 8.222971e-08, & 8.414509e-08, & 8.610508e-08, & 8.811073e-08, & 9.016309e-08, & 9.226326e-08, & 9.441235e-08, & 9.661149e-08, & 9.886186e-08, & 1.011647e-07, & 1.035211e-07, & 1.059324e-07, & 1.083999e-07, & 1.109248e-07, & 1.135086e-07, & 1.161526e-07, & 1.188581e-07, & 1.216267e-07, & 1.244597e-07, & 1.273587e-07, & 1.303253e-07, & 1.333610e-07, & 1.364674e-07, & 1.396461e-07, & 1.428989e-07, & 1.462274e-07, & 1.496335e-07, & 1.531189e-07, & 1.566855e-07, & 1.603352e-07, & 1.640698e-07, & 1.678915e-07, & 1.718022e-07, & 1.758040e-07, & 1.798990e-07, & 1.840894e-07, & 1.883774e-07, & 1.927653e-07, & 1.972554e-07, & 2.018500e-07, & 2.065517e-07, & 2.113629e-07, & 2.162862e-07, & 2.213241e-07, & 2.264794e-07, & 2.317548e-07, & 2.371531e-07, & 2.426771e-07, & 2.483298e-07, & 2.541141e-07, & 2.600332e-07, & 2.660901e-07, & 2.722882e-07, & 2.786306e-07, & 2.851207e-07, & 2.917620e-07, & 2.985580e-07, & 3.055124e-07, & 3.126287e-07, & 3.199107e-07, & 3.273624e-07, & 3.349876e-07, & 3.427905e-07, & 3.507751e-07, & 3.589457e-07, & 3.673066e-07, & 3.758623e-07, & 3.846173e-07, & 3.935762e-07, & 4.027437e-07, & 4.121248e-07, & 4.217244e-07, & 4.315477e-07, & 4.415997e-07, & 4.518859e-07, & 4.624117e-07, & 4.731826e-07, & 4.842045e-07, & 4.954830e-07, & 5.070243e-07, & 5.188344e-07, & 5.309196e-07, & 5.432863e-07, & 5.559411e-07, & 5.688906e-07, & 5.821418e-07, & 5.957016e-07, & 6.095773e-07, & 6.237762e-07, & 6.383058e-07, & 6.531738e-07, & 6.683882e-07, & 6.839570e-07, & 6.998884e-07, & 7.161909e-07, & 7.328731e-07, & 7.499439e-07, & 7.674124e-07, & 7.852877e-07, & 8.035794e-07, & 8.222971e-07, & 8.414509e-07, & 8.610508e-07, & 8.811072e-07, & 9.016309e-07, & 9.226326e-07, & 9.441234e-07, & 9.661149e-07, & 9.886186e-07, & 1.011647e-06, & 1.035211e-06, & 1.059324e-06, & 1.083999e-06, & 1.109248e-06, & 1.135086e-06, & 1.161526e-06, & 1.188581e-06, & 1.216267e-06, & 1.244597e-06, & 1.273587e-06, & 1.303253e-06, & 1.333610e-06, & 1.364674e-06, & 1.396461e-06, & 1.428989e-06, & 1.462274e-06, & 1.496335e-06, & 1.531189e-06, & 1.566855e-06, & 1.603352e-06, & 1.640698e-06, & 1.678915e-06, & 1.718022e-06, & 1.758040e-06, & 1.798990e-06, & 1.840894e-06, & 1.883774e-06, & 1.927653e-06, & 1.972553e-06, & 2.018500e-06, & 2.065517e-06, & 2.113629e-06, & 2.162862e-06, & 2.213241e-06, & 2.264794e-06, & 2.317548e-06, & 2.371531e-06, & 2.426771e-06, & 2.483298e-06, & 2.541141e-06, & 2.600332e-06, & 2.660901e-06, & 2.722882e-06, & 2.786306e-06, & 2.851207e-06, & 2.917620e-06, & 2.985581e-06, & 3.055124e-06, & 3.126287e-06, & 3.199107e-06, & 3.273624e-06, & 3.349876e-06, & 3.427905e-06, & 3.507751e-06, & 3.589457e-06, & 3.673066e-06, & 3.758623e-06, & 3.846173e-06, & 3.935762e-06, & 4.027437e-06, & 4.121248e-06, & 4.217245e-06, & 4.315477e-06, & 4.415997e-06, & 4.518859e-06, & 4.624117e-06, & 4.731826e-06, & 4.842045e-06, & 4.954830e-06, & 5.070243e-06, & 5.188344e-06, & 5.309196e-06, & 5.432863e-06, & 5.559411e-06, & 5.688906e-06, & 5.821418e-06, & 5.957016e-06, & 6.095773e-06, & 6.237762e-06, & 6.383058e-06, & 6.531739e-06, & 6.683882e-06, & 6.839570e-06, & 6.998884e-06, & 7.161909e-06, & 7.328731e-06, & 7.499439e-06, & 7.674124e-06, & 7.852877e-06, & 8.035794e-06, & 8.222971e-06, & 8.414509e-06, & 8.610508e-06, & 8.811072e-06, & 9.016309e-06, & 9.226326e-06, & 9.441235e-06, & 9.661149e-06, & 9.886186e-06, & 1.011647e-05, & 1.035211e-05, & 1.059324e-05, & 1.083999e-05, & 1.109248e-05, & 1.135086e-05, & 1.161526e-05, & 1.188581e-05, & 1.216267e-05, & 1.244597e-05, & 1.273588e-05, & 1.303253e-05, & 1.333610e-05, & 1.364674e-05, & 1.396461e-05, & 1.428989e-05, & 1.462274e-05, & 1.496335e-05, & 1.531189e-05, & 1.566855e-05, & 1.603352e-05, & 1.640699e-05, & 1.678915e-05, & 1.718022e-05, & 1.758040e-05, & 1.798990e-05, & 1.840894e-05, & 1.883774e-05, & 1.927653e-05, & 1.972553e-05, & 2.018500e-05, & 2.065517e-05, & 2.113629e-05, & 2.162862e-05, & 2.213241e-05, & 2.264794e-05, & 2.317548e-05, & 2.371531e-05, & 2.426771e-05, & 2.483298e-05, & 2.541141e-05, & 2.600332e-05, & 2.660901e-05, & 2.722882e-05, & 2.786306e-05, & 2.851207e-05, & 2.917620e-05, & 2.985580e-05, & 3.055124e-05, & 3.126287e-05, & 3.199107e-05, & 3.273624e-05, & 3.349876e-05, & 3.427905e-05, & 3.507751e-05, & 3.589457e-05, & 3.673066e-05, & 3.758623e-05, & 3.846173e-05, & 3.935762e-05, & 4.027437e-05, & 4.121248e-05, & 4.217245e-05, & 4.315477e-05, & 4.415997e-05, & 4.518859e-05, & 4.624117e-05, & 4.731826e-05, & 4.842045e-05, & 4.954830e-05, & 5.070243e-05, & 5.188344e-05, & 5.309196e-05, & 5.432863e-05, & 5.559411e-05, & 5.688906e-05, & 5.821418e-05, & 5.957016e-05, & 6.095773e-05, & 6.237762e-05, & 6.383058e-05, & 6.531738e-05, & 6.683882e-05, & 6.839570e-05, & 6.998884e-05, & 7.161909e-05, & 7.328731e-05, & 7.499439e-05, & 7.674123e-05, & 7.852877e-05, & 8.035794e-05, & 8.222971e-05, & 8.414509e-05, & 8.610508e-05, & 8.811073e-05, & 9.016309e-05, & 9.226326e-05, & 9.441235e-05, & 9.661149e-05, & 9.886186e-05, & 1.011646e-04, & 1.035211e-04, & 1.059324e-04, & 1.083999e-04, & 1.109248e-04, & 1.135086e-04, & 1.161526e-04, & 1.188581e-04, & 1.216267e-04, & 1.244597e-04, & 1.273587e-04, & 1.303253e-04, & 1.333610e-04, & 1.364674e-04, & 1.396461e-04, & 1.428989e-04, & 1.462274e-04, & 1.496335e-04, & 1.531189e-04, & 1.566855e-04, & 1.603352e-04, & 1.640699e-04, & 1.678915e-04, & 1.718022e-04, & 1.758040e-04, & 1.798990e-04, & 1.840894e-04, & 1.883774e-04, & 1.927653e-04, & 1.972553e-04, & 2.018500e-04, & 2.065517e-04, & 2.113629e-04, & 2.162862e-04, & 2.213241e-04, & 2.264794e-04, & 2.317548e-04, & 2.371531e-04, & 2.426771e-04, & 2.483298e-04, & 2.541141e-04, & 2.600332e-04, & 2.660901e-04, & 2.722882e-04, & 2.786306e-04, & 2.851207e-04, & 2.917620e-04, & 2.985580e-04, & 3.055124e-04, & 3.126286e-04, & 3.199107e-04, & 3.273624e-04, & 3.349876e-04, & 3.427905e-04, & 3.507751e-04, & 3.589457e-04, & 3.673066e-04, & 3.758623e-04, & 3.846173e-04, & 3.935762e-04, & 4.027437e-04, & 4.121248e-04, & 4.217245e-04, & 4.315477e-04, & 4.415997e-04, & 4.518859e-04, & 4.624117e-04, & 4.731826e-04, & 4.842045e-04, & 4.954830e-04, & 5.070243e-04, & 5.188344e-04, & 5.309196e-04, & 5.432863e-04, & 5.559411e-04, & 5.688906e-04, & 5.821418e-04, & 5.957016e-04, & 6.095773e-04, & 6.237762e-04, & 6.383058e-04, & 6.531738e-04, & 6.683882e-04, & 6.839570e-04, & 6.998884e-04, & 7.161909e-04, & 7.328731e-04, & 7.499439e-04, & 7.674124e-04, & 7.852877e-04, & 8.035794e-04, & 8.222971e-04, & 8.414509e-04, & 8.610508e-04, & 8.811072e-04, & 9.016309e-04, & 9.226326e-04, & 9.441234e-04, & 9.661149e-04, & 9.886186e-04, & 1.011646e-03, & 1.035211e-03, & 1.059324e-03, & 1.083999e-03, & 1.109248e-03, & 1.135086e-03, & 1.161526e-03, & 1.188581e-03, & 1.216267e-03, & 1.244597e-03, & 1.273587e-03, & 1.303253e-03, & 1.333610e-03, & 1.364674e-03, & 1.396461e-03, & 1.428989e-03, & 1.462274e-03, & 1.496335e-03, & 1.531189e-03, & 1.566855e-03, & 1.603352e-03, & 1.640698e-03, & 1.678915e-03, & 1.718022e-03, & 1.758040e-03, & 1.798990e-03, & 1.840894e-03, & 1.883774e-03, & 1.927653e-03, & 1.972553e-03, & 2.018500e-03, & 2.065517e-03, & 2.113629e-03, & 2.162862e-03, & 2.213241e-03, & 2.264794e-03, & 2.317548e-03, & 2.371531e-03, & 2.426771e-03, & 2.483298e-03, & 2.541141e-03, & 2.600332e-03, & 2.660901e-03, & 2.722882e-03, & 2.786306e-03, & 2.851207e-03, & 2.917620e-03, & 2.985580e-03, & 3.055124e-03, & 3.126286e-03, & 3.199107e-03, & 3.273624e-03, & 3.349876e-03, & 3.427905e-03, & 3.507751e-03, & 3.589457e-03, & 3.673066e-03, & 3.758623e-03, & 3.846173e-03, & 3.935762e-03, & 4.027437e-03, & 4.121248e-03, & 4.217245e-03, & 4.315477e-03, & 4.415997e-03, & 4.518859e-03, & 4.624117e-03, & 4.731826e-03, & 4.842044e-03, & 4.954830e-03, & 5.070243e-03, & 5.188344e-03, & 5.309196e-03, & 5.432863e-03, & 5.559411e-03, & 5.688906e-03, & 5.821418e-03, & 5.957016e-03, & 6.095773e-03, & 6.237762e-03, & 6.383058e-03, & 6.531739e-03, & 6.683882e-03, & 6.839570e-03, & 6.998884e-03, & 7.161909e-03, & 7.328731e-03, & 7.499439e-03, & 7.674124e-03, & 7.852877e-03, & 8.035794e-03, & 8.222971e-03, & 8.414509e-03, & 8.610508e-03, & 8.811073e-03, & 9.016309e-03, & 9.226326e-03, & 9.441234e-03, & 9.661149e-03, & 9.886186e-03, & 1.011646e-02, & 1.035211e-02, & 1.059324e-02, & 1.083999e-02, & 1.109248e-02, & 1.135086e-02, & 1.161526e-02, & 1.188581e-02, & 1.216267e-02, & 1.244597e-02, & 1.273587e-02, & 1.303253e-02, & 1.333610e-02, & 1.364674e-02, & 1.396461e-02, & 1.428989e-02, & 1.462274e-02, & 1.496335e-02, & 1.531189e-02, & 1.566855e-02, & 1.603352e-02, & 1.640698e-02, & 1.678915e-02, & 1.718022e-02, & 1.758040e-02, & 1.798990e-02, & 1.840894e-02, & 1.883774e-02, & 1.927653e-02, & 1.972553e-02, & 2.018500e-02, & 2.065517e-02, & 2.113629e-02, & 2.162862e-02, & 2.213241e-02, & 2.264794e-02, & 2.317548e-02, & 2.371531e-02, & 2.426771e-02, & 2.483298e-02, & 2.541141e-02, & 2.600332e-02, & 2.660901e-02, & 2.722882e-02, & 2.786306e-02, & 2.851207e-02, & 2.917620e-02, & 2.985580e-02, & 3.055124e-02, & 3.126287e-02, & 3.199107e-02, & 3.273624e-02, & 3.349876e-02, & 3.427905e-02, & 3.507751e-02, & 3.589457e-02, & 3.673066e-02, & 3.758623e-02, & 3.846173e-02, & 3.935762e-02, & 4.027437e-02, & 4.121248e-02, & 4.217245e-02, & 4.315477e-02, & 4.415997e-02, & 4.518859e-02, & 4.624117e-02, & 4.731826e-02, & 4.842045e-02, & 4.954830e-02, & 5.070243e-02, & 5.188344e-02, & 5.309196e-02, & 5.432863e-02, & 5.559411e-02, & 5.688906e-02, & 5.821418e-02, & 5.957016e-02, & 6.095773e-02, & 6.237762e-02, & 6.383058e-02, & 6.531738e-02, & 6.683882e-02, & 6.839570e-02, & 6.998884e-02, & 7.161909e-02, & 7.328731e-02, & 7.499439e-02, & 7.674123e-02, & 7.852877e-02, & 8.035794e-02, & 8.222971e-02, & 8.414509e-02, & 8.610508e-02, & 8.811073e-02, & 9.016309e-02, & 9.226326e-02, & 9.441234e-02, & 9.661149e-02, & 9.886186e-02, & 1.011646e-01, & 1.035211e-01, & 1.059324e-01, & 1.083999e-01, & 1.109248e-01, & 1.135086e-01, & 1.161526e-01, & 1.188581e-01, & 1.216267e-01, & 1.244597e-01, & 1.273587e-01, & 1.303253e-01, & 1.333610e-01, & 1.364674e-01, & 1.396461e-01, & 1.428989e-01, & 1.462274e-01, & 1.496335e-01, & 1.531189e-01, & 1.566855e-01, & 1.603352e-01, & 1.640699e-01, & 1.678915e-01, & 1.718022e-01, & 1.758040e-01, & 1.798990e-01, & 1.840894e-01, & 1.883774e-01, & 1.927653e-01, & 1.972553e-01, & 2.018500e-01, & 2.065517e-01, & 2.113629e-01, & 2.162862e-01, & 2.213241e-01, & 2.264794e-01, & 2.317548e-01, & 2.371531e-01, & 2.426771e-01, & 2.483298e-01, & 2.541141e-01, & 2.600332e-01, & 2.660901e-01, & 2.722882e-01, & 2.786306e-01, & 2.851207e-01, & 2.917620e-01, & 2.985581e-01, & 3.055124e-01, & 3.126287e-01, & 3.199107e-01, & 3.273624e-01, & 3.349876e-01, & 3.427905e-01, & 3.507751e-01, & 3.589457e-01, & 3.673066e-01, & 3.758623e-01, & 3.846173e-01, & 3.935762e-01, & 4.027437e-01, & 4.121248e-01, & 4.217245e-01, & 4.315477e-01, & 4.415997e-01, & 4.518859e-01, & 4.624117e-01, & 4.731826e-01, & 4.842045e-01, & 4.954830e-01, & 5.070243e-01, & 5.188344e-01, & 5.309196e-01, & 5.432863e-01, & 5.559411e-01, & 5.688906e-01, & 5.821418e-01, & 5.957016e-01, & 6.095773e-01, & 6.237762e-01, & 6.383058e-01, & 6.531738e-01, & 6.683882e-01, & 6.839570e-01, & 6.998884e-01, & 7.161909e-01, & 7.328731e-01, & 7.499439e-01, & 7.674123e-01, & 7.852877e-01, & 8.035794e-01, & 8.222971e-01, & 8.414509e-01, & 8.610508e-01, & 8.811073e-01, & 9.016309e-01, & 9.226326e-01, & 9.441234e-01, & 9.661149e-01, & 9.886186e-01/ * 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, & 1.298172e+03, & 1.904547e+03, & 1.979433e+03, & 2.053128e+03, & 2.113408e+03, & 2.183430e+03, & 2.253166e+03, & 6.107975e+04, & 6.033124e+04, & 6.161803e+04, & 6.285116e+04, & 6.384949e+04, & 6.496772e+04, & 6.627731e+04, & 6.745324e+04, & 6.897230e+04, & 7.015387e+04, & 7.110013e+04, & 7.205257e+04, & 7.316505e+04, & 7.423485e+04, & 7.533692e+04, & 7.585111e+04, & 7.666678e+04, & 7.745878e+04, & 7.832771e+04, & 7.897953e+04, & 7.974127e+04, & 8.006441e+04, & 8.081209e+04, & 8.117260e+04, & 8.137891e+04, & 8.155716e+04, & 8.182761e+04, & 8.171056e+04, & 8.172409e+04, & 8.170204e+04, & 8.121010e+04, & 8.089478e+04, & 8.018013e+04, & 7.977229e+04, & 7.913907e+04, & 7.835423e+04, & 7.740730e+04, & 7.639541e+04, & 7.548152e+04, & 7.436404e+04, & 7.315116e+04, & 7.167284e+04, & 7.050098e+04, & 6.906649e+04, & 6.761841e+04, & 6.626786e+04, & 6.451853e+04, & 6.294980e+04, & 6.115838e+04, & 5.947762e+04, & 5.765314e+04, & 5.589939e+04, & 5.418876e+04, & 5.237000e+04, & 5.068134e+04, & 4.881320e+04, & 4.709309e+04, & 4.533207e+04, & 4.368588e+04, & 4.192920e+04, & 4.019882e+04, & 3.864330e+04, & 3.711586e+04, & 3.543141e+04, & 3.395433e+04, & 3.254742e+04, & 3.102951e+04, & 2.967494e+04, & 2.833312e+04, & 2.706891e+04, & 2.589990e+04, & 2.467824e+04, & 2.364744e+04, & 2.259032e+04, & 2.172457e+04, & 2.082001e+04, & 1.994377e+04, & 1.917980e+04, & 1.847481e+04, & 1.784225e+04, & 1.734220e+04, & 1.673615e+04, & 1.623385e+04, & 1.568117e+04, & 1.525791e+04, & 1.492831e+04, & 1.458010e+04, & 1.425362e+04, & 1.409100e+04, & 1.386577e+04, & 1.373975e+04, & 1.342496e+04, & 1.328016e+04, & 1.320555e+04, & 1.314362e+04, & 1.293126e+04, & 1.288325e+04, & 1.280495e+04, & 1.278742e+04, & 1.270477e+04, & 1.271402e+04, & 1.268078e+04, & 1.263796e+04, & 1.259401e+04, & 1.261164e+04, & 1.247532e+04, & 1.255386e+04, & 1.260605e+04, & 1.262573e+04, & 1.248173e+04, & 1.256497e+04, & 1.267679e+04, & 1.260632e+04, & 1.257839e+04, & 1.247733e+04, & 1.241515e+04, & 1.244797e+04, & 1.243605e+04, & 1.233610e+04, & 1.234966e+04, & 1.251236e+04, & 1.220118e+04, & 1.226735e+04, & 1.213882e+04, & 1.217403e+04, & 1.203139e+04, & 1.204511e+04, & 1.201017e+04, & 1.200308e+04, & 1.183536e+04, & 1.191270e+04, & 1.184755e+04, & 1.186814e+04, & 1.188857e+04, & 1.188582e+04, & 1.185717e+04, & 1.173617e+04, & 1.186083e+04, & 1.180347e+04, & 1.185724e+04, & 1.181394e+04, & 1.177899e+04, & 1.176406e+04, & 1.170089e+04, & 1.180704e+04, & 1.191833e+04, & 1.180981e+04, & 1.170911e+04, & 1.171646e+04, & 1.162167e+04, & 1.164121e+04, & 1.168711e+04, & 1.160232e+04, & 1.166446e+04, & 1.167069e+04, & 1.167975e+04, & 1.156432e+04, & 1.165113e+04, & 1.160083e+04, & 1.170457e+04, & 1.161063e+04, & 1.166398e+04, & 1.175907e+04, & 1.164992e+04, & 1.160783e+04, & 1.155609e+04, & 1.159553e+04, & 1.152608e+04, & 1.167764e+04, & 1.150842e+04, & 1.160214e+04, & 1.158257e+04, & 1.165161e+04, & 1.150025e+04, & 1.153957e+04, & 1.154878e+04, & 1.156519e+04, & 1.152349e+04, & 1.162911e+04, & 1.173302e+04, & 1.146811e+04, & 1.159045e+04, & 1.156528e+04, & 1.157017e+04, & 1.164703e+04, & 1.151527e+04, & 1.165353e+04, & 1.145264e+04, & 1.157366e+04, & 1.138327e+04, & 1.167679e+04, & 1.166421e+04, & 1.145674e+04, & 1.166537e+04, & 1.161423e+04, & 1.165163e+04, & 1.167426e+04, & 1.154968e+04, & 1.159803e+04, & 1.164895e+04, & 1.163035e+04, & 1.170682e+04, & 1.146302e+04, & 1.157871e+04, & 1.159620e+04, & 1.150456e+04, & 1.162663e+04, & 1.175954e+04, & 1.150654e+04, & 1.158538e+04, & 1.153789e+04, & 1.144510e+04, & 1.160537e+04, & 1.147068e+04, & 1.164400e+04, & 1.144484e+04, & 1.153889e+04, & 1.158294e+04, & 1.159771e+04, & 1.164968e+04, & 1.158238e+04, & 1.174142e+04, & 1.181917e+04, & 1.154788e+04, & 1.165187e+04, & 1.170454e+04, & 1.183983e+04, & 1.161659e+04, & 1.178201e+04, & 1.175453e+04, & 1.175582e+04, & 1.169119e+04, & 1.173998e+04, & 1.192239e+04, & 1.174169e+04, & 1.160997e+04, & 1.164417e+04, & 1.176621e+04, & 1.172904e+04, & 1.155342e+04, & 1.190261e+04, & 1.162636e+04, & 1.152790e+04, & 1.174014e+04, & 1.173105e+04, & 1.184953e+04, & 1.172194e+04, & 1.181405e+04, & 1.194161e+04, & 1.179535e+04, & 1.168717e+04, & 1.183555e+04, & 1.174827e+04, & 1.185852e+04, & 1.153410e+04, & 1.191629e+04, & 1.179004e+04, & 1.194662e+04, & 1.176206e+04, & 1.181440e+04, & 1.179181e+04, & 1.184997e+04, & 1.178445e+04, & 1.190317e+04, & 1.201659e+04, & 1.214109e+04, & 1.203119e+04, & 1.194551e+04, & 1.222608e+04, & 1.200804e+04, & 1.213698e+04, & 1.193575e+04, & 1.204054e+04, & 1.198824e+04, & 1.205553e+04, & 1.213488e+04, & 1.189820e+04, & 1.209229e+04, & 1.174793e+04, & 1.225458e+04, & 1.204104e+04, & 1.190194e+04, & 1.218937e+04, & 1.209000e+04, & 1.211213e+04, & 1.202110e+04, & 1.229502e+04, & 1.190932e+04, & 1.239532e+04, & 1.193874e+04, & 1.191288e+04, & 1.238335e+04, & 1.217524e+04, & 1.229768e+04, & 1.173470e+04, & 1.216761e+04, & 1.243566e+04, & 1.225240e+04, & 1.225759e+04, & 1.227304e+04, & 1.205459e+04, & 1.248692e+04, & 1.220830e+04, & 1.220217e+04, & 1.232141e+04, & 1.211633e+04, & 1.246527e+04, & 1.234186e+04, & 1.225659e+04, & 1.213723e+04, & 1.250128e+04, & 1.225751e+04, & 1.233811e+04, & 1.221671e+04, & 1.252521e+04, & 1.246868e+04, & 1.229953e+04, & 1.218755e+04, & 1.216680e+04, & 1.208872e+04, & 1.269272e+04, & 1.221081e+04, & 1.238501e+04, & 1.239347e+04, & 1.227930e+04, & 1.212807e+04, & 1.242263e+04, & 1.260411e+04, & 1.237529e+04, & 1.251711e+04, & 1.264028e+04, & 1.240896e+04, & 1.244747e+04, & 1.244701e+04, & 1.230197e+04, & 1.238837e+04, & 1.288812e+04, & 1.282639e+04, & 1.262732e+04, & 1.265463e+04, & 1.237226e+04, & 1.201825e+04, & 1.287569e+04, & 1.239773e+04, & 1.266910e+04, & 1.262583e+04, & 1.262271e+04, & 1.245894e+04, & 1.268603e+04, & 1.278778e+04, & 1.261259e+04, & 1.274568e+04, & 1.236188e+04, & 1.291925e+04, & 1.305701e+04, & 1.286175e+04, & 1.281197e+04, & 1.249182e+04, & 1.237160e+04, & 1.303837e+04, & 1.228648e+04, & 1.270864e+04, & 1.307987e+04, & 1.254718e+04, & 1.278573e+04, & 1.269285e+04, & 1.276933e+04, & 1.300571e+04, & 1.294656e+04, & 1.321342e+04, & 1.288205e+04, & 1.282879e+04, & 1.244703e+04, & 1.258180e+04, & 1.229106e+04, & 1.279392e+04, & 1.335228e+04, & 1.258333e+04, & 1.270935e+04, & 1.287988e+04, & 1.275637e+04, & 1.292778e+04, & 1.252098e+04, & 1.322585e+04, & 1.252545e+04, & 1.309385e+04, & 1.286841e+04, & 1.301121e+04, & 1.316715e+04, & 1.278377e+04, & 1.244677e+04, & 1.359709e+04, & 1.284871e+04, & 1.288594e+04, & 1.323112e+04, & 1.308941e+04, & 1.324029e+04, & 1.268674e+04, & 1.269172e+04, & 1.334935e+04, & 1.386404e+04, & 1.295440e+04, & 1.228927e+04, & 1.067144e+04, & 9.542012e+03, & 1.028028e+04, & 1.079921e+04, & 1.157240e+04, & 1.224763e+04, & 1.268253e+04, & 1.273672e+04, & 1.219884e+04, & 1.276176e+04, & 1.277131e+04, & 1.321071e+04, & 1.220336e+04, & 1.340327e+04, & 1.310747e+04, & 1.351254e+04, & 1.339301e+04, & 1.348503e+04, & 1.353388e+04, & 1.353369e+04, & 1.338790e+04, & 1.284719e+04, & 1.305327e+04, & 1.302725e+04, & 1.306161e+04, & 1.340422e+04, & 1.358410e+04, & 1.334882e+04, & 1.398133e+04, & 1.297120e+04, & 1.362595e+04, & 1.317582e+04, & 1.308189e+04, & 1.346410e+04, & 1.356517e+04, & 1.338804e+04, & 1.345903e+04, & 1.329921e+04, & 1.376886e+04, & 1.305782e+04, & 1.355686e+04, & 1.322402e+04, & 1.320462e+04, & 1.375225e+04, & 1.348581e+04, & 1.420645e+04, & 1.364764e+04, & 1.344479e+04, & 1.336559e+04, & 1.344743e+04, & 1.345485e+04, & 1.321667e+04, & 1.289245e+04, & 1.166830e+04, & 1.413811e+04, & 1.395685e+04, & 1.346633e+04, & 1.407681e+04, & 1.328791e+04, & 1.380343e+04, & 1.382439e+04, & 1.378234e+04, & 1.403786e+04, & 1.372544e+04, & 1.426810e+04, & 1.397285e+04, & 1.399125e+04, & 1.374634e+04, & 1.274304e+04, & 1.322516e+04, & 1.374377e+04, & 1.368096e+04, & 1.357333e+04, & 1.382487e+04, & 1.435437e+04, & 1.497302e+04, & 1.396349e+04, & 1.393899e+04, & 1.401346e+04, & 1.329215e+04, & 1.459822e+04, & 1.372683e+04, & 1.371604e+04, & 1.391793e+04, & 1.305926e+04, & 1.263040e+04, & 1.326549e+04, & 1.312794e+04, & 1.292955e+04, & 1.375481e+04, & 1.336430e+04, & 1.453334e+04, & 1.427546e+04, & 1.350531e+04, & 1.486116e+04, & 1.354710e+04, & 1.467652e+04, & 1.441238e+04, & 1.426011e+04, & 1.464197e+04, & 1.382086e+04, & 1.426167e+04, & 1.363423e+04, & 1.508589e+04, & 1.502297e+04, & 1.551959e+04, & 1.428355e+04, & 1.428312e+04, & 1.374762e+04, & 1.410865e+04, & 1.477573e+04, & 1.507622e+04, & 1.425692e+04, & 1.478842e+04, & 1.442796e+04, & 1.510890e+04, & 1.525265e+04, & 1.467292e+04, & 1.494210e+04, & 1.475304e+04, & 1.521745e+04, & 1.435134e+04, & 1.509088e+04, & 1.546619e+04, & 1.539801e+04, & 1.259408e+04, & 1.186332e+04, & 1.478774e+04, & 1.507021e+04, & 1.528230e+04, & 1.497692e+04, & 1.529834e+04, & 1.499909e+04, & 1.510147e+04, & 1.476942e+04, & 1.501786e+04, & 1.445070e+04, & 1.480055e+04, & 1.451589e+04, & 1.510733e+04, & 1.508938e+04, & 1.558506e+04, & 1.567018e+04, & 1.558110e+04, & 1.530236e+04, & 1.601684e+04, & 1.650231e+04, & 1.591757e+04, & 1.603603e+04, & 1.654472e+04, & 1.662095e+04, & 1.597282e+04, & 1.716857e+04, & 1.602729e+04, & 1.625079e+04, & 1.658221e+04, & 1.705737e+04, & 1.708313e+04, & 1.617335e+04, & 1.654201e+04, & 1.656876e+04, & 1.635137e+04, & 1.694882e+04, & 1.705613e+04, & 1.683449e+04, & 1.700052e+04, & 1.702121e+04, & 1.694094e+04, & 1.699194e+04, & 1.692447e+04, & 1.705421e+04, & 1.642077e+04, & 1.732792e+04, & 1.799758e+04, & 1.727093e+04, & 1.796241e+04, & 1.740731e+04, & 1.801280e+04, & 1.837918e+04, & 1.842544e+04, & 1.833711e+04, & 1.844690e+04, & 1.809165e+04, & 1.854113e+04, & 2.012226e+04, & 1.912625e+04, & 2.027809e+04, & 2.047055e+04, & 2.000514e+04, & 2.004902e+04, & 2.036879e+04, & 1.983973e+04, & 1.857994e+04, & 1.948308e+04, & 2.063442e+04, & 2.181667e+04, & 2.271705e+04, & 2.382927e+04, & 2.184123e+04, & 2.056741e+04, & 1.741384e+04, & 1.358251e+04, & 1.026562e+04, & 8.894763e+03, & 1.012705e+04, & 1.252133e+04, & 1.402591e+04, & 1.515958e+04, & 1.647536e+04, & 1.702351e+04, & 1.716663e+04, & 1.944076e+04, & 1.846268e+04, & 1.949886e+04, & 1.947514e+04, & 2.040921e+04, & 1.979024e+04, & 2.173701e+04, & 2.091824e+04, & 2.117814e+04, & 2.179550e+04, & 2.118691e+04, & 2.316157e+04, & 2.322422e+04, & 2.298080e+04, & 2.278141e+04, & 2.531163e+04, & 2.392483e+04, & 2.550796e+04, & 2.497314e+04, & 2.503807e+04, & 2.678720e+04, & 2.659299e+04, & 2.737628e+04, & 2.815029e+04, & 2.794558e+04, & 2.794430e+04, & 2.803181e+04, & 2.725606e+04, & 2.347680e+04, & 2.173585e+04, & 1.741068e+04, & 1.676768e+04, & 1.654198e+04, & 1.769772e+04, & 1.925237e+04, & 2.177193e+04, & 2.302036e+04, & 2.500741e+04, & 2.617068e+04, & 2.465408e+04, & 2.467432e+04, & 2.613219e+04, & 2.829308e+04, & 2.866313e+04, & 2.867131e+04, & 2.046879e+04, & 2.797607e+04, & 3.334175e+04, & 3.648048e+04, & 3.729489e+04, & 3.510796e+04, & 2.987080e+04, & 2.596593e+04, & 2.427716e+04, & 2.510034e+04, & 2.551829e+04, & 2.546528e+04, & 2.328114e+04, & 2.545069e+04, & 2.947940e+04, & 3.267221e+04, & 3.319936e+04, & 3.528375e+04, & 3.682339e+04, & 3.711969e+04, & 3.844831e+04, & 3.972142e+04, & 3.706058e+04, & 3.246412e+04, & 3.143490e+04, & 3.558150e+04, & 4.128527e+04, & 3.725605e+04, & 3.802035e+04, & 4.105422e+04, & 4.522794e+04, & 4.358403e+04, & 4.650706e+04, & 4.425377e+04, & 4.958341e+04, & 5.234419e+04, & 4.850503e+04, & 4.123810e+04, & 4.228530e+04, & 4.354858e+04, & 4.715144e+04, & 4.888158e+04, & 4.900005e+04, & 5.063202e+04, & 5.663450e+04, & 5.707960e+04, & 5.704133e+04, & 5.842670e+04, & 5.194750e+04, & 5.077467e+04, & 5.949167e+04, & 6.043932e+04, & 5.870261e+04, & 5.697021e+04, & 4.526410e+04, & 4.216297e+04, & 4.091271e+04, & 4.547030e+04, & 6.013171e+04, & 6.667975e+04, & 7.825452e+04, & 8.459661e+04, & 9.069122e+04, & 8.935499e+04, & 7.261188e+04, & 8.370061e+04, & 9.023339e+04, & 8.429514e+04, & 8.916780e+04, & 9.088530e+04, & 9.585164e+04, & 9.141152e+04, & 8.929451e+04, & 8.414587e+04, & 9.945975e+04, & 1.050508e+05, & 9.932380e+04, & 9.761591e+04, & 9.859966e+04, & 1.028864e+05, & 1.043741e+05, & 9.946281e+04, & 8.945520e+04, & 9.106286e+04, & 1.016413e+05, & 9.427774e+04, & 1.023644e+05, & 1.042830e+05, & 1.017297e+05, & 9.540350e+04, & 8.965743e+04, & 8.311816e+04, & 7.984286e+04 / * ecited/(excited+ground state) cross section 10B(n,alpha)7Li * cross section ratio *------------------- DATA RATIO/ & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937102, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937101, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937100, & 0.937099, & 0.937099, & 0.937099, & 0.937099, & 0.937099, & 0.937099, & 0.937099, & 0.937099, & 0.937099, & 0.937099, & 0.937099, & 0.937098, & 0.937098, & 0.937098, & 0.937098, & 0.937098, & 0.937098, & 0.937097, & 0.937097, & 0.937097, & 0.937097, & 0.937096, & 0.937096, & 0.937096, & 0.937095, & 0.937095, & 0.937094, & 0.937094, & 0.937093, & 0.937092, & 0.937092, & 0.937092, & 0.937092, & 0.937091, & 0.937091, & 0.937091, & 0.937091, & 0.937090, & 0.937090, & 0.937090, & 0.937090, & 0.937089, & 0.937089, & 0.937089, & 0.937089, & 0.937088, & 0.937088, & 0.937088, & 0.937087, & 0.937087, & 0.937086, & 0.937086, & 0.937085, & 0.937085, & 0.937084, & 0.937084, & 0.937083, & 0.937083, & 0.937082, & 0.937081, & 0.937081, & 0.937080, & 0.937079, & 0.937078, & 0.937077, & 0.937077, & 0.937076, & 0.937075, & 0.937075, & 0.937074, & 0.937074, & 0.937073, & 0.937073, & 0.937072, & 0.937072, & 0.937071, & 0.937070, & 0.937070, & 0.937069, & 0.937069, & 0.937068, & 0.937068, & 0.937067, & 0.937067, & 0.937066, & 0.937066, & 0.937065, & 0.937064, & 0.937064, & 0.937063, & 0.937062, & 0.937061, & 0.937061, & 0.937060, & 0.937057, & 0.937055, & 0.937053, & 0.937051, & 0.937048, & 0.937045, & 0.937044, & 0.937043, & 0.937042, & 0.937040, & 0.937039, & 0.937038, & 0.937037, & 0.937035, & 0.937034, & 0.937033, & 0.937031, & 0.937030, & 0.937029, & 0.937027, & 0.937026, & 0.937024, & 0.937023, & 0.937021, & 0.937019, & 0.937017, & 0.937016, & 0.937014, & 0.937011, & 0.937009, & 0.937007, & 0.937004, & 0.937002, & 0.936999, & 0.936996, & 0.936993, & 0.936991, & 0.936989, & 0.936987, & 0.936984, & 0.936982, & 0.936979, & 0.936976, & 0.936974, & 0.936971, & 0.936968, & 0.936965, & 0.936961, & 0.936958, & 0.936954, & 0.936950, & 0.936946, & 0.936942, & 0.936938, & 0.936933, & 0.936928, & 0.936923, & 0.936917, & 0.936913, & 0.936908, & 0.936904, & 0.936900, & 0.936895, & 0.936890, & 0.936885, & 0.936879, & 0.936874, & 0.936868, & 0.936862, & 0.936855, & 0.936848, & 0.936841, & 0.936834, & 0.936828, & 0.936822, & 0.936816, & 0.936810, & 0.936803, & 0.936796, & 0.936789, & 0.936781, & 0.936773, & 0.936765, & 0.936756, & 0.936749, & 0.936741, & 0.936733, & 0.936724, & 0.936715, & 0.936706, & 0.936696, & 0.936685, & 0.936675, & 0.936666, & 0.936656, & 0.936645, & 0.936635, & 0.936623, & 0.936611, & 0.936599, & 0.936588, & 0.936577, & 0.936566, & 0.936553, & 0.936541, & 0.936527, & 0.936514, & 0.936501, & 0.936487, & 0.936472, & 0.936457, & 0.936442, & 0.936428, & 0.936413, & 0.936397, & 0.936380, & 0.936367, & 0.936353, & 0.936340, & 0.936326, & 0.936311, & 0.936295, & 0.936279, & 0.936263, & 0.936245, & 0.936227, & 0.936207, & 0.936187, & 0.936166, & 0.936144, & 0.936121, & 0.936096, & 0.936070, & 0.936043, & 0.936015, & 0.935985, & 0.935963, & 0.935943, & 0.935921, & 0.935898, & 0.935875, & 0.935850, & 0.935824, & 0.935797, & 0.935768, & 0.935739, & 0.935707, & 0.935674, & 0.935645, & 0.935622, & 0.935598, & 0.935572, & 0.935546, & 0.935518, & 0.935489, & 0.935459, & 0.935427, & 0.935395, & 0.935361, & 0.935325, & 0.935288, & 0.935249, & 0.935208, & 0.935165, & 0.935120, & 0.935073, & 0.935033, & 0.934994, & 0.934953, & 0.934910, & 0.934866, & 0.934819, & 0.934771, & 0.934720, & 0.934668, & 0.934613, & 0.934555, & 0.934494, & 0.934431, & 0.934365, & 0.934295, & 0.934222, & 0.934145, & 0.934065, & 0.933989, & 0.933910, & 0.933828, & 0.933742, & 0.933652, & 0.933558, & 0.933460, & 0.933357, & 0.933252, & 0.933148, & 0.933040, & 0.932927, & 0.932808, & 0.932685, & 0.932555, & 0.932420, & 0.932278, & 0.932130, & 0.931975, & 0.931813, & 0.931645, & 0.931468, & 0.931266, & 0.931054, & 0.930834, & 0.930603, & 0.930363, & 0.930107, & 0.929839, & 0.929559, & 0.929267, & 0.928962, & 0.928651, & 0.928330, & 0.927975, & 0.927585, & 0.927179, & 0.926757, & 0.926316, & 0.925857, & 0.925378, & 0.924878, & 0.924328, & 0.923737, & 0.923121, & 0.922481, & 0.921814, & 0.921119, & 0.920395, & 0.919639, & 0.918851, & 0.918027, & 0.917133, & 0.916196, & 0.915219, & 0.914200, & 0.913134, & 0.912011, & 0.910825, & 0.909584, & 0.908283, & 0.906920, & 0.905501, & 0.904020, & 0.902482, & 0.900894, & 0.899244, & 0.897562, & 0.895810, & 0.894025, & 0.892160, & 0.890259, & 0.888317, & 0.886339, & 0.884304, & 0.882153, & 0.879992, & 0.877716, & 0.875343, & 0.872877, & 0.870262, & 0.867462, & 0.864514, & 0.861385, & 0.858043, & 0.854445, & 0.850568, & 0.846459, & 0.842075, & 0.837248, & 0.832188, & 0.826852, & 0.821033, & 0.815019, & 0.808746, & 0.802179, & 0.795467, & 0.788815, & 0.782124, & 0.775475, & 0.769474, & 0.763239, & 0.758088, & 0.752865, & 0.748255, & 0.743886, & 0.740198, & 0.736471, & 0.733107, & 0.729269, & 0.725663, & 0.721613, & 0.717517, & 0.713134, & 0.708122, & 0.702588, & 0.697275, & 0.691266, & 0.684417, & 0.678544, & 0.671992, & 0.664620, & 0.660160, & 0.655260, & 0.649523, & 0.641827, & 0.633356, & 0.624880, & 0.616330, & 0.607074, & 0.598158, & 0.588621, & 0.579077, & 0.569336 / *======================================================================* * * * BASIC VERSION * * * *======================================================================* NOMORE = 0 * +-------------------------------------------------------------------* * | First call initializations: IF ( LFIRST ) THEN * | *** The following 3 cards are mandatory *** TKESUM = ZERZER LFIRST = .FALSE. LUSSRC = .TRUE. CALL myusrini() * | *** User initialization *** END IF * | * +-------------------------------------------------------------------* * WRITE (LNUOUT,*) ' SOURCE ', NPARMA * 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 * +-------------------------------------------------------------------* * Calculate the cumulative function TEGR(0)=FLUX(0) DO I = 1, 799 TEGR(I) = TEGR(I-1) + FLUX(I) END DO * Normalize the cumulative function to 1 SOMMA = TEGR(799) DO I =0, 799 TEGR(I)=TEGR(I)/SOMMA END DO * +-------------------------------------------------------------------* * The kinetic energy of the spallation neutron is sampled from the * neutron flux 5 XYZ= FLRNDM(XYYZ) DO 10 I=0, 798 IF (TEGR(I+1) .GT. XYZ .AND. TEGR(I+1) - TEGR(I).NE.0) THEN TNEUTENE = ENEDGE(I) + (ENEDGE(I+1) -ENEDGE(I))* & (XYZ-TEGR(I))/(TEGR(I+1) - TEGR(I)) INDICE = I GO TO 20 END IF IF (TEGR(I+1) .GT. XYZ. AND.TEGR(I+1) - TEGR(I). EQ.0) THEN WRITE(LNUOUT,*) 'BIZARRE ', NCASE GO TO 5 END IF 10 CONTINUE 20 CONTINUE **** * +-------------------------------------------------------------------* * | (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 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 neutrons in GeV * * * *------------------------------------------------------------* * * ENE = TNEUTENE *0.001 * Check on the neutron energy IF (TNEUTENE.LT.1.011646e-08) THEN WRITE(LNUOUT,*) ' SAMPLING ERROR: Neutron energy below limit ' END IF IF (TNEUTENE.GT.9.886186e-01) THEN WRITE(LNUOUT,*) ' SAMPLING ERROR: Neutron energy above limit ' END IF CALL treefill(NCASE, ENE) TKEFLK (NPFLKA) = ENE * Particle momentum * PMOFLK (NPFLKA) = PBEAM PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( & TKEFLK (NPFLKA) + TWOTWO * AM (IONID) ) ) * * * * Cosines (tx,ty,tz) * You want that neutrons are emitted along the z axis UBEAM = 0.0 VBEAM = 0.0 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 CALL FLNRR2 (GAUSS1, GAUSS2) IF (ENE.GE.1.0e-11.AND.ENE.LE.2.15443e-11) THEN X1 = GAUSS1*5.67161e-01 Y1 = GAUSS2*5.39079e-01 ELSE IF (ENE.GT.2.15443e-11.AND.ENE.LE.4.64159e-11) THEN X1= GAUSS1*5.27498e-01 Y1= GAUSS2*5.40903e-01 ELSE IF (ENE.GT.4.64159e-11.AND.ENE.LE.1.e-10) THEN X1= GAUSS1*5.42656e-01 Y1= GAUSS2*5.55597e-01 ELSE IF (ENE.GT.1.e-10.AND.ENE.LE.2.15443e-10) THEN X1= GAUSS1*5.47098e-01 Y1= GAUSS2*5.49595e-01 ELSE IF (ENE.GT.2.15443e-10.AND.ENE.LE.4.64159e-10) THEN X1= GAUSS1*5.71072e-01 Y1= GAUSS2*5.80309e-01 ELSE IF (ENE.GT.4.64159e-10.AND.ENE.LE.1.e-09) THEN X1= GAUSS1*5.76742e-01 Y1= GAUSS2*5.80790e-01 ELSE IF (ENE.GT.1.e-09.AND.ENE.LE.2.15443e-09) THEN X1= GAUSS1*5.64364e-01 Y1= GAUSS2*5.71395e-01 ELSE IF (ENE.GT.2.15443e-09.AND.ENE.LE.4.64159e-09) THEN X1= GAUSS1*5.68613e-01 Y1= GAUSS2*5.70629e-01 ELSE IF (ENE.GT.4.64159e-09.AND.ENE.LE.1.e-08) THEN X1= GAUSS1*5.71982e-01 Y1= GAUSS2*5.61733e-01 ELSE IF (ENE.GT.1.e-08.AND.ENE.LE.2.15443e-08) THEN X1= GAUSS1*5.69528e-01 Y1= GAUSS2*5.77204e-01 ELSE IF (ENE.GT.2.15443e-08.AND.ENE.LE.4.64159e-08) THEN X1= GAUSS1*5.73158e-01 Y1= GAUSS2*5.58874e-01 ELSE IF (ENE.GT.4.64159e-08.AND.ENE.LE.1.e-07) THEN X1= GAUSS1*5.72613e-01 Y1= GAUSS2*5.78091e-01 ELSE IF (ENE.GT.1.e-07.AND.ENE.LE.2.15443e-07) THEN X1= GAUSS1*5.58794e-01 Y1= GAUSS2*5.62267e-01 ELSE IF (ENE.GT.2.15443e-07.AND.ENE.LE.4.64159e-07) THEN X1= GAUSS1*5.67800e-01 Y1= GAUSS2*5.70722e-01 ELSE IF (ENE.GT.4.64159e-07.AND.ENE.LE.1.e-06) THEN X1= GAUSS1*5.78429e-01 Y1= GAUSS2*5.74007e-01 ELSE IF (ENE.GT.1.e-06.AND.ENE.LE.2.15443e-06) THEN X1= GAUSS1*5.71280e-01 Y1= GAUSS2*5.64139e-01 ELSE IF (ENE.GT.2.15443e-06.AND.ENE.LE.4.64159e-06) THEN X1= GAUSS1*5.69914e-01 Y1= GAUSS2*5.64641e-01 ELSE IF (ENE.GT.4.64159e-06.AND.ENE.LE.1.e-05) THEN X1= GAUSS1*5.71614e-01 Y1= GAUSS2*5.68463e-01 ELSE IF (ENE.GT.1.e-05.AND.ENE.LE.2.15443e-05) THEN X1= GAUSS1*5.68169e-01 Y1= GAUSS2*5.76375e-01 ELSE IF (ENE.GT.2.15443e-05.AND.ENE.LE.4.64159e-05) THEN X1= GAUSS1*5.71080e-01 Y1= GAUSS2*5.63680e-01 ELSE IF (ENE.GT.4.64159e-05.AND.ENE.LE.1.e-04) THEN X1= GAUSS1*5.72134e-01 Y1= GAUSS2*5.78628e-01 ELSE IF (ENE.GT.1.e-04.AND.ENE.LE.2.15443e-04) THEN X1= GAUSS1*5.69799e-01 Y1= GAUSS2*5.60572e-01 ELSE IF (ENE.GT.2.15443e-04.AND.ENE.LE.4.64159e-04) THEN X1= GAUSS1*5.56503e-01 Y1= GAUSS2*5.45938e-01 ELSE X1= GAUSS1*5.42785e-01 Y1= GAUSS2*5.32587e-01 END IF XFLK (NPFLKA) = X1 YFLK (NPFLKA) = Y1 ZFLK (NPFLKA) = 0.9 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