Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Update MF7 recipe to latest changes in ENDF-102 #3

Closed
wants to merge 2 commits into from
Closed
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 2 additions & 1 deletion endf_parserpy/endf_recipes/__init__.py
Original file line number Diff line number Diff line change
Expand Up @@ -10,7 +10,7 @@
from .endf_recipe_mf4 import ENDF_RECIPE_MF4
from .endf_recipe_mf5 import ENDF_RECIPE_MF5
from .endf_recipe_mf6 import ENDF_RECIPE_MF6
from .endf_recipe_mf7 import ENDF_RECIPE_MF7_MT2, ENDF_RECIPE_MF7_MT4
from .endf_recipe_mf7 import ENDF_RECIPE_MF7_MT2, ENDF_RECIPE_MF7_MT4, ENDF_RECIPE_MF7_MT451
from .endf_recipe_mf8 import ENDF_RECIPE_MF8
from .endf_recipe_mf8_mt454 import ENDF_RECIPE_MF8_MT454
from .endf_recipe_mf8_mt457 import ENDF_RECIPE_MF8_MT457
Expand Down Expand Up @@ -61,6 +61,7 @@ def __defdic(dic, index, recipe):
__defdic(endf_recipe_dictionary, (6,), ENDF_RECIPE_MF6)
__defdic(endf_recipe_dictionary, (7, 2), ENDF_RECIPE_MF7_MT2)
__defdic(endf_recipe_dictionary, (7, 4), ENDF_RECIPE_MF7_MT4)
__defdic(endf_recipe_dictionary, (7, 451), ENDF_RECIPE_MF7_MT451)
__defdic(endf_recipe_dictionary, (8, -1), ENDF_RECIPE_MF8)
__defdic(endf_recipe_dictionary, (8, 454), ENDF_RECIPE_MF8_MT454)
__defdic(endf_recipe_dictionary, (8, 457), ENDF_RECIPE_MF8_MT457)
Expand Down
20 changes: 20 additions & 0 deletions endf_parserpy/endf_recipes/endf_recipe_mf7.py
Original file line number Diff line number Diff line change
Expand Up @@ -22,6 +22,13 @@
# Incoherent Elastic Scattering
elif LTHR==2:
[MAT, 7, 2/ SB, 0.0, 0, 0, NR, NP/ Tint / Wp ]TAB1
elif LTHR==3:
[MAT, 7, 2/ T0, 0.0, LT, 0, NR, NP/ Eint / S ] TAB1 (S_T0_table)
for i=1 to LT:
[MAT, 7, 2/ T[i], 0.0, LI, 0, NP, 0/
{S[q,i]}{q=1 to NP} ] LIST
endfor
[MAT, 7, 2/ SB, 0.0, 0, 0, NR, NP/ Tint / Wp ]TAB1
endif
SEND
"""
Expand Down Expand Up @@ -52,3 +59,16 @@
endif
SEND
"""

ENDF_RECIPE_MF7_MT451 = """

# Generalized Information File
[MAT, 7, 451 / ZA, AWR, NA, 0, 0, 0]HEAD
for i=1 to NA:
[MAT, 7, 451 / 0, 0, NAS, 0, 6*NI[i], NI[i]/
{ZAI[i,j], LISI[i,j], AFI[i,j], AWRI[i,j], SFI[i,j], 0}{j=1 to NI[i]} ] LIST
endfor
SEND
"""