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

Updating c5g7 model

Browse source
This commit is contained in:
Adam Nelson 2018-04-27 09:33:05 -04:00
parent dfefc46019
commit 2147331fcd
3 changed files with 12 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

4
mg-criticality/make_library.py
View file

@ -552,10 +552,10 @@ def create_6g():
[0.000, 0.000, 2.000, 0.00, 0.00, 0.00],
[0.000, 0.000, 0.000, 2.00, 0.00, 0.00],
[0.000, 0.000, 0.000, 0.275, 0.60, 0.00],
[0.000, 0.000, 2.000, 0.033, 0.171, 0.024]]])
[0.000, 0.000, 0.000, 0.033, 0.171, 0.024]]])
scatter = np.rollaxis(scatter, 0, 3)
total = [0.240, 0.975, 3.10, 3.10, 0.975, 0.240]
chi = [0.48, 0.02, 0.0, 0.0, 0.02, 0.048]
chi = [0.48, 0.02, 0.0, 0.0, 0.02, 0.48]
URR = set_it('URR', groups, 0, fiss, nu, absorption, scatter, total, chi)
mg_cross_sections_file.add_xsdata(URR)

7
mg-criticality/repo.py
View file

@ -3,10 +3,11 @@
from make_model import Case
# GLOBAL DATA
batches = 5000
inactive = 500
particles = 10000
batches = 600
inactive = 100
particles = 100000
tab_leg = {'enable': True, 'num_points': 33}
# tab_leg = {'enable': False, 'num_points': 33}
def build_cases():

6
mg-criticality/run_cases.py
View file

@ -2,6 +2,7 @@
from __future__ import print_function
from optparse import OptionParser
import pickle
import matplotlib.pyplot as plt
import numpy as np
@ -32,6 +33,7 @@ if options.list_cases:
print('Case Name: {0}'.format(key))
exit()
results = []
# Run specific case, if requested
if options.case_name:
if options.case_name in repo.names:
@ -48,6 +50,7 @@ if options.case_name:
print(case.name + ' Failed Execution!')
else:
print_case(case)
results.append(case)
else:
print('Invalid Case Name: ' + options.case)
@ -66,3 +69,6 @@ else:
print(case.name + ' Failed Execution!')
else:
print_case(case)
results.append(case)
pickle.dump(results, open("save.pkl", "wb"))