!****************************************************************************
! 出力
!****************************************************************************
subroutine output
   use fdtd
   implicit none

   integer :: i,j,k

   open(10,file="field.dat")
   j=ny/2
   do k=2,nz-1
      write(10,*) (dabs(hy(i,j,k))**2,i=2,nx-1)
      write(10,"()")    ! 改行
      write(10,"()")    ! 改行
   end do
   close(10)

   open(10,file="field.dat")
   k=nz/2
   do j=2,ny-1
      write(10,*) (dabs(ez(i,j,k))**2,i=2,nx-1)
      write(10,"()")    ! 改行
      write(10,"()")    ! 改行
   end do
   close(10)

   return
end subroutine

!
! End of file
!