Merge branch 'master' into mcstas_sources
This commit is contained in:
+2
-1
@@ -7,4 +7,5 @@ out*
|
||||
slurm-*.out
|
||||
tmp
|
||||
runs
|
||||
|
||||
*.pyc
|
||||
*.pyo
|
||||
|
||||
+1
-1
@@ -2,7 +2,7 @@
|
||||
#SBATCH -J mcnpEstia
|
||||
#SBATCH -N 2
|
||||
#SBATCH --ntasks-per-node=24
|
||||
#SBATCH --time=48:00:00
|
||||
#SBATCH --time=7-00:00:00
|
||||
#SBATCH --mail-type=fail
|
||||
#SBATCH --mail-user=artur.glavic@psi.ch
|
||||
#SBATCH -o tmp/stdout.log
|
||||
|
||||
+8
-10
@@ -101,16 +101,14 @@ c read file=tallies/e1_e2_f4.i echo
|
||||
c read file=tallies/e1_e2_f5.i echo
|
||||
c read file=tallies/shutter_f4.i echo
|
||||
c read file=tallies/shutter_f5.i echo
|
||||
c read file=tallies/selene1_walls_f4.i echo
|
||||
c read file=tallies/optical_cave_f5.i echo
|
||||
c read file=tallies/experimental_cave_f5.i echo
|
||||
read file=tallies/cave_mesh.i echo
|
||||
c read file=tallies/cave_mesh.i echo
|
||||
c
|
||||
c ==========================================
|
||||
c
|
||||
c ==========================================
|
||||
c dxt:n 730 -1380 0 31 120 1020 -1170 0 31 120
|
||||
c 779.4 -1093.7 -2.6 13 80 0.1 1e-15
|
||||
dxc0:n 0 508r 0.25 64r 1 116r
|
||||
c
|
||||
c 1 2 3 4 5 6-14 15-22 23 24 25 26-34 35 36 37 38 39
|
||||
c 40-42 43-266 267 268-488 489-494 //
|
||||
@@ -122,7 +120,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 27r
|
||||
0 0.1V3 0 0.1V3 0 0.1V3 5r 0 30r
|
||||
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
|
||||
@@ -131,12 +129,13 @@ c -------------------------------------------------------
|
||||
c --------------- PHYSICS CARDS --------------------------
|
||||
c -------------------------------------------------------
|
||||
c dbcn 375642321 j j j j j j j j j j j 15291711
|
||||
nps 3e7
|
||||
nps 2e9
|
||||
c histp
|
||||
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 46r
|
||||
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
|
||||
cut:h,/,d,t,s,a 1e+08 1e-3
|
||||
@@ -152,11 +151,10 @@ lea 1 4 1 0 1 0 0 1
|
||||
c * Source for simulation from proton beam (needs different WW and ESplit)
|
||||
c read file=sources/proton.i echo
|
||||
c * Source Term from Valentina for 2m adapted to size
|
||||
c read file=sources/e2_2m.i echo
|
||||
read file=sources/e2_2m.i echo
|
||||
c * McStas source terms for gamma doeses
|
||||
c read file=sources/mcstas_source_vs.i echo
|
||||
read file=sources/mcstas_source_sample.i echo
|
||||
prdmp 5e6 -60 1 2 1e6
|
||||
prdmp 5e7 -60 1 2 5e6
|
||||
print 30 32 $ 120 128 150 160 161 162
|
||||
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
||||
|
||||
|
||||
+6
-66
@@ -79,35 +79,19 @@ c Main ESS plant geometry without the instrument
|
||||
( -330 : -326 : -333 : 334 : 325 : 64 : 3 )
|
||||
( -3 : -330 : -325 : -333 : 334 : 326 : 64 )
|
||||
( -565 : 349 : 566 : -1001 : 332 )
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
c (-1100 : 1101 : -1201 : -1202 : -1203 : -1204)
|
||||
c (-1101 : 1105 : -1211 : 1504 : 1503)
|
||||
c #590 #591
|
||||
43 4 -1.2E-3 ( -65 : -119 : 63 : 101 : 124 ) 3 -102 65 2 66 -67 64 #257
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
44 71 0.0843223 64 -67 66 69 -65 -106 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
45 71 0.0843223 64 70 -66 2 69 -106 3
|
||||
46 71 0.0843223 -67 66 -73 65 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
47 71 0.0843223 -67 66 -74 73 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
48 71 0.0843223 -67 66 -75 74 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
49 71 0.0843223 -67 66 -76 75 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
50 71 0.0843223 -67 66 -77 76 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
51 71 0.0843223 -67 66 -78 77 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
52 71 0.0843223 -67 66 -79 78 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
53 71 0.0843223 -67 66 -80 79 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
54 71 0.0843223 -67 66 -81 80 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
55 71 0.0843223 -67 66 2 81 -106 102 3
|
||||
dxc1:n=0.05 dxc2:n=0.05 dxc3:n=0.05 dxc4:n=0.05
|
||||
56 122 0.0824808 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -106 67 69
|
||||
-73 -82 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
57 122 0.0824808 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -106 69 -73
|
||||
@@ -509,38 +493,27 @@ c #590 #591
|
||||
255 122 0.0824808 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -106 -71 81
|
||||
2 100 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
256 4 -1.2E-3 ( -65 : -119 : 63 : 101 : 124 ) 3 -68 -2 -101 66 -67 64
|
||||
-112 dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
-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
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
259 49 0.074981 ( -3 : -69 : -2 : -70 : -64 : 71 : 72 ) 3 -72 103
|
||||
-104 -66 70 64
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
260 200 -3.9 -67 66 -107 -2 -72 68 3
|
||||
dxc1:n=0.1 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1
|
||||
261 200 -3.9 -67 66 -108 107 -72 68 3
|
||||
dxc1:n=0.1 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1
|
||||
262 200 -3.9 -67 66 -109 108 -72 68 3
|
||||
dxc1:n=0.1 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1
|
||||
263 200 -3.9 -67 66 -110 109 -72 68 3
|
||||
dxc1:n=0.3 dxc2:n=0.3 dxc3:n=1 dxc4:n=1
|
||||
264 200 -3.9 -67 66 -111 110 -72 68 3
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
265 200 -3.9 -67 66 -112 111 -72 68 3
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
c PSI Schwerbeton mit 5% B4C
|
||||
266 200 -3.9 -67 66 -113 112 -72 68 3
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
267 200 -3.9 -67 66 -114 113 -72 68 3 #(
|
||||
(68 -6100 -6201 -6202 -6203 -6204):(
|
||||
6100 -6101 -6301 -6302 -6303 -6304):(
|
||||
6101 -6102 -6401 -6402 -6403 -6404):(
|
||||
6102 -6103 -6501 -6502 -6503 -6504):(
|
||||
6103 -72 -6601 -6602 -6603 -6604 -6605))
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
268 200 -3.9 -67 66 -101 114 -72 68 3
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
269 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) ( -3 : -69 :
|
||||
-2 : -70 : -64 : 71 : 72 ) 64 -115 -107 103 67 -72 3 ( -65 :
|
||||
-119 : 63 : 101 : 124 )
|
||||
@@ -576,13 +549,10 @@ c PSI Schwerbeton mit 5% B4C
|
||||
-109 117 -118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
283 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 -105
|
||||
108 -109 118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
284 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 67
|
||||
109 -110 -115 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
285 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 109
|
||||
-110 115 -116 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
286 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 109
|
||||
-110 116 -117 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
287 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 109
|
||||
@@ -591,10 +561,8 @@ c PSI Schwerbeton mit 5% B4C
|
||||
109 -110 118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
289 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 67
|
||||
110 -111 -115 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
290 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 110
|
||||
-111 115 -116 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
291 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 110
|
||||
-111 116 -117 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
292 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 110
|
||||
@@ -603,61 +571,44 @@ c PSI Schwerbeton mit 5% B4C
|
||||
110 -111 118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
294 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 67
|
||||
111 -112 -115 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
295 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 111
|
||||
-112 115 -116 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
296 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 111
|
||||
-112 116 -117 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
297 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 111
|
||||
-112 117 -118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
298 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 -105
|
||||
111 -112 118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
299 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 67
|
||||
112 -113 -115 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
300 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 112
|
||||
-113 115 -116 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
301 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 112
|
||||
-113 116 -117 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
302 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 112
|
||||
-113 117 -118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
303 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 -105
|
||||
112 -113 118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
304 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 67
|
||||
113 -114 -115 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
305 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 113
|
||||
-114 115 -116 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
306 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 113
|
||||
-114 116 -117 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
6 dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
-114 116 -117 64 ( -65 : -119 : 63 : 101 : 124 ) 6
|
||||
307 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 113
|
||||
-114 117 -118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
308 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 -105
|
||||
113 -114 118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
309 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 67
|
||||
114 -104 -115 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
310 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 114
|
||||
-104 115 -116 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
6 dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
-104 115 -116 64 ( -65 : -119 : 63 : 101 : 124 ) 6
|
||||
311 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 114
|
||||
-104 116 -117 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
312 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 114
|
||||
-104 117 -118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
313 200 -3.9 ( -121 : -65 : -122 : 120 : 101 : 123 ) 3 -72 -105
|
||||
114 -104 118 64 ( -65 : -119 : 63 : 101 : 124 )
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
314 49 0.074981 63 -122 -101 65 -124 121
|
||||
315 49 0.074981 121 -123 65 -101 -125 122
|
||||
316 49 0.074981 121 -123 65 -101 -126 125
|
||||
@@ -902,26 +853,22 @@ c PSI Schwerbeton mit 5% B4C
|
||||
487 5 0.0582256 290 -320 33 -21 317 -319
|
||||
488 5 0.0582256 290 -320 33 -21 318 -316
|
||||
489 54 0.0833854 ( -153 : 154 : 159 ) -327 -142 -326 325 -330 61 3
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
490 54 0.0833854 -109 -334 325 -326 -64 330 3
|
||||
( -153 : 154 : 159 ) -329
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
491 54 0.0833854 ( -153 : 154 : 159 ) -327 -142 326 -325 -330 61 -3
|
||||
492 54 0.0833854 ( -153 : 154 : 159 ) -3 330 -64 -325 326 -329
|
||||
-334
|
||||
493 54 0.0833854 ( -153 : 154 : 159 ) -332 331 -325 326 -330 61 -3
|
||||
494 54 0.0833854 556 ( -153 : 154 : 159 )
|
||||
3 61 -330 -326 325 335 -332
|
||||
dxc1:n=0.01 dxc2:n=0.01 dxc3:n=0.01 dxc4:n=0.01
|
||||
c port
|
||||
585 0 ( 352 : -351 : 345 : 458 ) ( -554 : -575 : 457 ) -330 -457 575
|
||||
349 -350 61 348 dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
349 -350 61 348
|
||||
586 3 0.0876276 ( -554 : -575 : 457 ) -1103 #589
|
||||
-458 -345 351 -352 61 348
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
587 0 ( 356 : -355 : -461 : 462 ) -64 -460 459 353 -354 330 348
|
||||
588 3 0.0876276 -64 -332 459 355 -356 330 -1107
|
||||
348 #589 dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
348 #589
|
||||
c open beam 3 SKADI 30 x 30 mm2
|
||||
589 0 61 -64 -1597 1596 -1595 1594 1593
|
||||
590 0 ( 568 : -567 : -1000 : 458 ) ( -349 : -1001 : 457 ) -330 -457 1001
|
||||
@@ -930,23 +877,16 @@ c open beam 3 SKADI 30 x 30 mm2
|
||||
-458 1000 567 -568 61 -564 -2 -1103
|
||||
( 1101 : -1305 : -1306 : -1307 : -1308)
|
||||
(-1101 : -1213 : 1507 : 1508)
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
592 0 ( 572 : -571 : -461 : 462 ) -64 -460 459 569 -570 330 -564
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
593 3 0.0876276 -64 -332 459 571 -572 330 -564 -2
|
||||
( 1115 : -1213 : 1507 : 1508) -1107
|
||||
(-1115 : -1214 : 1510 : 1509)
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
c neutron guide
|
||||
595 0 ( 340 : -339 : 345 : 458 ) ( -565 : -575 : 457 ) -330 -457 575
|
||||
337 -338 61 336
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
596 3 0.0876276 ( -565 : -575 : 457 ) -1103
|
||||
-458 -345 339 -340 61 336
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
597 0 ( 344 : -343 : -461 : 462 ) -64 -460 459 341 -342 330 336
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
598 3 0.0876276 -64 -332 459 343 -344 330 336 -1107
|
||||
dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
c closed beam 1
|
||||
600 0 -600 dxc1:n=0.25 dxc2:n=0.25 dxc3:n=0.25 dxc4:n=0.25
|
||||
600 0 -600
|
||||
|
||||
@@ -53,7 +53,7 @@ c concrete floor
|
||||
c ceiling outside
|
||||
8002 49 0.074981 ( 72 -8101 8207 8203 8201 8202)
|
||||
c left outer wall (E01)
|
||||
8003 49 0.074981 ( 72 -8100 8201 -8205 59 -8207 3)
|
||||
8003 4 -1.2E-3 ( 72 -8100 8201 -8205 59 -8207 3)
|
||||
c Selene 1/2 separator
|
||||
8004 200 -3.9 (8100 -8101 8201 8202 59 -8207 3)
|
||||
#(8100 -8101 -8221 -8222 -8223 -8224)
|
||||
@@ -69,12 +69,13 @@ c air
|
||||
#(-8121 -8931)
|
||||
#707 #710
|
||||
c B4C concrete inner layers
|
||||
8011 16 -2.4 (8211 8215 -8213 -8214 72 -8100 3)
|
||||
8011 4 -1.2E-3 (8211 8215 -8213 -8214 72 -8100 3)
|
||||
8012 16 -2.4 (8205 8202 8213 -8218 72 -8100 3)
|
||||
8013 16 -2.4 (8212 -8216 -8213 -8214 8111 -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 (8212 8202 -8213 -8214 72 -8111 3)
|
||||
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)
|
||||
|
||||
@@ -51,13 +51,10 @@ c extended cut-out in reflector
|
||||
c split of wedge of monolith components for increased importance
|
||||
1236 5 0.0582256 ( 142 -332 -556 325 ) ( -1201 : -1202 : -1203 : -1204 : -58)
|
||||
( -39 : 40 : -41 : 42 ) -61 3 21 -46 45 -44 43
|
||||
dxc1:n=0.001
|
||||
1240 121 0.0862739 ( 142 -332 -556 325 )
|
||||
( -43 : -45 : 46 : 44 ) -58 ( -53 : 54 : 55 )
|
||||
dxc1:n=0.001 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1
|
||||
1241 121 0.0862739 ( 142 -332 -556 325 ) ( -1201 : -1202 : -1203 : -1204 )
|
||||
-61 58 ( -43 : -45 : 46 : 44 )
|
||||
dxc1:n=0.1 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1
|
||||
c partitioned air in bunker
|
||||
1255 4 -1.2E-3 3 -68 -2 -101 66 -67 64 112
|
||||
( 329 : 4615 : -342 : 353 : 4100 )
|
||||
|
||||
@@ -505,8 +505,8 @@ c Wedge
|
||||
-1080.498 -3.724 769.345 -1083.823 -3.525
|
||||
c ----------- Ouside Bunker --------------
|
||||
c Wall radii
|
||||
8100 cz 2102.01
|
||||
8101 cz 2187.17
|
||||
8100 cz 2112.01
|
||||
8101 cz 2157.17
|
||||
c 8102 cz 3023.41
|
||||
8103 cz 1355.64
|
||||
8104 cz 2075.58
|
||||
@@ -522,6 +522,7 @@ 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
|
||||
8204 pz -200.000
|
||||
c steel/concrete interface
|
||||
|
||||
@@ -0,0 +1,47 @@
|
||||
#!/bin/bash
|
||||
# script to submit interactive job to mcc05 or mcc06 with 64 cores
|
||||
[[ -z "$1" ]] && { echo "You have to supply mcc05 or mcc06" ; exit 1; }
|
||||
# make sure the tag for the generated files is correct
|
||||
# by checking that all changes have been commited to GIT
|
||||
# (this does not include newly created files)
|
||||
if ( ! git diff-index --quiet --cached HEAD -- ) || ( ! git diff-files --quiet )
|
||||
then
|
||||
echo "Uncommitted changes in the GIT repository:"
|
||||
git status --untracked-files=no
|
||||
|
||||
read -p "Add changes and commit, amend last commit or exit [c/a/e]: " res
|
||||
case $res in
|
||||
[Cc]* ) git commit -a;;
|
||||
[Aa]* ) git commit -a --amend;;
|
||||
* ) exit 1;;
|
||||
esac
|
||||
fi
|
||||
|
||||
CURRENT_VERSION=`git describe --tags --long`
|
||||
NCORES=${2:-64}
|
||||
# version to continue, by default this or the previous version
|
||||
CONTINUE_VERSION=${3:-$CURRENT_VERSION}
|
||||
if [ ! -f "runs/e2-estia_$CONTINUE_VERSION.run" ]; then
|
||||
# get describe from previous version
|
||||
echo "No run for revision $CONTINUE_VERSION, try to use previous version"
|
||||
PREV_ID=`git log -n 2 --oneline --pretty=format:"%h" | tail -1`
|
||||
CONTINUE_VERSION=`git describe --tags --long $PREV_ID`
|
||||
fi
|
||||
|
||||
echo "Starting at `date`"
|
||||
echo "Running on hosts: $1"
|
||||
echo "Running on $NCORES processors."
|
||||
echo "Current working directory is `pwd`"
|
||||
echo "Estia GIT revision is $CURRENT_VERSION"
|
||||
echo "Continue run of revision $CONTINUE_VERSION"
|
||||
|
||||
mpirun -np $NCORES --host $1 mcnp6.mpi c \
|
||||
inp=crun.i outp=results/e2-estia_$CURRENT_VERSION.out \
|
||||
runtpe=runs/e2-estia_$CONTINUE_VERSION.run \
|
||||
mctal=results/e2-estia_$CURRENT_VERSION.mct \
|
||||
mdata=results/e2-estia_$CURRENT_VERSION.mdt \
|
||||
wwout=weight_windows/e2-estia_$CURRENT_VERSION.wwg \
|
||||
wwone=weight_windows/e2-estia_$CURRENT_VERSION.wwo \
|
||||
wwinp=weight_windows/current.wwg
|
||||
|
||||
echo "Program finished with exit code $? at: `date`"
|
||||
+137
@@ -0,0 +1,137 @@
|
||||
#!/usr/bin/env python
|
||||
#-*- coding: utf-8 -*-
|
||||
"""
|
||||
Generate nice graphs from mesh-tally 1 mctal files.
|
||||
|
||||
At the current stage only simple text processing is done that doesn't work with
|
||||
multiple tallies or energy dependent meshes, yet.
|
||||
"""
|
||||
|
||||
import gzip
|
||||
from numpy import *
|
||||
from pylab import *
|
||||
from matplotlib.colors import LogNorm
|
||||
|
||||
COLORMAP='tab20c'
|
||||
SCALING=1.56e16/8.7374e15
|
||||
|
||||
def read_tally(fname, use_tally=None):
|
||||
if fname.endswith('.gz'):
|
||||
txt=gzip.open(fname, 'r').read()
|
||||
else:
|
||||
txt=open(fname, 'r').read()
|
||||
|
||||
tally_ranges=[]
|
||||
tidx=0
|
||||
while 'tally' in txt[tidx+1:]:
|
||||
tidx=txt.index('tally ', tidx+1)
|
||||
tnr=int(txt[tidx:].split(None,3)[1])
|
||||
if len(tally_ranges)>0:
|
||||
tally_ranges[-1][2]=tidx-1
|
||||
tally_ranges.append([tnr, tidx, -1])
|
||||
if use_tally is None:
|
||||
txt=txt[tally_ranges[0][1]:tally_ranges[0][2]]
|
||||
else:
|
||||
tidx=[ti[0] for ti in tally_ranges].index(use_tally)
|
||||
txt=txt[tally_ranges[tidx][1]:tally_ranges[tidx][2]]
|
||||
header, data=txt.split('vals')
|
||||
|
||||
|
||||
# get dimensions of data
|
||||
fstart=header.index('\nf ')
|
||||
hi=header[fstart+3:].splitlines()[0]
|
||||
N,nE,nX,nY,nZ=map(int, hi.strip().split())
|
||||
|
||||
grid=header[fstart+3+len(hi):header.index('\nd ')]
|
||||
grid=array(grid.replace('\n', '').strip().split(), dtype=float)
|
||||
|
||||
x=grid[:nX+1]
|
||||
y=grid[nX+1:nX+nY+2]
|
||||
z=grid[nX+nY+2:nX+nY+nZ+3]
|
||||
|
||||
data=array(data.replace('\n','').strip().split(), dtype=float)
|
||||
I,dI=data.reshape(nZ, nY, nX, 2).transpose(3,2,1,0)
|
||||
|
||||
return x,y,z,I,dI
|
||||
|
||||
def plot_xy(outfile, res, zidx=5, Imin=1e-6, Imax=1e4, cticks=None):
|
||||
x,y,z,I,dI=res
|
||||
fig=gcf()
|
||||
ax=gca()
|
||||
p=ax.pcolormesh(x/100., y/100., I[:,:,zidx].T, norm=LogNorm(Imin, Imax), cmap=COLORMAP)
|
||||
ax.set_aspect('equal')
|
||||
fig.colorbar(p, orientation='horizontal', shrink=0.9, ticks=cticks, label='Dose rate [Sv/h]')
|
||||
xlabel('x [m]')
|
||||
ylabel('y [m]')
|
||||
title('XY plane at z=%.1fm'%((z[zidx]+z[zidx+1])/200.))
|
||||
fig.subplots_adjust()
|
||||
#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))
|
||||
if outfile is not None:
|
||||
fig.savefig(outfile, transparent=False, dpi=300)
|
||||
|
||||
def plot_xz(outfile, res, yidx=5, Imin=1e-6, Imax=1e4, cticks=None):
|
||||
x,y,z,I,dI=res
|
||||
fig=gcf()
|
||||
ax=gca()
|
||||
p=ax.pcolormesh(x/100., z/100., I[:,yidx].T, norm=LogNorm(Imin, Imax), cmap=COLORMAP)
|
||||
ax.set_aspect('equal')
|
||||
fig.colorbar(p, orientation='horizontal', shrink=0.9, ticks=cticks, label='Dose rate [Sv/h]')
|
||||
xlabel('x [m]')
|
||||
ylabel('z [m]')
|
||||
ypos=(y[yidx]+y[yidx+1])/200.
|
||||
title('XZ plane at y=%.1fm'%ypos)
|
||||
fig.subplots_adjust()
|
||||
#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)
|
||||
if outfile is not None:
|
||||
fig.savefig(outfile, transparent=False, dpi=300)
|
||||
|
||||
def plot_yz(outfile, res, xidx=5, Imin=1e-6, Imax=1e4, cticks=None):
|
||||
x,y,z,I,dI=res
|
||||
fig=gcf()
|
||||
ax=gca()
|
||||
p=ax.pcolormesh(y/100., z/100., I[xidx].T, norm=LogNorm(Imin, Imax), cmap=COLORMAP)
|
||||
ax.set_aspect('equal')
|
||||
fig.colorbar(p, orientation='horizontal', shrink=0.9, ticks=cticks, label='Dose rate [Sv/h]')
|
||||
xlabel('y [m]')
|
||||
ylabel('z [m]')
|
||||
title('YZ plane at x=%.1fm'%((x[xidx]+x[xidx+1])/200.))
|
||||
fig.subplots_adjust()
|
||||
#fig.savefig(outfile)
|
||||
#pc=p.get_facecolors()
|
||||
#pc[:,3]=where(I[xidx].T.flatten()>0, 0.9-0.25*dI[xidx].T.flatten(), 0.)
|
||||
if outfile is not None:
|
||||
fig.savefig(outfile, transparent=False, dpi=300)
|
||||
|
||||
if __name__=='__main__':
|
||||
import sys, os
|
||||
res=read_tally(sys.argv[1], use_tally=1)
|
||||
x,y,z,I,dI=res
|
||||
I*=SCALING
|
||||
fpre=os.path.join('.', 'images', os.path.basename(sys.argv[1])[:-4])
|
||||
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])
|
||||
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])
|
||||
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])
|
||||
|
||||
+1
-1
@@ -2,7 +2,7 @@
|
||||
#SBATCH -J mcnpEstia
|
||||
#SBATCH -N 2
|
||||
#SBATCH --ntasks-per-node=24
|
||||
#SBATCH --time=48:00:00
|
||||
#SBATCH --time=7-00:00:00
|
||||
#SBATCH --mail-type=fail
|
||||
#SBATCH --mail-user=artur.glavic@psi.ch
|
||||
#SBATCH -o tmp/stdout.log
|
||||
|
||||
+1
-1
@@ -1,7 +1,7 @@
|
||||
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,p j j 4 j -1 j 6.55305E14 $ WW was generated for other source scaling
|
||||
WWP:n j j 5 j -1 j 1.0
|
||||
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
|
||||
|
||||
@@ -0,0 +1,18 @@
|
||||
c f4 tally sphere at BBG
|
||||
fc14 neutron flux at BBG
|
||||
f14:n 701
|
||||
e14 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
|
||||
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
|
||||
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
|
||||
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
|
||||
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
|
||||
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
|
||||
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
|
||||
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
|
||||
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
|
||||
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
|
||||
1e-0 1.58 2.51 3.98 6.31
|
||||
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
|
||||
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
|
||||
1e+3 1.58e+3 2.51e+3
|
||||
fq14 e f
|
||||
@@ -0,0 +1,18 @@
|
||||
c f5 tally at BBG
|
||||
fc15 neutron flux at BBG
|
||||
f15:n 331.2 -448.0 13.3 0
|
||||
e15 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
|
||||
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
|
||||
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
|
||||
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
|
||||
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
|
||||
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
|
||||
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
|
||||
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
|
||||
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
|
||||
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
|
||||
1e-0 1.58 2.51 3.98 6.31
|
||||
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
|
||||
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
|
||||
1e+3 1.58e+3 2.51e+3
|
||||
fq15 e f
|
||||
@@ -0,0 +1,11 @@
|
||||
c f4 tally for all selene guide 1 side walls for WWG runs
|
||||
fc904 neutron flux on Selene 1 walls (not properly normalized)
|
||||
f904:n (8011 8013 8004 8401) 8011 8013 (8004 8401)
|
||||
e904 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
|
||||
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
|
||||
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
|
||||
1e-0 1.58 2.51 3.98 6.31
|
||||
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
|
||||
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
|
||||
1e+3 1.58e+3 2.51e+3
|
||||
fq904 e f
|
||||
@@ -1,6 +1,6 @@
|
||||
c f5 detector tally within the beam just outside the selene 1/2 separator wall
|
||||
fc55 neutron flux behind instrument shutter
|
||||
f55:n 1282 -1794 -13 0.0
|
||||
f55:n 1236.1 -1730.5 -12.1 0.0
|
||||
e55 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
|
||||
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
|
||||
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
|
||||
@@ -17,4 +17,5 @@ e55 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
|
||||
1e+3 1.58e+3 2.51e+3
|
||||
fq55 e f
|
||||
c 1-255 256 257 268 259-265 266-268 269-508 509-end
|
||||
pd55 0.001 254r 0.1 0.001 1 0.001 6r 1 2r 0.001 239r 0.1 110r 1 71r
|
||||
pd55 0.0001 254r 0.01 0.0001 0.1 0.0001 6r 1 2r 0.0001 239r
|
||||
0.01 141r 1 40r
|
||||
|
||||
@@ -1 +1 @@
|
||||
e2-estia_st2m_v02b-1-gc3eae61.wwg
|
||||
e2-estia_bioS_01-3-g2b98cef.wwg
|
||||
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
Reference in New Issue
Block a user