=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