Implement hole in cabinet wall and run on mcc200

This commit is contained in:
2025-01-06 16:59:12 +01:00
parent fc4b141804
commit 443c1f4dc5
5 changed files with 22 additions and 15 deletions

View File

@@ -220,7 +220,7 @@ c -------------------------------------------------------
c dbcn 375642321 j j j j j j j j j j j 15291711
c stop calculation either at XeY histories, 6.5 days or tally 15 error <0.001
c nps 39242646 $ works with rssa file from McStas, original tracks 39242646
stop nps 1e8 CTME 4e5 $f31 0.01
stop nps 1e6 CTME 4e5 $f31 0.01
c nps 28e11 $ works with rssa file, value is scaled to 20x source particles
c write particle tracks
c ssw 206505 pty=n $ record all neutron tracks crossing outward
@@ -272,8 +272,8 @@ c read file=sources/mcstas_rssa.i echo
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
read file=sources/mcstas_gammas.i echo
c read file=sources/mcstas_gammas_2.i echo
c read file=sources/mcstas_gammas.i echo
read file=sources/mcstas_gammas_2.i echo
c read file=sources/mcstas_source_shutter.i echo
prdmp -120 -30 1 2 5e5
print 30 32 40

View File

@@ -231,15 +231,16 @@ c Selene 2 Vacuum tank (stainless steel)
-208932) #(208626 208636 -208932)
c #(8626 8636 -8932)
c Air outside the guide shielding
208601 4 -1.2E-3 100610 -208709 -100104
208601 4 -1.2E-3 (100610 -208709 -100104
((208232 -208223):(100101 208223))
208203 208251
208203 208251 ):(100059 -208235 -208236 208237 208231 -208232)
208602 49 0.074981 100003 100072 -100610 100101 -100104 100067 208203
c Thin wall shielding the cabinet area
208603 48 -2.25 100610 -208709 208231 -208232 100059 -208233 308234
208604 4 -1.2E-3 100610 -208709 208231 -208232 208233 -208223 308234
208603 48 -2.25 100610 -208709 208231 -208232 100059 -208233 208234
#(-208235 -208236 208237)
208604 4 -1.2E-3 100610 -208709 208231 -208232 208233 -208223 208234
208605 4 -1.2E-3 (100610:(-100101 208250)) -208709
208231 -208232 100059 -208223 -308234
208231 -208232 100059 -208223 -208234
c Tube after Selene 2 with sapphire window
208611 0 208636 208641 -208932
208612 5 0.0582256 -208626 208642 -208931 208932

View File

@@ -468,8 +468,12 @@ c Wall E01 side for cabinets
208231 py 289.4
208232 py 309.4
c diagonal limiting element from door
208233 pz 13.7
308234 p 0.813101 0.582123 0 2169.0
208233 pz 90.0
208234 p 0.813101 0.582123 0 2169.0
c cut-out for cables
208235 pz -180.0
208236 p -1.0000002116649709 -1.148348833801549e-07 0.0 -2845.6615
208237 p -1.0000002116649709 -1.148348833801549e-07 0.0 -3003.6615
c New Gemoetry
c SKADI wall
c Sides

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
/home/l_glavic/shielding/mcnp6/MCNP_SM/bin/mcnp6 p runtpe=$1 \
/mnt/nfs/home/glavic_a/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
/home/l_glavic/shielding/mcnp6/MCNP_SM/bin/mcnp6 ip inp=e2-estia.i \
/mnt/nfs/home/glavic_a/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 \

View File

@@ -1,13 +1,13 @@
#!/usr/bin/env tcsh
#SBATCH -J mcnpEstia
#SBATCH -N 5
#SBATCH -N 6
#SBATCH --ntasks-per-node=24
#SBATCH --time=7-00:00:00
#SBATCH --mail-type=fail
#SBATCH --mail-user=artur.glavic@psi.ch
#SBATCH -o tmp/stdout.log
#SBATCH -e tmp/stderr.log
#SBATCH --partition=long
#SBATCH --partition=week
setenv CURRENT_VERSION $1
#"`/usr/bin/git describe --tags --long`"
@@ -19,7 +19,9 @@ echo "Running on $SLURM_NPROCS processors."
echo "Current working directory is `pwd`"
echo "Estia GIT revision is $CURRENT_VERSION"
mpirun -np $SLURM_NPROCS /home/l_glavic/shielding/mcnp6/MCNP_SM/bin/mcnp6.mpi \
module load intel/22.2
mpiexec -n $SLURM_NPROCS /mnt/nfs/home/glavic_a/shielding/mcnp6/MCNP_SM/bin/mcnp6.mpi \
inp=e2-estia.i outp=results/e2-estia_$CURRENT_VERSION.out \
runtpe=runs/e2-estia_$CURRENT_VERSION.run \
mctal=results/e2-estia_$CURRENT_VERSION.mct \