23.12.2012 Views

sp-fispact2001 - ENEA AFS Cell

sp-fispact2001 - ENEA AFS Cell

sp-fispact2001 - ENEA AFS Cell

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!