ESS-Estia/Estia-MCNP
1
0
Fork 0
Code Actions Packages Releases Activity

Merge branch 'full_open' geometry and model changes

Browse Source
This commit is contained in:
glavic_a
2019-04-16 09:42:28 +02:00
parent e175b19c59 e0ec1162b1
commit 6fbf034af8
18 changed files with 718 additions and 155 deletions
Download Patch File Download Diff File Expand all files Collapse all files
e2-estia.i
+13 -7
View File
@@ -17,6 +17,7 @@ c tallies
read file=geometry/estia_inbunker.i echo
read file=geometry/estia_guides.i echo
read file=geometry/estia_cave.i echo
read file=geometry/estia_universes.i echo
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
c -------------------------------------------------------
@@ -31,10 +32,10 @@ c 700 s 117.0 -155.00 13.5 3.0
703 1 s 1082 -25.0 0 5.0
704 1 s 1082 -0.25 0 5.0
705 1 s 1165 0.7 0 5.0
706 1 s 1328 5 0 12.0
706 1 s 1320 5 0 12.0
707 1 s 2190 3.5 0 12.0
708 s 730 -1380 0 30.0
709 s 1020 -1170 0 30.0
708 s 746 -1370 0 30.0
709 s 1030 -1155 0 30.0
710 s 1328 -1756 -13 12.0
c Estia geometry
read file=geometry/estia_surfaces.i echo
@@ -95,6 +96,10 @@ c ------------------- TALLY CARDS ---------------------------
c
c read file=tallies/bunker_mesh.i echo
c read file=tallies/bunker_mesh_p.i echo
c read file=tallies/bunker_mesh_np.i echo
c read file=tallies/full_mesh.i echo
c read file=tallies/full_mesh_p.i echo
c read file=tallies/full_mesh_np.i echo
c read file=tallies/inbunker_f4.i echo
c read file=tallies/chopper_pit_doses.i echo
c read file=tallies/bunker_wall_f4.i echo
@@ -122,7 +127,7 @@ ext:n 0 0 0.5V2 0 0.5V2 0 8r 0.5V2 7r 0 0.5V2 0 0.5V2 8r 0 0.5V2 0 0.5V2 0
0 0.5V2 0 0.5V2 0 0 0.5V2 0 0.5V2 0 0.5V2 0 0.5V2 0 22r
0.5V2 10r 0 4r 0.5V2 3r 0 0.5V2 13r 0 48r
-0.5V1 2r 0 4r -0.1V1 9r 0 0.1V3 3r 0 11r 0.1V3 2r
0 0.1V3 0 0.1V3 0 0.1V3 5r 0 30r
0 0.1V3 0 0 0.1V3 0 0 0.1V3 5r 0 42r
c exponential transform reference vectors: origin, selene 1, inst. shutter
vect v1 0 0 0 v2 1000 -1394 -6.5 v3 1282 -1793 -13
c
@@ -137,9 +142,9 @@ mode n p n h / d t s a
c volume calculated for tally cells (spheres)
c and set for chopper disk/spindle and stages x1/x2/z1/z2/rot
vol 1 507r 10j 1 107r 3848.5 10104 1 13r 479.5 479.5 913 913 1734
1 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 2.05e5 1 12r $ wall segments
imp:h,/,d,t,s,a 0 0 1 691r
imp:n,p 0 0 1 691r
1 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 2.05e5 1 26r $ wall segments
imp:h,/,d,t,s,a 0 0 1 705r
imp:n,p 0 0 1 705r
cut:h,/,d,t,s,a 1e+08 1e-3
c cut:n j 0.99E-6
cut:n j 1e-11
@@ -158,6 +163,7 @@ c * McStas source terms for gamma doeses
c read file=sources/mcstas_source_vs.i echo
c read file=sources/mcstas_source_sample.i echo
c read file=sources/mcstas_gammas.i echo
c read file=sources/mcstas_gammas_2.i echo
prdmp 5e7 -60 1 2 5e6
print 30 32 $ 120 128 150 160 161 162
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
geometry/ess_cells.i
+4 -9
View File
@@ -5,15 +5,10 @@ c Main ESS plant geometry without the instrument
159 ) -1 ( -154 : 161 : 166 ) ( -62 : 64 : 63 ) ( -3 : -103 :
-70 : -64 : 105 : 104 : 72 ) ( -121 : -65 : -63 : 122 : 101 :
124 )
#((68 -72 8209 -70 -104 -2)
:(72 -8101 8201 8202 8203 8209 -2)
:(8709 8710 8701 8703 8706 8707 8708 8709 8710 (8704 : 8705))
:(8709 8710 8702 8703 8711 8712)
:(8709 8809 -8311 -8712 -8701 -8313)
:(8101 -8709 8301 8302 8303 8209) 3) 600 #707
#(72 -8709 -104 -8201 8203 8204)
#(((112 -8202):(8202 -8302 8101)) 72 -8709 8203 8204)
#(8709 8720 112 -104 8703 8701)
#(-104 112 72 8710 8209 8203)
#( 104 8708 8707 8710 8703 8701)
#(68 -72 8209 -70 -104 -2)
3 118 0.122325 24 4 -19
4 0 159 -137 -21 153 -154
5 118 0.122325 37 -5 -19
@@ -496,7 +491,7 @@ c Main ESS plant geometry without the instrument
-112
257 200 -3.9 ( -3 : -69 : -2 : -70 : -64 : 71 : 68 ) 3 -102 2 103
66 -67 64
258 49 0.074981 64 -67 66 -104 101 -72 3
258 49 0.074981 -104 101 -67 ((64 66 -72 3):(72 -610 59))
259 49 0.074981 ( -3 : -69 : -2 : -70 : -64 : 71 : 72 ) 3 -72 103
-104 -66 70 64
260 200 -3.9 -67 66 -107 -2 -72 68 3
geometry/ess_surfaces.i
+8 -2
View File
@@ -100,10 +100,13 @@ c 70 pz -334.4
98 pz 341.389474
99 pz 351.126316
100 pz 360.863158
101 p -0.906307787 -0.422618262 0.0 -5.52451529e-14
c Building Wall E01
c 101 p -0.906307787 -0.422618262 0.0 -5.52451529e-14
101 p -0.866025 -0.5 0.0 80.345
102 cz 2450
103 p 0.0 -1 0.0 -80
104 p -0.906307787 -0.422618262 0.0 80
c 104 p -0.906307787 -0.422618262 0.0 80
104 p -0.866025 -0.5 0.0 180.345
105 pz 375.6
106 cz 2770
107 p -0.125718017 -0.992066016 0.0 0.0
@@ -379,3 +382,6 @@ c 463 pz 13.5
572 p 0.809016994 0.587785252 0.0 17.5262
575 pz 7.4
600 s 1603.5 2265.0 13.0 3.0
c additional surfaces from CATIA
c D01/D02 boundary
610 px 1090
geometry/estia_cave.i
+9 -4
View File
@@ -1,7 +1,7 @@
c Experimental cave
c First largest area: up-stream/down-stream/roof/floor/walls
c then cut corner and add extended roof
8701 49 0.074981 ((8709 8710 8701 8703 8706 8707 8708 8709 8710 (8704 : 8705))
8701 4 -1.2E-3 ((8709 8710 8701 8703 8706 8707 8708 8709 8710 (8704 : 8705))
:(8709 8710 8702 8703 8711 8712))
#((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810
(-8804 : -8805))
@@ -13,10 +13,15 @@ c then cut corner and add extended roof
(-8804 : -8805))
:(-8809 -8810 -8802 -8803 -8811 -8812)
:(-8311 -8312 8709 8809 -8313 -8803)) 3
#(8621 8622 8623 8624 8626)
#6111 #8501 #8511 #9001
c sample position
9001 19 -8.65 -9001
c 9001 4 -1.2E-3 -9001
c area around cave
9900 0 8709 8710 -104 (-8708:-8707) 8703 8701
9901 0 8709 8710 112 ((-8704 -8705):-8706) 8703 8701
9902 0 8720 -8710 112 -104 8703 8701
9900 4 -1.2E-3 8709 8710 -104 (-8708:-8707) 8703 8701
9901 4 -1.2E-3 8709 8710 112 ((-8704 -8705):-8706) 8703 8701
9902 4 -1.2E-3 8720 -8710 112 -104 8703 8701
c visualize Selene gamma sources (comment for real runs)
c 50001 0 -50001
c 50002 0 -50002
geometry/estia_guides.i
+52 -33
View File
@@ -13,19 +13,23 @@ c -----Wall Cut-Out---
6003 4 -1.2E-3 (6103 -72 -6601 -6602 -6603 -6604 -6605) 3
#(7103 8405 -8921)
#(7102 -7103 7205 7206 7207 7208)
#(-8405 8401 8402 8403 8404)
#(8421 8422 8423 8424 8425 8426)
c #(-8405 8401 8402 8403 8404)
c Vacuum before and inside Selene Guide 1
6101 0 ((7103 8405 -8922 3)
:(-8405 8406 -8411 8412 8413 8402 3))
#706 #8311 #8411 #8412 #8413
6102 5 0.0582256 (7103 8405 8922 -8921 3)
c Vaccum after Selene 1, until and inside Selen 2
6111 0 ((-8406 -8100 -8932 3)
6101 0 (7103 8435 -8922 3) #706
6102 5 0.0582256 (7103 -8425 8922 -8921 3)
6103 0 (-8431 -8432 -8433 -8434 -8435 -8436 (-8441:-8442) (-8443:-8444))
#8301 #8311 #8411 #8412 #8413
c Vacuum after Selene 1, until Selen 2
6111 0 ((8436 -8100 -8932 3)
:(8100 -8101 -8501 -8502 -8503 -8504 3)
:(8101 -8121 -8932 3)
:(8121 -8122 8601 -8611 -8612 -8613 3))
#707 #8511
6112 5 0.0582256 (-8406 -8121 8932 -8931 3 #(8100 -8101))
:(8101 8635 -8932 3))
#707
6112 5 0.0582256 (-8426 -8625 8932 -8931 3 #(8100 -8101))
c Vacuum inside Selene 2
6113 0 ((-8631 -8632 -8633 -8634 -8635 -8636
#(8121 -8122 8601 8602 8603 8604))
:(8121 -8122 8601 -8611 -8612 -8613 3)) #8511
c -----Bunker Wall Collimator (BWC)-----
7001 0 (7100 -7103 -7301 -7302 -7303 -7304 3)
#(7101 -7401 7402)
@@ -48,38 +52,35 @@ c ---- Outside of Bunker
c -- Shielding Walls --
c concrete floor
8001 49 0.074981 ((68 -72 8209 -70 -104 -2)
:(72 -8101 8209 -59 8201 8202)
:(8101 -8709 8209 -59 8301 8302)) 3
:(72 8710 8707 112 8209 (-8304:(-8305 -610 -101))))
3 -104
c ceiling outside
8002 49 0.074981 ( 72 -8101 8207 8203 8201 8202)
8002 4 -1.2E-3 ( 72 -8101 8207 8203 8201 8202)
c left outer wall (E01)
8003 4 -1.2E-3 ( 72 -8100 8201 -8205 59 -8207 3)
8003 4 -1.2E-3 ( 72 -8100 8201 -8205 59 -8207 3) #(-8305 -610)
c Selene 1/2 separator
8004 200 -3.9 (8100 -8101 8201 8202 59 -8207 3)
#(8100 -8101 -8221 -8222 -8223 -8224)
c Optical cave
8005 200 -3.9 (8101 -8709 59 8301 8302 8303 3)
8005 4 -1.2E-3 (8101 -8709 59 8201 8302 8203 3)
#(8101 -8709 -8311 -8312 -8313 -8314)
c air
8008 4 -1.2E-3 (72 -8100 -8211 -8212 -8213 -8214 3)
#(-8406 -8931)
#(-8405 8406 8401 8402 8403 8404)
c #8370 #8371 #8372
#(8421 8422 8423 8424 8425 8426)
#(-8426 -8931) #(-8305 -610)
8009 4 -1.2E-3 (8101 -8709 -8311 -8312 -8313 -8314 3)
#(8121 -8122 8601 8602 8603 8604)
#(-8121 -8931)
#(8621 8622 8623 8624 8625 8626)
#707 #710
c B4C concrete inner layers
8011 4 -1.2E-3 (8211 8215 -8213 -8214 72 -8100 3)
8012 16 -2.4 (8205 8202 8213 -8218 72 -8100 3)
8013 4 -1.2E-3 ((-8216 8111 -8100):(72 -8111 18202))
8212 -8213 -8214 3
8014 16 -2.4 (8202 -8206 8217 -8213 -8214 8111 -8100 3)
c Steel middle layers
8020 71 0.0843223 (-18202 8202 -8213 -8214 72 -8111 3)
8021 71 0.0843223 (-8215 8205 -8213 -8214 72 -8100 3)
8022 71 0.0843223 (8205 8202 8218 -8207 72 -8100 3)
8023 71 0.0843223 (8216 8202 (-8217:8206) -8213 -8214 8111 -8100 3)
c Selene 1 shielding
8011 4 -1.2E-3 (8211 8205 -8213 -8214 72 -8100 3) #(-8305 -610)
8012 16 -2.400 (8205 8202 8213 -8218 72 -8100 3)
8013 4 -1.2E-3 72 8202 8212 8217 -8213 -8214 3 #(-8305 -610)
c
8020 71 0.0843223 (-8217 -8100 8215 8202 -8213 -8214 3)
8021 71 0.0843223 (-8217 -8100 -8215 8212 -8213 -8214 3)
8022 71 0.0843223 (8205 8202 8218 -8207 72 -8100 3)
c
c Selene 1 carrier
8301 18 -2.4 -8405 8406 ((8401 8402 8403 8404)
@@ -87,6 +88,14 @@ c Selene 1 carrier
c Selene 1 glass (change material!!)
8311 12 -2.51 -8405 8406 -8903 -8904 -8518
#(-8901 -8902) #(-8515 -8517) 3
c Selene 1 Vacuum tank (stainless steel)
8321 3 0.0876276 ( 8421 8422 8423 8424 8425 8426)
#(-8431 -8432 -8433 -8434 -8435 -8436)
#(8425 8435 -8922)
#(8426 8436 -8932)
c Lead shielding layer at top-right corner of tank
c 13 -11.34
8422 0 -8435 -8436 -8432 -8433 ((8441 8442):(8443 8444))
c Reference sensor material for Selene 1 dose simulations
c 8370 71 0.0843223 -8370
c 8371 71 0.0843223 -8371
@@ -110,9 +119,19 @@ c Selene 2 carrier
c Selene 2 glass (change material!!)
8511 12 -2.51 8121 -8122 -8913 -8914 8915
#(-8911 -8912) #(-8515 -8517) 3
c Selene 2 Vacuum tank (stainless steel)
8521 3 0.0876276 ( 8621 8622 8623 8624 8625 8626)
#(-8631 -8632 -8633 -8634 -8635 -8636)
#(8625 8635 -8932)
c #(8626 8636 -8932)
c Air outside the guide shielding
8601 0 72 -8709 -104 -8201 8203 8204
8601 4 -1.2E-3 ((72 -8709 -104 -8201 8203 8204)
:(-8701 8203 8710 8709 -104 8302))
#(101 -610 -67) #(-8305 -610)
#(8302 8702 (8712:(8809 -8311 -8313)))
#708
8602 0 ((112 -8202)
:(8202 -8302 8101)) 72 -8709 8203 8204
8602 4 -1.2E-3 ((112 -8202 -8709)
:(8202 -8302 8101 -8709)
:(-8701 8203 8710 8709 -8302 112))
72 8203 8204 #(-8305 -610)
#709
geometry/estia_surfaces.i
+83 -22
View File
@@ -471,7 +471,7 @@ c Sides 4
c Sides 5
6601 p -0.809017 -0.587785 0.000000 85.500
6602 p 0.809017 0.587785 0.000000 41.800
6603 pz 76.000
6603 pz 85.0
6604 p 0.000000 0.000000 -1.000000 81.000
6605 p -0.987688 0.156434 0.000000 -895.551
c Cu-collimator
@@ -508,8 +508,8 @@ c Wall radii
8100 cz 2112.01
8101 cz 2157.17
c 8102 cz 3023.41
8103 cz 1355.64
8104 cz 2075.58
c 8103 cz 1355.64
c 8104 cz 2075.58
8111 cz 1655.0
8112 cz 1665.0
8113 cz 1710.0
@@ -521,9 +521,9 @@ c 8102 cz 3023.41
c Shielding wall geometry
c first section outside
8201 p 0.813101 0.582123 0.000000 -148.964
8202 p -0.776046 -0.630676 0.000000 1.917
18202 p -0.776046 -0.630676 0.000000 16.917
8203 p 0.000000 0.000000 -1.000000 -166.400
8202 p -0.771625 -0.636078 0.000000 -21.127
c 18202 p -0.776046 -0.630676 0.000000 16.917
8203 p 0.000000 0.000000 -1.000000 -666.400
8204 pz -200.000
c steel/concrete interface
8205 p 0.813101 0.582123 0.000000 -98.964
@@ -533,14 +533,14 @@ c 8208 p -0.776046 -0.630676 0.000000 31.917
8209 pz -300.000
c inside
8211 p -0.813101 -0.582123 0.000000 68.794
8212 p 0.813101 0.582123 0.000000 71.037
8212 p 0.790690 0.612217 0.000000 0.000
8213 pz 85.900
8214 p 0.000000 0.000000 -1.000000 200.000
c B4C concrete/steel interface
8215 p 0.813101 0.582123 0.000000 -88.964
8216 p 0.813101 0.582123 0.000000 91.037
8217 p 0.813101 0.582123 0.000000 101.037
8218 pz 105.900
8215 p 0.739506 0.673150 0.000000 -102.437
8216 pz 50.000
8217 p -0.612217 0.790690 0.000000 -1958.974
8218 pz 105.900
c 8219 p -0.776046 -0.630676 0.000000 111.917
c wall breakthrough
8221 p -0.813101 -0.582123 0.000000 32.866
@@ -551,8 +551,10 @@ c Optical Cave Part
c Outside
8301 p 0.813101 0.582123 0.000000 -147.866
8302 p -0.813101 -0.582123 0.000000 -102.134
8303 p 0.000000 0.000000 -1.000000 -132.500
8303 p 0.000000 0.000000 -1.000000 -666
c -132.500
8304 pz -200.000
8305 pz -150.000
c Inside
8311 p -0.813101 -0.582123 0.000000 117.866
8312 p 0.813101 0.582123 0.000000 72.134
@@ -578,6 +580,26 @@ c Inside
8411 p -0.813101 -0.582123 0.000000 14.866
8412 p -0.007112 0.009934 -0.999925 -38.883
8413 p 0.007112 -0.009934 0.999925 -11.117
c Selene 1 Vacuum
c Outside
8421 p 0.813101 0.582123 0.000000 -66.866
8422 p -0.813101 -0.582123 0.000000 -35.134
8423 p -0.007112 0.009934 -0.999925 -96.883
8424 p 0.007112 -0.009934 0.999925 -52.617
8425 p 0.582080 -0.813040 -0.012217 1345.888
8426 p -0.582080 0.813040 0.012217 -2097.888
c Inside
8431 p -0.813101 -0.582123 0.000000 64.166
8432 p 0.813101 0.582123 0.000000 31.934
8433 p 0.007112 -0.009934 0.999925 94.183
8434 p -0.007112 0.009934 -0.999925 50.917
8435 p -0.582080 0.813040 0.012217 -1347.388
8436 p 0.582080 -0.813040 -0.012217 2096.388
c Anti-gamma lead shield concept
8441 p 0.813101 0.582123 0.000000 28.934
8442 p 0.007112 -0.009934 0.999925 44.183
8443 p 0.007112 -0.009934 0.999925 91.183
8444 p 0.813101 0.582123 0.000000 -9.934
c Beam through wall
c Inside
8501 p -0.792291 -0.610143 -0.000426 79.405
@@ -597,19 +619,39 @@ c 8516 p -0.007112 0.009934 -0.999925 -5.400
8519 p 0.813101 0.582123 0.000000 -9.866
c Selene 2 Carrier
c Outside
8601 p 0.813101 0.582123 0.000000 -1.866
8602 p -0.813101 -0.582123 0.000000 -58.134
8603 p -0.007112 0.009934 -0.999925 -68.883
8604 p 0.007112 -0.009934 0.999925 -41.117
8601 p 0.813101 0.582123 0.000000 -2.7658
8602 p -0.813101 -0.582123 0.000000 -42.9342
8603 p -0.007112 0.009934 -0.999925 -58.883
8604 p 0.007112 -0.009934 0.999925 -31.117
c Inside
8611 p 0.813101 0.582123 0.000000 19.634
8612 p 0.007112 -0.009934 0.999925 31.383
8613 p -0.007112 0.009934 -0.999925 1.617
c Experimental Cave
8611 p 0.813101 0.582123 0.000000 19.1342
8612 p 0.007112 -0.009934 0.999925 38.8829
8613 p -0.007112 0.009934 -0.999925 11.1171
c Guide field plate
8616 p -0.813101 -0.582123 0.000000 -0.1342
8617 p 0.813101 0.582123 0.000000 0.8658
8618 p -0.007112 0.009934 -0.999925 -30.5829
8619 p 0.007112 -0.009934 0.999925 -2.8171
c Selene 2 Vacuum
c Outside
8621 p 0.813101 0.582123 0.000000 -30.866
8622 p -0.813101 -0.582123 0.000000 -71.134
8623 p -0.007112 0.009934 -0.999925 -96.883
8624 p 0.007112 -0.009934 0.999925 -52.617
8625 p 0.582080 -0.813040 -0.012217 2545.388
8626 p -0.582080 0.813040 0.012217 -3297.388
c Inside
8631 p -0.813101 -0.582123 0.000000 28.166
8632 p 0.813101 0.582123 0.000000 67.934
8633 p 0.007112 -0.009934 0.999925 94.183
8634 p -0.007112 0.009934 -0.999925 50.117
8635 p -0.582080 0.813040 0.012217 -2546.888
8636 p 0.582080 -0.813040 -0.012217 3295.888
c Experimental Cave
c Outside
8701 p 0.000000 0.000000 -1.000000 -87.940
8702 p 0.000000 0.000000 -1.000000 -136.940
8703 pz -300.000
8703 pz -200.000
8704 p -0.813101 -0.582123 0.000000 -142.134
8705 p 0.582123 -0.813101 0.000000 3364.795
8706 p -0.736097 -0.676876 0.000000 196.513
@@ -620,7 +662,7 @@ c Outside
8711 p -0.813101 -0.582123 0.000000 -102.134
8712 p 0.813101 0.582123 0.000000 -77.866
8720 p -0.582123 0.813101 0.000000 -4464.795
c Inside
c Inside
8801 pz 57.941
8802 pz 108.941
8803 p 0.000000 0.000000 -1.000000 166.960
@@ -652,4 +694,23 @@ c Vacuum vessels before S1/S1-S2
8932 3 c/x 0 -2 15.0
c Sphere at sample position to simulate gamma generation
9001 15 s 0 0 0 12.0
c UNIVERSES
c Selene 1
18301 3 px -295
18302 3 px -247
18311 3 c/y -291.25 1.75 1.8
18312 3 c/y -291.25 12.35 1.8
18313 3 c/y -291.25 18.75 1.8
18314 3 c/y -291.25 -1.75 1.8
18315 3 c/y -291.25 -12.35 1.8
18316 3 c/y -291.25 -18.75 1.8
18321 3 c/y -251.75 22.25 1.8
18322 3 c/y -251.75 12.35 1.8
18323 3 c/y -251.75 18.75 1.8
18324 3 c/y -251.75 -22.25 1.8
18325 3 c/y -251.75 -12.35 1.8
18326 3 c/y -251.75 -18.75 1.8
c visualize Selene gamma sources
50001 1 rpp 1339.25 2060.75 2 8 2 8
50002 1 rpp 2539.25 3260.75 -8 -2 -8 -2
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
geometry/estia_universes.i
+17
View File
@@ -0,0 +1,17 @@
c UNIVERSES for Selene geometry lattices
c Selene 1 holes for adjusters
18300 0 18301 -18302 U=1 fill=2 lat=1
18301 18 -2.4 18311 18312 18313 18314 18315 18316
18321 18322 18323 18324 18325 18326 U=2
18311 0 -18311:-18314 U=2
18312 0 -18312 U=2
18313 0 -18313 U=2
c 18314 0 -18314 U=2
18315 0 -18315 U=2
18316 0 -18316 U=2
18321 0 -18321:-18323 U=2
18322 0 -18322 U=2
c 18323 0 -18323 U=2
18324 0 -18324:-18326 U=2
18325 0 -18325 U=2
c 18326 0 -18326 U=2
plot_tally.py
+63 -12
View File
@@ -13,7 +13,7 @@ from pylab import *
from matplotlib.colors import LogNorm
COLORMAP='tab20c'
SCALING=1.56e16/8.7374e15
SCALING=1.0
def read_tally(fname, use_tally=None):
if fname.endswith('.gz'):
@@ -68,9 +68,16 @@ def plot_xy(outfile, res, zidx=5, Imin=1e-6, Imax=1e4, cticks=None):
#fig.savefig(outfile)
#pc=p.get_facecolors()
#pc[:,3]=where(I[:,:,zidx].T.flatten()>0, 0.9-0.25*dI[:,:,zidx].T.flatten(), 0.)
ax.add_patch(Circle((0,0), radius=5.50, color='black', fill=False, lw=2))
ax.add_patch(Circle((0,0), radius=11.50, color='black', fill=False, lw=2))
ax.add_patch(Circle((0,0), radius=15.00, color='black', fill=False, lw=2))
ax.add_patch(Circle((-0.151,0.021), radius=5.50, color='black', fill=False, lw=2))
ax.add_patch(Circle((-0.151,0.021), radius=11.50, color='black', fill=False, lw=2))
ax.add_patch(Circle((-0.151,0.021), radius=15.00, color='black', fill=False, lw=2))
# shielding walls outside bunker
ax.add_line(Line2D([21.12-0.151, 21.12-0.151], [-1.5,1.5], color='black', lw=2))
ax.add_line(Line2D([14.845,19.426], [-0.541,-0.713], color='black', lw=2))
ax.add_line(Line2D([19.426,19.405], [-0.713,-1.248], color='black', lw=2))
ax.add_line(Line2D([19.405,21.205], [-1.248,-1.460], color='black', lw=2))
ax.add_line(Line2D([14.723, 22.740], [2.232, 3.084], color='black', lw=2))
ax.add_line(Line2D([14.801, 20.894], [-1.210, -1.724], color='black', lw=2))
if outfile is not None:
fig.savefig(outfile, transparent=False, dpi=300)
@@ -89,9 +96,12 @@ def plot_xz(outfile, res, yidx=5, Imin=1e-6, Imax=1e4, cticks=None):
#fig.savefig(outfile)
#pc=p.get_facecolors()
#pc[:,3]=where(I[:,yidx].T.flatten()>0, 0.9-0.25*dI[:,yidx].T.flatten(), 0.)
ax.plot([sqrt(5.5**2-ypos**2), sqrt(5.5**2-ypos**2)], [z.min()/100., z.max()/100.], color='black', lw=2)
ax.plot([sqrt(11.5**2-ypos**2), sqrt(11.5**2-ypos**2)], [z.min()/100., z.max()/100.], color='black', lw=2)
ax.plot([sqrt(15.0**2-ypos**2), sqrt(15.0**2-ypos**2)], [z.min()/100., z.max()/100.], color='black', lw=2)
ax.plot([sqrt(5.5**2-ypos**2)-0.151, sqrt(5.5**2-ypos**2)-0.151],
[z.min()/100., z.max()/100.], color='black', lw=2)
ax.plot([sqrt(11.5**2-ypos**2)-0.151, sqrt(11.5**2-ypos**2)-0.151],
[z.min()/100., z.max()/100.], color='black', lw=2)
ax.plot([sqrt(15.0**2-ypos**2)-0.151, sqrt(15.0**2-ypos**2)-0.151],
[z.min()/100., z.max()/100.], color='black', lw=2)
if outfile is not None:
fig.savefig(outfile, transparent=False, dpi=300)
@@ -112,26 +122,67 @@ def plot_yz(outfile, res, xidx=5, Imin=1e-6, Imax=1e4, cticks=None):
if outfile is not None:
fig.savefig(outfile, transparent=False, dpi=300)
def get_name(fname):
# extract tag-name and version info from filename
stripped=os.path.basename(fname)[9:].split('.')[0]
tag, version=stripped.split('-', 1)
return tag, version
if __name__=='__main__':
import sys, os
res=read_tally(sys.argv[1], use_tally=1)
if '-o' in sys.argv:
idx=sys.argv.index('-o')
sys.argv.pop(idx)
prefix=sys.argv.pop(idx)
else:
prefix='estia_dosemap'
print sys.argv[1]
if ':' in sys.argv[1]:
fname, tally=sys.argv[1].split(':',2)
tally=int(tally)
else:
fname=sys.argv[1]
tally=None
res=read_tally(fname, use_tally=tally)
x,y,z,I,dI=res
I*=SCALING
fpre=os.path.join('.', 'images', os.path.basename(sys.argv[1])[:-4])
fpath=os.path.join('.', 'images', get_name(fname)[0]+'-'+get_name(fname)[1])
if tally is not None:
fpath+='[%s]'%tally
if len(sys.argv)>2:
for fi in sys.argv[2:]:
print fi
if ':' in fi:
fi, tally=fi.split(':',2)
tally=int(tally)
else:
tally=None
resi=read_tally(fi, use_tally=tally)
xi,yi,zi,Ii,dIi=resi
Ii*=SCALING
I+=Ii
dI=sqrt(dI**2+dIi**2)
res=(x,y,z,I,dI)
fpath+='+'+get_name(fi)[1]
if tally is not None:
fpath+='[%s]'%tally
os.mkdir(fpath)
fpre=os.path.join(fpath, prefix)
xi, yi, zi=res[3].shape
for zidx in range(zi):
print fpre+'_xy_%02i.png'%zidx
figure(figsize=(12,8))
plot_xy(fpre+'_xy_%02i.png'%zidx, res, zidx=zidx,
Imin=1e-8, Imax=1e2, cticks=[1e-8, 1e-6, 1e-4, 1e-2, 1, 100])
Imin=1e-7, Imax=1e-2, cticks=[1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2])
for yidx in range(yi):
print fpre+'_xz_%02i.png'%yidx
figure(figsize=(12,8))
plot_xz(fpre+'_xz_%02i.png'%yidx, res, yidx=yidx,
Imin=1e-8, Imax=1e2, cticks=[1e-8, 1e-6, 1e-4, 1e-2, 1, 100])
Imin=1e-7, Imax=1e-2, cticks=[1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2])
for xidx in range(xi):
print fpre+'_yz_%02i.png'%xidx
figure(figsize=(8,10))
plot_yz(fpre+'_yz_%02i.png'%xidx, res, xidx=xidx,
Imin=1e-8, Imax=1e2, cticks=[1e-8, 1e-6, 1e-4, 1e-2, 1, 100])
Imin=1e-7, Imax=1e-2, cticks=[1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2])
sources/e2_2m.i
+2 -2
View File
@@ -1,10 +1,10 @@
c source definition based on tally_181.txt
c but adapted to model geometry and actual size of Estia beam extraction
c make sure to supply the correct WWINP file
WWP:n j j 5 j -1 j 1.0
WWP:n j j 5 j -1 j 1.79
c the source definition with neutrons/s
sdef dir=d1 erg=fdir d2 x=199.98 y=d23 z=d24 vec=1 0 0 tr=101
par=n ara=165.0 wgt=3.6041775E13
par=n ara=165.0 wgt=6.436556E13
si1 -1 0 0.5 0.866 0.9659 0.9848
0.99619 0.9975640 0.998629 0.999390 0.999847 1
sp1 0 5.01272E+12 1.04029E+12 2.01432E+12 1.39792E+12
sources/mcstas_gammas.i
+58 -59
View File
@@ -1,145 +1,144 @@
c source definition based on tally_181.txt
c but adapted to model geometry and actual size of Estia beam extraction
c source definition based on McStas
c make sure to supply the correct WWINP file
WWP:n j j 5 j -1 j 1.0
c the source definition with neutrons/s
c the source definition with gammas/s
sdef erg=fx d1 x=d31 y=d32 z=d33 tr=1 par=p wgt=9.59833e+09
ds1 q 1339.25 100
1387.35 101 1435.45 102 1483.55 103 1531.65 104 1579.75 105
1627.85 106 1675.95 107 1724.05 108 1772.15 109 1820.25 110
1387.35 101 1435.45 102 1483.55 103 1531.65 104 1579.75 105
1627.85 106 1675.95 107 1724.05 108 1772.15 109 1820.25 110
1868.35 111 1916.45 112 1964.55 113 2012.65 114 2060.75 115
c Energy distribution for all positions along the guide
si100 0 1
sp100 0 1e-20
si101 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si101 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp101 0.00000E+00 1.56557E+05 2.76591E+05 8.82267E+05 3.79148E+06
7.22100E+09 7.15782E+04 2.36360E+05 1.64987E+06 1.20477E+07
5.83592E+06 5.13734E+06 1.05131E+07 9.13771E+06 2.23406E+06
5.55688E+06 3.73177E+06 9.84371E+06 3.67345E+04 5.62062E+04
C ------------------------------------------------------------------------
si102 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si102 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp102 0.00000E+00 1.75066E+05 3.14271E+05 1.00187E+06 4.19715E+06
4.29576E+08 6.49313E+04 8.35318E+04 1.59407E+06 1.13391E+07
3.24825E+06 2.40414E+06 3.23346E+06 2.35587E+06 1.81032E+06
4.22018E+06 2.81568E+06 1.08595E+07 3.50283E+04 1.16665E+04
C ------------------------------------------------------------------------
si103 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si103 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp103 0.00000E+00 1.14282E+05 2.09752E+05 6.68141E+05 2.71818E+06
2.88507E+08 4.31570E+04 5.55428E+04 1.05765E+06 7.30837E+06
2.10495E+06 1.58164E+06 2.11537E+06 1.54343E+06 1.20473E+06
2.79310E+06 1.87709E+06 7.24741E+06 2.27791E+04 7.70741E+03
C ------------------------------------------------------------------------
si104 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si104 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp104 0.00000E+00 1.19942E+05 2.24594E+05 7.14916E+05 2.83168E+06
2.69546E+08 4.59551E+04 5.82997E+04 1.12511E+06 7.56863E+06
2.17669E+06 1.65382E+06 2.17880E+06 1.58300E+06 1.28282E+06
2.95750E+06 2.00035E+06 7.75812E+06 2.37907E+04 7.90258E+03
C ------------------------------------------------------------------------
si105 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si105 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp105 0.00000E+00 1.02808E+05 1.94757E+05 6.19690E+05 2.41655E+06
2.29099E+08 3.97544E+04 5.03268E+04 9.72514E+05 6.44006E+06
1.85566E+06 1.42096E+06 1.86323E+06 1.35351E+06 1.11024E+06
2.55211E+06 1.73259E+06 6.72701E+06 2.03449E+04 6.77894E+03
C ------------------------------------------------------------------------
si106 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si106 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp106 0.00000E+00 8.86042E+04 1.69642E+05 5.39580E+05 2.07421E+06
1.99597E+08 3.45601E+04 4.37451E+04 8.44765E+05 5.51349E+06
1.59281E+06 1.22897E+06 1.60647E+06 1.16764E+06 9.65712E+05
2.21413E+06 1.50825E+06 5.85933E+06 1.74992E+04 5.87135E+03
C ------------------------------------------------------------------------
si107 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si107 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp107 0.00000E+00 7.87877E+04 1.51840E+05 4.82849E+05 1.83971E+06
1.73871E+08 3.08863E+04 3.89901E+04 7.54666E+05 4.88100E+06
1.41072E+06 1.09310E+06 1.42351E+06 1.03396E+06 8.63179E+05
1.97565E+06 1.34859E+06 5.24414E+06 1.55375E+04 5.20473E+03
C ------------------------------------------------------------------------
si108 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si108 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp108 0.00000E+00 5.01840E+04 9.68964E+04 3.08110E+05 1.17096E+06
1.12987E+08 1.97075E+04 2.49206E+04 4.81428E+05 3.10577E+06
8.98765E+05 6.97579E+05 9.09029E+05 6.60776E+05 5.50891E+05
1.26041E+06 8.60883E+05 3.34661E+06 9.89468E+03 3.33164E+03
C ------------------------------------------------------------------------
si109 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si109 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp109 0.00000E+00 4.11446E+04 7.95509E+04 2.52944E+05 9.59525E+05
9.23981E+07 1.61749E+04 2.04454E+04 3.95097E+05 2.54402E+06
7.36328E+05 5.72028E+05 7.44921E+05 5.41445E+05 4.52162E+05
1.03416E+06 7.06655E+05 2.74751E+06 8.11001E+03 2.73082E+03
C ------------------------------------------------------------------------
si110 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si110 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp110 0.00000E+00 2.51204E+04 4.88245E+04 1.55217E+05 5.84610E+05
5.28251E+07 9.90982E+03 1.24431E+04 2.42028E+05 1.54694E+06
4.46975E+05 3.48156E+05 4.50560E+05 3.26731E+05 2.76963E+05
6.32447E+05 4.32872E+05 1.68610E+06 4.94348E+03 1.64531E+03
C ------------------------------------------------------------------------
si111 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si111 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp111 0.00000E+00 2.36608E+04 4.56443E+04 1.45144E+05 5.52264E+05
4.88854E+07 9.27597E+03 1.16381E+04 2.26683E+05 1.46399E+06
4.22049E+05 3.26897E+05 4.23743E+05 3.07060E+05 2.59131E+05
5.92778E+05 4.04758E+05 1.57629E+06 4.66257E+03 1.54109E+03
C ------------------------------------------------------------------------
si112 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si112 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp112 0.00000E+00 2.34986E+04 4.48276E+04 1.42601E+05 5.50860E+05
4.89730E+07 9.13069E+03 1.14783E+04 2.23305E+05 1.46457E+06
4.21412E+05 3.23890E+05 4.21794E+05 3.05694E+05 2.54959E+05
5.84877E+05 3.97945E+05 1.54818E+06 4.64119E+03 1.52926E+03
C ------------------------------------------------------------------------
si113 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si113 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp113 0.00000E+00 1.87764E+04 3.52661E+04 1.12245E+05 4.42775E+05
3.94292E+07 7.20582E+03 9.07992E+03 1.76424E+05 1.18187E+06
3.39132E+05 2.57888E+05 3.37896E+05 2.44903E+05 2.01079E+05
4.63087E+05 3.13513E+05 1.21806E+06 3.72001E+03 1.21944E+03
C ------------------------------------------------------------------------
si114 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si114 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp114 0.00000E+00 1.82897E+04 3.37346E+04 1.07439E+05 4.34219E+05
4.04983E+07 6.92221E+03 8.78390E+03 1.69672E+05 1.16466E+06
3.33775E+05 2.50952E+05 3.32061E+05 2.41077E+05 1.93078E+05
4.46799E+05 3.00726E+05 1.16534E+06 3.63777E+03 1.19682E+03
C ------------------------------------------------------------------------
si115 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
si115 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp115 0.00000E+00 1.08301E+04 1.94946E+04 6.21412E+04 2.59394E+05
2.48374E+07 4.02189E+03 5.13569E+03 9.87454E+04 6.99903E+05
@@ -147,17 +146,17 @@ sp115 0.00000E+00 1.08301E+04 1.94946E+04 6.21412E+04 2.59394E+05
2.61025E+05 1.74300E+05 6.73544E+05 2.16453E+03 7.10515E+02
C ------------------------------------------------------------------------
c x-distribution
si31 1339.25
1387.35 1435.45 1483.55 1531.65 1579.75
1627.85 1675.95 1724.05 1772.15 1820.25
si31 1339.25
1387.35 1435.45 1483.55 1531.65 1579.75
1627.85 1675.95 1724.05 1772.15 1820.25
1868.35 1916.45 1964.55 2012.65 2060.75
sp31 0
7.29220E+09 4.79340E+08 3.21180E+08 3.03858E+08 2.58578E+08
2.25073E+08 1.96544E+08 1.27443E+08 1.04253E+08 6.00576E+07
5.56826E+07 5.57082E+07 4.47933E+07 4.57106E+07 2.79109E+07
sp31 0
7.29220E+09 4.79340E+08 3.21180E+08 3.03858E+08 2.58578E+08
2.25073E+08 1.96544E+08 1.27443E+08 1.04253E+08 6.00576E+07
5.56826E+07 5.57082E+07 4.47933E+07 4.57106E+07 2.79109E+07
c y-distribution
si32 -6 -1.0
si32 2 8
sp32 0 1
c z-distribution
si33 -6 6
si33 2 8
sp33 0 1
sources/mcstas_gammas_2.i
+162
View File
@@ -0,0 +1,162 @@
c source definition based on McStas
c make sure to supply the correct WWINP file
WWP:n j j 5 j -1 j 1.0
c the source definition with gammas/s
sdef erg=fx d1 x=d31 y=d32 z=d33 tr=1 par=p wgt=1.04957e+09
ds1 q 2539.25 100
2587.35 101 2635.45 102 2683.55 103 2731.65 104 2779.75 105
2827.85 106 2875.95 107 2924.05 108 2972.15 109 3020.25 110
3068.35 111 3116.45 112 3164.55 113 3212.65 114 3260.75 115
c Energy distribution for all positions along the guide
si100 0 1
sp100 0 1e-20
si101 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp101 0.00000E+00 6.74946E+03 1.26349E+04 4.02190E+04 1.59434E+05
4.37259E+07 2.64491E+03 3.95862E+03 6.43047E+04 4.33561E+05
1.34669E+05 1.05571E+05 1.54428E+05 1.18439E+05 7.48297E+04
1.73991E+05 1.17757E+05 4.37637E+05 1.36106E+03 6.36018E+02
C ------------------------------------------------------------------------
si102 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp102 0.00000E+00 1.17151E+04 2.23209E+04 7.10081E+04 2.74744E+05
1.91981E+07 4.53666E+03 5.59356E+03 1.11043E+05 7.29343E+05
2.07983E+05 1.59135E+05 2.04496E+05 1.47055E+05 1.26489E+05
2.89983E+05 1.97214E+05 7.70669E+05 2.31033E+03 7.27223E+02
C ------------------------------------------------------------------------
si103 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp103 0.00000E+00 1.20232E+04 2.32095E+04 7.38021E+04 2.80546E+05
2.04130E+07 4.70687E+03 5.81143E+03 1.15089E+05 7.42422E+05
2.12512E+05 1.64201E+05 2.10375E+05 1.51470E+05 1.31339E+05
3.00162E+05 2.04993E+05 8.01335E+05 2.36550E+03 7.53558E+02
C ------------------------------------------------------------------------
si104 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp104 0.00000E+00 1.75885E+04 3.43282E+04 1.09117E+05 4.08635E+05
3.05406E+07 6.94852E+03 8.58466E+03 1.69754E+05 1.07843E+06
3.09623E+05 2.41213E+05 3.08149E+05 2.22054E+05 1.94011E+05
4.42225E+05 3.03071E+05 1.18520E+06 3.45334E+03 1.10997E+03
C ------------------------------------------------------------------------
si105 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp105 0.00000E+00 1.85638E+04 3.67221E+04 1.16674E+05 4.28978E+05
3.21759E+07 7.41406E+03 9.14695E+03 1.80954E+05 1.12799E+06
3.24748E+05 2.55490E+05 3.24699E+05 2.34011E+05 2.07133E+05
4.70575E+05 3.23871E+05 1.26777E+06 3.63481E+03 1.17512E+03
C ------------------------------------------------------------------------
si106 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp106 0.00000E+00 2.21857E+04 4.43306E+04 1.40802E+05 5.10575E+05
3.77016E+07 8.93175E+03 1.09930E+04 2.17852E+05 1.33861E+06
3.85960E+05 3.05846E+05 3.86785E+05 2.78627E+05 2.49620E+05
5.65666E+05 3.90547E+05 1.53035E+06 4.33437E+03 1.40296E+03
C ------------------------------------------------------------------------
si107 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp107 0.00000E+00 3.02306E+04 6.08550E+04 1.93239E+05 6.93589E+05
4.94114E+07 1.22401E+04 1.50130E+04 2.98419E+05 1.81412E+06
5.23227E+05 4.16740E+05 5.24421E+05 3.77374E+05 3.42124E+05
7.73787E+05 5.35477E+05 2.10065E+06 5.89540E+03 1.90221E+03
C ------------------------------------------------------------------------
si108 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp108 0.00000E+00 2.94446E+04 5.92364E+04 1.88103E+05 6.75727E+05
4.79016E+07 1.19154E+04 1.46108E+04 2.90519E+05 1.76765E+06
5.09678E+05 4.05737E+05 5.10570E+05 3.67353E+05 3.33032E+05
7.53323E+05 5.21216E+05 2.04477E+06 5.74269E+03 1.85095E+03
C ------------------------------------------------------------------------
si109 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp109 0.00000E+00 4.35594E+04 8.73170E+04 2.77304E+05 1.00114E+06
7.16969E+07 1.75773E+04 2.15781E+04 4.28663E+05 2.62181E+06
7.55655E+05 6.00004E+05 7.56553E+05 5.44496E+05 4.91231E+05
1.11219E+06 7.68647E+05 3.01416E+06 8.50262E+03 2.74125E+03
C ------------------------------------------------------------------------
si110 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp110 0.00000E+00 5.08122E+04 1.01261E+05 3.21649E+05 1.17065E+06
8.44071E+07 2.04082E+04 2.50830E+04 4.97895E+05 3.07096E+06
8.84254E+05 6.99122E+05 8.83955E+05 6.36310E+05 5.70220E+05
1.29292E+06 8.91913E+05 3.49561E+06 9.93105E+03 3.19801E+03
C ------------------------------------------------------------------------
si111 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp111 0.00000E+00 5.54824E+04 1.09734E+05 3.48652E+05 1.28218E+06
9.19204E+07 2.21468E+04 2.72337E+04 5.40609E+05 3.37053E+06
9.68856E+05 7.61698E+05 9.65731E+05 6.95046E+05 6.18581E+05
1.40515E+06 9.67039E+05 3.78823E+06 1.08606E+04 3.48353E+03
C ------------------------------------------------------------------------
si112 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp112 0.00000E+00 5.97040E+04 1.16814E+05 3.71278E+05 1.38574E+06
9.88383E+07 2.36237E+04 2.90770E+04 5.77107E+05 3.65346E+06
1.04777E+06 8.17221E+05 1.04035E+06 7.48620E+05 6.59504E+05
1.50209E+06 1.03023E+06 4.03280E+06 1.17128E+04 3.73774E+03
C ------------------------------------------------------------------------
si113 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp113 0.00000E+00 4.85803E+04 9.36078E+04 2.97674E+05 1.13437E+06
7.99568E+07 1.89850E+04 2.33908E+04 4.64308E+05 3.00271E+06
8.58300E+05 6.62033E+05 8.47386E+05 6.09533E+05 5.29619E+05
1.21081E+06 8.26424E+05 3.23184E+06 9.55949E+03 3.02649E+03
C ------------------------------------------------------------------------
si114 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp114 0.00000E+00 5.88344E+04 1.11186E+05 3.53808E+05 1.38411E+06
9.77202E+07 2.26369E+04 2.79595E+04 5.54391E+05 3.68221E+06
1.04882E+06 7.98041E+05 1.02934E+06 7.40419E+05 6.30958E+05
1.44953E+06 9.83220E+05 3.83908E+06 1.16223E+04 3.65352E+03
C ------------------------------------------------------------------------
si115 1e-12 0.15 0.2 0.3 0.4
0.5 0.6 0.8 1.0 1.5
2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 11.0
sp115 0.00000E+00 5.61276E+04 1.03236E+05 3.28823E+05 1.33382E+06
9.44561E+07 2.11336E+04 2.61958E+04 5.18590E+05 3.57223E+06
1.01277E+06 7.56511E+05 9.86081E+05 7.09333E+05 5.88350E+05
1.36095E+06 9.15062E+05 3.56505E+06 1.11463E+04 3.47063E+03
C ------------------------------------------------------------------------
c x-distribution
si31 2539.25
2587.35 2635.45 2683.55 2731.65 2779.75
2827.85 2875.95 2924.05 2972.15 3020.25
3068.35 3116.45 3164.55 3212.65 3260.75
sp31 0
4.57688E+07 2.25345E+07 2.38501E+07 3.55841E+07 3.75155E+07
4.40951E+07 5.81307E+07 5.63921E+07 8.42501E+07 9.90333E+07
1.07862E+08 1.15949E+08 9.38290E+07 1.14450E+08 1.10325E+08
c y-distribution
si32 -8 -2
sp32 0 1
c z-distribution
si33 -8 -2
sp33 0 1
sources/proton.i
+1 -1
View File
@@ -5,7 +5,7 @@ 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 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=8.7374e15
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
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8
sp1 D 0 0.12109 0.23438 0.33984 0.4375 0.52734 0.60938 0.68359 0.75
tallies/bunker_mesh.i
+2 -2
View File
@@ -2,8 +2,8 @@ c large mesh tally covering the beamline from monolith to S1/2 separator
tmesh
rmesh1:n dose 1 trans 1
cora1 500 10i 1000 20i 1500 16i 2300
corb1 -305 5i -55 4i 45 5i 295
corc1 -200 1i -100 5i 100 2i 300
corb1 -305 5i -55 4i 45 9i 495
corc1 -200 1i -100 5i 100 5i 600
c Official ESS dose conversion factors from n/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf1 0.000000000 1.00000e-08
tallies/bunker_mesh_np.i
+80
View File
@@ -0,0 +1,80 @@
c large mesh tally covering the beamline from monolith to S1/2 separator
tmesh
rmesh1:n dose 1 trans 1
cora1 500 10i 1000 20i 1500 16i 2300
corb1 -305 5i -55 4i 45 9i 495
corc1 -200 1i -100 5i 100 5i 600
c Official ESS dose conversion factors from n/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf1 0.000000000 1.00000e-08
1.00000e-09 1.11000e-08 1.00000e-08 1.28000e-08
2.50000e-08 1.44000e-08 1.00000e-07 1.87000e-08
2.00000e-07 2.11000e-08 5.00000e-07 2.37000e-08
1.00000e-06 2.53000e-08 2.00000e-06 2.66000e-08
5.00000e-06 2.78000e-08 1.00000e-05 2.82000e-08
2.00000e-05 2.82000e-08 5.00000e-05 2.82000e-08
1.00000e-04 2.80000e-08 2.00000e-04 2.78000e-08
5.00000e-04 2.71000e-08 1.00000e-03 2.71000e-08
2.00000e-03 2.74000e-08 5.00000e-03 2.87000e-08
1.00000e-02 3.28000e-08 2.00000e-02 4.39000e-08
3.00000e-02 5.65000e-08 5.00000e-02 8.28000e-08
7.00000e-02 1.10000e-07 1.00000e-01 1.51000e-07
1.50000e-01 2.18000e-07 2.00000e-01 2.84000e-07
3.00000e-01 4.10000e-07 5.00000e-01 6.37000e-07
7.00000e-01 8.35000e-07 9.00000e-01 1.00000e-06
1.00000e+00 1.08000e-06 1.20000e+00 1.19000e-06
1.50000e+00 1.31000e-06 2.00000e+00 1.47000e-06
3.00000e+00 1.65000e-06 4.00000e+00 1.74000e-06
5.00000e+00 1.78000e-06 6.00000e+00 1.79000e-06
7.00000e+00 1.80000e-06 8.00000e+00 1.80000e-06
9.00000e+00 1.80000e-06 1.00000e+01 1.80000e-06
1.20000e+01 1.80000e-06 1.40000e+01 1.78000e-06
1.50000e+01 1.77000e-06 1.60000e+01 1.76000e-06
1.80000e+01 1.74000e-06 2.00000e+01 1.72000e-06
2.10000e+01 1.71000e-06 3.00000e+01 1.63000e-06
5.00000e+01 1.56000e-06 7.50000e+01 1.58000e-06
1.00000e+02 1.60000e-06 1.30000e+02 1.61000e-06
1.50000e+02 1.61000e-06 1.80000e+02 1.61000e-06
2.00000e+02 1.61000e-06 3.00000e+02 1.67000e-06
4.00000e+02 1.79000e-06 5.00000e+02 1.92000e-06
6.00000e+02 2.05000e-06 7.00000e+02 2.16000e-06
8.00000e+02 2.24000e-06 9.00000e+02 2.30000e-06
1.00000e+03 2.35000e-06 2.00000e+03 2.76000e-06
5.00000e+03 3.64000e-06 1.00000e+04 4.75000e-06
c
rmesh11:p dose 2 trans 1
cora11 500 10i 1000 20i 1500 16i 2300
corb11 -305 5i -55 4i 45 9i 495
corc11 -200 1i -100 5i 100 5i 600
c Official ESS dose conversion factors from gamma/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf2 0.000000 1.00000e-04
1.00E-02 2.47E-10 1.50E-02 5.62E-10
2.00E-02 8.1E-10 3.00E-02 1.13E-09
4.00E-02 1.26E-09 5.00E-02 1.33E-09
6.00E-02 1.4E-09 7.00E-02 1.49E-09
8.00E-02 1.6E-09 1.00E-01 1.87E-09
1.50E-01 2.69E-09 2.00E-01 3.6E-09
3.00E-01 5.44E-09 4.00E-01 7.2E-09
5.00E-01 8.89E-09 5.11E-01 9.07E-09
6.00E-01 1.05E-08 6.62E-01 1.14E-08
8.00E-01 1.34E-08 1.00E+00 1.62E-08
1.12E+00 1.76E-08 1.33E+00 2.01E-08
1.50E+00 2.2E-08 2.00E+00 2.69E-08
3.00E+00 3.51E-08 4.00E+00 4.21E-08
5.00E+00 4.82E-08 6.00E+00 5.4E-08
6.13E+00 5.47E-08 8.00E+00 6.7E-08
1.00E+01 7.92E-08 1.50E+01 1.09E-07
2.00E+01 1.38E-07 3.00E+01 1.85E-07
4.00E+01 2.23E-07 5.00E+01 2.6E-07
6.00E+01 2.95E-07 8.00E+01 3.52E-07
1.00E+02 3.96E-07 1.50E+02 4.68E-07
2.00E+02 5.15E-07 3.00E+02 5.8E-07
4.00E+02 6.19E-07 5.00E+02 6.48E-07
6.00E+02 6.7E-07 8.00E+02 7.02E-07
1.00E+03 7.24E-07 1.50E+03 7.63E-07
2.00E+03 7.92E-07 3.00E+03 8.35E-07
4.00E+03 8.75E-07 5.00E+03 9.04E-07
6.00E+03 9.29E-07 8.00E+03 9.65E-07
1.00E+04 9.94E-07
endmd
tallies/bunker_mesh_p.i
+2 -2
View File
@@ -2,8 +2,8 @@ c large mesh tally covering the beamline from monolith to S1/2 separator
tmesh
rmesh1:p dose 1 trans 1
cora1 500 10i 1000 20i 1500 16i 2300
corb1 -305 5i -55 4i 45 5i 295
corc1 -200 1i -100 5i 100 2i 300
corb1 -305 5i -55 4i 45 9i 495
corc1 -200 1i -100 5i 100 5i 600
c Official ESS dose conversion factors from gamma/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf1 0.000000 1.00000e-04
tallies/full_mesh.i
+44
View File
@@ -0,0 +1,44 @@
c large and detailed mesh tally covering the beamline from monolith to full cave
tmesh
rmesh1:n dose 1 trans 1
cora1 500 10i 1000 20i 1500 60i 4500
corb1 -705 13i -55 20i 55 13i 705
corc1 -200 2i -50 30i 100 2i 300
c Official ESS dose conversion factors from n/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf1 0.000000000 1.00000e-08
1.00000e-09 1.11000e-08 1.00000e-08 1.28000e-08
2.50000e-08 1.44000e-08 1.00000e-07 1.87000e-08
2.00000e-07 2.11000e-08 5.00000e-07 2.37000e-08
1.00000e-06 2.53000e-08 2.00000e-06 2.66000e-08
5.00000e-06 2.78000e-08 1.00000e-05 2.82000e-08
2.00000e-05 2.82000e-08 5.00000e-05 2.82000e-08
1.00000e-04 2.80000e-08 2.00000e-04 2.78000e-08
5.00000e-04 2.71000e-08 1.00000e-03 2.71000e-08
2.00000e-03 2.74000e-08 5.00000e-03 2.87000e-08
1.00000e-02 3.28000e-08 2.00000e-02 4.39000e-08
3.00000e-02 5.65000e-08 5.00000e-02 8.28000e-08
7.00000e-02 1.10000e-07 1.00000e-01 1.51000e-07
1.50000e-01 2.18000e-07 2.00000e-01 2.84000e-07
3.00000e-01 4.10000e-07 5.00000e-01 6.37000e-07
7.00000e-01 8.35000e-07 9.00000e-01 1.00000e-06
1.00000e+00 1.08000e-06 1.20000e+00 1.19000e-06
1.50000e+00 1.31000e-06 2.00000e+00 1.47000e-06
3.00000e+00 1.65000e-06 4.00000e+00 1.74000e-06
5.00000e+00 1.78000e-06 6.00000e+00 1.79000e-06
7.00000e+00 1.80000e-06 8.00000e+00 1.80000e-06
9.00000e+00 1.80000e-06 1.00000e+01 1.80000e-06
1.20000e+01 1.80000e-06 1.40000e+01 1.78000e-06
1.50000e+01 1.77000e-06 1.60000e+01 1.76000e-06
1.80000e+01 1.74000e-06 2.00000e+01 1.72000e-06
2.10000e+01 1.71000e-06 3.00000e+01 1.63000e-06
5.00000e+01 1.56000e-06 7.50000e+01 1.58000e-06
1.00000e+02 1.60000e-06 1.30000e+02 1.61000e-06
1.50000e+02 1.61000e-06 1.80000e+02 1.61000e-06
2.00000e+02 1.61000e-06 3.00000e+02 1.67000e-06
4.00000e+02 1.79000e-06 5.00000e+02 1.92000e-06
6.00000e+02 2.05000e-06 7.00000e+02 2.16000e-06
8.00000e+02 2.24000e-06 9.00000e+02 2.30000e-06
1.00000e+03 2.35000e-06 2.00000e+03 2.76000e-06
5.00000e+03 3.64000e-06 1.00000e+04 4.75000e-06
endmd
tallies/full_mesh_np.i
+80
View File
@@ -0,0 +1,80 @@
c large and detailed mesh tally covering the beamline from monolith to full cave
tmesh
rmesh1:n dose 1 trans 1
cora1 500 10i 1000 20i 1500 60i 4500
corb1 -705 13i -55 20i 55 13i 705
corc1 -200 2i -50 30i 100 2i 300
c Official ESS dose conversion factors from n/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf1 0.000000000 1.00000e-08
1.00000e-09 1.11000e-08 1.00000e-08 1.28000e-08
2.50000e-08 1.44000e-08 1.00000e-07 1.87000e-08
2.00000e-07 2.11000e-08 5.00000e-07 2.37000e-08
1.00000e-06 2.53000e-08 2.00000e-06 2.66000e-08
5.00000e-06 2.78000e-08 1.00000e-05 2.82000e-08
2.00000e-05 2.82000e-08 5.00000e-05 2.82000e-08
1.00000e-04 2.80000e-08 2.00000e-04 2.78000e-08
5.00000e-04 2.71000e-08 1.00000e-03 2.71000e-08
2.00000e-03 2.74000e-08 5.00000e-03 2.87000e-08
1.00000e-02 3.28000e-08 2.00000e-02 4.39000e-08
3.00000e-02 5.65000e-08 5.00000e-02 8.28000e-08
7.00000e-02 1.10000e-07 1.00000e-01 1.51000e-07
1.50000e-01 2.18000e-07 2.00000e-01 2.84000e-07
3.00000e-01 4.10000e-07 5.00000e-01 6.37000e-07
7.00000e-01 8.35000e-07 9.00000e-01 1.00000e-06
1.00000e+00 1.08000e-06 1.20000e+00 1.19000e-06
1.50000e+00 1.31000e-06 2.00000e+00 1.47000e-06
3.00000e+00 1.65000e-06 4.00000e+00 1.74000e-06
5.00000e+00 1.78000e-06 6.00000e+00 1.79000e-06
7.00000e+00 1.80000e-06 8.00000e+00 1.80000e-06
9.00000e+00 1.80000e-06 1.00000e+01 1.80000e-06
1.20000e+01 1.80000e-06 1.40000e+01 1.78000e-06
1.50000e+01 1.77000e-06 1.60000e+01 1.76000e-06
1.80000e+01 1.74000e-06 2.00000e+01 1.72000e-06
2.10000e+01 1.71000e-06 3.00000e+01 1.63000e-06
5.00000e+01 1.56000e-06 7.50000e+01 1.58000e-06
1.00000e+02 1.60000e-06 1.30000e+02 1.61000e-06
1.50000e+02 1.61000e-06 1.80000e+02 1.61000e-06
2.00000e+02 1.61000e-06 3.00000e+02 1.67000e-06
4.00000e+02 1.79000e-06 5.00000e+02 1.92000e-06
6.00000e+02 2.05000e-06 7.00000e+02 2.16000e-06
8.00000e+02 2.24000e-06 9.00000e+02 2.30000e-06
1.00000e+03 2.35000e-06 2.00000e+03 2.76000e-06
5.00000e+03 3.64000e-06 1.00000e+04 4.75000e-06
c
rmesh11:p dose 2 trans 1
cora11 500 10i 1000 20i 1500 60i 4500
corb11 -705 13i -55 20i 55 13i 705
corc11 -200 2i -50 30i 100 2i 300
c Official ESS dose conversion factors from gamma/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf2 0.000000 1.00000e-04
1.00E-02 2.47E-10 1.50E-02 5.62E-10
2.00E-02 8.1E-10 3.00E-02 1.13E-09
4.00E-02 1.26E-09 5.00E-02 1.33E-09
6.00E-02 1.4E-09 7.00E-02 1.49E-09
8.00E-02 1.6E-09 1.00E-01 1.87E-09
1.50E-01 2.69E-09 2.00E-01 3.6E-09
3.00E-01 5.44E-09 4.00E-01 7.2E-09
5.00E-01 8.89E-09 5.11E-01 9.07E-09
6.00E-01 1.05E-08 6.62E-01 1.14E-08
8.00E-01 1.34E-08 1.00E+00 1.62E-08
1.12E+00 1.76E-08 1.33E+00 2.01E-08
1.50E+00 2.2E-08 2.00E+00 2.69E-08
3.00E+00 3.51E-08 4.00E+00 4.21E-08
5.00E+00 4.82E-08 6.00E+00 5.4E-08
6.13E+00 5.47E-08 8.00E+00 6.7E-08
1.00E+01 7.92E-08 1.50E+01 1.09E-07
2.00E+01 1.38E-07 3.00E+01 1.85E-07
4.00E+01 2.23E-07 5.00E+01 2.6E-07
6.00E+01 2.95E-07 8.00E+01 3.52E-07
1.00E+02 3.96E-07 1.50E+02 4.68E-07
2.00E+02 5.15E-07 3.00E+02 5.8E-07
4.00E+02 6.19E-07 5.00E+02 6.48E-07
6.00E+02 6.7E-07 8.00E+02 7.02E-07
1.00E+03 7.24E-07 1.50E+03 7.63E-07
2.00E+03 7.92E-07 3.00E+03 8.35E-07
4.00E+03 8.75E-07 5.00E+03 9.04E-07
6.00E+03 9.29E-07 8.00E+03 9.65E-07
1.00E+04 9.94E-07
endmd
tallies/full_mesh_p.i
+38
View File
@@ -0,0 +1,38 @@
c large and detailed mesh tally covering the beamline from monolith to full cave
tmesh
rmesh1:p dose 1 trans 1
cora1 500 10i 1000 20i 1500 60i 4500
corb1 -705 13i -55 20i 55 13i 705
corc1 -200 2i -50 30i 100 2i 300
c Official ESS dose conversion factors from gamma/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf1 0.000000 1.00000e-04
1.00E-02 2.47E-10 1.50E-02 5.62E-10
2.00E-02 8.1E-10 3.00E-02 1.13E-09
4.00E-02 1.26E-09 5.00E-02 1.33E-09
6.00E-02 1.4E-09 7.00E-02 1.49E-09
8.00E-02 1.6E-09 1.00E-01 1.87E-09
1.50E-01 2.69E-09 2.00E-01 3.6E-09
3.00E-01 5.44E-09 4.00E-01 7.2E-09
5.00E-01 8.89E-09 5.11E-01 9.07E-09
6.00E-01 1.05E-08 6.62E-01 1.14E-08
8.00E-01 1.34E-08 1.00E+00 1.62E-08
1.12E+00 1.76E-08 1.33E+00 2.01E-08
1.50E+00 2.2E-08 2.00E+00 2.69E-08
3.00E+00 3.51E-08 4.00E+00 4.21E-08
5.00E+00 4.82E-08 6.00E+00 5.4E-08
6.13E+00 5.47E-08 8.00E+00 6.7E-08
1.00E+01 7.92E-08 1.50E+01 1.09E-07
2.00E+01 1.38E-07 3.00E+01 1.85E-07
4.00E+01 2.23E-07 5.00E+01 2.6E-07
6.00E+01 2.95E-07 8.00E+01 3.52E-07
1.00E+02 3.96E-07 1.50E+02 4.68E-07
2.00E+02 5.15E-07 3.00E+02 5.8E-07
4.00E+02 6.19E-07 5.00E+02 6.48E-07
6.00E+02 6.7E-07 8.00E+02 7.02E-07
1.00E+03 7.24E-07 1.50E+03 7.63E-07
2.00E+03 7.92E-07 3.00E+03 8.35E-07
4.00E+03 8.75E-07 5.00E+03 9.04E-07
6.00E+03 9.29E-07 8.00E+03 9.65E-07
1.00E+04 9.94E-07
endmd