crp/openmc-designs
1
0
Fork 1
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

fixed c5g7

Browse source
This commit is contained in:
Adam Nelson 2019-10-19 19:15:58 -05:00
parent 955a075e37
commit 448e2ffd12
3 changed files with 7 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

10
mg-criticality/make_model.py
View file

@ -1,5 +1,3 @@
#!/usr/bin/env python
import numpy as np
import openmc
@ -225,10 +223,10 @@ class Case(object):
surfs = []
for r, rad in enumerate(self.rad):
if r == len(self.rad) - 1:
surfs.append(openmc.ZCylinder(R=rad,
surfs.append(openmc.ZCylinder(r=rad,
boundary_type='vacuum'))
else:
surfs.append(openmc.ZCylinder(R=rad))
surfs.append(openmc.ZCylinder(r=rad))
# Instantiate Cells
cells = []
@ -244,9 +242,9 @@ class Case(object):
surfs = []
for r, rad in enumerate(self.rad):
if r == len(self.rad) - 1:
surfs.append(openmc.Sphere(R=rad, boundary_type='vacuum'))
surfs.append(openmc.Sphere(r=rad, boundary_type='vacuum'))
else:
surfs.append(openmc.Sphere(R=rad))
surfs.append(openmc.Sphere(r=rad))
# Instantiate Cells
cells = []

2
mg-criticality/repo.py
View file

@ -1,5 +1,3 @@
#!/usr/bin/env python
from make_model import Case
# GLOBAL OPTIONS

5
mg-criticality/run_cases.py
View file

@ -83,8 +83,9 @@ with open("results.csv", mode="w") as csv_file:
# write the header
writer.writerow(['Case ID', 'Case Name', 'Eigenvalue Bias [pcm]',
'Error Code'])
'Eigenvalue Bias Std. Dev. [pcm]', 'Error Code'])
# Now run each case
for i in range(len(case_nums)):
writer.writerow([case_nums[i], case_names[i], biases[i], codes[i]])
writer.writerow([case_nums[i], case_names[i], biases[i].n,
biases[i].s, codes[i]])