Floating Point Exception Error when using Source.f

From: <Lewis.Macfarlane_at_nuclear.co.uk>
Date: Sat, 11 May 2013 10:20:50 +0200

Dear FLUKA experts,

I have some issues trying to run a simple (test) proton case with a source.f subroutine.
The same subroutine has been used on earlier versions of Fluka/Flair without any issues. Current versions are:

Fluka: 2011.2 (Feb. 2013)
Flair: 1.1-1 (March 2013)

The input file runs fine when source.f file is not attached, i.e. with a mono-energetic source.
The source.f file allows the use of a simple rectangular spectrum with a conical beam profile defined by a divergence half-angle.
    When used with the source.f file, the code terminates with the error:

*** Vacuum stopping: Ij, Pla, Ekin

Within the log file, is:

Program received signal 8 (SIGFPE): Floating-point exception.

Backtrace for this error:
     + /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fbcdb71a4a0]
     + function kaskad_ (0x501D2C)
       at line 4139 of file kaskad.f
     + function feeder_ (0x4F67B6)
       at line 257 of file feeder.f
     + function flukam_ (0x40D055)
       at line 3095 of file flukam.f
     + function fluka (0x40785C)
       at line 289 of file fluka.f
     + /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7fbcdb70576d]

I have attached my input file and source.f for your consideration. Any help
    would be greatly appreciated,

Many thanks,
Lewis MacFarlane.

Received on Sat May 11 2013 - 10:57:36 CEST

This archive was generated by hypermail 2.3.0 : Sat May 11 2013 - 10:57:40 CEST