!****************************************************************************
! 磁界の計算
!****************************************************************************
subroutine magnetic_field
   use fdtd
   implicit none

   integer :: i,j,id

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

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

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

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

   return
end subroutine

!
! End of file
!