Hello everyone!
When I tried to read the statepoint.h5 file, I got the error message quoted in the subject line of this post. Why did I receive this message, and how do I fix the problem? Any assistance will be appreciated. thanks.
The complete error message goes as:
495/1 1.21424 1.21479 +/- 0.00005
496/1 1.21424 1.21479 +/- 0.00005
497/1 1.21666 1.21480 +/- 0.00005
498/1 1.21424 1.21479 +/- 0.00005
499/1 1.21594 1.21480 +/- 0.00005
500/1 1.21349 1.21479 +/- 0.00005
Creating state point statepoint.500.h5...
HDF5-DIAG: Error detected in HDF5 (1.10.6) thread 140047653668672:
#000: H5A.c line 284 in H5Acreate2(): unable to create attribute
major: Attribute
minor: Unable to initialize object
#001: H5Aint.c line 275 in H5A__create(): unable to create attribute in object header
major: Attribute
minor: Unable to insert object
#002: H5Oattribute.c line 335 in H5O__attr_create(): unable to create new attribute in header
major: Attribute
minor: Unable to insert object
#003: H5Omessage.c line 222 in H5O__msg_append_real(): unable to create new message
major: Object header
minor: No space available for allocation
#004: H5Omessage.c line 1933 in H5O__msg_alloc(): unable to allocate space for message
major: Object header
minor: Unable to initialize object
#005: H5Oalloc.c line 1312 in H5O__alloc(): object header message is too large
major: Object header
minor: Unable to initialize object
HDF5-DIAG: Error detected in HDF5 (1.10.6) thread 140047653668672:
#000: H5A.c line 620 in H5Awrite(): not an attribute
major: Invalid arguments to routine
minor: Inappropriate type
HDF5-DIAG: Error detected in HDF5 (1.10.6) thread 140047653668672:
#000: H5A.c line 1528 in H5Aclose(): not an attribute
major: Invalid arguments to routine
minor: Inappropriate type
=======================> TIMING STATISTICS <=======================
Total time for initialization = 7.5945e+00 seconds
Reading cross sections = 4.9610e+00 seconds
Total time in simulation = 1.0470e+06 seconds
Time in transport only = 1.0465e+06 seconds
Time in inactive batches = 5.9320e+04 seconds
Time in active batches = 9.8765e+05 seconds
Time synchronizing fission bank = 1.7903e+02 seconds
Sampling source sites = 1.5099e+02 seconds
SEND/RECV source sites = 2.8036e+01 seconds
Time accumulating tallies = 2.0102e+02 seconds
Time writing statepoints = 7.6662e+00 seconds
Total time for finalization = 2.5537e-01 seconds
Total time elapsed = 1.0470e+06 seconds
Calculation Rate (inactive) = 6743.06 particles/second
Calculation Rate (active) = 607.501 particles/second
============================> RESULTS <============================
k-effective (Collision) = 1.21479 +/- 0.00005
k-effective (Track-length) = 1.21479 +/- 0.00005
k-effective (Absorption) = 1.21486 +/- 0.00004
Combined k-effective = 1.21483 +/- 0.00003
Leakage Fraction = 0.11252 +/- 0.00001
---------------------------------------------------------------------------
KeyError Traceback (most recent call last)
###############################################################################
547 # Load the last statepoint file
--> 548 with openmc.StatePoint(statepoint_newfile) as sp:
549 #-----------------------------------------------------------------------------#
550 # get keff form .h5 file
551 #-----------------------------------------------------------------------------#
552 if Current_Version_Newer_flag == True:
553 MC_keff = sp.keff
File ~/anaconda3/envs/fenicsproject/lib/python3.8/site-packages/openmc/statepoint.py:144, in StatePoint.__init__(self, filepath, autolink)
142 if os.path.exists(path_summary):
143 su = openmc.Summary(path_summary)
--> 144 self.link_with_summary(su)
146 path_volume = os.path.join(os.path.dirname(filename), 'volume_*.h5')
147 for path_i in glob.glob(path_volume):
File ~/anaconda3/envs/fenicsproject/lib/python3.8/site-packages/openmc/statepoint.py:664, in StatePoint.link_with_summary(self, summary)
660 raise ValueError(msg)
662 cells = summary.geometry.get_all_cells()
--> 664 for tally in self.tallies.values():
665 tally.with_summary = True
667 for tally_filter in tally.filters:
File ~/anaconda3/envs/fenicsproject/lib/python3.8/site-packages/openmc/statepoint.py:366, in StatePoint.tallies(self)
363 # Read a list of the IDs for each Tally
364 if n_tallies > 0:
365 # Tally user-defined IDs
--> 366 tally_ids = tallies_group.attrs['ids']
367 else:
368 tally_ids = []
File h5py/_objects.pyx:54, in h5py._objects.with_phil.wrapper()
File h5py/_objects.pyx:55, in h5py._objects.with_phil.wrapper()
File ~/anaconda3/envs/fenicsproject/lib/python3.8/site-packages/h5py/_hl/attrs.py:56, in AttributeManager.__getitem__(self, name)
52 @with_phil
53 def __getitem__(self, name):
54 """ Read the value of an attribute.
55 """
---> 56 attr = h5a.open(self._id, self._e(name))
58 if is_empty_dataspace(attr):
59 return Empty(attr.dtype)
File h5py/_objects.pyx:54, in h5py._objects.with_phil.wrapper()
File h5py/_objects.pyx:55, in h5py._objects.with_phil.wrapper()
File h5py/h5a.pyx:80, in h5py.h5a.open()
KeyError: "Can't open attribute (can't locate attribute: 'ids')"