Skip to content

Warn when ACE average heating numbers are written as zeros #27

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.

Sign up for GitHub

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

Open
jon-proximafusion wants to merge 1 commit into
base: master
Choose a base branch
from
Open

Warn when ACE average heating numbers are written as zeros #27

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
4 changes: 4 additions & 0 deletions brownies/LANL/toACE/gndsToACE.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,9 +10,12 @@
-) When delayed neutrons are present, their multiplicities are added to prompt, and prompt and total nu_bar
are written. However, their decay constants, etc are currently not written.
-) Only neutrons distribution data are outputted.
-) The average heating numbers (KERMA) in the ESZ block are not calculated and are written as zeros. There is
currently no equivalent of NJOY's HEATR module in FUDGE. A warning is emitted when the zeros are written.
"""

import time
import warnings

from xData import enums as xDataEnumsModule
from xData import axes as axesModule
Expand Down Expand Up @@ -217,6 +220,7 @@ def toACE(self, styleLabel, cdf_style, fileName, evaluationId, productData, dela
else:
updateXSSInfo( 'absorption cross section', annotates, XSS, mapCrossSectionToGrid(absorptionXSec) )
updateXSSInfo( 'elastic cross section', annotates, XSS, elasticXSec )
warnings.warn('Average heating numbers (KERMA) are not calculated and are written to the ESZ block as zeros.')
averageHeating = len( energyGrid ) * [ 0. ]
updateXSSInfo( 'average heating', annotates, XSS, averageHeating )

Expand Down