subroutine old_neighbors(child_igrid,child_ipe,igrid,ipe)
use mod_global_parameters
integer, dimension(2,2,2), intent(in) :: child_igrid, child_ipe
integer, intent(in) :: igrid, ipe
integer :: iside, i1,i2,i3, ic1,ic2,ic3
do ic3=1,2
do ic2=1,2
do ic1=1,2
old_neighbor(:,:,:,:,child_igrid(ic1,ic2,ic3))=-1
do iside=1,2
if (ic1==iside) then
i1=kr(1,1)*(2*iside-3);i2=kr(2,1)*(2*iside-3)
i3=kr(3,1)*(2*iside-3);
old_neighbor(1,i1,i2,i3,child_igrid(ic1,ic2,ic3))=fine_neighbors(ic1,&
ic2,ic3,1,igrid)%igrid
old_neighbor(2,i1,i2,i3,child_igrid(ic1,ic2,ic3))=fine_neighbors(ic1,&
ic2,ic3,1,igrid)%ipe
end if
end do
do iside=1,2
if (ic2==iside) then
i1=kr(1,2)*(2*iside-3);i2=kr(2,2)*(2*iside-3)
i3=kr(3,2)*(2*iside-3);
old_neighbor(1,i1,i2,i3,child_igrid(ic1,ic2,ic3))=fine_neighbors(ic1,&
ic2,ic3,2,igrid)%igrid
old_neighbor(2,i1,i2,i3,child_igrid(ic1,ic2,ic3))=fine_neighbors(ic1,&
ic2,ic3,2,igrid)%ipe
end if
end do
do iside=1,2
if (ic3==iside) then
i1=kr(1,3)*(2*iside-3);i2=kr(2,3)*(2*iside-3)
i3=kr(3,3)*(2*iside-3);
old_neighbor(1,i1,i2,i3,child_igrid(ic1,ic2,ic3))=fine_neighbors(ic1,&
ic2,ic3,3,igrid)%igrid
old_neighbor(2,i1,i2,i3,child_igrid(ic1,ic2,ic3))=fine_neighbors(ic1,&
ic2,ic3,3,igrid)%ipe
end if
end do
end do
end do
end do
end subroutine old_neighbors