You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2. htape3x int=int2 outt=resid<br />
(„int2‟ is the input file for HTAPE3X with the options IOPT=108 and KOPT=0. „resid‟ is the<br />
output containing the residual nuclei).<br />
Here is an example for int2:<br />
Residual nuclei<br />
in 25 cells<br />
108,,,,0,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 />
Until now this is just a standard MCNPX procedure.<br />
SP-FISPACT2001 will read the values of „resid‟ (stable and metastable isotopes) and subtract<br />
these values from „elast‟ to obtain the residual nuclei per source proton; finally it will create<br />
the preliminary file COLLAPX with the new pseudo cross-sections.<br />
3.2 Preparation of the preliminary input COLLAPX<br />
The file FLUXES for FISPACT needs to be prepared with the values of the flux <strong>sp</strong>ectrum<br />
under 20 MeV as provided by MCNPX.<br />
In order to prepare properly the preliminary input „collapx.i‟ for SP-FISPACT, the user must<br />
<strong>sp</strong>ecify before the code word COLLAPSE: the value of the flux under 20MeV irradiating the<br />
material, the source proton current used, the material irradiated and the number of atoms of<br />
this material. Five new code words have been introduced for this purpose: CURRENT,<br />
TARGET, ATTARGET, HCOLLAPSE, TCOLLAPSE.<br />
The file collapx.i has to be prepared following these steps:<br />
1. Input the integrated flux under 20MeV using the code word FLUX (this is a standard<br />
FISPACT code word, but it has to be <strong>sp</strong>ecified also in the preliminary input and not only in<br />
the main input). If activation calculations are needed only above 20MeV then enter in any<br />
case a value of the flux as well as a <strong>sp</strong>ectrum of the flux in the file FLUXES. These will be<br />
fictitious values without any consequence for the results. Use a realistic value, such as FLUX<br />
1E15 (see WARNINGS n. 5 at page 14).<br />
Example: if there is a flux under 20MeV of 10 14 neutrons/cm 2 sec, then write: FLUX 1E14.<br />
2. Input the proton current (in mA) used as source particles in MCNPX. Since the results in<br />
MCNPX are, as stated above, provided as residual nuclei per source particle, the value of the<br />
current has to be fixed in order to provide the correct number of residual nuclei per second.<br />
The new code word is CURRENT.<br />
Example: if the proton current of the accelerator is 3mA, then write: CURRENT 3.<br />
11