Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
SP-FISPACT, together with these cross-sections, stores the „pseudo cross-sections‟ in this file<br />
and, thus, the different utilization of SP-FISPACT with re<strong>sp</strong>ect to FISPACT concerns only<br />
the way the file COLLAPX is prepared by the input file „collapx.i‟.<br />
See APPENDIX A for the detailed modifications to the subroutines of FISPACT2001 and<br />
APPENDIX B for the definition of the „pseudo cross-section‟.<br />
3. Calculations with SP-FISPACT2001<br />
The procedure for carrying out the calculations is shown in the following page and is<br />
explained here step by step.<br />
3.1 Preparation of the outputs of MCNPX<br />
The MCNPX input file requires the following cards :<br />
- „f4:n‟. This card calculates the volumetric neutron flux. Choose one of the energy group<br />
structures of FISPACT2001 (XMAS, VITAMIN-J, TRIPOLI, see [4]) 2 and, using the card<br />
„e4‟, define the energy bins accordingly.<br />
- „histp‟. This card creates the file „histp‟ (it can be quite a lengthy file). This file is postprocessed<br />
by the code 3 HTAPE3X [1] for calculation of residual nuclei per source proton<br />
after interactions >20MeV.<br />
Two output files are needed from HTAPE3X: the first providing the elastic collisions (this<br />
output file must be called „elast‟), the second providing the residual nuclei after elastic<br />
and inelastic collisions above 20MeV (this ouput file must be called „resid‟):<br />
1. htape3x int=int1 outt=elast<br />
(„int1‟ is the input file for HTAPE3X with the options IOPT=115 and KOPT=2. „elast‟ is the<br />
output containing the elastic collisions).<br />
Here is an example for int1:<br />
Elastic collision<br />
in 25 cells<br />
115,,,,2,25,,,,/<br />
15,20,25,30,44,46,48,52,54,56,58,60,62,64,92,94,96,98,104,<br />
106,108,110,112,114,116<br />
2 The group structures GAM-II and WIMS are not suitable since their highest energy values are re<strong>sp</strong>ectively 14.9<br />
MeV and 10 MeV.<br />
3 This code belongs to the MCNPX package.<br />
9