165 lines
6.1 KiB
OpenEdge ABL
165 lines
6.1 KiB
OpenEdge ABL
Input File ESTIA
|
|
c -------------------------------------------------------
|
|
c --------------- CELL CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
read file=geometry/ess_cells.i echo
|
|
c tallies
|
|
701 9 2.45e-05 -701
|
|
702 0 -702
|
|
703 0 -703
|
|
704 0 -704
|
|
705 0 -705
|
|
706 0 -706
|
|
707 0 -707
|
|
708 0 -708
|
|
709 0 -709
|
|
710 0 -710
|
|
read file=geometry/estia_inbunker.i echo
|
|
read file=geometry/estia_guides.i echo
|
|
read file=geometry/estia_cave.i echo
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
|
|
c -------------------------------------------------------
|
|
c --------------- SURFACE CARDS -------------------------
|
|
c -------------------------------------------------------
|
|
c modified ESS geometry
|
|
read file=geometry/ess_surfaces.i echo
|
|
c tallies
|
|
c 700 s 117.0 -155.00 13.5 3.0
|
|
701 1 s 542.0 -6.3 6.3 3.0
|
|
702 1 s 1010 -0.8 0 5.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
|
|
707 1 s 2190 3.5 0 12.0
|
|
708 s 730 -1380 0 30.0
|
|
709 s 1020 -1170 0 30.0
|
|
710 s 1328 -1756 -13 12.0
|
|
8370 s 1211 -1660 -11 10.0
|
|
8371 s 1211 -1660 20 10.0
|
|
8372 s 1153 -1700 -11 10.0
|
|
c Estia geometry
|
|
read file=geometry/estia_surfaces.i echo
|
|
|
|
read file=geometry/materials.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 1352 -1889 -12 0.590799503 -0.806726225 -0.0121961161 &
|
|
-0.80678623 -0.590843447 -8.67361738E-19 &
|
|
0.00720599525 -0.0098396585 0.999925625 1
|
|
c Transformation from SKADI model coordinate system here
|
|
TR4 0 5.0 1.0 0.66913061 -0.74314483 0 0.74314483 0.66913061 0 0 0 1 1
|
|
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 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 -------------------------------------------------------
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
c -------------------------------------------------------
|
|
c --------------- WEIGHT CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
c weight window generator, can be run with previous weight window files
|
|
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=133.7 -176.7 13.5 AXS=0 0 1
|
|
c IMESH=200 550 1150 1500 5000
|
|
c IINTS= 2 5 1 10 5
|
|
c JMESH=100 300 460 500 540 700 900 1000
|
|
c JINTS= 2 4 4 8 8 4 4 2
|
|
c KMESH=298 315 360
|
|
c KINTS= 14 3 2
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
c -----------------------------------------------------------
|
|
c ------------------- TALLY CARDS ---------------------------
|
|
c
|
|
c read file=tallies/bunker_mesh.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/selene1_walls_f4.i echo
|
|
c read file=tallies/optical_cave_f5.i echo
|
|
c read file=tallies/experimental_cave_f5.i echo
|
|
c read file=tallies/cave_mesh.i echo
|
|
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//
|
|
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 33r
|
|
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
|
|
c -----------------------------------------------------------
|
|
c -------------------------------------------------------
|
|
c --------------- PHYSICS CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
c dbcn 375642321 j j j j j j j j j j j 15291711
|
|
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 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 3j 2.05e5 1 12r $ wall segments
|
|
imp:h,/,d,t,s,a 0 0 1 694r
|
|
imp:n,p 0 0 1 694r
|
|
cut:h,/,d,t,s,a 1e+08 1e-3
|
|
c cut:n j 0.99E-6
|
|
cut:n j 1e-11
|
|
cut:p 1e+08 0.5
|
|
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 * 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
|
|
read file=sources/e2_2m.i echo
|
|
c * McStas source terms for gamma doeses
|
|
c read file=sources/mcstas_source_vs.i echo
|
|
prdmp 5e7 -60 1 2 5e6
|
|
print 30 32 $ 120 128 150 160 161 162
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
|