Re: Index '100001' of dimension 1 of array 'tx' above upper bound of 100000

From: Stefan Roesler <sroesler_at_mail.cern.ch>
Date: Wed, 6 Feb 2013 16:46:56 +0100

Dear Sam

I am afraid that you have to optimize your parenthesis. The code arrived
at 100k terms in the expansion that would slow down a lot the calculation
if the array dimension was increased.

There is a protection missing in the code that may need to be fixed. But
for your case it would only give you a 'stop' instead of an overflow...

Cheers
Stefan


On Tue, 5 Feb 2013, Sam Tygier wrote:

> Hi,
>
> I am hitting the error:
> At line 84 of file rpnorm.f
> Fortran runtime error: Index '100001' of dimension 1 of array 'tx' above
> upper bound of 100000
> With a fairly simple geometry in FLUKA2011 Version 2.17 Dec-12. I have
> striped out most of my geometry to just the problematic elements.
>
> I have 7 rings/coils, defined using 14 planes, and 8 cylinders (ZCC). A
> vacuum region surrounds them. This surrounding region needs to have the
> 7 rings cut out of it, so it has a couple of lines in its definition.
> This still seems to me like a simple region, and the 100001 vs 100000
> would make me suspicious of an out-by-one error.
>
> The 'Handling of errors and crashes' tutorial suggests that infinite
> bodies are are problem, but the fluka documentation recommends infinite
> bodies over finite ones.
>
> What should I do? Thanks,
>
> Sam
>
> ===============
>
> TITLE
> piontarget
> FREE
> GLOBAL, 10000, , , 3 , 1 ,
> DEFAULTS , , , , , , ,NEW-DEFA
> BEAM, -8, 0, 0, 0, 0, 0, PROTON
> BEAMPOS, -15.643446504, 0, 98.7688340595, 0.15643446504, 0, , NEGATIVE
> FREE
> GEOBEGIN , , , , , , , COMBNAME
> 0, 0, pt
> XYP r0pzs -95.0
> XYP r0pze -77.5
> XYP r1pzs -72.5
> XYP r1pze 62.5
> XYP r2pzs 67.5
> XYP r2pze 247.5
> XYP r3pzs 254.5
> XYP r3pze 292.5
> XYP r4pzs 284.5
> XYP r4pze 305.5
> XYP r5pzs 318.0
> XYP r5pze 357.0
> XYP r6pzs 367.0
> XYP r6pze 398.0
> ZCC zcc74_2 0 0 74.2
> ZCC zcc66_2 0 0 66.2
> ZCC zcc29_8 0 0 29.8
> ZCC zcc26_6 0 0 26.6
> ZCC zcc25_0 0 0 25.0
> ZCC zcc80_6 0 0 80.6
> ZCC zcc28_2 0 0 28.2
> ZCC zcc79_0 0 0 79.0
> ZCC sur 0 0 100
> XYP surps -105
> XYP surpe 408.0
> SPH bbs 0.0 0.0 0.0 1000
> END
> R0 5 -r0pzs +r0pze -zcc66_2 +zcc80_6
> R1 5 -r1pzs +r1pze -zcc66_2 +zcc80_6
> R2 5 -r2pzs +r2pze -zcc66_2 +zcc74_2
> R3 5 -r3pzs +r3pze -zcc66_2 +zcc79_0
> R4 5 -r4pzs +r4pze -zcc25_0 +zcc26_6
> R5 5 -r5pzs +r5pze -zcc25_0 +zcc28_2
> R6 5 -r6pzs +r6pze -zcc25_0 +zcc29_8
> SUR 5 - surps + surpe + sur -( -r0pzs +r0pze -zcc66_2 +zcc80_6) -(
> -r1pzs +r1pze -zcc66_2 +zcc80_6) -( -r2pzs +r2pze -zcc66_2 +zcc74_2) -(
> -r3pzs +r3pze -zcc66_2 +zcc79_0) -( -r4pzs +r4pze -zcc25_0 +zcc26_6) -(
> -r5pzs +r5pze -zcc25_0 +zcc28_2) -( -r6pzs +r6pze -zcc25_0 +zcc29_8)
> BB 5 +bbs - ( - surps + surpe + sur)
> END
> GEOEND
> ASSIGNMAT, COPPER, R0, , ,1
> ASSIGNMAT, COPPER, R1, , ,1
> ASSIGNMAT, COPPER, R2, , ,1
> ASSIGNMAT, COPPER, R3, , ,1
> ASSIGNMAT, COPPER, R4, , ,1
> ASSIGNMAT, COPPER, R5, , ,1
> ASSIGNMAT, COPPER, R6, , ,1
> ASSIGNMAT, VACUUM, SUR, , ,1
> ASSIGNMAT, BLCKHOLE, BB, , ,0
> SCORE, 208
> RANDOMIZE, 1, 0
> START , 10.0
> STOP
>
>
>
Received on Wed Feb 06 2013 - 18:00:56 CET