Re: Joint calculation with source file

Dear Ogawa,

You have defined HEAVYION as primary particle. Thus IJBEAM=-2 and the
source will never enter the branch where you set

         ILOFLK (NPFLKA) = jpart(sample)

Furthermore, you have modified in the default source the index of AM in

your source:
      PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
     & + TWOTWO * AM (ILOFLK(NPFLKA)) ) )

default:
      PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
     & + TWOTWO * AM (IONID) ) )

In case of a HEAVYION primary ILOFLK(NPFLKA) contains the identity of the
ion, in your case carbon (601230) and IONID=-2. Thus, using ILOFLK(NPFLKA)
to index AM causes a crash, while IONID would not.

I suggest to set your beam particle to a hadron of the highest energy
of any particle in your data file for your simulation.

Best regards
Stefan

On Mon, 11 Jan 2010, Ogawa Tatsuhiko wrote:

> Dear FLUKA experst,
>
> Please find the attached input file. This is the input file I used actually. I
> also attached the source.f file again for convenience.
>
> I would appreciate if you could point out what is wrong in this input.
>
> Sincerely,
> Ogawa
>
>

-- 
___________________________________
Stefan Roesler
CERN, DGS-RP
CH-1211 Geneva 23
Switzerland
Phone:  +41-22-7679891
Fax:    +41-22-7669639
E-mail: Stefan.Roesler_at_cern.ch