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

Generate weight window mesh for NBPI window

Browse Source
This commit is contained in:
glavic_a
2019-05-31 20:20:38 +02:00
parent 44c6e53750
commit 4cefa58ab3
2 changed files with 15 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
e2-estia.i
View File

@@ -93,16 +93,17 @@ c -------------------------------------------------------
c --------------- WEIGHT CARDS --------------------------
c -------------------------------------------------------
c weight window generator, can be run with previous weight window files
c WWGE:n 1e-3 1e-1 1 10 100 3e3
WWGE:n 1e-3 1e-1 1 10 100 3e3
c mesh based weight windows, can be reused after change of cell structure
c WWG 55 0 0
WWG 514 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 2112.01 2157.17 4600
c IINTS= 2 5 1 10 1 1 4
c JMESH=100 300 460 500 540 700 900 1200
c JINTS= 2 4 4 8 8 4 4 6
c KMESH=298 303 309 315 360
c KINTS= 14 1 20 1 2
MESH GEOM=CYL ORIGIN=0 0 -500 REF=0 0 13.4 AXS=0 0 1
IMESH=200 550 1150 1500 2112.01 2157.17 4600
IINTS= 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 ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
@@ -160,7 +161,7 @@ c -------------------------------------------------------
c --------------- PHYSICS CARDS --------------------------
c -------------------------------------------------------
c dbcn 375642321 j j j j j j j j j j j 15291711
nps 1e5
nps 1e7
c histp
mode n p h /
d $ deuteron
@@ -204,7 +205,7 @@ 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 5e8 -60 1 2 5e7
prdmp 1e6 -60 1 2 1e5
print 30 32 40 $ 120 128 150 160 161 162
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++

7
sources/proton.i
View File

@@ -1,11 +1,12 @@
c Original definition of source from proton beam
cut:n j 0.99e-06
c cut:n j 0.99e-06
cut:n j 0
cut:p 1e+08 0.5
c weight windows from external file generated by WWG
c make sure to supply the correct WWINP file
c 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
c esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000
c energy splitting for higher neutron energies and roulette below 1keV
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=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