Tally data over multiple burnup steps

Hello,

I am currently running some comparative studies between OpenMC, SERPENT, and Polaris, and one of the parameters I wanted to compare was the pin peaking factor with respect to burnup.

I believe I have setup the tallies I want correctly, but I am currently unable to extract these same tallies for subsequent burnup steps. I tried using the openmc ‘get_tally’ command to extract the data, but I am only able to extract the values at BOC. Opening my ‘tallies.out’ file also indicates that no tallies are recorded beyond those at BOC.

So my question is, how exactly do I extract my tally data for subsequent burnup steps?

I’ve attached the build script (ASSMO.py) and the post-processing script for convenience.

Thanks in advance.

ASSMO.py (32.4 KB)
PostProcess.py (8.7 KB)

There should be multiple data files present of the form openmc_simulation_n<N>.h5 present in your working directory. Each file contains the tally data for step N, e.g., openmc_simulation_n0.h5, openmc_simulation_n1.h5, …

Each of these files contains the tally data for user declared tallies for that step.

Drew