crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Added 3d capability to hex_lattice

Browse source
This commit is contained in:
Sterling Harper 2014-08-30 21:23:14 -04:00
parent b454fec351
commit 6438c667dc
1 changed files with 4 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

9
src/geometry.F90
View file

@ -306,10 +306,9 @@ contains
xyz_t(2) = xyz(2) - (lat % center(2) + &
&(2 * (i_xyz(2) - lat % n_rings)) * lat % pitch(1) + &
&(i_xyz(1) - lat % n_rings) * lat % pitch(1))
! TODO transition to center from lower_left
if (lat % is_3d) then
xyz_t(3) = xyz(3) - (lat % center(3) + &
&(i_xyz(3) - 0.5_8) * lat % pitch(3))
xyz_t(3) = xyz(3) - lat % center(3) &
&+ (lat % n_axial/2 - i_xyz(3) + 1) * lat % pitch(2)
else
xyz_t(3) = xyz(3)
end if
@ -404,9 +403,9 @@ contains
type is (HexLattice)
! Index z direction.
! TODO transition to center from lower_left
if (lat % is_3d) then
i_xyz(3) = ceiling((xyz(3) - lat % center(3)) / lat % pitch(2))
i_xyz(3) = ceiling((xyz(3) - lat % center(3))/lat % pitch(2) + 0.5_8) &
&+ lat % n_axial/2
else
i_xyz(3) = 1
end if