!****************************************************************************
! モデリング
!****************************************************************************

!----------------------------------------------------------------------------
! 物体のモデリング
!----------------------------------------------------------------------------
subroutine modeling
   use fdtd
   implicit none

   integer :: i,j

   ! すべて真空 (=1) に初期化する
   do i=1,nx
      do j=1,ny
         media_id(i,j)=1
      end do
   end do

   ! 金属板
   i=nx/2
   do j=1,ny/2
      media_id(i,j)=2
   end do

   return
end subroutine

!----------------------------------------------------------------------------
! モデリングの確認
!----------------------------------------------------------------------------
subroutine test_modeling
   use fdtd
   implicit none

   integer :: fp    ! ファイル識別子
   integer :: i,j

   open(fp,file="test_modeling.dat")
   do j=1,ny
      write(fp,*) (media_id(i,j),i=1,nx)
      write(fp,'()')    ! 改行
      write(fp,'()')    ! 改行
   end do
   close(fp)

   return
end subroutine

!
! End of file
!