Hi,
I would like to make a collision tape for particle fragments in order to
calculate differential and double differential cross sections.
Using
ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )
IF ( .NOT. LFCOPE ) THEN
LFCOPE = .TRUE.
OPEN ( UNIT = 28, FILE = "Li-6_production.dat",
$ STATUS = 'unknown', FORM = 'FORMATTED' )
end if
call geor2n ( mreg, mrgnam, ierr1 )
*
if (mrgnam .eq. "regC4" .and. Icres .eq. 3 .and.
$ Ibres .eq. 6 .and. jtrack .ge. -6) then
write(28,'(i3,5(1p,1g25.15))' ) jtrack, etrack -am (jtrack),
$ Tkeflk, Txflk, Tyflk, Tzflk
endif
*
do iii = 1, Np
if (mrgnam .eq. "region" .and. Icheav(iii) .eq. 3 .and.
$ Ibheav(iii) .eq. 6) then
call USRDCI(jtrack, IONA, IONZ, IONM)
write(28,'(i3,1p,1g25.15,1p,2i3,4(1p,1g25.15))') jtrack,
$ etrack - am (jtrack), IONZ, IONA, Tkeflk, Txflk, Tyflk, Tzflk
endif
enddo
for Li-6 production, the output file always gets very large (~ 150 MB)
and except the first four lines the data do not make any sense (-1E+30,
2.0 or -1.0 in every column and line). What can I do in order to solve
this problem?
Best regards,
David Ellmeyer
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Fri May 13 2016 - 01:09:45 CEST