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

   integer :: fp    ! ファイル識別子
   integer :: i,j
   real(8) :: vecabs,xx,yy
   vecabs(xx,yy)=dsqrt(xx**2+yy**2)

   open(fp,file="field.dat")
!   do j=1,ny
!      !write(fp,*) (vecabs(ex(i,j),ey(i,j)),i=1,nx)    ! 磁界の強さ
!      write(fp,*) (dabs(hz(i,j)),i=1,nx)              ! 電界の強さ
!      write(fp,'()')    ! 改行
!      write(fp,'()')    ! 改行
!   end do
   do j=(ny/4),ny-(ny/4)
      write(fp,*) (dabs(hz(i,j)),i=nx/4,nx-nx/4)
      write(fp,'()')    ! 改行
      write(fp,'()')    ! 改行
   end do
   close(fp)

   return
end subroutine

!
! End of file
!