search

openmc-dev / openmc

OpenMC Monte Carlo Code
https://docs.openmc.org
Other
699 stars 444 forks source link

Track file not generated on particle restart run #2956

Closed Edgar-21 closed 1 week ago

Edgar-21 commented 2 weeks ago

Bug Description

When restarting openmc from a particle.h5 restart file, with tracks=True, no track file is written.

Steps to Reproduce

The following code will lose some particles, then attempt to restart from a particle.h5 with tracks=True, and only a summary.h5 file will be outputted.

import openmc
model = openmc.examples.pwr_pin_cell()
model.settings.max_tracks = 20

# change the boundary condition of some surfaces to induce lost particles
for surf in model.geometry.get_all_surfaces().values():
    if isinstance(surf, openmc.XPlane):
        surf.boundary_type = 'transmission'

model.run()

# try to restart and write a track
openmc.run(restart_file='particle_1_11.h5', tracks=True)

Environment

I have tried this both on the following conda install of openmc, and on the latest version of the develop-dagmc docker image. openmc -v

OpenMC version 0.14.0
Git SHA1: e1a8ee7794b441c992426f17fafe216391cbba83
Copyright (c) 2011-2023 MIT, UChicago Argonne LLC, and contributors
MIT/X license at <https://docs.openmc.org/en/latest/license.html>
Build type:            Release
Compiler ID:           GNU 12.3.0
MPI enabled:           yes
Parallel HDF5 enabled: yes
PNG support:           yes
DAGMC support:         yes
libMesh support:       no
MCPL support:          no
NCrystal support:      no
Coverage testing:      no
Profiling flags:       no