crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 06:05:58 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

fix energy unit conversion from Mev to eV

Browse source
This commit is contained in:
Jingang Liang 2020-09-14 11:12:36 +08:00
parent 364422fea8
commit 0c95b93d5c
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