subroutine identifyphysbound_connectivity(igrid)
use mod_global_parameters
integer, intent(in) :: igrid
if (neighbor_type(-1,0,0,igrid) == neighbor_boundary .and. neighbor_type(+1,0,0,igrid) == neighbor_boundary ) then
idphyb(1,igrid) = 2
else if (neighbor_type(-1,0,0,igrid) == neighbor_boundary) then
idphyb(1,igrid) = -1
else if (neighbor_type(+1,0,0,igrid) == neighbor_boundary) then
idphyb(1,igrid) = +1
else
idphyb(1,igrid) = 0
end if
if (neighbor_type(0,-1,0,igrid) == neighbor_boundary .and. neighbor_type(0,+1,0,igrid) == neighbor_boundary ) then
idphyb(2,igrid) = 2
else if (neighbor_type(0,-1,0,igrid) == neighbor_boundary) then
idphyb(2,igrid) = -1
else if (neighbor_type(0,+1,0,igrid) == neighbor_boundary) then
idphyb(2,igrid) = +1
else
idphyb(2,igrid) = 0
end if
if (neighbor_type(0,0,-1,igrid) == neighbor_boundary .and. neighbor_type(0,0,+1,igrid) == neighbor_boundary ) then
idphyb(3,igrid) = 2
else if (neighbor_type(0,0,-1,igrid) == neighbor_boundary) then
idphyb(3,igrid) = -1
else if (neighbor_type(0,0,+1,igrid) == neighbor_boundary) then
idphyb(3,igrid) = +1
else
idphyb(3,igrid) = 0
end if
end subroutine identifyphysbound_connectivity