crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 13:45:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Merge pull request #1659 from liangjg/fix_dose_energy

Browse source 
fix typo in energy unit in dose coefficients
This commit is contained in:
Paul Romano 2020-09-14 11:44:46 -05:00 • committed by GitHub
commit 87f2b0e928
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
2 changed files with 7 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

2
openmc/data/effective_dose/dose.py
View file

@ -24,7 +24,7 @@ def _load_dose_icrp116():
for particle, filename in _FILES:
path = Path(__file__).parent / filename
data = np.loadtxt(path, skiprows=3)
data[:, 0] *= 1e-6 # Change energies to eV
data[:, 0] *= 1e6 # Change energies to eV
_DOSE_ICRP116[particle] = data

12
tests/unit_tests/test_data_dose.py
View file

@ -5,21 +5,21 @@ from pytest import approx, raises
def test_dose_coefficients():
# Spot checks on values from ICRP tables
energy, dose = dose_coefficients('photon', 'AP')
assert energy[0] == approx(0.01e-6)
assert energy[0] == approx(0.01e6)
assert dose[0] == approx(0.0685)
assert energy[-1] == approx(10e-3)
assert energy[-1] == approx(10e9)
assert dose[-1] == approx(90.4) # updated in corrigendum
energy, dose = dose_coefficients('neutron', 'LLAT')
assert energy[0] == approx(1e-15)
assert energy[0] == approx(1e-3)
assert dose[0] == approx(1.04)
assert energy[-1] == approx(10e-3)
assert energy[-1] == approx(10e9)
assert dose[-1] == approx(1.23e3)
energy, dose = dose_coefficients('electron', 'ISO')
assert energy[0] == approx(0.01e-6)
assert energy[0] == approx(0.01e6)
assert dose[0] == approx(0.0188)
assert energy[-1] == approx(10e-3)
assert energy[-1] == approx(10e9)
assert dose[-1] == approx(699.0)
# Invalid particle/geometry should raise an exception