Re: [fluka-discuss]: Collision tape for fragments

From: Luigi Salvatore Esposito <luigi.salvatore.esposito_at_cern.ch>
Date: Fri, 13 May 2016 09:30:54 +0000

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

This archive was generated by hypermail 2.3.0 : Fri May 13 2016 - 13:10:34 CEST