Add supermirror reflection for NBOA and NFGA

This commit is contained in:
2020-02-11 15:08:43 +01:00
parent c965b26e75
commit b491a58a21
3 changed files with 20 additions and 5 deletions

View File

@@ -43,6 +43,7 @@ c 700 s 117.0 -155.00 13.5 3.0
c Estia geometry
read file=geometry/estia_surfaces.i echo
read file=geometry/sm_reflectivity.i echo
read file=geometry/materials.i echo
read file=geometry/materials_monolith.i echo
c
@@ -143,17 +144,17 @@ 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
read file=tallies/full_mesh_np.i echo
c read file=tallies/full_mesh_np.i echo
c read file=tallies/boundary_meshes.i echo
c read file=tallies/boundary_wall_meshes.i echo
c read file=tallies/selene_source_mesh_p.i echo
c read file=tallies/NBPI_f4.i echo
c read file=tallies/BBG_f4.i echo
c read file=tallies/BBG_f5.i echo
c read file=tallies/inbunker_f4.i echo
read file=tallies/inbunker_f4.i echo
c read file=tallies/chopper_pit_doses.i echo
c read file=tallies/bunker_feedthrough_f1.i echo
c read file=tallies/bunker_wall_f4.i echo
read file=tallies/bunker_wall_f4.i echo
c read file=tallies/bunker_wall_f5.i echo
c read file=tallies/pre_selene1_mesh.i
c read file=tallies/e1_e2_f4.i echo

View File

@@ -0,0 +1,14 @@
c Specify the reflecitive surfaces for neutron supermirrors
c NBOA
REFLE1 201501 36 -201001 -201002 -201003 -201004
-201005 -201006 -201007 -201008
REFLE2 201502 36 -201001 -201002 -201003 -201004
-201005 -201006 -201007 -201008
c NFGA
REFLE3 203501 36 -203001 -203002 -203003 -203004 -203005 -203006 -203007
-203011 -203012 -203013 -203014 -203015 -203016 -203017
REFLE4 203502 36 -203001 -203002 -203003 -203004 -203005 -203006 -203007
-203011 -203012 -203013 -203014 -203015 -203016 -203017
c Selene 1
c Selene 2
REFF36 0.99 2.19e-2 3.6 6.07 3.0e-3

4
plot
View File

@@ -3,13 +3,13 @@ PLOT_NAME=`git describe --tags --long`
# if the first argument is an runtpe file, use different plot command
if [[ $1 == runs/* ]]; then
mcnp6 p runtpe=$1 \
/home/ll_glavic/shielding/mcnp6/MCNP_SM/bin/mcnp6 p runtpe=$1 \
outp=tmp/plota.out comout=tmp/plota.com \
wwout=tmp/plota.wwg wwone=tmp/plota.wwo \
wwinp=weight_windows/current.wwg \
plotm=results/e2-estia_$PLOT_NAME ${@:2}
else
mcnp6 ip inp=e2-estia.i \
/home/ll_glavic/shielding/mcnp6/MCNP_SM/bin/mcnp6 ip inp=e2-estia.i \
outp=tmp/plota.out comout=tmp/plota.com \
wwout=tmp/plota.wwg wwone=tmp/plota.wwo \
mctal=tmp/plota.mct mdata=tmp/plota.mdt \