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

From: Stefan Roesler <sroesler_at_mail.cern.ch>
Date: Fri, 8 Feb 2013 18:49:58 +0100

Hi Sam

Yes, perhaps redefine region SUR and avoid (too many) parenthesis..

Cheers
Stefan


On Fri, 8 Feb 2013, Sam Tygier wrote:

> Hi Stefan,
>
> Do you mean I need to find a simpler expression of the SUR region into my
> inputfile, or that I should be using the WHAT(5) option in GEOBEGIN? or is it
> the GEOBEGIN optimisation that has expanded my expression from 31 to 100k
> terms?
>
> Thanks,
> Sam
>
> On 06/02/13 15:46, Stefan Roesler wrote:
>> 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 Mon Feb 11 2013 - 10:02:15 CET