From 14c728ca2a33e98e2a2a10ecb1c1bbb4dbb86048 Mon Sep 17 00:00:00 2001 From: glavic_a Date: Fri, 31 May 2019 12:56:11 +0200 Subject: [PATCH] Fix fatal errors and many warnings --- e2-estia.i | 9 +- geometry/ess_monolith_surfaces.i | 22 +-- geometry/ess_outside_surfaces.i | 9 +- geometry/estia_cave.i | 4 +- geometry/estia_inbunker.i | 10 +- geometry/materials.i | 272 +++++++++++++++---------------- geometry/materials_monolith.i | 2 +- sources/proton.i | 4 +- 8 files changed, 168 insertions(+), 164 deletions(-) diff --git a/e2-estia.i b/e2-estia.i index c5a3172..f741986 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -165,8 +165,9 @@ c volume calculated for tally cells (spheres) c and set for chopper disk/spindle and stages x1/x2/z1/z2/rot c vol 1 506r 11j 1 183r 3848.5 10104 1j 1 7r 479.5 479.5 913 913 1734 c 1 96r $ wall segments -imp:h,/,d,t,s,a 1 290r 0 0 1 812r -imp:n,p 1 290r 0 0 1 812r +c -234 +imp:h,/,d,t,s,a 1 290r 0 0 1 578r +imp:n,p 1 290r 0 0 1 578r cut:h,/,d,t,s,a 1e+08 1e-3 phys:n 3000 0 phys:p 100 @@ -177,9 +178,9 @@ lea 1 4 1 0 1 0 0 1 c * For testing for geometry errors c read file=sources/geometry_test.i echo c * Source for simulation from proton beam (needs different WW and ESplit) -c read file=sources/proton.i echo +read file=sources/proton.i echo c * Source Term from Valentina for 2m adapted to size -read file=sources/e2_2m.i echo +c read file=sources/e2_2m.i echo c read file=sources/e2_2m_act.i echo c * McStas source terms for gamma doeses c read file=sources/mcstas_source_vs.i echo diff --git a/geometry/ess_monolith_surfaces.i b/geometry/ess_monolith_surfaces.i index 4b682d3..4be4647 100644 --- a/geometry/ess_monolith_surfaces.i +++ b/geometry/ess_monolith_surfaces.i @@ -992,7 +992,7 @@ C ------------------------------------------------------------------------------ 33172 301 P 9.973595E-01 0.000000E+00 7.262183E-02 0.000000E+00 $ !!! C ------------------------------------------------------------------------------ 34001 301 P -2.880816E-01 0.000000E+00 0.000000E+00 1.000000E+00 $ !!! -34002 301 P 2.880816E-01 0.000000E+00 0.000000E+00 1.000000E+00 $ !!! +34002 301 PX 3.4712387 $ !!! 34003 301 P -2.837050E-01 0.000000E+00 5.002484E-02 1.000000E+00 $ !!! 34004 301 P 2.837050E-01 0.000000E+00 -5.002484E-02 1.000000E+00 $ !!! 34005 301 P -2.707081E-01 0.000000E+00 9.852970E-02 1.000000E+00 $ !!! @@ -1009,7 +1009,7 @@ C ------------------------------------------------------------------------------ 34016 301 P 9.852970E-02 0.000000E+00 -2.707081E-01 1.000000E+00 $ !!! 34017 301 P -5.002484E-02 0.000000E+00 2.837050E-01 1.000000E+00 $ !!! 34018 301 P 5.002484E-02 0.000000E+00 -2.837050E-01 1.000000E+00 $ !!! -34019 301 P 0.000000E+00 0.000000E+00 2.880816E-01 1.000000E+00 $ !!! +34019 301 PZ 3.4712387 $ !!! 34020 301 P 0.000000E+00 0.000000E+00 -2.880816E-01 1.000000E+00 $ !!! 34021 301 P 5.002484E-02 0.000000E+00 2.837050E-01 1.000000E+00 $ !!! 34022 301 P -5.002484E-02 0.000000E+00 -2.837050E-01 1.000000E+00 $ !!! @@ -1029,7 +1029,7 @@ C ------------------------------------------------------------------------------ 34036 301 P -2.837050E-01 0.000000E+00 -5.002484E-02 1.000000E+00 $ !!! C ------------------------------------------------------------------------------ 34101 301 P -3.743581E-01 0.000000E+00 0.000000E+00 1.000000E+00 $ !!! -34102 301 P 3.743581E-01 0.000000E+00 0.000000E+00 1.000000E+00 $ !!! +34102 301 PX 2.6712391 $ !!! 34103 301 P -3.686708E-01 0.000000E+00 6.500661E-02 1.000000E+00 $ !!! 34104 301 P 3.686708E-01 0.000000E+00 -6.500661E-02 1.000000E+00 $ !!! 34105 301 P -3.517816E-01 0.000000E+00 1.280380E-01 1.000000E+00 $ !!! @@ -1046,7 +1046,7 @@ C ------------------------------------------------------------------------------ 34116 301 P 1.280380E-01 0.000000E+00 -3.517816E-01 1.000000E+00 $ !!! 34117 301 P -6.500661E-02 0.000000E+00 3.686708E-01 1.000000E+00 $ !!! 34118 301 P 6.500661E-02 0.000000E+00 -3.686708E-01 1.000000E+00 $ !!! -34119 301 P 0.000000E+00 0.000000E+00 3.743581E-01 1.000000E+00 $ !!! +34119 301 PZ 2.6712391 $ !!! 34120 301 P 0.000000E+00 0.000000E+00 -3.743581E-01 1.000000E+00 $ !!! 34121 301 P 6.500661E-02 0.000000E+00 3.686708E-01 1.000000E+00 $ !!! 34122 301 P -6.500661E-02 0.000000E+00 -3.686708E-01 1.000000E+00 $ !!! @@ -1066,7 +1066,7 @@ C ------------------------------------------------------------------------------ 34136 301 P -3.686708E-01 0.000000E+00 -6.500661E-02 1.000000E+00 $ !!! C ------------------------------------------------------------------------------ 34201 301 P -4.103197E-01 0.000000E+00 0.000000E+00 1.000000E+00 $ !!! -34202 301 P 4.103197E-01 0.000000E+00 0.000000E+00 1.000000E+00 $ !!! +34202 301 PX 2.437124 $ !!! 34203 301 P -4.040860E-01 0.000000E+00 7.125127E-02 1.000000E+00 $ !!! 34204 301 P 4.040860E-01 0.000000E+00 -7.125127E-02 1.000000E+00 $ !!! 34205 301 P -3.855744E-01 0.000000E+00 1.403376E-01 1.000000E+00 $ !!! @@ -1083,7 +1083,7 @@ C ------------------------------------------------------------------------------ 34216 301 P 1.403376E-01 0.000000E+00 -3.855744E-01 1.000000E+00 $ !!! 34217 301 P -7.125127E-02 0.000000E+00 4.040860E-01 1.000000E+00 $ !!! 34218 301 P 7.125127E-02 0.000000E+00 -4.040860E-01 1.000000E+00 $ !!! -34219 301 P 0.000000E+00 0.000000E+00 4.103197E-01 1.000000E+00 $ !!! +34219 301 PZ 2.437124 $ !!! 34220 301 P 0.000000E+00 0.000000E+00 -4.103197E-01 1.000000E+00 $ !!! 34221 301 P 7.125127E-02 0.000000E+00 4.040860E-01 1.000000E+00 $ !!! 34222 301 P -7.125127E-02 0.000000E+00 -4.040860E-01 1.000000E+00 $ !!! @@ -1103,7 +1103,7 @@ C ------------------------------------------------------------------------------ 34236 301 P -4.040860E-01 0.000000E+00 -7.125127E-02 1.000000E+00 $ !!! C ------------------------------------------------------------------------------ 34301 301 P -1.000000E+00 0.000000E+00 0.000000E+00 1.000000E+00 $ B04 -34302 301 P 1.000000E+00 0.000000E+00 0.000000E+00 1.000000E+00 $ B04 +34302 301 PX 1 $ B04 34303 301 P -9.396926E-01 0.000000E+00 3.420201E-01 1.000000E+00 $ B04 34304 301 P 9.396926E-01 0.000000E+00 -3.420201E-01 1.000000E+00 $ B04 34305 301 P -7.660444E-01 0.000000E+00 6.427876E-01 1.000000E+00 $ B04 @@ -1259,10 +1259,10 @@ C ------------------------------------------------------------------------------ C ------------------------------------------------------------------------------ 36201 620 P 0.000000E+00 -2.044649E+01 -1.000000E+00 -2.486579E+02 $ RHG 36202 620 P -1.000000E+00 0.000000E+00 0.000000E+00 6.900000E+00 $ RHG -36203 620 P 1.000000E+00 0.000000E+00 0.000000E+00 6.900000E+00 $ RHG +36203 620 PX 6.900000E+00 $ RHG 36204 626 P 0.000000E+00 -1.683191E+01 -1.000000E+00 -2.071391E+02 $ RHG 36205 626 P -1.000000E+00 0.000000E+00 0.000000E+00 6.900000E+00 $ RHG -36206 626 P 1.000000E+00 0.000000E+00 0.000000E+00 6.900000E+00 $ RHG +36206 626 PX 6.900000E+00 $ RHG 36207 632 P 0.000000E+00 -1.430067E+01 -1.000000E+00 -1.834000E+02 $ RHG 36208 632 P 0.000000E+00 1.430067E+01 -1.000000E+00 2.067821E+02 $ RHG 36209 632 P 0.000000E+00 1.430067E+01 -1.000000E+00 -1.834000E+02 $ RHG @@ -1583,8 +1583,8 @@ C ------------------------------------------------------------------------------ 61005 CZ 1.550000E+01 $ B01 61006 CZ 2.100000E+01 $ B01 C ------------------------------------------------------------------------------ -65201 520 C/Z 0.000000E+00 0.000000E+00 9.550000E+00 $ RHG -65202 520 C/Z 0.000000E+00 0.000000E+00 8.550000E+00 $ RHG +65201 520 CZ 9.550000E+00 $ RHG +65202 520 CZ 8.550000E+00 $ RHG 65203 520 C/Z 0.000000E+00 -2.000000E+01 9.550000E+00 $ RHG 65204 520 C/Z 0.000000E+00 -2.000000E+01 8.550000E+00 $ RHG 65205 520 C/Z 0.000000E+00 -1.000000E+01 2.640000E+01 $ RHG diff --git a/geometry/ess_outside_surfaces.i b/geometry/ess_outside_surfaces.i index bc3eca7..8e2dfaa 100644 --- a/geometry/ess_outside_surfaces.i +++ b/geometry/ess_outside_surfaces.i @@ -159,9 +159,12 @@ c 104 p -0.906307787 -0.422618262 0.0 80 100154 pz 7 100155 pz -5.8 100156 pz 5.8 -100157 gq 6.01392e-05 6.01392e-05 0.00138 0 0 0 0 -0.0134711808 0 & - -0.2456138752 -100158 gq 1 1 23.01803079 0 0 0 0 -224 0 -4135.732457 +c 100157 gq 6.01392e-05 6.01392e-05 0.00138 0 0 0 0 -0.0134711808 0 & +c -0.2456138752 +100157 sq 6.01392e-05 6.01392e-05 0.00138 0 -0.0134711808 0 -0.2456138752 + 0 0 0 +c 100158 gq 1 1 23.01803079 0 0 0 0 -224 0 -4135.732457 +100158 sq 1 1 23.01803079 0 -224 0 -4135.732457 0 0 0 100159 c/z 0 112 131.15 100160 c/z 0 112 64.07 100161 pz 435 diff --git a/geometry/estia_cave.i b/geometry/estia_cave.i index 98aa5b9..92d6339 100644 --- a/geometry/estia_cave.i +++ b/geometry/estia_cave.i @@ -16,8 +16,8 @@ c floor c walls 208712 49 0.074981 -208803 -208801 208709 -208705 -208312 -208302 208713 49 0.074981 -208803 -208801 -208302 208806 208705 208805 -208714 200 -3.9 -208803 -208801 -208805 -208810 208806 208706 $ heavy - concrete +c heavy concrete +208714 200 -3.9 -208803 -208801 -208805 -208810 208806 208706 208715 49 0.074981 -208803 -208801 208706 208707 208810 208710 208716 49 0.074981 -208803 -208801 -208810 208707 208807 -208808 208717 49 0.074981 -208803 -208801 208707 208808 208708 208709 diff --git a/geometry/estia_inbunker.i b/geometry/estia_inbunker.i index c48341a..f6ffd6f 100644 --- a/geometry/estia_inbunker.i +++ b/geometry/estia_inbunker.i @@ -1,9 +1,9 @@ c Estia geometry inside monolith and bunker c elliptic guide -c NBOA universe to fill the insert cells -201000 2636 -7.78E+00 -100001 #((100061 -201101 201307 201308 201306 +c NBOA universe to fill the insert cells|universe boundary instead 100061/4 +201000 2636 -7.78E+00 -100001 #(( -201101 201307 201308 201306 201213): (201101 -201115 -201507 -201508 201213): (201115 -201117 - -201509 -201510 201214): (201117 -100064 -201507 -201508 201213)) + -201509 -201510 201214): (201117 -201507 -201508 201213)) U=10 c NBOA beam 201001 9 2.45e-05 201100 -201101 -201301 -201302 -201303 -201304 U=10 @@ -51,11 +51,11 @@ c NBOA wedges 201372)) -201502 -201311 U=10 201036 14 -8.92 201107 -201108 201381 201382 -201502 -201311 U=10 c Al window and He in front of NBOA -201041 5 0.0582256 100061 -201110 201305 201306 201307 201308 U=10 +201041 5 0.0582256 -201110 201305 201306 201307 201308 U=10 201042 9 2.45e-05 201110 -201100 201305 201306 201307 201308 U=10 c Al window and He behind NBOA 201051 9 2.45e-05 201108 -201118 201213 -201508 -201507 U=10 -201052 5 0.0582256 201118 -100064 201213 -201508 -201507 U=10 +201052 5 0.0582256 201118 201213 -201508 -201507 U=10 c extended cut-out in reflector c 1135 0 -61 3 58 1205 1206 1207 1208 c 1136 5 0.0582256 -61 3 #(41 -42) 58 1201 1202 1203 1204 #1135 diff --git a/geometry/materials.i b/geometry/materials.i index b286526..4b2bad6 100644 --- a/geometry/materials.i +++ b/geometry/materials.i @@ -3,65 +3,65 @@ c --------------- MATERIAL CARDS ------------------------ c ------------------------------------------------------- c c Gd3Ga5O12 -7.08 (Gadolinium Gallium Garnet - GGG) -m2 8016.80c 0.599760 8017.80c 0.000240 $ no data for 8018 - 31000.42c 0.250000 64152.80c 0.000300 64154.80c 0.003270 - 64155.80c 0.022200 64156.80c 0.030705 64157.80c 0.023475 - 64158.80c 0.037260 64160.80c 0.032790 +m2 8016 0.599760 8017 0.000240 $ no data for 8018 + 31000.42c 0.250000 64152 0.000300 64154 0.003270 + 64155 0.022200 64156 0.030705 64157 0.023475 + 64158 0.037260 64160 0.032790 mx2:h j model 8r c c Material : Stainless304 rho=0.0876276 -m3 6000.70c 0.00031864 14028.70c 0.00170336 15031.70c 6.95038e-05 - 16032.70c 4.4752e-05 24000.50c 0.0174813 25055.70c 0.00174159 - 26054.70c 0.00338 26056.70c 0.053455 26057.70c 0.00128218 +m3 6000 0.00031864 14028 0.00170336 15031 6.95038e-05 + 16032 4.4752e-05 24000.50c 0.0174813 25055 0.00174159 + 26054 0.00338 26056 0.053455 26057 0.00128218 28000.50c 0.00815128 pnlib=70u mx3:h model j j model 2r j j j model c C NIST STANDARD DRY AIR rho=1.2E-3 -m4 6000.70C 0.000150190 7014.70C 0.781575139 7015.70C 0.002855327 - 8016.80C 0.210748199 18036.70C 0.000015718 18038.70C 0.000002952 - 18040.70C 0.004652475 +m4 6000 0.000150190 7014 0.781575139 7015 0.002855327 + 8016 0.210748199 18036 0.000015718 18038 0.000002952 + 18040 0.004652475 mx4:h model j model j model 2r c Material : Aluminium rho=0.0582256 -m5 13027.70c 0.054381 12024.70c 0.00207 12025.70c 0.0002621 - 12026.70c 0.0002885 14028.70c 0.000217 22048.70c 8.7e-05 - 24052.24c 7.8e-05 24053.24c 1e-05 25055.70c 0.000408 - 26056.70c 0.000214 26057.70c 5e-06 29063.70c 4e-05 - 29065.70c 2e-05 30000.70c 0.000145 pnlib=70u -mt5 al27.10t +m5 13027 0.054381 12024 0.00207 12025 0.0002621 + 12026 0.0002885 14028 0.000217 22048 8.7e-05 + 24052.24c 7.8e-05 24053.24c 1e-05 25055 0.000408 + 26056 0.000214 26057 5e-06 29063 4e-05 + 29065 2e-05 30000 0.000145 pnlib=70u +mt5 al27.12t mx5:h j model 2r j model j j model 4j model c c polyethylene PSI -0.94 g/cm3 -c m6 1001.70c -0.1437 6000.70c -0.8563 -c mt6 poly.01t +c m6 1001 -0.1437 6000 -0.8563 +c mt6 poly.10t c mx6:h j model c borated polyethylene PSI (5% B4C) -0.95 g/cm3 -c m7 1001.70c -0.1183 5010.80c -0.0268 5011.80c -0.1077 6000.70c -0.7472 +c m7 1001 -0.1183 5010 -0.0268 5011 -0.1077 6000 -0.7472 c mx7:h j j j model -c mt7 poly.01t +c mt7 poly.10t c Material : helium rho=2.45e-05 c PSI Borax for HE-shutter c EP + B4C 50%-50% (1.8 g/cc) -c m8 5010.80c -0.0782 5011.80c -0.3133 6000.70c -0.4481 7015.70c -0.132 & -c 1001.70c -0.0284 +c m8 5010 -0.0782 5011 -0.3133 6000 -0.4481 7015 -0.132 & +c 1001 -0.0284 c mx8:h j j model j j c mt8 lwtr.01t c Helium -m9 2004.70c 2.45e-05 pnlib=70u +m9 2004 2.45e-05 pnlib=70u c c BorAL 1100 alloy (2.7 g/cc) -m10 5010.80c -0.04 5011.80c -0.01 13027.70c -0.9288 - 14028.70c -0.01 25055.70c -5.0E-3 26000.50c -5.0E-3 - 29000.50c -1.0E-3 30000.70c -2.0E-4 +m10 5010 -0.04 5011 -0.01 13027 -0.9288 + 14028 -0.01 25055 -5.0E-3 26000.50c -5.0E-3 + 29000.50c -1.0E-3 30000 -2.0E-4 mx10:h j model j j model 3r c c Material : H2O rho=0.100283 -m11 1001.70c 0.0668456 1002.70c 1.003e-05 8016.80c 0.0334278 - pnlib=70u -mt11 lwtr.01t +c m11 1001 0.0668456 1002 1.003e-05 8016 0.0334278 +c pnlib=70u +c mt11 lwtr.01t c c bor-float glass 13% B2O3 2.51 g/cm3 -m12 8016.80c -0.4078 11023.70c -0.096441 14028.70c -0.336563 - 20000.60c -0.107205 5010.80c -0.0104 5011.80c -0.0416 +m12 8016 -0.4078 11023 -0.096441 14028 -0.336563 + 20000.60c -0.107205 5010 -0.0104 5011 -0.0416 mx12:h j model j model j model c lead (11.34 g/cm3) m13 82000.42c 1.0 @@ -70,137 +70,137 @@ c Copper density=8.92 g/cm3 m14 29000.50c 1.0 mx14:h model c PSI Beton Ba/Fe/Bor 4.9 (4.9% B4C) (4.99 g/cm3) -c m15 1001.70c -3.62e-3 -c 6000.70c -5.23e-3 -c 8016.80c -0.08511 +c m15 1001 -3.62e-3 +c 6000 -5.23e-3 +c 8016 -0.08511 c 12000.60c -7e-4 -c 13027.70c -4.71e-3 -c 14028.70c -1.983e-2 -c 16000.70c -1.368e-2 +c 13027 -4.71e-3 +c 14028 -1.983e-2 +c 16000 -1.368e-2 c 20000.60c -2.425e-2 c 24000.50c -0.13896 -c 25055.70c -1.544e-2 +c 25055 -1.544e-2 c 26000.50c -0.52608 c 27059.50c -0.007 c 28000.50c -7.72e-2 c 56138.50c -3.985e-2 -c 5010.80c -7.66e-3 -c 5011.80c -3.068e-2 +c 5010 -7.66e-3 +c 5011 -3.068e-2 c mx15:h j model j j j j j j j j j j j j j j -c mt15 poly.01t +c mt15 poly.10t c Boron Concrete (5% B4C)(2.4 g/cc) -m16 1001.70c -6.09e-3 6000.70c -0.1726 8016.80c -0.404 - 12000.60c -0.0322 13027.70c -0.0115 14028.70c -3.3906e-2 - 20000.60c -0.3222 26054.70c -4.22e-4 26056.70c -6.92e-3 - 26057.70c -1.62e-4 5010.80c -2e-3 5011.80c -8e-3 +m16 1001 -6.09e-3 6000 -0.1726 8016 -0.404 + 12000.60c -0.0322 13027 -0.0115 14028 -3.3906e-2 + 20000.60c -0.3222 26054 -4.22e-4 26056 -6.92e-3 + 26057 -1.62e-4 5010 -2e-3 5011 -8e-3 mx16:h j model j model j j model 3r j model -mt16 lwtr.01t +mt16 lwtr.10t c Mirrobor (1.36 g/cm3) -m17 1001.70c -0.01 6000.70c -0.3 8016.80c -0.05 - 5010.80c -0.127 5011.80c -0.513 +m17 1001 -0.01 6000 -0.3 8016 -0.05 + 5010 -0.127 5011 -0.513 mx17:h j model j j model -mt17 poly.01t +mt17 poly.10t c Mineral cast 2% Boron (2.4 g/cm3) -m18 1001.70c -0.05 6000.70c -0.26 8016.80c -0.2325 - 12000.60c -0.007 13027.70c -0.023 14028.70c -0.38 - 26000.50c -0.02 28000.50c -0.004 5010.80c -0.0047 - 5011.80c -0.0188 +m18 1001 -0.05 6000 -0.26 8016 -0.2325 + 12000.60c -0.007 13027 -0.023 14028 -0.38 + 26000.50c -0.02 28000.50c -0.004 5010 -0.0047 + 5011 -0.0188 mx18:h j model j model j j model 1r j model -mt18 poly.01t +mt18 poly.10t c c Cadmium density=8.65 g/cm3 m19 48000.42c 1.0 mx19:h model c B4C sintered (2.3 g/cc) -m20 5010.80c -0.1841 5011.80c -0.74345 8016.80c -0.07245 +m20 5010 -0.1841 5011 -0.74345 8016 -0.07245 mx20:h j model j c c Material : Concrete rho=0.074981 -m49 1001.70c 0.00776555 1002.70c 1.16501e-06 8016.80c 0.0438499 - 11023.70c 0.00104778 12000.60c 0.000148662 13027.70c 0.00238815 - 14028.70c 0.0158026 16032.70c 5.63433e-05 19000.60c 0.000693104 - 20000.60c 0.00291501 26054.70c 1.84504e-05 26056.70c 0.000286826 - 26057.70c 6.5671e-06 26058.70c 8.75613e-07 +m49 1001 0.00776555 1002 1.16501e-06 8016 0.0438499 + 11023 0.00104778 12000.60c 0.000148662 13027 0.00238815 + 14028 0.0158026 16032 5.63433e-05 19000.60c 0.000693104 + 20000.60c 0.00291501 26054 1.84504e-05 26056 0.000286826 + 26057 6.5671e-06 26058 8.75613e-07 pnlib=70u -mt49 lwtr.01t +mt49 lwtr.10t mx49:h 3j model 1r j j model 2r 3j model c c Material : CastIron rho=0.0833854 -m54 6000.70c 0.00636 14028.70c 0.0051404 25055.70c 0.000386 - 26054.70c 0.00419 26056.70c 0.065789 26057.70c 0.00152 +m54 6000 0.00636 14028 0.0051404 25055 0.000386 + 26054 0.00419 26056 0.065789 26057 0.00152 pnlib=70u mx54:h model j model j j j c c Material : Steel71 rho=0.0843223 -m71 6000.70c 0.000101187 14028.70c 0.000295128 24000.50c 0.000168645 - 25055.70c 0.000505934 26054.70c 0.0048414 26056.70c 0.0759996 - 26057.70c 0.00175516 26058.70c 0.00023358 28058.70c 0.000172218 - 28060.70c 6.6335e-05 28061.70c 2.88126e-06 28062.70c 9.19272e-06 - 28064.70c 2.33992e-06 29063.70c - 0.000116651 29065.70c 5.19931e-05 pnlib=70u +m71 6000 0.000101187 14028 0.000295128 24000.50c 0.000168645 + 25055 0.000505934 26054 0.0048414 26056 0.0759996 + 26057 0.00175516 26058 0.00023358 28058 0.000172218 + 28060 6.6335e-05 28061 2.88126e-06 28062 9.19272e-06 + 28064 2.33992e-06 29063 + 0.000116651 29065 5.19931e-05 pnlib=70u mx71:h model j model 1r j j j model j j j j j j j c c Material : Copper rho=0.0844385 c (rho=8.91001 g/cc) -m73 29063.70c 0.0583892 29065.70c 0.0260493 pnlib=70u +m73 29063 0.0583892 29065 0.0260493 pnlib=70u c c Material : Be5H2O rho=0.122325 -m118 1001.70c 0.00334228 4009.70c 0.117311 8016.80c 0.00167114 - pnlib=70u -mx118:h j model j -mt118 be.60t lwtr.01t -c -c Material : Be10H2O rho=0.121164 -m119 1001.70c 0.00668456 4009.70c 0.111137 8016.80c 0.00334278 - pnlib=70u -c -mx119:h j model j -mt119 be.60t lwtr.01t +c m118 1001 0.00334228 4009 0.117311 8016 0.00167114 +c pnlib=70u +c mx118:h j model j +c mt118 be.60t lwtr.01t +c c +c c Material : Be10H2O rho=0.121164 +c m119 1001 0.00668456 4009 0.111137 8016 0.00334278 +c pnlib=70u +c c +c mx119:h j model j +c mt119 be.60t lwtr.01t c c Material : Tungsten_15.1g rho=0.0494621 -m120 74182.70c 0.0131232 74183.70c 0.00708653 74184.70c 0.0151734 - 74186.70c 0.014079 pnlib=70u +c m120 74182 0.0131232 74183 0.00708653 74184 0.0151734 +c 74186 0.014079 pnlib=70u c c Material : Iron_10H2O rho=0.0862739 -m121 1001.70c 0.00668924 8016.80c 0.00334462 26054.70c 0.00445623 - 26056.70c 0.0699533 26057.70c 0.00161553 26058.70c 0.000214997 - pnlib=70u -mt121 lwtr.01t -mx121:h 5j model +c m121 1001 0.00668924 8016 0.00334462 26054 0.00445623 +c 26056 0.0699533 26057 0.00161553 26058 0.000214997 +c pnlib=70u +c mt121 lwtr.01t +c mx121:h 5j model c c Material : CarstonConc rho=0.0824808 -m122 1001.70c 0.0229186 1002.70c 1.23628e-06 5010.80c 0.0001187 - 5011.80c 0.00047781 6000.70c 0.00623398 8016.80c 0.0406466 - 14028.70c 0.0049816 14029.70c 0.0002529 14030.70c 0.00016675 - 20040.70c 0.0062902 20042.70c 4.1981e-05 20043.70c 8.759e-06 - 20044.70c 0.0001353 20046.70c 2.5956e-07 20048.70c 1.21339e-05 - 26054.70c 1.1341e-05 26056.70c 0.00017803 26057.70c 4.1115e-06 - 26058.70c 5.47166e-07 pnlib=70u +m122 1001 0.0229186 1002 1.23628e-06 5010 0.0001187 + 5011 0.00047781 6000 0.00623398 8016 0.0406466 + 14028 0.0049816 14029 0.0002529 14030 0.00016675 + 20040 0.0062902 20042 4.1981e-05 20043 8.759e-06 + 20044 0.0001353 20046 2.5956e-07 20048 1.21339e-05 + 26054 1.1341e-05 26056 0.00017803 26057 4.1115e-06 + 26058 5.47166e-07 pnlib=70u mx122:h j j j model model 5j model 4r 3j model -mt122 poly.01t +mt122 poly.10t c c Material : ParaOrtho%0.5 rho=0.041975 -m123 1001.70c 0.000209875 1005.60c 0.0417651 pnlib=70u -mx123:h j 1001.70h -mt123 hpara.10t hortho.10t +c m123 1001 0.000209875 1005.60c 0.0417651 pnlib=70u +c mx123:h j 1001.70h +c mt123 hpara.10t hortho.10t c c Material : ESS Magnadense Heavy Concrete (density -3.9g/cm³) -m200 1001.70c 1.3431e-01 1002.70c 1.5448e-05 8016.80c 5.2962e-01 - 8017.80c 2.0174e-04 11023.70c 5.1116e-03 13027.70c 6.0596e-03 - 14028.70c 3.9343e-02 14029.70c 1.9987e-03 14030.70c 1.3190e-03 - 15031.70c 3.6290e-03 16032.70c 2.2703e-04 16033.70c 1.7929e-06 - 16034.70c 1.0158e-05 16036.70c 2.4026e-08 19039.70c 9.1400e-04 - 19040.70c 1.1512e-07 19041.70c 6.5961e-05 20040.70c 1.2235e-02 - 20042.70c 8.1657e-05 20043.70c 1.7038e-05 20044.70c 2.6327e-04 - 20046.70c 5.0454e-07 20048.70c 2.3602e-05 26054.70c 1.5462e-02 - 26056.70c 2.4273e-01 26057.70c 5.6057e-03 26058.70c 7.4601e-04 -mt200 lwtr.01t +m200 1001 1.3431e-01 1002 1.5448e-05 8016 5.2962e-01 + 8017 2.0174e-04 11023 5.1116e-03 13027 6.0596e-03 + 14028 3.9343e-02 14029 1.9987e-03 14030 1.3190e-03 + 15031 3.6290e-03 16032 2.2703e-04 16033 1.7929e-06 + 16034 1.0158e-05 16036 2.4026e-08 19039 9.1400e-04 + 19040 1.1512e-07 19041 6.5961e-05 20040 1.2235e-02 + 20042 8.1657e-05 20043 1.7038e-05 20044 2.6327e-04 + 20046 5.0454e-07 20048 2.3602e-05 26054 1.5462e-02 + 26056 2.4273e-01 26057 5.6057e-03 26058 7.4601e-04 +mt200 lwtr.10t mx200:h 3j model model 5j model 6r j model 4r 3j model c c Ar 2.7bar, CO2 0.1bar, 3He 0.2bar -m300 2003.70c 0.02381 18036.70C 0.00004046 - 18038.70C 0.0000075208 18040.70C 0.0118524 - 6000.70c 0.32143 8016.80c 0.64286 +m300 2003 0.02381 18036 0.00004046 + 18038 0.0000075208 18040 0.0118524 + 6000 0.32143 8016 0.64286 mx300:h j model 3r j c suppress warnings for particles with no available tables mx2:d model 9r @@ -227,10 +227,10 @@ mx10:d model 7r mx10:t model 7r mx10:s model 7r mx10:a model 7r -mx11:d model 2r -mx11:t model 2r -mx11:s model 2r -mx11:a model 2r +c mx11:d model 2r +c mx11:t model 2r +c mx11:s model 2r +c mx11:a model 2r mx12:d model 5r mx12:t model 5r mx12:s model 5r @@ -279,30 +279,30 @@ mx73:d model 1r mx73:t model 1r mx73:s model 1r mx73:a model 1r -mx118:d model 2r -mx118:t model 2r -mx118:s model 2r -mx118:a model 2r -mx119:d model 2r -mx119:t model 2r -mx119:s model 2r -mx119:a model 2r -mx120:d model 3r -mx120:t model 3r -mx120:s model 3r -mx120:a model 3r -mx121:d model 5r -mx121:t model 5r -mx121:s model 5r -mx121:a model 5r +c mx118:d model 2r +c mx118:t model 2r +c mx118:s model 2r +c mx118:a model 2r +c mx119:d model 2r +c mx119:t model 2r +c mx119:s model 2r +c mx119:a model 2r +c mx120:d model 3r +c mx120:t model 3r +c mx120:s model 3r +c mx120:a model 3r +c mx121:d model 5r +c mx121:t model 5r +c mx121:s model 5r +c mx121:a model 5r mx122:d model 18r mx122:t model 18r mx122:s model 18r mx122:a model 18r -mx123:d model 1r -mx123:t model 1r -mx123:s model 1r -mx123:a model 1r +c mx123:d model 1r +c mx123:t model 1r +c mx123:s model 1r +c mx123:a model 1r mx200:d model 26r mx200:t model 26r mx200:s model 26r diff --git a/geometry/materials_monolith.i b/geometry/materials_monolith.i index ae0feee..272a7f4 100644 --- a/geometry/materials_monolith.i +++ b/geometry/materials_monolith.i @@ -1401,7 +1401,7 @@ MT6140 LWTR.10T C ------------------------------------------------------------------------------ MT6141 LWTR.10T FE56.12T - C============================================================================== +C ============================================================================== C ======================================== REPLACEMENT NUCLIDES AND/OR MODELS == C ============================================================================== MX0400:H MODEL diff --git a/sources/proton.i b/sources/proton.i index 509caa4..650f280 100644 --- a/sources/proton.i +++ b/sources/proton.i @@ -3,9 +3,9 @@ cut:n j 0.99e-06 cut:p 1e+08 0.5 c weight windows from external file generated by WWG c make sure to supply the correct WWINP file -WWP:n,p j j 10 j -1 j 8.7374e15 $ WW was generated for source scaling 1 +c WWP:n,p j j 10 j -1 j 8.7374e15 $ WW was generated for source scaling 1 c energy splitting for higher neutron energies and roulette below 1MeV -esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000 +c esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000 c the source definition with protons/s sdef dir=1 erg=2000 par=9 vec=0.0 1 0.0 x=d1 y=-30 z=d2 wgt=1.560377E16 si1 A -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0