Problem with Heavy-Ion-Fluka

From: Gutermuth Frank (F.Gutermuth@gsi.de)
Date: Mon Dec 01 2003 - 13:20:58 CET

  • Next message: Gutermuth Frank: "Problem with Heavy-Ion-Fluka"

    Hello everybody,

    especially to the organizers and participants of the Fluka-Course in
    Frascati! In the meantime we tried to use Fluka with the heavy ion option
    (DPMJET) to simulate the neutron production by a uranium beam (at GSI we
    have really "heavy ions" ;-)). For the planned extension of the GSI
    accelerator we would like to use Fluka for the shielding design. Uranium
    shall be used as primary beam with energies up to approximately 35 GeV/u.

    We started with a very simple geometry consisting of an iron target (box
    shape 10cm X 10cm X 20cm). Everything worked smoothly as long as the energy
    of the primary uranium particles stays below approximately 5 GeV/u. With
    higher energies Fluka crashes and we get the following message in the log
    file:

    >Subscript out of range on file line 284, procedure ficonf.f/ficonf.
    >Attempt to access the 211-th element of variable nexc[subscript-2-of-2].

    We used "NEW-DEFAULTS" both with EMF switched on and switched off. Maybe
    some of the specialists/authors of the RQMD/DPMJET-extension of Fluka can
    comment on this?

    Thanks in advance

    Frank

    P.S. See our input file in the attachment.

    ------_=_NextPart_000_01C3B805.93872240
    Content-Type: application/octet-stream;
            name="uran10.inp"
    Content-Transfer-Encoding: quoted-printable
    Content-Disposition: attachment;
            filename="uran10.inp"

    TITLE=0A=
     28-11-2003: 10X10X20CM**3 IRON TARGET=0A=
    DEFAULTS =
    NEW-DEFA=0A=
    EMF =
    EMF-OFF=0A=
    *23456789 123456789 123456789 123456789 123456789 123456789 123456789 =
    123456789=0A=
    BEAM -10.0 0.0 0.0 0.0 0.0 1.0 =
    HEAVYION=0A=
    HI-PROPE 92.0 238.0=0A=
    BEAMPOS 0.0 0.0 -100.0=0A=
    EVENTYPE 2.0 =
    DPMJET=0A=
    GEOBEGIN =
    COMBINAT=0A=
              target box=0A=
    *23456789 123456789 123456789 123456789 123456789 123456789 123456789 =
    123456789=0A=
      SPH 1 0.0 0.0 0.0 1000.0=0A=
      SPH 2 0.0 0.0 0.0 990.0=0A=
      RPP 3 -5.0 5.0 -5.0 5.0 -10.0 =
    10.0=0A=
      END=0A=
    *234x6789x123456x89 123x56789 x234567x9 1234x6789 1x345678x 12345x789 =
    12x stop at 73=0A=
    * black hole=0A=
        1 5 1 -2=0A=
    * vacuum=0A=
        2 5 2 -3=0A=
    * target=0A=
        3 5 3=0A=
      END=0A=
    GEOEND=0A=
    *=0A=
    *234x6789x123456x89 123x56789 x234567x9 1234x6789 1x345678x 12345x789 =
    12x stop at 73=0A=
    * black hole=0A=
    ASSIGNMAT 1.0 1.0=0A=
    * vacuum=0A=
    ASSIGNMAT 2.0 2.0=0A=
    * target=0A=
    ASSIGNMAT 11.0 3.0=0A=
    *=0A=
    *23456789 123456789 123456789 123456789 123456789 123456789 123456789 =
    123456789=0A=
    *PLOTGEOM 0.0 1.0 0.0 2000.0 0.0 0.0 =
    FORMAT=0A=
    * sth=0A=
    * -20.0 -20.0 0.0 20.0 20.0 0.0=0A=
    * 0.0 1.0 0.0 1.0 0.0 0.0=0A=
    * 1.0 1.0 1.0=0A=
    *=0A=
    *=0A=
    *=0A=
    *=0A=
    *scoring of neutron yield=0A=
    USRYIELD 2399.0 8.0 50.0 3.0 2.0 1.0 =
    highyld=0A=
    USRYIELD 2.0196 0.0196 50.0 1.0 0.0 3.0 =
     &=0A=
    USRYIELD 2399.0 8.0 51.0 3.0 2.0 1.0 =
    highyld=0A=
    USRYIELD 2.0196 0.0196 50.0 8.0 7.0 3.0 =
     &=0A=
    USRYIELD 2399.0 8.0 52.0 3.0 2.0 1.0 =
    highyld=0A=
    USRYIELD 2.0196 0.0196 50.0 15.5 14.5 3.0 =
     &=0A=
    USRYIELD 2399.0 8.0 53.0 3.0 2.0 1.0 =
    highyld=0A=
    USRYIELD 2.0196 0.0196 50.0 31.0 29.0 3.0 =
     &=0A=
    USRYIELD 2399.0 8.0 54.0 3.0 2.0 1.0 =
    highyld=0A=
    USRYIELD 2.0196 0.0196 50.0 51.0 49.0 3.0 =
     &=0A=
    USRYIELD 2399.0 8.0 55.0 3.0 2.0 1.0 =
    highyld=0A=
    USRYIELD 2.0196 0.0196 50.0 92.0 88.0 3.0 =
     &=0A=
    *=0A=
    RANDOMIZE 1.0=0A=
    *=0A=
    *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 =
    *23456789 =0A=
    START 5.E+1=0A=
    STOP=0A=

    ------_=_NextPart_000_01C3B805.93872240--


  • Next message: Gutermuth Frank: "Problem with Heavy-Ion-Fluka"

    This archive was generated by hypermail 2.1.6 : Mon Dec 01 2003 - 13:38:53 CET