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

Merge with source_term_2m, insert mesh tally and test without WW

Browse Source
This commit is contained in:
glavic_a
2018-10-02 16:05:01 +02:00
parent 185ac9cb44
commit dcba60afa4
8 changed files with 18269 additions and 24 deletions
Download Patch File Download Diff File Expand all files Collapse all files

52
e2-estia.i
View File

@@ -695,16 +695,24 @@ c Parallel to c-axis but horizontal and centered at VS
TR5 650.7875 -905.3441 0.261 0.58207953 -0.81304008 0.0 &
-0.81310076 -0.58212297 0 &
0.0 0.0 1.0 1
c Coordinate system of the Estia beamport focal point
TR100 8.9 -5.4 13.7 0.58778525 -0.80901699 0.0 &
-0.80901699 -0.58778525 0 &
0.0 0.0 1.0 1
c Same coordinates but with origin at TCS
TR101 0.0 0.0 0.0 0.58778525 -0.80901699 0.0 &
-0.80901699 -0.58778525 0 &
0.0 0.0 1.0 1
c -------------------------------------------------------
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
c -------------------------------------------------------
c --------------- WEIGHT CARDS --------------------------
c -------------------------------------------------------
WWP:n,p j j 100 j -1
c weight windows from external file generated by WWG
c WWP:n,p j j 100 j -1
c weight window generator, can be run with previous weight window files
c WWGE:n,p 1e-6 1e-3 1 3e3
c WWG:n,p 4 0 0
c c WWG:n,p 5 0 0
c WWGE:n,p 1e-3 1e-1 1e1 3e3
c WWG:n,p 5 0 0
c MESH GEOM=CYL ORIGIN=0 0 -500 REF=0 -8 0.5 AXS=0 0 1
c IMESH=200 550 1150 1500 5000
c IINTS= 2 5 1 10 5
@@ -723,6 +731,7 @@ tmesh
corc1 -200 1i -100 5i 100 2i 300
endmd
fm1 8.7374e15
c
c
c ==========================================
c
@@ -730,18 +739,21 @@ 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 115r
c pd0 1e-5 254r 0.1 73r 1e-5 179r 1 64r 1 113r
c has to be for specific detector as pd0 does not work properly
c 1-255 256 257 268 259-265 266-268 269-508 509-end
c pd5 0.001 254r 0.1 0.001 1 0.001 6r 1 2r 0.001 239r 0.1 110r 1 70r
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 //
c 495 496 497 498 499 500 501 502 503 504 505 506 507 508-537 //
c 538 539 540-543 544-554 555-585 //
c 586 587-591 592 593 594-601 602 603-end//
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.5V2 2r 0 222r 0.02V3 2r 0 219r 0.25V2 5r
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 26r
c 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
c 0.5V2 2r 0 222r 0.02V3 2r 0 219r 0.25V2 5r
c 0 0.5V2 0 0.5V2 0 0 0.5V2 0 0.5V2 0 0.5V2 0 0.5V2 0 22r
c 0.5V2 10r 0 4r 0.5V2 3r 0 0.5V2 13r 0 48r
c -0.5V1 2r 0 4r -0.1V1 9r 0 0.1V3 3r 0 11r 0.1V3 2r
c 0 0.1V3 0 0.1V3 0 0.1V3 5r 0 26r
vect v1 0 0 0 v2 1352 -1889 0 v3 1282 -1793 -13
c
c -----------------------------------------------------------
@@ -749,7 +761,7 @@ c -------------------------------------------------------
c --------------- PHYSICS CARDS --------------------------
c -------------------------------------------------------
c dbcn 375642321 j j j j j j j j j j j 15291711
nps 4e6
nps 2e7
c histp
mode n p n h / d t s a
c calculate volume for tally cells (spheres)
@@ -759,7 +771,7 @@ imp:n,p 0 0 1 687r
cut:h,/,d,t,s,a 1e+08 1e-3
c cut:n j 0.99E-6
cut:n j 1e-11
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
cut:p 1e+08 0.5
phys:n 3000 0
phys:p 100
@@ -767,17 +779,9 @@ phys:/,d,t,s,a 3000
phys:h 3000
lca 2 1 1 0023 1 1 0 1 1 0
lea 1 4 1 0 1 0 0 1
sdef dir=1 erg=2000 par=9 vec=0.0 1 0.0 x=d1 y=-30 z=d2
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
0.80859 0.85938 0.90234 0.9375 0.96484 0.98438 0.99609 1
0.99609 0.98438 0.96484 0.9375 0.90234 0.85938 0.80859 0.75
0.68359 0.60938 0.52734 0.4375 0.33984 0.23438 0.12109 0
si2 A -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3
sp2 D 0 0.30556 0.55556 0.75 0.88889 0.97222 1 0.97222 0.88889 0.75
0.55556 0.30556 0
prdmp 1e6 -60 1 2 1e5
c read file=source_proton.i echo
read file=e2-source_2m.i echo
prdmp 5e6 -60 1 2 1e6
print 30 32 $ 120 128 150 160 161 162
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++

1055
e2-source_2m.i Normal file
View File

File diff suppressed because it is too large Load Diff

3435
ess_reference/tally_181.txt Normal file
View File

File diff suppressed because it is too large Load Diff

3435
ess_reference/tally_191.txt Normal file
View File

File diff suppressed because it is too large Load Diff

3435
ess_reference/tally_201.txt Normal file
View File

File diff suppressed because it is too large Load Diff

3435
ess_reference/tally_211.txt Normal file
View File

File diff suppressed because it is too large Load Diff

3435
ess_reference/tally_221.txt Normal file
View File

File diff suppressed because it is too large Load Diff

11
source_proton.i Normal file
View File

@@ -0,0 +1,11 @@
c Original definition of source from proton beam
sdef dir=1 erg=2000 par=9 vec=0.0 1 0.0 x=d1 y=-30 z=d2
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
0.80859 0.85938 0.90234 0.9375 0.96484 0.98438 0.99609 1
0.99609 0.98438 0.96484 0.9375 0.90234 0.85938 0.80859 0.75
0.68359 0.60938 0.52734 0.4375 0.33984 0.23438 0.12109 0
si2 A -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3
sp2 D 0 0.30556 0.55556 0.75 0.88889 0.97222 1 0.97222 0.88889 0.75
0.55556 0.30556 0