Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
APPENDIX A – Fortran Modifications to FISPACT2001<br />
Only one subroutine has been added to FISPACT2001. The subroutine is called <strong>sp</strong>all.f. The<br />
other modifications to the original code only concern the following subroutines:<br />
fmain.f<br />
col069.f<br />
col100.f<br />
col172.f<br />
col175.f<br />
col315.f<br />
The subroutine fmain.f has been modified in order to add five new code words: CURRENT,<br />
TARGET, ATTARGET, HCOLLAPSE, TCOLLAPSE. All of them concern the<br />
preliminary input.<br />
In FISPACT2001 the subroutines col*.f are called by fmain.f (code word COLLAPSE) in<br />
order to create the file COLLAPX according to the desired energy group structure. The<br />
subroutine <strong>sp</strong>all.f is called by the subroutines col*.f (code words TCOLLAPSE and<br />
HCOLLAPSE) for calculations of the pseudo cross-sections for interactions over 20MeV.<br />
Subroutine <strong>sp</strong>all.f<br />
SUBROUTINE SPALL<br />
C C.Petrovich - <strong>ENEA</strong><br />
C***************************************************************<br />
C SUBROUTINE FOR CALCULATIONS >20MeV.<br />
C IT USES OUTPUTS PRODUCED BY MCNPX 2.1.5 AND HTAPE3X<br />
C***************************************************************<br />
C<br />
CHARACTER*80 FILES<br />
COMMON /FILES1/ FILES(40)<br />
COMMON /SPAL1/ CPFLUX,CURR,IDTARG,ATTARG<br />
COMMON /SPAL3/ NP2(15000),ND2(15000),ZSEKT(15000),MCP2<br />
DOUBLE PRECISION HISOT(0:300,0:300),ELAS(0:300,0:300),EL,HIS,<br />
* HISOTM(0:300,0:300)<br />
PARAMETER (IS=3000)<br />
INTEGER IZ(IS),IN(IS),IA(IS),IZE(IS),INE(IS),<br />
* IZF(IS),INF(IS),IAF(IS),IZM(IS),INM(IS),IAM(IS),<br />
* ID(IS),IDF(0:300,0:300),IDFM(0:300,0:300)<br />
CHARACTER B*6,C*5,D*3<br />
c Protons per second in 1mA<br />
PROTSC=6.24146E15<br />
17