fdtd1.f90

求目标散射体建模

!-----------建模--------------
Program main
real ob(-700:700,-700:700)
real Media(1:4,1:15)
character*10 Filename
real WaveLength,WL
real Ismin,Ismax,Jsmin,Jsmax
real MediaNo
integer Record,i,j

!------------Initiate------------
Filename='rect.sct'
WaveLength=1.
WL=40.
Ismin=-60
Ismax=60
Jsmin=-60
Jsmax=60
MediaNo=1	                !只有一种介质
Media(1,1)=1.			    !epsilon_r
Media(2,1)=1.		    	!mu_r
Media(3,1)=3.72e+7		    !sigma
Media(4,1)=0.			    !sigma_m

!------------建立散射体------------
do i=Ismin,Ismax			!空间所有网格初始值设为背景
   do j=Jsmin,Jsmax
     ob(i,j)=0.0
   end do
end do
do i=Ismin,Ismax
   do j=Jsmin,Jsmax
    if((i.ge.-WL).and.(i.le.WL-1).and.(j.ge.-WL).and.(j.le.WL-1)) then	 !判断网格是否在目标上
	 ob(i,j)=1.0
	end if
   end do
end do
!------------输出--------------
open(1,File=FileName,access='direct',recl=4)
write(1,rec=1)WaveLength
write(1,rec=2)WL
write(1,rec=3)Ismin
write(1,rec=4)Ismax
write(1,rec=5)Jsmin
write(1,rec=6)Jsmax
write(1,rec=7)MediaNo
Record=7
do j=1,MediaNo
  do i=1,4
   Record=Record+1
   write(1,rec=Record)Media(i,j)
  end do
end do
do i=Ismin,Ismax
  do j=Jsmin,Jsmax
      Record=Record+1
	  write(1,rec=Record)ob(i,j)
  end do
end do
close(1)
write(*,*)Ismin,Ismax,Jsmin,Jsmax
end