!****************************************************************************
! �o��
!****************************************************************************
subroutine output
   use fdtd
   implicit none

   integer :: i,j,k

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

1  format(/)        ! ���s
   return
end subroutine

!
! End of file
!