=shell
rm *.data
end
=csas26       parm=size=1000000
heu-sol-therm-037-009
238g infhommedium
solnuo2(no3)2 1 83.0 0.500 den=1.129 1 293.0 92234 0.90 92235 89.08
            92236 0.21 92238 9.81 end  
h2o           2 den=0.9983 end  
fe            3 0 5.9088e-2 end  
cr            3 0 1.6532e-2 end  
ni            3 0 8.1369e-3 end  
mn            3 0 1.3039e-3 end  
si            3 0 1.3603e-3 end  
ti            3 0 5.9844e-4 end 
b4c           4 den=1.25 end  
end comp  
read param npg=1000 gen=1020 nsk=20 far=yes flx=yes fdn=yes end param  
read geometry
unit 1
com='holes in grid plate'
cylinder 10 1.625 15.5 15.0
hexprism 20 3.0 15.5 15.0
media 1 1 10
media 3 1 20 -10
boundary 20
unit 2
com='blank grid location'
hexprism 20 3.0 15.5 15.0
media 3 1 20
boundary 20
unit 3
com='absorber rods through grid plate'
cylinder 10 1.6 15.5 15.0
cylinder 15 1.65 15.5 15.0
cylinder 20 1.35 15.5 15.0
hexprism 30 3.0 15.5 15.0
media 4 1 20
media 3 1 10 -15
media 1 1 15 -20
media 3 1 30 -15
boundary 30
unit 4
com='solution below grid plate'
hexprism 10 3.0 15.0 0.0
media 1 1 10
boundary 10
unit 5
com='void above solution'
hexprism 10 3.0 198.35 35.3848
media 0 1 10
boundary 10
unit 6
com='absorber rod below grid plate'
cylinder 10 1.6 15.0 0.0
cylinder 20 1.35 15.0 0.4
hexprism 30 3.0 15.0 0.0
media 4 1 20
media 3 1 10 -20
media 1 1 30 -10
boundary 30
unit 7
com='absorber rod above grid plate, through solution'
cylinder 10 1.6 35.3848 15.5
cylinder 20 1.35 35.3848 15.5
hexprism 30 3.0 35.3848 15.5
media 4 1 20
media 3 1 10 -20
media 1 1 30 -10
boundary 30
unit 8
com='absorber rod above solution'
cylinder 10 1.6 198.35 35.3848
cylinder 20 1.35 198.35 35.3848
hexprism 30 3.0 198.35 35.3848
media 4 1 20
media 3 1 10 -20
media 0 1 30 -10
boundary 30
unit 9
com='solution above grid plate'
hexprism 10 3.0 35.3848 15.5
media 1 1 10
boundary 10
global unit 11
cylinder  5 79.7 198.35 0.0
cylinder 10 80.0 198.35 35.3848
cylinder 15 80.0 35.3848 0.0
cylinder 20 80.65 198.35 -0.65
cylinder 30 84.65 150.0 -4.65
cylinder 40 85.3 150.65 -5.3
cylinder 50 99.7 198.35 -101.15
cylinder 60 100.5 198.35 -102.15
array 1 5 place 17 17 1 0.0 0.0 0.0
media 0 1 10 -5
media 1 1 15 -5
media 3 1 20 -10 -15
media 2 1 30 -20
media 3 1 40 -30 -20
media 0 1 50 -40 -20
media 3 1 60 -50
boundary 60
end geom
read array ara=1 typ=hexagonal nux=33 nuy=33 nuz=4
fill
33r4
33r4
33r4
33r4
18r4                     6 6 6 6 6 6 6 6 6                   6r4
16r4                  6 6 6 6 6 6 6 6 6 6 6 6                5r4
14r4               6 6 6 6 6 6 6 6 6 6 6 6 6 6 6             4r4
13r4              6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6            4r4
12r4             6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6           4r4
11r4            6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6          4r4
10r4           6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6         4r4
9r4           6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6        4r4
8r4          6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6       4r4
7r4         6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6      4r4
6r4        6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6     4r4
6r4         6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6      5r4
5r4        6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6     5r4
5r4         6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6      6r4
4r4        6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6     6r4
4r4         6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6      7r4
4r4          6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6       8r4
4r4           6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6        9r4
4r4            6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6         10r4
4r4             6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6          11r4
4r4              6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6           12r4
4r4               6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6            13r4
4r4                6 6 6 6 6 6 6 6 6 6 6 6 6 6 6             14r4
5r4                   6 6 6 6 6 6 6 6 6 6 6 6                16r4
6r4                      6 6 6 6 6 6 6 6 6                   18r4
33r4
33r4
33r4
33r4

33r2
33r2
33r2
33r2
18r2                     3 3 3 3 3 3 3 3 3                   6r2
16r2                  3 3 3 3 3 3 3 3 3 3 3 3                5r2
14r2               3 3 3 3 3 3 3 3 3 3 3 3 3 3 3             4r2
13r2              3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3            4r2
12r2             3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3           4r2
11r2            3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3          4r2
10r2           3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3         4r2
9r2           3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3        4r2
8r2          3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3       4r2
7r2         3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3      4r2
6r2        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3     4r2
6r2         3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3      5r2
5r2        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3     5r2
5r2         3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3      6r2
4r2        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3     6r2
4r2         3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3      7r2
4r2          3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3       8r2
4r2           3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3        9r2
4r2            3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3         10r2
4r2             3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3          11r2
4r2              3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3           12r2
4r2               3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3            13r2
4r2                3 3 3 3 3 3 3 3 3 3 3 3 3 3 3             14r2
5r2                   3 3 3 3 3 3 3 3 3 3 3 3                16r2
6r2                      3 3 3 3 3 3 3 3 3                   18r2
33r2
33r2
33r2
33r2

33r9
33r9
33r9
33r9
18r9                     7 7 7 7 7 7 7 7 7                   6r9
16r9                  7 7 7 7 7 7 7 7 7 7 7 7                5r9
14r9               7 7 7 7 7 7 7 7 7 7 7 7 7 7 7             4r9
13r9              7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7            4r9
12r9             7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7           4r9
11r9            7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7          4r9
10r9           7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7         4r9
9r9           7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7        4r9
8r9          7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7       4r9
7r9         7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7      4r9
6r9        7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7     4r9
6r9         7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7      5r9
5r9        7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7     5r9
5r9         7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7      6r9
4r9        7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7     6r9
4r9         7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7      7r9
4r9          7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7       8r9
4r9           7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7        9r9
4r9            7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7         10r9
4r9             7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7          11r9
4r9              7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7           12r9
4r9               7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7            13r9
4r9                7 7 7 7 7 7 7 7 7 7 7 7 7 7 7             14r9
5r9                   7 7 7 7 7 7 7 7 7 7 7 7                16r9
6r9                      7 7 7 7 7 7 7 7 7                   18r9
33r9
33r9
33r9
33r9

33r5
33r5
33r5
33r5
18r5                     8 8 8 8 8 8 8 8 8                   6r5
16r5                  8 8 8 8 8 8 8 8 8 8 8 8                5r5
14r5               8 8 8 8 8 8 8 8 8 8 8 8 8 8 8             4r5
13r5              8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8            4r5
12r5             8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8           4r5
11r5            8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8          4r5
10r5           8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8         4r5
9r5           8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8        4r5
8r5          8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8       4r5
7r5         8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8      4r5
6r5        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8     4r5
6r5         8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8      5r5
5r5        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8     5r5
5r5         8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8      6r5
4r5        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8     6r5
4r5         8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8      7r5
4r5          8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8       8r5
4r5           8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8        9r5
4r5            8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8         10r5
4r5             8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8          11r5
4r5              8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8           12r5
4r5               8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8            13r5
4r5                8 8 8 8 8 8 8 8 8 8 8 8 8 8 8             14r5
5r5                   8 8 8 8 8 8 8 8 8 8 8 8                16r5
6r5                      8 8 8 8 8 8 8 8 8                   18r5
33r5
33r5
33r5
33r5

end fill
end array
read plot  scr=yes  lpi=10
ttl='y-z slice at x=0.0 '
xul=0.0  yul=-100.5  zul=190 xlr=0.0  ylr=100.5  zlr=-102.15
vax=1 wdn=-1 nax=400  end plt0
ttl='x-y slice at z=10  below grid plate'
xul=-100.5 yul=100.5  zul=10.0 xlr=100.5 ylr=-100.5  zlr=10.0
uax=1 vdn=-1 nax=400 end plt1
ttl='x-y slice at z=15.2  through grid plate'
xul=-100.5 yul=100.5  zul=15.2 xlr=100.5 ylr=-100.5  zlr=15.2
uax=1 vdn=-1 nax=400 end plt2
ttl='x-y slice at z=100  above solution'
xul=-100.5 yul=100.5  zul=100 xlr=100.5 ylr=-100.5  zlr=100
uax=1 vdn=-1 nax=400 end plt4
end plot   
end data
end