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

   integer :: i,j
   real(8) :: vecabs,xx,yy
   vecabs(xx,yy)=dsqrt(xx**2+yy**2)

   open(10,file="field.dat")
   do j=1,ny
      write(10,*) (vecabs(ex(i,j),ey(i,j)),i=1,nx)
      write(10,1)
      write(10,1)
   end do
   close(10)

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

!
! End of file
!