abc.f90

来自「laser_damage_te用fdtd方法研究激光损伤」· F90 代码 · 共 23 行

!============================================================================================
!                        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