Added basic flux tally over coreHEADmaster

1 files changed, 8 insertions, 0 deletions

diff --git a/hpmr.py b/hpmr.py
index 8713518..91b4b2d 100755
--- a/hpmr.py
+++ b/hpmr.py
@@ -824,6 +824,14 @@ tal_spectrum = openmc.Tally(name="Global Energy Spectrum")
 tal_spectrum.filters = [openmc.EnergyFilter(np.logspace(np.log10(1.0),np.log10(20e6),500))]
 tal_spectrum.scores = ["flux"]
 
+tal_flux_mesh = openmc.RegularMesh()
+tal_flux_mesh.lower_left=(-set_dyn_radial,-set_dyn_radial,-set_dyn_axial)
+tal_flux_mesh.upper_right=(+set_dyn_radial,+set_dyn_radial,+set_dyn_axial)
+tal_flux_mesh.dimension=(300,300,300)
+tal_flux = openmc.Tally(name="2D flux tally")
+tal_flux.filters = [openmc.MeshFilter(tal_flux_mesh)]
+tal_flux.scores = ["flux"]
+
 tallies = openmc.Tallies([tal_global, tal_spectrum])
 tallies.export_to_xml()