Hello,
I am using TRITON (t6-depl currently) to deplete an assembly containing TRISO fuel. I had believed to have debugged the model successfully, but for a certain input combination I have a model that appears to work fine, but for one set of input parameters it runs fine, and for another if fails.
The parameters that it runs for are as follows:
read parameter
gen=350
npg=1000
nsk=35
htm=yes
nub=no
cfx=yes
end parameter
The parameters if fails for is as follows:
read parameter
gen=350
npg=5000
nsk=35
htm=yes
nub=no
cfx=yes
end parameter
I am running on WINDOWS (and will check whether it occurs with a LINUX install) and the error message is as follows:
"attempt to access nonexistant record, unit 90"
It appears the line of code is line 145 of direct_access_M.f90 called from rerite_I.f90 called from concupdate.f90 called from triton6.f90. Which leads me to believe this occurs during update of nuclide concentrations. Is there some input causing this difficulty? I have the input below.
'Input generated by GeeWiz SCALE 6.0.13.04 Compiled on January 4, 2010
=t6-depl parm=(centrm)
oakridge deep burn
v7-238
read composition
wtptuo2-4pct 1 10.176 3
8016 11.85
92235 3.526
92238 84.624
1 900 end
wtptth 2 9.6 2
90232 86.21
8016 13.79
1 900 end
zr 4 den=6.55 1 900 end
h2o 5 den=0.7 1 600 end
b 5 0 4.012417e-05 600 end
wtpttru 6 9.461 10
8016 10.78666
92235 0.001242
92238 0.17433
93237 4.40183
94238 2.669427
94239 51.85101
94240 19.60122
94241 3.235999
94242 5.89453
95241 1.383744
1 900 end
c-graphite 7 den=1.05 1 900 end
wtptsic 8 3.18 2
14000 70.05
6000 29.95
1 900 end
wtptsic 9 3.18 2
14000 70.05
6000 29.95
1 900 end
wtptipyc 10 1.9 2
6000 86.87672
1001 13.12328
1 900 end
wtptopyc 11 1.9 2
6000 86.87672
1001 13.12328
1 900 end
h2o 13 den=0.7 1 600 end
b 13 0 4.012417e-05 600 end
zr 14 den=6.55 1 900 end
c 16 0 0.0952621 900 end
c 17 0 0.0952621 900 end
end composition
read celldata
latticecell squarepitch fuelr=0.4095 1 gapr=0.417 0 cladr=0.475 4 hpitch=0.643 5 end
doublehet fuelmix=18 end
gfr=0.0175 6 coatr=0.0275 7 coatr=0.031 16 coatr=0.0345 9 coatr=0.0385 17 vf=0.33 matrix=8 end grain
rod
squarepitch fuelr=0.4095 cladr=0.475 14 hpitch=0.643 13 fuelh=381 right_bdy=white left_bdy=reflected end
end celldata
read depletion
-1 -6
end depletion
read burndata
power=44.44 burn=2 end
power=44.44 burn=10 end
power=44.44 burn=15 end
power=44.44 burn=30 end
power=44.44 burn=1045 nlib=26 end
end burndata
read model
read parameter
gen=350
npg=5000
nsk=35
htm=yes
nub=no
cfx=yes
end parameter
read geometry
global unit 1
com='assembly box'
cuboid 30 20.576 0 20.576 0 427.02 0
array 3 30 place 1 1 1 0.643 0.643 427.02
boundary 30
unit 2
com='fuel pin'
cylinder 1 0.418 0 -427.02
cylinder 2 0.475 0 -427.02
cuboid 3 0.643 -0.643 0.643 -0.643 0 -427.02
media 1 1 1 vol=234.2316
media 4 1 -1 2 vol=28.23685
media 5 1 3 -2 vol=50.38546
boundary 3
unit 4
com='water'
cylinder 2 0.418 0 -427.02
cylinder 3 0.475 0 -427.02
cuboid 4 0.643 -0.643 0.643 -0.643 0 -427.02
media 5 1 2 vol=234.2316
media 4 1 3 -2 vol=68.23685
media 5 1 4 -3 vol=50.38546
boundary 4
unit 5
com='sic-tru'
cylinder 1 0.418 0 -427.02
cylinder 2 0.475 0 -427.02
cuboid 3 0.643 -0.643 0.643 -0.643 0 -427.02
media 18 1 1 vol=234.2316
media 14 1 -1 2 vol=28.23685
media 13 1 3 -2 vol=50.38546
boundary 3
end geometry
read array
ara=2 nux=16 nuy=16 nuz=1 typ=square
com='config-b'
fill
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 5 4 4 5 2 2 2 2 5 4 4 5 2 2
2 2 5 4 4 5 2 2 2 2 5 4 4 5 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 2 2 2 2 2 5 5 2 2 2 2 2 2 2
2 2 2 2 2 2 5 4 4 5 2 2 2 2 2 2
2 2 2 2 2 2 5 4 4 5 2 2 2 2 2 2
2 2 2 2 2 2 2 5 5 2 2 2 2 2 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 5 4 4 5 2 2 2 2 5 4 4 5 2 2
2 2 5 4 4 5 2 2 2 2 5 4 4 5 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 end fill
ara=3 nux=16 nuy=16 nuz=1 typ=square
com='cfg-c'
fill
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 5 4 4 5 2 2 2 2 5 4 4 5 2 2
2 2 5 4 4 5 2 2 2 2 5 4 4 5 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 2 2 2 2 5 5 5 5 2 2 2 2 2 2
2 2 2 2 2 2 5 4 4 5 2 2 2 2 2 2
2 2 2 2 2 2 5 4 4 5 2 2 2 2 2 2
2 2 2 2 2 2 5 5 5 5 2 2 2 2 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 5 4 4 5 2 2 2 2 5 4 4 5 2 2
2 2 5 4 4 5 2 2 2 2 5 4 4 5 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 end fill
ara=4 nux=16 nuy=16 nuz=1 typ=square
com='cfg-d'
fill
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 5 2 5 5 2 5 2 2 5 2 5 5 2 5 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 5 2 4 4 2 5 2 2 5 2 4 4 2 5 2
2 5 2 4 4 2 5 2 2 5 2 4 4 2 5 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 5 2 5 5 2 5 2 2 5 2 5 5 2 5 2
2 2 2 2 2 2 2 4 4 2 2 2 2 2 2 2
2 2 2 2 2 2 2 4 4 2 2 2 2 2 2 2
2 5 2 5 5 2 5 2 2 5 2 5 5 2 5 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 5 2 4 4 2 5 2 2 5 2 4 4 2 5 2
2 5 2 4 4 2 5 2 2 5 2 4 4 2 5 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 5 2 5 5 2 5 2 2 5 2 5 5 2 5 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 end fill
ara=5 nux=16 nuy=16 nuz=1 typ=square
com='cfg-e'
fill
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 5 2 4 4 2 5 2 2 5 2 4 4 2 5 2
2 5 2 4 4 2 5 2 2 5 2 4 4 2 5 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 2 2 2 2 2 4 4 2 2 2 2 2 2 2
2 2 2 2 2 2 2 4 4 2 2 2 2 2 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 5 2 4 4 2 5 2 2 5 2 4 4 2 5 2
2 5 2 4 4 2 5 2 2 5 2 4 4 2 5 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 5 5 2 2 2 2 2 2 5 5 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 end fill
ara=1 nux=16 nuy=16 nuz=1 typ=square gbl=1
com='cfg-a'
fill
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 4 4 2 2 2 2 2 2 4 4 2 2 2
2 2 2 4 4 2 2 2 2 2 2 4 4 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 4 4 2 2 2 2 2 2 2
2 2 2 2 2 2 2 4 4 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 4 4 2 2 2 2 2 2 4 4 2 2 2
2 2 2 4 4 2 2 2 2 2 2 4 4 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 end fill
end array
read plot
scr=yes
ttl='plot2'
pic=mixtures
xul=-10.4
yul=10.4
zul=0
xlr=10.4
ylr=-10.4
zlr=0
nax=1000
clr=1 255 0 0
2 0 0 205
4 0 238 0
5 205 205 0
6 238 0 0
7 145 44 238
8 150 150 150
9 240 200 220
10 0 191 255
11 224 255 255
13 255 255 224
14 255 0 0
16 67 110 238
17 174 238 238
end color
uax=1 vdn=-1
end
end plot
read bnds
body=30
all=mirror
end bnds
read start
nst=0
xsm=0
xsp=0
ysm=0
ysp=0
zsm=0
zsp=0
end start
end data
end model
end
Thank you
John Lyons