From 9bfce4ee1c65496f06b718a1fb9fb54d16f3a1a2 Mon Sep 17 00:00:00 2001 From: Stefano Segantin <92783079+SteSeg@users.noreply.github.com> Date: Thu, 6 Mar 2025 16:33:26 -0500 Subject: [PATCH] Determine nuclides correctly for DAGMC models in d1s.get_radionuclides (#3335) Co-authored-by: Paul Romano --- openmc/deplete/d1s.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/openmc/deplete/d1s.py b/openmc/deplete/d1s.py index 22fc6ec23d..8f6bffa865 100644 --- a/openmc/deplete/d1s.py +++ b/openmc/deplete/d1s.py @@ -33,8 +33,10 @@ def get_radionuclides(model: openmc.Model, chain_file: str | None = None) -> lis List of nuclide names """ - # Determine what nuclides appear in model - model_nuclides = set(model.geometry.get_all_nuclides()) + + # Determine what nuclides appear in the model + model_nuclides = {nuc for mat in model._materials_by_id.values() + for nuc in mat.get_nuclides()} # Load chain file if chain_file is None: