From 4ede491f4d9ada7b9b4c4dc1795c3c8f47281a0a Mon Sep 17 00:00:00 2001
From: Jonathan Shimwell <mail@jshimwell.com>
Date: Thu, 1 Feb 2024 17:59:30 +0000
Subject: [PATCH] neutron activation draft looks ok

---
 ..._mesh_based_shut_down_dose_rate_example.py | 33 ++++++++++---------
 1 file changed, 17 insertions(+), 16 deletions(-)

diff --git a/tasks/task_11_CSG_shut_down_dose_tallies/3_mesh_based_shut_down_dose_rate_example.py b/tasks/task_11_CSG_shut_down_dose_tallies/3_mesh_based_shut_down_dose_rate_example.py
index 8f8d4e8f..7b48c08b 100644
--- a/tasks/task_11_CSG_shut_down_dose_tallies/3_mesh_based_shut_down_dose_rate_example.py
+++ b/tasks/task_11_CSG_shut_down_dose_tallies/3_mesh_based_shut_down_dose_rate_example.py
@@ -192,6 +192,8 @@
         voxel_mat = openmc.Material.mix_materials(materials_in_voxel, norm_fracs)
     else:
         voxel_mat = openmc.Material()
+    
+    voxel_mat.volume = sum(volumes_in_voxel)
     material_in_voxel.append(voxel_mat)
 
 openmc.lib.finalize()
@@ -199,26 +201,25 @@
 
 
 # # # # constructing the operator, note we pass in the flux and micro xs
-# operator = openmc.deplete.IndependentOperator().from_nuclides(
-#     volume=volume_of_material_in_voxel
-#     materials=
-#     fluxes=flux_in_each_group,
-#     micros=micro_xs,
-#     reduce_chain=True,  # reduced to only the isotopes present in depletable materials and their possible progeny
-#     reduce_chain_level=5,
-#     normalization_mode="source-rate"
-# )
+operator = openmc.deplete.IndependentOperator(
+    materials=openmc.Materials(material_in_voxel),
+    fluxes=[flux[0] for flux in flux_in_each_group],
+    micros=micro_xs,
+    reduce_chain=True,  # reduced to only the isotopes present in depletable materials and their possible progeny
+    reduce_chain_level=5,
+    normalization_mode="source-rate"
+)
 
-# integrator = openmc.deplete.PredictorIntegrator(
-#     operator=operator,
-#     timesteps=timesteps,
-#     source_rates=source_rates,
-#     timestep_units='s'
-# )
+integrator = openmc.deplete.PredictorIntegrator(
+    operator=operator,
+    timesteps=timesteps,
+    source_rates=source_rates,
+    timestep_units='s'
+)
 
 # # this runs the depletion calculations for the timesteps
 # # this does the neutron activation simulations and produces a depletion_results.h5 file
-# integrator.integrate()
+integrator.integrate()
 # # TODO add output dir to integrate command so we don't have to move the file like this
 # # integrator.integrate(path=statepoints_folder / "neutrons" / "depletion_results.h5")
 # # PR on openmc is open