Re: source error

Hi Mojdeh,

there were a couple of small problems in the input
1. BEAMPOS, unfortunately the sdum for POSITIVE should be empty. flair
was wrongly treating that option
2. MAT-PROP for DPA-ENE doesn't work with the present release version.
It will work with the next FLUKA release

however the real problem causing the crash was coming from the source
routine.

In the CUMENE array which has the normalized to one cumulative
probability is starting from the value of 2.33568E-1 up to 1.0. However
once you select a random number and you scan the array, you search for
the interval between CUMENE(K-1) and CUMENE(K), therefore if you get a
random number XI less than 2.33568E-1 then it exits the loop without
finding the interval. This was happening exactly on the first particle
and the ENERGY variable after the loop had the value 0.0 which was
latter creating a crash in the KASKAD routine.

Please find attached the modified routine that works. I added an extra
value in the CUMENE (the zero) and perform the search from K=0,99
checking the interval CUMENE(K) to CUMENE(K+1)

Best Regards
Vasilis

On 08/06/10 08:27, mojdeh najm wrote:
> Dear fluka experts
> I have a problem when run inp with a source and i don't know it's
> reason. I attached inp and source file. please help me.
> Regards
> Mojdeh Anjom

-- 
-------------------------------------------------------
Vasilis Vlachoudis
Dep. EN-STI-EET
Web: home.cern.ch/bnv
Tel: +41-22-7679851
Fax: +41-22-7669644
-------------------------------------------------------