boundary.f90

laser_damage_te用fdtd方法研究激光损伤

!===================================================================================
!                         electric_boundary_condition
!===================================================================================
subroutine electric_boundary_condition
use fdtd
implicit none

integer i,j

j=1
do i=1,nx
   ex(i,j)=0.0d0
end do 

j=ny
do i=1,nx
   ex(i,j)=0.0d0
end do

i=nx/2
do j=1,ny/2
   ey(i,j)=0.0d0
end do

i=nx/2+1
do j=1,ny/2
   ey(i,j)=0.0d0
end do

i=nx/2
j=ny/2+1
ey(i,j)=0.0d0

return 
end subroutine