=shell rm *.data end =csas26 parm=size=1000000 heu-sol-therm-037-002 238g infhommedium solnuo2(no3)2 1 41.9 0.25 den=1.063 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 31.3895 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 31.3895 15.5 cylinder 20 1.35 31.3895 15.5 hexprism 30 3.0 31.3895 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 31.3895 cylinder 20 1.35 198.35 31.3895 hexprism 30 3.0 198.35 31.3895 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 31.3895 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 31.3895 cylinder 15 80.0 31.3895 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 33r4 33r4 18r4 6 6 6 6 6 6 6 8r4 16r4 6 6 6 6 6 6 6 6 6 6 7r4 14r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6r4 13r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6r4 12r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6r4 11r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6r4 10r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6r4 9r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6r4 8r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6r4 8r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 7r4 7r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 7r4 7r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 8r4 6r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 8r4 6r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 9r4 6r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 10r4 6r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 11r4 6r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 12r4 6r4 6 6 6 6 6 6 6 6 6 6 6 6 6 6 13r4 6r4 6 6 6 6 6 6 6 6 6 6 6 6 6 14r4 7r4 6 6 6 6 6 6 6 6 6 6 16r4 8r4 6 6 6 6 6 6 6 18r4 33r4 33r4 33r4 33r4 33r4 33r4 33r2 33r2 33r2 33r2 18r2 1 1 1 1 1 1 1 1 1 6r2 16r2 1 1 1 1 1 1 1 1 1 1 1 1 5r2 14r2 1 1 1 1 3 3 3 3 3 3 3 1 1 1 1 4r2 13r2 1 1 1 3 3 3 3 3 3 3 3 3 3 1 1 1 4r2 12r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 4r2 11r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 4r2 10r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 4r2 9r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 4r2 8r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 4r2 7r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 4r2 6r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 4r2 6r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 5r2 5r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 5r2 5r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 6r2 4r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 6r2 4r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 7r2 4r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 8r2 4r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 9r2 4r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 10r2 4r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 11r2 4r2 1 1 3 3 3 3 3 3 3 3 3 3 3 3 3 1 1 12r2 4r2 1 1 1 3 3 3 3 3 3 3 3 3 3 1 1 1 13r2 4r2 1 1 1 1 3 3 3 3 3 3 3 1 1 1 1 14r2 5r2 1 1 1 1 1 1 1 1 1 1 1 1 16r2 6r2 1 1 1 1 1 1 1 1 1 18r2 33r2 33r2 33r2 33r2 33r9 33r9 33r9 33r9 33r9 33r9 18r9 7 7 7 7 7 7 7 8r9 16r9 7 7 7 7 7 7 7 7 7 7 7r9 14r9 7 7 7 7 7 7 7 7 7 7 7 7 7 6r9 13r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 6r9 12r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 6r9 11r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 6r9 10r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 6r9 9r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 6r9 8r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 6r9 8r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7r9 7r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7r9 7r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 8r9 6r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 8r9 6r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 9r9 6r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 10r9 6r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 11r9 6r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 12r9 6r9 7 7 7 7 7 7 7 7 7 7 7 7 7 7 13r9 6r9 7 7 7 7 7 7 7 7 7 7 7 7 7 14r9 7r9 7 7 7 7 7 7 7 7 7 7 16r9 8r9 7 7 7 7 7 7 7 18r9 33r9 33r9 33r9 33r9 33r9 33r9 33r5 33r5 33r5 33r5 33r5 33r5 18r5 8 8 8 8 8 8 8 8r5 16r5 8 8 8 8 8 8 8 8 8 8 7r5 14r5 8 8 8 8 8 8 8 8 8 8 8 8 8 6r5 13r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 6r5 12r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 6r5 11r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 6r5 10r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 6r5 9r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 6r5 8r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 6r5 8r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7r5 7r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7r5 7r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8r5 6r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8r5 6r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 9r5 6r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 10r5 6r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 11r5 6r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 12r5 6r5 8 8 8 8 8 8 8 8 8 8 8 8 8 8 13r5 6r5 8 8 8 8 8 8 8 8 8 8 8 8 8 14r5 7r5 8 8 8 8 8 8 8 8 8 8 16r5 8r5 8 8 8 8 8 8 8 18r5 33r5 33r5 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