Dear David,
from just your piece of code I try to make a couple of comments
(it would be more efficient if you had provided the input file — even a simplified version — and the user routines).
1. Did you read what it is reported in the manual section about USDRAW entry and
how to access secondary particle info and, in particular, the exceptions?
2. I suggest you to check the do-loop: the loop is over NP (from GENSTK common) but then
there is an if-statement against Icheav(iii) Ibheav(iii) (from FHEAVY common), and finally
you dump TKEFLK and following variables from FLKSTK common.
I would not be surprised if you are accessing meaningless memory arrays.
Hope this could help.
Regards, luigi
On 12 May 2016, at 23:25, David Ellmeyer <e1126053_at_student.tuwien.ac.at<mailto:e1126053_at_student.tuwien.ac.at>> wrote:
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
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Fri May 13 2016 - 13:10:32 CEST