220 lines
8.6 KiB
OpenEdge ABL
220 lines
8.6 KiB
OpenEdge ABL
Input File ESTIA
|
|
c -------------------------------------------------------
|
|
c --------------- CELL CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
read file=geometry/ess_monolith_cells.i echo
|
|
read file=geometry/ess_outside_cells.i echo
|
|
c tallies
|
|
100701 0 -100701:-101701
|
|
100702 0 -100702
|
|
100703 0 -100703
|
|
100704 0 -100704
|
|
100705 0 -100705
|
|
100706 0 -100706
|
|
100707 0 -100707
|
|
100708 0 -100708
|
|
100709 0 -100709
|
|
100710 0 -100710
|
|
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 -------------------------------------------------------
|
|
c --------------- SURFACE CARDS -------------------------
|
|
c -------------------------------------------------------
|
|
c modified ESS geometry
|
|
read file=geometry/ess_monolith_surfaces.i echo
|
|
read file=geometry/ess_outside_surfaces.i echo
|
|
c tallies
|
|
c 700 s 117.0 -155.00 13.5 3.0
|
|
100701 1 s 542.0 -6.3 6.3 3.0
|
|
101701 1 s 542.0 -6.3 -6.3 3.0
|
|
100702 1 s 1010 -0.8 0 5.0
|
|
100703 1 s 1082 -25.0 0 5.0
|
|
100704 1 s 1082 -0.25 0 5.0
|
|
100705 1 s 1165 0.7 0 5.0
|
|
100706 s 778 -1085 -3 12.0
|
|
100707 s 1267 -1769 -13 12.0
|
|
100708 s 746 -1370 0 30.0
|
|
100709 s 1030 -1155 0 30.0
|
|
100710 s 1341 -1717 -13 12.0
|
|
c Estia geometry
|
|
read file=geometry/estia_surfaces.i echo
|
|
|
|
read file=geometry/materials.i echo
|
|
read file=geometry/materials_monolith.i echo
|
|
c
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
c -------------------------------------------------------
|
|
c --------------- TRANSFORM CARDS -----------------------
|
|
c Estia ISCS with origin at first focal point of feeder
|
|
TR1 10.5 -11 13.7 0.58207953 -0.81304008 -0.012217 &
|
|
-0.81310076 -0.58212297 0 &
|
|
0.0071118 -0.00993365 0.99992537 1
|
|
C center of feeder vessel with x-axis parallel to cylinder
|
|
TR2 457.7511 -626.2452 4.3862 0.582079399 -0.813040169 -0.01221704 &
|
|
-0.813101788 -0.582121521 -0.000132635413 &
|
|
0.00700396398 -0.0100109014 0.99992536 1
|
|
c Beam parallel x-axis and origin at middle focus
|
|
TR3 1349.2829 -1880.9922 -14.3991 0.590799503 -0.806726225 -0.0121961161 &
|
|
-0.80678623 -0.590843447 -8.67361738E-19 &
|
|
0.00720599525 -0.0098396585 0.999925625 1
|
|
c Beam parallel x-axis and origin at virtual source
|
|
TR4 650.7875 -905.3441 0.261 0.590799503 -0.806726225 -0.0121961161 &
|
|
-0.80678623 -0.590843447 -8.67361738E-19 &
|
|
0.00720599525 -0.0098396585 0.999925625 1
|
|
c Parallel to c-axis but horizontal and centered at MF
|
|
TR13 1349.2829 -1880.9922 -14.3991 0.58207953 -0.81304008 0.0 &
|
|
-0.81310076 -0.58212297 0 &
|
|
0.0 0.0 1.0 1
|
|
c Parallel to c-axis but horizontal and centered at VS
|
|
TR14 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 Parallel to c-axis but horizontal and centered at sample
|
|
TR15 2047.78 -2856.64 -29.06 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 Transformation from SKADI model coordinate system here
|
|
TR102 0 5.0 1.0 0.66913061 -0.74314483 0 0.74314483 0.66913061 0 0 0 1 1
|
|
read file=geometry/ess_monolith_transformations.i echo
|
|
c -------------------------------------------------------
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
c -------------------------------------------------------
|
|
c --------------- WEIGHT CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
c weight window generator, can be run with previous weight window files
|
|
WWGE:n 1e-3 1e-1 1 10 100 3e3
|
|
c mesh based weight windows, can be reused after change of cell structure
|
|
WWG 15 0 0
|
|
c MESH GEOM=CYL ORIGIN=0 0 -500 REF=133.7 -176.7 13.5 AXS=0 0 1
|
|
MESH GEOM=CYL ORIGIN=0 0 -500 REF=0 -1 2.5 AXS=0 0 1
|
|
IMESH=20 200 550 1150 1500 2112.01 2157.17 4600
|
|
IINTS=2 2 5 1 10 1 1 4
|
|
JMESH=100 300 460 500 540 700 900 1200
|
|
JINTS= 2 4 4 8 8 4 4 6
|
|
KMESH=298 303 309 315 360
|
|
KINTS= 14 1 20 1 2
|
|
c generate cell based weight windows, needs to be updated on geometry change
|
|
c WWG 55 1042 0
|
|
c simple weight window game for other particles, high inside monolith
|
|
c source weight is 1e16, give large window with low split/roulette factor
|
|
c check only on collisions
|
|
c WWE:/,d,t,s,a 3e3
|
|
c WWP:/,d,t,s,a 20 10 2 -1 0
|
|
c WWN1:/,d,t,s,a 1e15 276r -1 -1 1e10 578r
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
c -----------------------------------------------------------
|
|
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/bunker_mesh_act.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/NBPI_f4.i echo
|
|
c read file=tallies/BBG_f4.i echo
|
|
read file=tallies/BBG_f5.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
|
|
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/shutter_wall_f4.i echo
|
|
c read file=tallies/shutter_sample_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/optical_cave_gammas.i echo
|
|
c read file=tallies/experimental_cave_f5.i echo
|
|
c read file=tallies/cave_mesh.i echo
|
|
c read file=tallies/selene1_f6.i echo
|
|
c
|
|
c ==========================================
|
|
c
|
|
c ==========================================
|
|
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//
|
|
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 1r 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 28r
|
|
c 0.5V2 4r 0 4r 0.5V2 3r 0 0.5V2 13r 0 82r
|
|
c -0.5V1 49r 0 24r -0.1V1 4r 0 13r 0.1V3 9r $ chopper pit to bunker wall
|
|
c 0 2r 0.1V3 1r 0 58r
|
|
c forced collisions in neutron guides
|
|
c fcl:n 0 528r -0.05 -0.1 4r 0 39r -0.1 0 13r -0.25 5r 0 171r
|
|
c -1 0 10r -1 -1 0 22r -1 0 8r $218r
|
|
c exponential transform reference vectors: origin, selene 1, inst. shutter
|
|
c vect v1 0 0 0 v2 1000 -1394 -6.5 v3 1282 -1793 -13
|
|
c
|
|
c -----------------------------------------------------------
|
|
c -------------------------------------------------------
|
|
c --------------- PHYSICS CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
c dbcn 375642321 j j j j j j j j j j j 15291711
|
|
nps 12e7
|
|
c histp
|
|
mode n p h /
|
|
d $ deuteron
|
|
t $ triton
|
|
s $ helion
|
|
a $ alpha
|
|
c in ESS model from target:
|
|
c MODE N $ neutron
|
|
c P $ photon
|
|
c H $ proton
|
|
c / $ pion-+
|
|
c | $ mu--
|
|
c U $ nu-e
|
|
c V $ nu-mu
|
|
c Z $ pion-0
|
|
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
|
|
c -234
|
|
imp:h,/,d,t,s,a 1 276r 0 0 1 578r
|
|
imp:n,p 1 276r 0 0 1 578r
|
|
cut:h,/,d,t,s,a 1e+08 1e-3
|
|
phys:n 3000 0
|
|
phys:p 100
|
|
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
|
|
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)
|
|
read file=sources/proton.i echo
|
|
c * Source Term from Valentina for 2m adapted to size
|
|
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
|
|
c read file=sources/mcstas_source_sample.i echo
|
|
c read file=sources/mcstas_source_mf.i echo
|
|
c read file=sources/mcstas_gammas.i echo
|
|
c read file=sources/mcstas_gammas_2.i echo
|
|
c read file=sources/mcstas_source_shutter.i echo
|
|
prdmp 2e7 -60 1 2 1e5
|
|
print 30 32 40 $ 120 128 150 160 161 162
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
|