diff --git a/docs/source/usersguide/output/statepoint.rst b/docs/source/usersguide/output/statepoint.rst index 02f4fdcf3..d1ebc7231 100644 --- a/docs/source/usersguide/output/statepoint.rst +++ b/docs/source/usersguide/output/statepoint.rst @@ -134,11 +134,7 @@ if run_mode == 'k-eigenvalue': **/tally/meshes/keys** (*int[]*) - User-identified unique ID of each mesh - -**/tallies/meshes/mesh /id** (*int*) - - Unique identifier of the mesh. + User-identified unique ID of each mesh. **/tallies/meshes/mesh /type** (*char[]*) diff --git a/docs/source/usersguide/output/summary.rst b/docs/source/usersguide/output/summary.rst index 0623693c5..6bf0cf2ad 100644 --- a/docs/source/usersguide/output/summary.rst +++ b/docs/source/usersguide/output/summary.rst @@ -243,7 +243,7 @@ The current revision of the summary file format is 1. **/tallies/mesh /index** (*int*) - Index in the meshes array used internally in OpenMC + Index in the meshes array used internally in OpenMC. **/tallies/mesh /type** (*char[]*) diff --git a/openmc/statepoint.py b/openmc/statepoint.py index cdb82d95e..133bd766c 100644 --- a/openmc/statepoint.py +++ b/openmc/statepoint.py @@ -265,8 +265,7 @@ class StatePoint(object): # Iterate over all Meshes for mesh_key in mesh_keys: - # Read the user-specified Mesh ID and type - mesh_id = self._f['{0}{1}/id'.format(base, mesh_key)].value + # Read the mesh type mesh_type = self._f['{0}{1}/type'.format(base, mesh_key)].value.decode() # Read the mesh dimensions, lower-left coordinates, @@ -277,7 +276,7 @@ class StatePoint(object): width = self._f['{0}{1}/width'.format(base, mesh_key)].value # Create the Mesh and assign properties to it - mesh = openmc.Mesh(mesh_id) + mesh = openmc.Mesh(mesh_key) mesh.dimension = dimension mesh.width = width mesh.lower_left = lower_left @@ -285,7 +284,7 @@ class StatePoint(object): mesh.type = mesh_type # Add mesh to the global dictionary of all Meshes - self._meshes[mesh_id] = mesh + self._meshes[mesh_key] = mesh self._meshes_read = True diff --git a/src/state_point.F90 b/src/state_point.F90 index a5c89a8a2..64ba7ef55 100644 --- a/src/state_point.F90 +++ b/src/state_point.F90 @@ -176,7 +176,6 @@ contains meshp => meshes(id_array(i)) mesh_group = create_group(meshes_group, "mesh " // trim(to_str(meshp%id))) - call write_dataset(mesh_group, "id", meshp%id) select case (meshp%type) case (MESH_REGULAR) call write_dataset(mesh_group, "type", "regular")