!****************************************************************************
! 電界の計算
!****************************************************************************
subroutine electric_field
   use fdtd
   implicit none

   integer :: i,j,id

   ! Ez
      do j=2,ny
         do i=2,nx
            id=media_id(i,j)

            if(id.eq.1) then
               ! 1: 自由空間
               ez(i,j)=cez0*ez(i,j) &
                      +cezrx0*(hy(i,j)-hy(i-1,j)) &
                      -cezry0*(hx(i,j)-hx(i,j-1))
            else if(id.eq.2) then
               ! 2: 完全導体
               ez(i,j)=0.0d0
            else
               ! 3以上: 任意媒質
               ez(i,j)=cez(id)*ez(i,j) &
                      +cezrx(id)*(hy(i,j)-hy(i-1,j)) &
                      -cezry(id)*(hx(i,j)-hx(i,j-1))
            end if
         end do
      end do

   return
end subroutine

!
! End of file
!