!============================================================================================
!                        absorbing_boundary_condition
!=============================================================================================
subroutine absorbing_boundary_condition
use fdtd
implicit none


integer i,j
!mur 1
do j=1,ny
   hz(1,j)=hz(2,j)+mur1_cx*(hz(2,j)-hz(1,j))             !左邊界
   hz(nx,j)=hz(nx+1,j)+mur1_cx*(hz(nx+1,j)-hz(nx,j))     !右邊界
end do
do i=1,nx
   hz(i,1)=hz(i,2)+mur1_cy*(hz(i,2)-hz(i,1))             !下邊界
   hz(i,ny)=hz(i,ny+1)+mur1_cy*(hz(i,ny+1)-hz(i,ny))     !上邊界
end do 

return
end subroutine