3d hex_lattice bugfixes

This commit is contained in:
Sterling Harper 2014-09-05 14:04:17 -04:00
parent 8aec6c9f6f
commit c3d4d243a9
2 changed files with 4 additions and 4 deletions

View file

@ -1569,7 +1569,7 @@ contains
! Top and bottom sides.
if (lat % is_3d) then
z0 = sign(lat % pitch(3) * 0.5_8, w)
z0 = sign(lat % pitch(2) * 0.5_8, w)
if (abs(z - z0) < FP_PRECISION) then
d = INFINITY

View file

@ -1301,8 +1301,8 @@ contains
end do RECT_LATTICES
HEX_LATTICES: do i = 1, n_hlats
allocate(HexLattice::lattices(i) % obj)
lat => lattices(i) % obj
allocate(HexLattice::lattices(n_rlats + i) % obj)
lat => lattices(n_rlats + i) % obj
select type (lat)
type is (HexLattice)
@ -1477,7 +1477,7 @@ contains
end if
! Add lattice to dictionary
call lattice_dict % add_key(lat % id, i)
call lattice_dict % add_key(lat % id, n_rlats + i)
end select
end do HEX_LATTICES