#/bin/bash
echo ''
echo '======================================================================='
echo ' Script to run NJOY to get elemental cross section from isotope ones'
echo ' and than group'
echo '======================================================================='
echo
end=jendl
Tape1=n-Fe-054.$end
Tape2=n-Fe-056.$end
Tape3=n-Fe-057.$end
Tape4=n-Fe-058.$end
end=endf
Tape1=n-026_Fe_054.$end
Tape2=n-026_Fe_056.$end
Tape3=n-026_Fe_057.$end
Tape4=n-026_Fe_058.$end
end=jeff32
Tape1=n-26-Fe-054.$end
Tape2=n-26-Fe-056.$end
Tape3=n-26-Fe-057.$end
Tape4=n-26-Fe-058.$end
end=tendl
Tape1=Fe054-n.$end
Tape2=Fe056-n.$end
Tape3=Fe057-n.$end
Tape4=Fe058-n.$end
cp $Tape1 tape20
cp $Tape2 tape30
cp $Tape3 tape40
cp $Tape4 tape50
MAT1=2625
MAT2=2631
MAT3=2634
MAT4=2637
Abund1=0.059
Abund2=0.9172
Abund3=0.021
Abund4=0.0028
MAT=2600
TEMP=300
echo ---------------------------------------------------------------------
echo to construct Element MAT = $MAT at Temperature = $TEMP
echo from Evaluations: $Tape1, $Tape2, $Tape3, $Tape4
echo and Abundancies: $Abund1, $Abund2, $Abund3, $Abund4
echo ---------------------------------------------------------------------
echo ' Creating Input for MODER ...'
cat > input << EOF
moder / Extract Files
1 21 / combining Tapes1-4 on nout #21
'$Tape1 + $Tape2 + $Tape3 + $Tape4' / nout tape description
20 $MAT1 / nin
30 $MAT2 / nin
40 $MAT3 / nin
50 $MAT4 / nin
0 / termination
stop
EOF
echo ' Running NJOY ...'
../njoy2012/xnjoy input << EOF
reconr / Reconstruct XS
21 -22 / input endf tape, output pendf tape
'pendf tape for $Tape1 MAT = $MAT1'/
$MAT1 1 / MAT, No of descriptive cards
.001 0 0.005 / Reconstruction Tolerance, Temp, Reconstruction Tolerance
'$Tape1 processed by NJOY-2012.50'/ Descriptive card
0/
broadr / Doppler Broadening
21 -22 -23 / input endf tape, input pendf tape, output pendf tape
$MAT1 1 0 0 0. / MAT, No of Temp, restart, bootstrap, starting Temp
.001 2.e+6 / Tolerance for thinning, max. energy for broadening and thinning (default 1MeV)
$TEMP / Temp
0/
purr / Process Unresolved Resonance Range if any
21 -23 -24 / input endf tape, input pendf tape, output pendf tape
$MAT1 1 1 20 100/ MAT, print option, ?, number of probability bins, number of resonance ladders
$TEMP / Temp
1.E+10 / sigma zero values
0/
heatr / Adding of Heating (Kerma) and Damage Energy
21 -24 -25 / input endf tape, input pendf tape, output pendf tape, output for graph
$MAT1 3 0 0 0 1/ MAT, No of Kerma, nqa, No of Temp, iprint
303 443 444 / MT Kerma (301-tot, 304-inel, 318-fis, 443-high_limit) or Damage Energy (444-total,445-el,447-disappear)
gaspr / Adding of gas production
21 -25 -26 / input endf tape, input pendf tape, output pendf tape
stop
EOF
echo ' Running NJOY for '$MAT1 ...
../njoy2012/xnjoy input << EOF
reconr / Reconstruct XS
21 -32 / input endf tape, output pendf tape
'pendf tape for $Tape2 MAT = $MAT2'/
$MAT2 1 / MAT, No of descriptive cards
.001 0 0.005 / Reconstruction Tolerance, Temp, Reconstruction Tolerance
'$Tape2 processed by NJOY-2012.50'/ Descriptive card
0/
broadr / Doppler Broadening
21 -32 -33 / input endf tape, input pendf tape, output pendf tape
$MAT2 1 0 0 0. / MAT, No of Temp, restart, bootstrap, starting Temp
.001 2.e+6 / Tolerance for thinning, max. energy for broadening and thinning (default 1MeV)
$TEMP / Temp
0/
purr / Process Unresolved Resonance Range if any
21 -33 -34 / input endf tape, input pendf tape, output pendf tape
$MAT2 1 1 20 100/ MAT, print option, ?, number of probability bins, number of resonance ladders
$TEMP / Temp
1.E+10 / sigma zero values
0/
heatr / Adding of Heating (Kerma) and Damage Energy
21 -34 -35 / input endf tape, input pendf tape, output pendf tape, output for graph
$MAT2 3 0 0 0 1 / MAT, No of Kerma, nqa, No of Temp, iprint
303 443 444 / MT Kerma (301-tot, 304-inel, 318-fis, 443-high_limit) or Damage Energy (444-total,445-el,447-disappear)
gaspr / Adding of gas production
21 -35 -36 / input endf tape, input pendf tape, output pendf tape
stop
EOF
echo ' Running NJOY for '$MAT2 ...
../njoy2012/xnjoy input << EOF
reconr / Reconstruct XS
21 -42 / input endf tape, output pendf tape
'pendf tape for $Tape3 MAT = $MAT3'/
$MAT3 1 / MAT, No of descriptive cards
.001 0 0.005 / Reconstruction Tolerance, Temp, Reconstruction Tolerance
'$Tape3 processed by NJOY-2012.50'/ Descriptive card
0/
broadr / Doppler Broadening
21 -42 -43 / input endf tape, input pendf tape, output pendf tape
$MAT3 1 0 0 0. / MAT, No of Temp, restart, bootstrap, starting Temp
.001 2.e+6 / Tolerance for thinning, max. energy for broadening and thinning (default 1MeV)
$TEMP / Temp
0/
purr / Process Unresolved Resonance Range if any
21 -43 -44 / input endf tape, input pendf tape, output pendf tape
$MAT3 1 1 20 100/ MAT, print option, ?, number of probability bins, number of resonance ladders
$TEMP / Temp
1.E+10 / sigma zero values
0/
heatr / Adding of Heating (Kerma) and Damage Energy
21 -44 -45 / input endf tape, input pendf tape, output pendf tape, output for graph
$MAT3 3 0 0 0 1 / MAT, No of Kerma, nqa, No of Temp, iprint
303 443 444 / MT Kerma (301-tot, 304-inel, 318-fis, 443-high_limit) or Damage Energy (444-total,445-el,447-disappear)
gaspr / Adding of gas production
21 -45 -46 / input endf tape, input pendf tape, output pendf tape
stop
EOF
echo ' Running NJOY for' $MAT3 ...
../njoy2012/xnjoy input << EOF
reconr / Reconstruct XS
21 -52 / input endf tape, output pendf tape
'pendf tape for $Tape4 MAT = $MAT4'/
$MAT4 1 / MAT, No of descriptive cards
.001 0 0.005 / Reconstruction Tolerance, Temp, Reconstruction Tolerance
'$Tape4 processed by NJOY-2012.50'/ Descriptive card
0/
broadr / Doppler Broadening
21 -52 -53 / input endf tape, input pendf tape, output pendf tape
$MAT4 1 0 0 0. / MAT, No of Temp, restart, bootstrap, starting Temp
.001 2.e+6 / Tolerance for thinning, max. energy for broadening and thinning (default 1MeV)
$TEMP / Temp
0/
purr / Process Unresolved Resonance Range if any
21 -53 -54 / input endf tape, input pendf tape, output pendf tape
$MAT4 1 1 20 100 / MAT, print option, ?, number of probability bins, number of resonance ladders
$TEMP / Temp
1.E+10 / sigma zero values
0/
heatr / Adding of Heating (Kerma) and Damage Energy
21 -54 -55 / input endf tape, input pendf tape, output pendf tape, output for graph
$MAT4 3 0 0 0 1 / MAT, No of Kerma, nqa, No of Temp, iprint
303 443 444 / MT Kerma (301-tot, 304-inel, 318-fis, 443-high_limit) or Damage Energy (444-total,445-el,447-disappear)
gaspr / Adding of gas production
21 -55 -56 / input endf tape, input pendf tape, output pendf tape
stop
EOF
echo ' Running NJOY for' $MAT4 ...
../njoy2012/xnjoy input << EOF
mixr / mixing the isotopic cross sections
90 -26 -36 -46 -56 / nout, nin1 (endf or pendf), nin2
203 204 205 206 207 301 444 / MTs for output
$MAT1 $Abund1
$MAT2 $Abund2
$MAT3 $Abund3
$MAT4 $Abund4 / MAT, abundance
$TEMP / Temp
$MAT 26000. 55.3670672 / output material: MAT, za, awr
'26-Fe- 0 from $end by NJOY2012.50/MIXR (update XX Sep 2015)' / Description
groupr / Prepare multigroup data for neutrons
90 90 0 92 /nendf, pendf, ngout1 - input gout (default=0), ngout2 - output gout (default=0)
$MAT 1 0 6 1 1 1 / matb, ign, igg, iwt, lord, ntemp, nsigz, iprint
'MAT = $MAT by NJOY2012.50 (update XX Sep 2015)' / title
$TEMP / Temp
1.E+10 / Card5: sig zero
700 / Card 6a: ngn(if ign=1) =640(<20MeV)+40(20-60)+20(60-100,step=2MeV)
1.00000E-4 1.05000E-4 1.10000E-4 1.15000E-4 1.20000E-4 1.27500E-4
1.35000E-4 1.42500E-4 1.50000E-4 1.60000E-4 1.70000E-4 1.80000E-4
1.90000E-4 2.00000E-4 2.10000E-4 2.20000E-4 2.30000E-4 2.40000E-4
2.55000E-4 2.70000E-4 2.80000E-4 3.00000E-4 3.20000E-4 3.40000E-4
3.60000E-4 3.80000E-4 4.00000E-4 4.25000E-4 4.50000E-4 4.75000E-4
5.00000E-4 5.25000E-4 5.50000E-4 5.75000E-4 6.00000E-4 6.30000E-4
6.60000E-4 6.90000E-4 7.20000E-4 7.60000E-4 8.00000E-4 8.40000E-4
8.80000E-4 9.20000E-4 9.60000E-4 .001000000 .001050000 .001100000
.001150000 .001200000 .001275000 .001350000 .001425000 .001500000
.001600000 .001700000 .001800000 .001900000 .002000000 .002100000
.002200000 .002300000 .002400000 .002550000 .002700000 .002800000
.003000000 .003200000 .003400000 .003600000 .003800000 .004000000
.004250000 .004500000 .004750000 .005000000 .005250000 .005500000
.005750000 .006000000 .006300000 .006600000 .006900000 .007200000
.007600000 .008000000 .008400000 .008800000 .009200000 .009600000
.010000000 .010500000 .011000000 .011500000 .012000000 .012750000
.013500000 .014250000 .015000000 .016000001 .017000001 .017999999
.018999999 .020000000 .021000000 .022000000 .023000000 .024000000
.025500000 .027000001 .028000001 .029999999 .032000002 .034000002
.035999998 .037999999 .039999999 .042500000 .045000002 .047499999
.050000001 .052499998 .055000000 .057500001 .059999999 .063000001
.066000000 .068999998 .071999997 .075999998 .079999998 .083999999
.088000000 .092000000 .096000001 .100000001 .104999997 .109999999
.115000002 .119999997 .127499998 .135000005 .142499998 .150000006
.159999996 .170000002 .180000007 .189999998 .200000003 .209999993
.219999999 .230000004 .239999995 .254999995 .270000011 .280000001
.300000012 .319999993 .340000004 .360000014 .379999995 .400000006
.425000012 .449999988 .474999994 .500000000 .524999976 .550000012
.574999988 .600000024 .629999995 .660000026 .689999998 .720000029
.759999990 .800000012 .839999974 .879999995 .920000017 .959999979
1.00000000 1.04999995 1.10000002 1.14999998 1.20000005 1.27499998
1.35000002 1.42499995 1.50000000 1.60000002 1.70000005 1.79999995
1.89999998 2.00000000 2.09999990 2.20000005 2.29999995 2.40000010
2.54999995 2.70000005 2.79999995 3.00000000 3.20000005 3.40000010
3.59999990 3.79999995 4.00000000 4.25000000 4.50000000 4.75000000
5.00000000 5.25000000 5.50000000 5.75000000 6.00000000 6.30000019
6.59999990 6.90000010 7.19999981 7.59999990 8.00000000 8.39999962
8.80000019 9.19999981 9.60000038 10.0000000 10.5000000 11.0000000
11.5000000 12.0000000 12.7500000 13.5000000 14.2500000 15.0000000
16.0000000 17.0000000 18.0000000 19.0000000 20.0000000 21.0000000
22.0000000 23.0000000 24.0000000 25.5000000 27.0000000 28.0000000
30.0000000 32.0000000 34.0000000 36.0000000 38.0000000 40.0000000
42.5000000 45.0000000 47.5000000 50.0000000 52.5000000 55.0000000
57.5000000 60.0000000 63.0000000 66.0000000 69.0000000 72.0000000
76.0000000 80.0000000 84.0000000 88.0000000 92.0000000 96.0000000
100.000000 105.000000 110.000000 115.000000 120.000000 127.500000
135.000000 142.500000 150.000000 160.000000 170.000000 180.000000
190.000000 200.000000 210.000000 220.000000 230.000000 240.000000
255.000000 270.000000 280.000000 300.000000 320.000000 340.000000
360.000000 380.000000 400.000000 425.000000 450.000000 475.000000
500.000000 525.000000 550.000000 575.000000 600.000000 630.000000
660.000000 690.000000 720.000000 760.000000 800.000000 840.000000
880.000000 920.000000 960.000000 1000.00000 1050.00000 1100.00000
1150.00000 1200.00000 1275.00000 1350.00000 1425.00000 1500.00000
1600.00000 1700.00000 1800.00000 1900.00000 2000.00000 2100.00000
2200.00000 2300.00000 2400.00000 2550.00000 2700.00000 2800.00000
3000.00000 3200.00000 3400.00000 3600.00000 3800.00000 4000.00000
4250.00000 4500.00000 4750.00000 5000.00000 5250.00000 5500.00000
5750.00000 6000.00000 6300.00000 6600.00000 6900.00000 7200.00000
7600.00000 8000.00000 8400.00000 8800.00000 9200.00000 9600.00000
10000.0000 10500.0000 11000.0000 11500.0000 12000.0000 12750.0000
13500.0000 14250.0000 15000.0000 16000.0000 17000.0000 18000.0000
19000.0000 20000.0000 21000.0000 22000.0000 23000.0000 24000.0000
25500.0000 27000.0000 28000.0000 30000.0000 32000.0000 34000.0000
36000.0000 38000.0000 40000.0000 42500.0000 45000.0000 47500.0000
50000.0000 52500.0000 55000.0000 57500.0000 60000.0000 63000.0000
66000.0000 69000.0000 72000.0000 76000.0000 80000.0000 84000.0000
88000.0000 92000.0000 96000.0000 100000.000 105000.000 110000.000
115000.000 120000.000 127500.000 135000.000 142500.000 150000.000
160000.000 170000.000 180000.000 190000.000 200000.000 210000.000
220000.000 230000.000 240000.000 255000.000 270000.000 280000.000
300000.000 320000.000 340000.000 360000.000 380000.000 400000.000
425000.000 450000.000 475000.000 500000.000 525000.000 550000.000
575000.000 600000.000 630000.000 660000.000 690000.000 720000.000
760000.000 800000.000 840000.000 880000.000 920000.000 960000.000
1000000.00 1100000.00 1200000.00 1300000.00 1400000.00 1500000.00
1600000.00 1700000.00 1800000.00 1900000.00 2000000.00 2100000.00
2200000.00 2300000.00 2400000.00 2500000.00 2600000.00 2700000.00
2800000.00 2900000.00 3000000.00 3100000.00 3200000.00 3300000.00
3400000.00 3500000.00 3600000.00 3700000.00 3800000.00 3900000.00
4000000.00 4100000.00 4200000.00 4300000.00 4400000.00 4500000.00
4600000.00 4700000.00 4800000.00 4900000.00 5000000.00 5100000.00
5200000.00 5300000.00 5400000.00 5500000.00 5600000.00 5700000.00
5800000.00 5900000.00 6000000.00 6100000.00 6200000.00 6300000.00
6400000.00 6500000.00 6600000.00 6700000.00 6800000.00 6900000.00
7000000.00 7100000.00 7200000.00 7300000.00 7400000.00 7500000.00
7600000.00 7700000.00 7800000.00 7900000.00 8000000.00 8100000.00
8200000.00 8300000.00 8400000.00 8500000.00 8600000.00 8700000.00
8800000.00 8900000.00 9000000.00 9100000.00 9200000.00 9300000.00
9400000.00 9500000.00 9600000.00 9700000.00 9800000.00 9900000.00
10000000.0 10100000.0 10200000.0 10300000.0 10400000.0 10500000.0
10600000.0 10700000.0 10800000.0 10900000.0 11000000.0 11100000.0
11200000.0 11300000.0 11400000.0 11500000.0 11600000.0 11700000.0
11800000.0 11900000.0 12000000.0 12100000.0 12200000.0 12300000.0
12400000.0 12500000.0 12600000.0 12700000.0 12800000.0 12900000.0
13000000.0 13100000.0 13200000.0 13300000.0 13400000.0 13500000.0
13600000.0 13700000.0 13800000.0 13900000.0 14000000.0 14100000.0
14200000.0 14300000.0 14400000.0 14500000.0 14600000.0 14700000.0
14800000.0 14900000.0 15000000.0 15100000.0 15200000.0 15300000.0
15400000.0 15500000.0 15600000.0 15700000.0 15800000.0 15900000.0
16000000.0 16100000.0 16200000.0 16300000.0 16400000.0 16500000.0
16600000.0 16700000.0 16800000.0 16900000.0 17000000.0 17100000.0
17200000.0 17300000.0 17400000.0 17500000.0 17600000.0 17700000.0
17800000.0 17900000.0 18000000.0 18100000.0 18200000.0 18300000.0
18400000.0 18500000.0 18600000.0 18700000.0 18800000.0 18900000.0
19000000.0 19100000.0 19200000.0 19300000.0 19400000.0 19500000.0
19600000.0 19700000.0 19800000.0 19900000.0 20000000.0 20500000.0
21000000.0 21500000.0 22000000.0 22500000.0 23000000.0 23500000.0
24000000.0 24500000.0 25000000.0 25500000.0 26000000.0 26500000.0
27000000.0 27500000.0 28000000.0 28500000.0 29000000.0 29500000.0
30000000.0 31000000.0 32000000.0 33000000.0 34000000.0 35000000.0
36000000.0 37000000.0 38000000.0 39000000.0 40000000.0 42000000.0
44000000.0 46000000.0 48000000.0 50000000.0 52000000.0 54000000.0
56000000.0 58000000.0 60000000.0 62000000.0 64000000.0 66000000.0
68000000.0 70000000.0 72000000.0 74000000.0 76000000.0 78000000.0
80000000.0 82000000.0 84000000.0 86000000.0 88000000.0 90000000.0
92000000.0 94000000.0 96000000.0 98000000.0 100000000.0 / ngn+1 group breaks
3 203 '(n,xp)' / Card9: mfd mtd mtname to process
3 204 '(n,xd)' /
3 205 '(n,xt)' /
3 205 '(n,xHe3)'/
3 207 '(n,xHe4)'/
3 301 'energy-balance total KERMA'/
3 444 'damage energy'/
0/ terminates this temperature or material
0/ Card10: next material to process (matd = 0 - terminate groupr run)
stop
EOF
echo ' Running NJOY ...'
../njoy2012/xnjoy