Hi
Is this the right approach for building the mesh tally to calculate and plot the axial power distribution?
mesh1 = openmc.RegularMesh()
mesh1.dimension = [1, 1, 1]
mesh1.lower_left = [-75.2626, -75.2626, 11.365]
mesh1.upper_right = [75.2626, 75.2626, 14.92]
mesh2 = openmc.RegularMesh()
mesh2.dimension = [1, 1, 1]
mesh2.lower_left = [-75.2626, -75.2626, 14.92]
mesh2.upper_right = [75.2626, 75.2626, 19.365]
mesh3 = openmc.RegularMesh()
mesh3.dimension = [1, 1, 4]
mesh3.lower_left = [-75.2626, -75.2626, 19.365]
mesh3.upper_right = [75.2626, 75.2626, 65.974]
mesh4 = openmc.RegularMesh()
mesh4.dimension = [1, 1, 1]
mesh4.lower_left = [-75.2626, -75.2626, 65.974]
mesh4.upper_right = [75.2626, 75.2626, 70.419]
mesh5 = openmc.RegularMesh()
mesh5.dimension = [1, 1, 4]
mesh5.lower_left = [-75.2626, -75.2626, 70.419]
mesh5.upper_right = [75.2626, 75.2626, 117.028]
mesh6 = openmc.RegularMesh()
mesh6.dimension = [1, 1, 1]
mesh6.lower_left = [-75.2626, -75.2626, 117.028]
mesh6.upper_right = [75.2626, 75.2626, 121.473]
mesh7 = openmc.RegularMesh()
mesh7.dimension = [1, 1, 4]
mesh7.lower_left = [-75.2626, -75.2626, 121.473]
mesh7.upper_right = [75.2626, 75.2626, 168.082]
mesh8 = openmc.RegularMesh()
mesh8.dimension = [1, 1, 1]
mesh8.lower_left = [-75.2626, -75.2626, 168.082]
mesh8.upper_right = [75.2626, 75.2626, 172.527]
mesh9 = openmc.RegularMesh()
mesh9.dimension = [1, 1, 4]
mesh9.lower_left = [-75.2626, -75.2626, 172.527]
mesh9.upper_right = [75.2626, 75.2626, 211.365]
Instantiate tally filters for each mesh
mesh_filter1 = openmc.MeshFilter(mesh1)
mesh_filter2 = openmc.MeshFilter(mesh2)
mesh_filter3 = openmc.MeshFilter(mesh3)
mesh_filter4 = openmc.MeshFilter(mesh4)
mesh_filter5 = openmc.MeshFilter(mesh5)
mesh_filter6 = openmc.MeshFilter(mesh6)
mesh_filter7 = openmc.MeshFilter(mesh7)
mesh_filter8 = openmc.MeshFilter(mesh8)
mesh_filter9 = openmc.MeshFilter(mesh9)
Create a tally to score the axial power distribution
tally = openmc.Tally(name=‘mesh tally’)
tally.filters = [mesh_filter1, mesh_filter2,
mesh_filter3, mesh_filter4,
mesh_filter5, mesh_filter6,
mesh_filter7, mesh_filter8, mesh_filter9]
tally.scores = [‘kappa-fission’]
