Fixed stacking error for 3D calculation

1 files changed, 2 insertions, 2 deletions

diff --git a/hpmr.py b/hpmr.py
index 477830a..ef818cc 100755
--- a/hpmr.py
+++ b/hpmr.py
@@ -713,8 +713,8 @@ geo_core_refl_bot_cell = openmc.Cell(name="Bottom Axial Reflector Cell",fill=geo
 # 3D Stacking Lattice
 geo_core_axial_lattice = openmc.RectLattice()
 geo_core_axial_lattice.lower_left = (-sys.float_info.max/2, -sys.float_info.max/2, -int(160//geo_const_height)/2*geo_const_height)
-geo_core_axial_lattice.pitch=(sys.float_info.max, sys.float_info.max, geo_const_height)
-geo_core_axial_lattice.universes = [[[geo_core_rad_universe]] for _ in range(int(160//geo_const_height))]
+geo_core_axial_lattice.pitch=(sys.float_info.max, sys.float_info.max, 2*geo_const_height)
+geo_core_axial_lattice.universes = [[[geo_core_rad_universe]] for _ in range(int(160//(2*geo_const_height)))]
 geo_core_axial_lattice.outer=geo_none_universe
 
 geo_core_cell = openmc.Cell(name="Core 3D Cell",fill=geo_core_axial_lattice, region=-geo_const_core_cyl & +geo_const_core_z_min & -geo_const_core_z_max)