[fluka-discuss]: Collision tape for fragments

From: David Ellmeyer <e1126053_at_student.tuwien.ac.at>
Date: Thu, 12 May 2016 23:25:29 +0200

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

This archive was generated by hypermail 2.3.0 : Fri May 13 2016 - 01:09:49 CEST