!**************************************************************************** ! ␳dEE !**************************************************************************** subroutine electric_boundary_condition use fdtd implicit none integer :: i,j i=nx/2+1 do j=1,ny/2 ez(i,j)=0.0d0 end do i=nx j=ny/2 ez(i,j)=0.0d0 i=nx+1 j=ny/2+1 ez(i,j)=0.0d0 return end subroutine !**************************************************************************** ! ␳EE !**************************************************************************** subroutine magnetic_boundary_condition use fdtd implicit none integer :: i,j j=1 do i=1,nx hx(i,j)=0.0d0 end do j=ny do i=1,nx hx(i,j)=0.0d0 end do return end subroutine ! ! End of file !