RE: More FLUKA Crashes

From: Harine, Katherine (Katherine.Harine@msfc.nasa.gov)
Date: Fri Sep 23 2005 - 20:40:12 CEST

  • Next message: Giuseppe Battistoni: "corrections to fluka2005.6"

    This is a multi-part message in MIME format.

    ------_=_NextPart_001_01C5C06E.3B6BA2F4
    Content-Type: text/plain;
            charset="iso-8859-1"
    Content-Transfer-Encoding: quoted-printable

    Alfredo,

    The problem appears to be in the index of the term: AM.

    The source.f file that is distributed gives the following lines:

    * Kinetic energy of the particle (GeV)
          TKEFLK (NPFLKA) =3D SQRT ( PBEAM**2 + AM (IONID)**2 ) - AM (IONID)
    * Particle momentum
          PMOFLK (NPFLKA) =3D PBEAM
    * PMOFLK (NPFLKA) =3D SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
    * & + TWOTWO * AM (ILOFLK(NPFLKA)) ) )

    I read in the energies in GeV per nucleaon so I do this for heavy ions:

    * Kinetic energy of the particle (GeV)
    * TKEFLK (NPFLKA) =3D SQRT ( PBEAM**2 + AM (IJBEAM)**2 ) - AM =
    (IJBEAM)
          I =3D 0 =
     KJH
          XFLRNDM =3D FLRNDM(XXX) =
     KJH
      301 I =3D I + 1 =
     KJH
          IF (XNRMINFL(I) - XFLRNDM) 201, 202, 202 =
    KJH
      201 GOTO 301 =
    KJH
      202 TKEFLK (NPFLKA) =3D XENERMID(I) * AM (IONID) =
     KJH
    * Particle momentum
    * PMOFLK (NPFLKA) =3D PBEAM
          PMOFLK (NPFLKA) =3D SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
         & + TWOTWO * AM (IONID) ) )

    Please excuse the commented lines. I left them in from the distributed =
    file.

    Evidently this does not work for particles other than heavy ions. =20

    To get successful runs with PROTONS and 4-HELIUM et cetera I put the =
    following in:

    * Kinetic energy of the particle (GeV)
    * TKEFLK (NPFLKA) =3D SQRT ( PBEAM**2 + AM (IJBEAM)**2 ) - AM =
    (IJBEAM)
          I =3D 0 =
     KJH
          XFLRNDM =3D FLRNDM(XXX) =
     KJH
      301 I =3D I + 1 =
     KJH
          IF (XNRMINFL(I) - XFLRNDM) 201, 202, 202 =
    KJH
      201 GOTO 301 =
    KJH
      202 TKEFLK (NPFLKA) =3D XENERMID(I) * AM (IJBEAM) =
     KJH
    * Particle momentum
    * PMOFLK (NPFLKA) =3D PBEAM
          PMOFLK (NPFLKA) =3D SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
         & + TWOTWO * AM (ILOFLK(NPFLKA)) ) )

    then it runs for PROTON et cetera, but crashes for heavy ions.

    I sent an email about this a couple of days ago and perhaps didn't =
    explain the problem. Exactly what do IJBEAM, ILOFLK(NPFLKA), and IONID =
    mean and when should each be used as an index of AM?

    - Katherine

    -----Original Message-----
    From: Alfredo Ferrari [mailto:alfredo.ferrari@cern.ch]
    Sent: Fri 9/23/2005 12:19 PM
    To: Harine, Katherine
    Cc: giuseppe.battistoni@mi.infn.it; alfredo.ferrari@cern.ch; =
    thomas.l.wilson@nasa.gov; brandon.d.reddell@boeing.com
    Subject: Re: More FLUKA Crashes
    =20
    Hi Katherine

    please send the messages to fluka-discuss, so that they are helpful for=20
    everybody. From the err file it appears that a proton is generated with=20
    zero energy in vacuum. Do you use a source routine? And in case are your =

    sure you are not generating protons with zero energy??

                Ciao
               Alfredo

    On Fri, 23 Sep 2005, Harine, Katherine wrote:

    >
    > Attached are the .err and .log files for two runs that appear to crash =
    on the first particle transported. This is strange because most runs =
    ran fine for 10 cycles of 1000 particles with heavy ions. However, when =
    I ran PROTON as the particle or 4-HELIUM, the runs crashed at what =
    appears to be the first particle transported.
    >
    > These runs were on a 32-bit Red Hat Linux box.
    >
    > - Katherine
    >
    >
    >
    >
    >

    --=20

    +------------------------------------------------------------------------=
    ----+
    | Alfredo Ferrari || Tel.: +41.22.767.6119 =
        |
    | CERN-AB || Fax.: +41.22.767.7555 =
        |
    | 1211 Geneva 23 || e-mail: Alfredo.Ferrari@cern.ch =
        |
    | Switzerland || Alfredo.Ferrari@mi.infn.it =
        |
    +------------------------------------------------------------------------=
    ----+

    ------_=_NextPart_001_01C5C06E.3B6BA2F4
    Content-Type: text/html;
            charset="iso-8859-1"
    Content-Transfer-Encoding: quoted-printable

    <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
    <HTML>
    <HEAD>
    <META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; =
    charset=3Diso-8859-1">
    <META NAME=3D"Generator" CONTENT=3D"MS Exchange Server version =
    6.5.7232.7">
    <TITLE>RE: More FLUKA Crashes</TITLE>
    </HEAD>
    <BODY>
    <!-- Converted from text/plain format -->

    <P><FONT SIZE=3D2>Alfredo,<BR>
    <BR>
    The problem appears to be in the index of the term:&nbsp; AM.<BR>
    <BR>
    The source.f file that is distributed gives the following lines:<BR>
    <BR>
    <BR>
    * Kinetic energy of the particle (GeV)<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; TKEFLK (NPFLKA) =3D SQRT ( PBEAM**2 + AM =
    (IONID)**2 ) - AM (IONID)<BR>
    * Particle momentum<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; PMOFLK (NPFLKA) =3D PBEAM<BR>
    *&nbsp;&nbsp;&nbsp;&nbsp; PMOFLK (NPFLKA) =3D SQRT ( TKEFLK (NPFLKA) * ( =
    TKEFLK (NPFLKA)<BR>
    *&nbsp;&nbsp;&nbsp; =
    &amp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; + =
    TWOTWO * AM (ILOFLK(NPFLKA)) ) )<BR>
    <BR>
    I read in the energies in GeV per nucleaon so I do this for heavy =
    ions:<BR>
    <BR>
    * Kinetic energy of the particle (GeV)<BR>
    *&nbsp;&nbsp;&nbsp;&nbsp; TKEFLK (NPFLKA) =3D SQRT ( PBEAM**2 + AM =
    (IJBEAM)**2 ) - AM (IJBEAM)<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I =3D =
    0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
    nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
    sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
    KJH<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; XFLRNDM =3D =
    FLRNDM(XXX)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
    sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
    p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; KJH<BR>
    &nbsp; 301 I =3D I + =
    1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
    nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
    sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; KJH<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; IF (XNRMINFL(I) - XFLRNDM)&nbsp; 201, =
    202, =
    202&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
    p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
    ; KJH<BR>
    &nbsp; 201 GOTO =
    301&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
    p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
    ;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
    nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; KJH<BR>
    &nbsp; 202 TKEFLK (NPFLKA) =3D XENERMID(I) * AM =
    (IONID)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
    KJH<BR>
    * Particle momentum<BR>
    *&nbsp;&nbsp;&nbsp;&nbsp; PMOFLK (NPFLKA) =3D PBEAM<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; PMOFLK (NPFLKA) =3D SQRT ( TKEFLK =
    (NPFLKA) * ( TKEFLK (NPFLKA)<BR>
    &nbsp;&nbsp;&nbsp;&nbsp; =
    &amp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; + =
    TWOTWO * AM (IONID) ) )<BR>
    <BR>
    Please excuse the commented lines.&nbsp; I left them in from the =
    distributed file.<BR>
    <BR>
    Evidently this does not work for particles other than heavy =
    ions.&nbsp;<BR>
    <BR>
    To get successful runs with PROTONS and 4-HELIUM et cetera I put the =
    following in:<BR>
    <BR>
    * Kinetic energy of the particle (GeV)<BR>
    *&nbsp;&nbsp;&nbsp;&nbsp; TKEFLK (NPFLKA) =3D SQRT ( PBEAM**2 + AM =
    (IJBEAM)**2 ) - AM (IJBEAM)<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I =3D =
    0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
    nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
    sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
    KJH<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; XFLRNDM =3D =
    FLRNDM(XXX)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
    sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
    p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; KJH<BR>
    &nbsp; 301 I =3D I + =
    1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
    nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
    sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; KJH<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; IF (XNRMINFL(I) - XFLRNDM)&nbsp; 201, =
    202, =
    202&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
    p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
    ; KJH<BR>
    &nbsp; 201 GOTO =
    301&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
    p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
    ;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
    nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; KJH<BR>
    &nbsp; 202 TKEFLK (NPFLKA) =3D XENERMID(I) * AM =
    (IJBEAM)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
    ;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
    KJH<BR>
    * Particle momentum<BR>
    *&nbsp;&nbsp;&nbsp;&nbsp; PMOFLK (NPFLKA) =3D PBEAM<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; PMOFLK (NPFLKA) =3D SQRT ( TKEFLK =
    (NPFLKA) * ( TKEFLK (NPFLKA)<BR>
    &nbsp;&nbsp;&nbsp;&nbsp; =
    &amp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; + =
    TWOTWO * AM (ILOFLK(NPFLKA)) ) )<BR>
    <BR>
    then it runs for PROTON et cetera, but crashes for heavy ions.<BR>
    <BR>
    I sent an email about this a couple of days ago and perhaps didn't =
    explain the problem.&nbsp; Exactly what do IJBEAM, ILOFLK(NPFLKA), and =
    IONID mean and when should each be used as an index of AM?<BR>
    <BR>
    - Katherine<BR>
    <BR>
    <BR>
    <BR>
    -----Original Message-----<BR>
    From: Alfredo Ferrari [<A =
    HREF=3D"mailto:alfredo.ferrari@cern.ch">mailto:alfredo.ferrari@cern.ch</A=
    >]<BR>
    Sent: Fri 9/23/2005 12:19 PM<BR>
    To: Harine, Katherine<BR>
    Cc: giuseppe.battistoni@mi.infn.it; alfredo.ferrari@cern.ch; =
    thomas.l.wilson@nasa.gov; brandon.d.reddell@boeing.com<BR>
    Subject: Re: More FLUKA Crashes<BR>
    <BR>
    Hi Katherine<BR>
    <BR>
    please send the messages to fluka-discuss, so that they are helpful =
    for<BR>
    everybody. From the err file it appears that a proton is generated =
    with<BR>
    zero energy in vacuum. Do you use a source routine? And in case are =
    your<BR>
    sure you are not generating protons with zero energy??<BR>
    <BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
    Ciao<BR>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Alfredo<BR>
    <BR>
    On Fri, 23 Sep 2005, Harine, Katherine wrote:<BR>
    <BR>
    &gt;<BR>
    &gt; Attached are the .err and .log files for two runs that appear to =
    crash on the first particle transported.&nbsp; This is strange because =
    most runs ran fine for 10 cycles of 1000 particles with heavy =
    ions.&nbsp; However, when I ran PROTON as the particle or 4-HELIUM, the =
    runs crashed at what appears to be the first particle transported.<BR>
    &gt;<BR>
    &gt; These runs were on a 32-bit Red Hat Linux box.<BR>
    &gt;<BR>
    &gt; - Katherine<BR>
    &gt;<BR>
    &gt;<BR>
    &gt;<BR>
    &gt;<BR>
    &gt;<BR>
    <BR>
    --<BR>
    <BR>
    +------------------------------------------------------------------------=
    ----+<BR>
    |&nbsp; Alfredo =
    Ferrari&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp; ||&nbsp; Tel.: =
    +41.22.767.6119&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
    p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; |<BR>
    |&nbsp; =
    CERN-AB&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
    ||&nbsp; Fax.: =
    +41.22.767.7555&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
    p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; |<BR>
    |&nbsp; 1211 Geneva =
    23&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
    ;&nbsp;&nbsp;&nbsp;&nbsp; ||&nbsp; e-mail: =
    Alfredo.Ferrari@cern.ch&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; |<BR>
    |&nbsp; =
    Switzerland&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
    bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
    ||&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
    Alfredo.Ferrari@mi.infn.it&nbsp;&nbsp;&nbsp;&nbsp; |<BR>
    +------------------------------------------------------------------------=
    ----+<BR>
    <BR>
    <BR>
    </FONT>
    </P>

    </BODY>
    </HTML>
    ------_=_NextPart_001_01C5C06E.3B6BA2F4--

  • Next message: Giuseppe Battistoni: "corrections to fluka2005.6"

    This archive was generated by hypermail 2.1.6 : Sat Sep 24 2005 - 09:06:50 CEST