!====================================================================================
!                           model
!===================================================================================
 subroutine model
   use fdtd
   implicit none

   integer  i,j

   !設(shè)置真空背景 
   do i=1,nx
      do j=1,ny
         media_id(i,j)=1
      end do
   end do

   !設(shè)置光柵
   do i=153,nx
      do j=100,400
         media_id(i,j)=2
      end do
   end do

   do i=150,153
      do j=178,250
         media_id(i,j)=2
      end do
   end do

   do i=150,153
      do j=328,400
         media_id(i,j)=2
	  end do 
   end do

   return
end subroutine

!=============================================================================
!                         輸出model
!=============================================================================
subroutine test_model
   use fdtd
   implicit none

   integer :: fp     
   integer :: i,j

   open(fp,file="test_model.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