boundary.f90

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

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