#/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=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 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 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 ...' 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 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 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 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 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 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