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

Operator aborts in unpacking, not call, if no energy produced

Browse source
This commit is contained in:
Andrew Johnson 2019-09-17 12:17:49 -05:00
parent 6b386f8447
commit cfbd4f2fd6
No known key found for this signature in database
GPG key ID: 253418E91B7F6FEB
1 changed files with 7 additions and 11 deletions
Download patch file Download diff file Expand all files Collapse all files

18
openmc/deplete/operator.py
View file

@ -277,16 +277,7 @@ class Operator(TransportOperator):
time_openmc = time.time()
# Extract results
try:
op_result = self._unpack_tallies_and_normalize(power)
except ZeroDivisionError:
if comm.rank == 0:
sys.stderr.flush()
print(" No energy reported from openmc tallies. Do you have "
"MT901 data?\n", file=sys.stderr, flush=True)
comm.barrier()
comm.Abort(1)
op_result = self._unpack_tallies_and_normalize(power)
return copy.deepcopy(op_result)
@ -671,7 +662,12 @@ class Operator(TransportOperator):
# Guard against divide by zero
if energy == 0:
raise ZeroDivisionError
if comm.rank == 0:
sys.stderr.flush()
print(" No energy reported from OpenMC tallies. Do your HDF5 "
"files have heating data?\n", file=sys.stderr, flush=True)
comm.barrier()
comm.Abort(1)
# Scale reaction rates to obtain units of reactions/sec
rates *= power / energy