!**************************************************************************** ! 補正電界境界条件 !**************************************************************************** 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 ! ! End of file !