BLOCKS: Nullify new neighbor fields in deallocate_metablock().

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2014-07-05 13:50:48 -03:00
parent 5262f29e3a
commit 1b3eb9c933

View File

@ -986,7 +986,7 @@ module blocks
! local variables
!
integer :: i, j, k
integer :: n, i, j, k
!
!-------------------------------------------------------------------------------
!
@ -1015,6 +1015,32 @@ module blocks
nullify(pmeta%child(i)%ptr)
end do
! nullify fields pointing to face, edge, and corner neighbors
!
#if NDIMS == 2
do i = 1, nsides
do j = 1, nsides
do n = 1, ndims
nullify(pmeta%edges(n,i,j)%ptr)
end do ! ndims
nullify(pmeta%corners(i,j)%ptr)
end do ! nsides
end do ! nsides
#endif /* NDIMS == 2 */
#if NDIMS == 3
do i = 1, nsides
do j = 1, nsides
do k = 1, nsides
do n = 1, ndims
nullify(pmeta%faces(n,i,j,k)%ptr)
nullify(pmeta%edges(n,i,j,k)%ptr)
end do ! ndims
nullify(pmeta%corners(i,j,k)%ptr)
end do ! nsides
end do ! nsides
end do ! nsides
#endif /* NDIMS == 3 */
! nullify fields pointing to neighbors
!
do i = 1, ndims