diff --git a/src/geometry.F90 b/src/geometry.F90 index c85c4241e2..b661bcbcd7 100644 --- a/src/geometry.F90 +++ b/src/geometry.F90 @@ -220,7 +220,7 @@ contains p % coord % lattice_x = i_x p % coord % lattice_y = i_y p % coord % lattice_z = i_z - p % coord % universe = lat % element(i_x,i_y,i_z) + p % coord % universe = lat % universes(i_x,i_y,i_z) end if call find_cell(found) @@ -613,7 +613,7 @@ contains end if else ! Find universe for next lattice element - p % coord % universe = lat % element(i_x, i_y, i_z) + p % coord % universe = lat % universes(i_x, i_y, i_z) ! Find cell in next lattice element call find_cell(found) diff --git a/src/geometry_header.F90 b/src/geometry_header.F90 index 5dda5e4d85..abbad38628 100644 --- a/src/geometry_header.F90 +++ b/src/geometry_header.F90 @@ -26,10 +26,10 @@ module geometry_header integer :: type ! Type of lattice (rectangular, hex, etc) integer :: level ! Level of lattice integer :: n_dimension ! Number of dimensions - integer, allocatable :: dimension(:) ! number of cells in each direction - real(8), allocatable :: lower_left(:) ! lower-left corner of lattice - real(8), allocatable :: width(:) ! width of each lattice cell - integer, allocatable :: element(:,:,:) ! specified universes + integer, allocatable :: dimension(:) ! number of cells in each direction + real(8), allocatable :: lower_left(:) ! lower-left corner of lattice + real(8), allocatable :: width(:) ! width of each lattice cell + integer, allocatable :: universes(:,:,:) ! specified universes end type Lattice !=============================================================================== diff --git a/src/hdf5_interface.F90 b/src/hdf5_interface.F90 index 67b6005d33..1785971e4f 100644 --- a/src/hdf5_interface.F90 +++ b/src/hdf5_interface.F90 @@ -388,12 +388,12 @@ contains n_z = 1 end if - ! Write lattice elements + ! Write lattice universes allocate(lattice_universes(n_x, n_y, n_z)) do j = 1, n_x do k = 1, n_y do m = 1, n_z - lattice_universes(j,k,m) = universes(lat % element(j,k,m)) % id + lattice_universes(j,k,m) = universes(lat % universes(j,k,m)) % id end do end do end do diff --git a/src/initialize.F90 b/src/initialize.F90 index fdca8549e5..ff3a2090f0 100644 --- a/src/initialize.F90 +++ b/src/initialize.F90 @@ -502,9 +502,9 @@ contains do m = 1, n_z do k = 1, n_y do j = 1, n_x - id = lat % element(j,k,m) + id = lat % universes(j,k,m) if (universe_dict % has_key(id)) then - lat % element(j,k,m) = universe_dict % get_key(id) + lat % universes(j,k,m) = universe_dict % get_key(id) else message = "Invalid universe number " // trim(to_str(id)) & // " specified on lattice " // trim(to_str(lat % id)) diff --git a/src/input_xml.F90 b/src/input_xml.F90 index f0a8acd5da..34b8268694 100644 --- a/src/input_xml.F90 +++ b/src/input_xml.F90 @@ -928,11 +928,11 @@ contains else n_z = 1 end if - allocate(lat % element(n_x, n_y, n_z)) + allocate(lat % universes(n_x, n_y, n_z)) ! Check that number of universes matches size if (size(lattice_(i) % universes) /= n_x*n_y*n_z) then - message = "Number of universes on does not match size of & + message = "Number of universes on does not match size of & &lattice " // trim(to_str(lat % id)) // "." call fatal_error() end if @@ -941,7 +941,7 @@ contains do m = 1, n_z do k = 0, n_y - 1 do j = 1, n_x - lat % element(j, n_y - k, m) = lattice_(i) % & + lat % universes(j, n_y - k, m) = lattice_(i) % & universes(j + n_x*k + n_x*n_y*(m-1)) end do end do