!**************************************************************************** ! 電界の計算 !**************************************************************************** subroutine electric_field use fdtd implicit none integer :: i,j,id ! Ex do j=2,ny do i=1,nx id=media_id(i,j) if(id.eq.1) then ! 1: 自由空間 ex(i,j)=cex0*ex(i,j) & +cexry0*(hz(i,j)-hz(i,j-1)) else if(id.eq.2) then ! 2: 完全導体 ex(i,j)=0.0d0 else ! 3以上: 任意媒質 ex(i,j)=cex(id)*ex(i,j) & +cexry(id)*(hz(i,j)-hz(i,j-1)) end if end do end do ! Ey do j=1,ny do i=2,nx id=media_id(i,j) if(id.eq.1) then ! 1: 自由空間 ey(i,j)=cey0*ey(i,j) & -ceyrx0*(hz(i,j)-hz(i-1,j)) else if(id.eq.2) then ! 2: 完全導体 ey(i,j)=0.0d0 else ! 3以上: 任意媒質 ey(i,j)=cey(id)*ey(i,j) & -ceyrx(id)*(hz(i,j)-hz(i-1,j)) end if end do end do return end subroutine ! ! End of file !